st: lrtest of nested logistic models with vce(cluster) specified

[Lily Yor <lilyyor1@gmail.com>]

Hi, I have two nested binary logistic regression models, both of which have
vce(cluster) specified to account for intraclass clustering at the classroom
level.  For example, my Stata commands look like this in simplified form:

logit [dep var] [ind_var1], vce(cluster class_id)
est store A

logit [dep var] [ind_var1 ind_var2], vce(cluster class_id)
est store B

I would like to test whether model B fits the data significantly better than
model A, so I ran the lrtest, but Stata gives me this error:

. lrtest A
LR test likely invalid for models with robust vce
How do I test the fit of these nested models if I cannot use lrtest?  I
found a similar question posted by another user earlier in the archives but
it was never answered.  I would much appreciate any and all tips.

Best,
Lily

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