Displaying 1 result from an estimated 1 matches for "anxietyscore".
2012 May 03
2
Very small random effect estimation in lmer but not in stata xtmixed
...nd 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....?
LMER MODEL
summary(lmer(AnxietyScore ~ cc + (1|SetID), data=mydf))
Linear mixed model fit by REML
Formula: AnxietyScore ~ cc + (1 | SetID)
Data: mydf
AIC BIC logLik deviance REMLdev
493.4 503.4 -242.7 486.6 485.4
Random effects:
Groups Name Variance Std.Dev.
SetID (Intercept) 1.2819e-09 3.5803e-05
Residual...