[R-sig-ME] Very small random effect estimation in lmer but not in stata xtmixed
Ben Bolker
bbolker at gmail.com
Fri May 4 14:39:58 CEST 2012
Mohammed Mohammed <M.A.MOHAMMED at ...> writes:
> I am using the lmer function (in the lme4 library) to analyse some
> data where individuals are clustered into sets (using the SetID
> variable) with a single fixed effect (cc - 0 or 1). The lmer model
> and output is shown below.
> Whilst the fixed effects are consistent with stata (using xtmixed,
> see below), the std dev of the random effect for SetID is very very
> small (3.5803e-05)compared to stata's (see below 1.002). Any ideas
> why this should be happening please....?
> I am happy to share my data if that helps.
>
Nothing really obvious pops out. Have you tried
nlme (i.e library(nlme); lme(AnxietyScore ~cc, random=~1|SetID,data=mydf))
as a cross-check?
Can you send/post data?
Ben Bolker
More information about the R-sig-mixed-models
mailing list