st: Using mi estimate with xtmixed and accounting for sampling weights

[Justina Kamiel Grayman <justina.grayman@nyu.edu>]

Hi everyone,

I am attempting to estimate a two-level model with a random intercept.
I am using xtmixed. My data was also imputed, using MI ICE. I now have
10 imputed datasets. Lastly, I am trying to include sampling weights
at both levels of the model. So basically I need to include sampling
weights in a two-level model estimated on imputed data.

The problem: When I perform the analysis (syntax below), it runs fine
on the original data.  It also runs when I use an if statement to run
the analysis on only one of the imputed datasets. However, when I try
to use mi estimate as a prefix, the model does not run. Is there some
incompatibility between mi estimate and using pweights in xtmixed?
Does anyone know why this is occurring?

Thanks!


Syntax That Does Not Work:

mi estimate: ///
xtmixed studeff ///
classclimmeanc classclimvarc ///
pctnonwhitec ///
educc female usborn ///
black latino asian other multi ///
btgtaugc wrkshpcntc someba  ///
bcna9c urban rural northeast southeast west privateschl [pw = totwgt] ///
|| idschool: , covariance(unstructured) pweight(schwgt) , mle pwscale(size)


Syntax That Does Work

xtmixed studeff ///
classclimmeanc classclimvarc ///
pctnonwhitec ///
educc female usborn ///
black latino asian other multi ///
btgtaugc wrkshpcntc someba  ///
bcna9c urban rural northeast southeast west privateschl [pw = totwgt] ///
if  _mi_m==1 ///
|| idschool:  , covariance(unstructured) pweight(schwgt) , mle pwscale(size)
--
Justina Kamiel Grayman
Doctoral Student, Psychology & Social Intervention
New York University
*
*   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/