Hi,

I use Stata 16.1 and I need to run a 3-level Zero Inflated Negative Binomial regression. I have seen some previous answers suggest using "gllamm"

Based on the "gllamm" manual in Stata I wrote the following code:

meglm y x1, x2, x3....xn || lev3: || lev2:, family(nbinomial)

After 11 iterations the log likelihood stops changing.

Does anyone know if I am doing something wrong with the code and how I can get around this.