st: different standard errors with gllamm vs. xtmelogit

["de Vries, Robert" <r.de-vries08@imperial.ac.uk>]

Hello everyone. I'm having a weird problem with gllamm and xtmelogit when r= unning a fairly simple 2-level random intercepts model.

The model is predicting a binary health outcome from 1 level 2 variable (me=
angini) and several level 1 variables. It is a sample if 5,4410 people in 1=
6 countries, with 'country' as the level 2 cluster.

The xtmelogit model looks like this:

xi: xtmelogit poorhealth meangini age47 gndr i.education if poorhealth_samp=
le=3D=3D1 || country1 :

(note that the number of integration points is left at the default of 7)

It converges fine on iteration 3 with a log likelihood if -15046.989. The r= esult I am interested in is for meangini and in this model it is -0.030 (SE=  =3D 0.035)

The gllamm model is identical (as far as I can tell)

xi: gllamm poorhealth meangini age47 gndr i.education if poorhealth_sample= =3D=3D1, i(country1) nip(7) link(logit) f(binomial)

However the coefficient for the same variable is different (-0.041). And th= e Standard Error (0.0043) is over 8 times smaller.

This is obviously extremely important in interpreting the statistical signi= ficance of the results so I'd appreciate any help anyone might be able to o= ffer as to what's going on.

Cheers
Rob


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