Iím trying to specify a multilevel model for count data in Mplus that I have run with MLwiN.
Because Mplus does not permit estimation of variances for count variables at the within level, I specified single-indicator latent variables and correlated the latent factors instead. Below is my model:
It sounds like you are using a random intercept Poisson model in Mplus. Are you sure you are using the same model in MLwiN and not for example negbin? Make sure the programs show the same number of parameters and ompare loglikelihood value and if Mplus has a less high value you can sharpen the convergence criterion (mconv).
Thank you so much for your response. I checked MLwiN and I am indeed specifying a random intercept Poisson with the same number of parameters (although I am using first-order marginal quasi-likelihood for estimation). Whereas the log-likelihood value for Mplus is -5611.421, MLwiN provides -2*log(lh) = 18554.5. I tried sharpening the convergence criterion (MCONVERGENCE = .00001), but still got the same results.
I should also note that when I treat the variables as continuous, I get the same exact results across programs.
Should I try lowering the convergence criterion even more? If so, what would you recommend?