st: comparing two mim gllamm logistic models

[jl591164@albany.edu]

Dear Stata listers,
I use -mim, cat(fit):gllamm- to fit three-level logistic regression
models. Now i need to compare the model fit of two nested models, e.g.
randome intercept models with and without a quadratic term of time.
Another example of two nested models is a randome intercept model and a
randome intercept and random slope model. Do i compare the -2logL, AIC and
BIC of the two models of concern? If so, how do i obtain these model fit
statistics after fitting mim:gllamm models?

If they are simple gllamm models without multiple imputation, i can use
estimates store after fitting each model, then run the -lrtest gllamm
gllamm2- to test which model fit better. However, the estimates store does
not work after fitting mim: gllamm models. I got the error message of
command estimates not recognised by mim; try specifying category option. I
tired -mim, cat(combine) est(e(MIM_Q)): regress Y-, it ran and the output
was Combined estimate = . (standard errors not supplied). But I am not
sure this -combine- is what i need for model comparision.

Anyone knows what is an appropirate way to do model comparision? Thanks a
lot.

Junqing

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/