st: Predicting Random Effects from a Crossed-Level Model using xtmelogit

["Lee, Albert" <Albert.Lee@hud.gov>]

Hi, 

This may be obvious.  If so, I apologize in advance.

I'm estimating a logistic model with crossed-level random effects.  The outcome binary variable is fail with one fixed-effect independent binary variable refinance.  The three crossed-level random effects are VAR1, VAR2 and VAR3.  The STATA command I used is as below:

xtmelogit fail refinance || _all:R.VAR1 || _all:R.VAR2 || _all:R.VAR3 

After estimating this model, I tried to recover the random effects of VAR3 using:

predict u, reffects level(VAR3)

However, u only contains missing values.

I would very much appreciate if someone can help me recovering these predicted random effects.

Thanks,

Albert Lee

*
*   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/