[R-sig-ME] Problematic correlations between random effects
David Duffy
David.Duffy at qimr.edu.au
Thu Feb 23 23:42:11 CET 2012
On Thu, 23 Feb 2012, Leite,Walter wrote:
> Folks,
>
> I fitted a model with several random effects to data with a
> cross-classification of students and schools using lmer. Two random
> effects have correlation of -1 and others near 1, as shown below. I
> thought that this would only occur if any random effects were not
> necessary in the model, but I built this model adding the random effects
> one at a time and comparing the fit with BIC and AIC. Also, none of the
> variances are very small. Could someone advise me on why these
> correlations occurred and ways I can solve this problem?
>
> Random effects:
> Groups Name Variance Std.Dev. Corr
> CHILDID (Intercept) 333.63787 18.26576
> time.period 7.99838 2.82814 -1.000
> I(time.period^2) 11.16830 3.34190 -0.514 0.514
> I(time.period^3) 0.31869 0.56453 -0.189 0.189 0.938
> School_ID (Intercept) 87.59255 9.35909
> time.period 33.30370 5.77094 0.093
> I(time.period^2) 1.48008 1.21659 -0.171 0.960
> Residual 39.94649 6.32032
> Number of obs: 11960, groups: CHILDID, 2990; School_ID, 996
Have you centred time.period? My simple minded understanding is that
polynomial terms like that will always be highly correlated unless you
orthogonalize them, but the LR based criteria will still guide you
correctly.
--
| David Duffy (MBBS PhD) ,-_|\
| email: davidD at qimr.edu.au ph: INT+61+7+3362-0217 fax: -0101 / *
| Epidemiology Unit, Queensland Institute of Medical Research \_,-._/
| 300 Herston Rd, Brisbane, Queensland 4029, Australia GPG 4D0B994A v
More information about the R-sig-mixed-models
mailing list