[R-sig-ME] random effect NaN corr interpretation?
Mollie Brooks
mbrooks at ufl.edu
Tue Oct 18 20:39:08 CEST 2011
Hello,
I fit an LMM and need help understanding the results. They may indicate that I'm overfitting, but I have many observations (37553). I think that different Tags (random) may respond differently to the fixed effect I(SHTS*HT) so I fit a random intercept and slope for each Tag. In the model with only a random intercept, the Variance is estimated to be 0 by lmer. In this more complex model with both a random intercept and slope, the variance of the intercept is still 0 and the variance of the slope is 5.7392e-03. The part I'm not sure about is that their correlation is NaN. How do I interpret this?
thanks,
Mollie
> sm1.3=lmer(grow1~I(SHTS*HT)+SIZE+(I(SHTS*HT)|Tag_PLOT)+(1|year)+(1|PLOT), data=sHlong)
> summary(sm1.3)
Linear mixed model fit by REML
Formula: grow1 ~ I(SHTS * HT) + SIZE + (I(SHTS * HT) | Tag_PLOT) + (1 | year) + (1 | PLOT)
Data: sHlong
AIC BIC logLik deviance REMLdev
448411 448488 -224197 448394 448393
Random effects:
Groups Name Variance Std.Dev. Corr
Tag_PLOT (Intercept) 0.0000e+00 0.000000
I(SHTS * HT) 5.7392e-03 0.075758 NaN
year (Intercept) 2.6760e+02 16.358366
PLOT (Intercept) 3.1679e+01 5.628452
Residual 8.7827e+03 93.716081
Number of obs: 37553, groups: Tag_PLOT, 5672; year, 11; PLOT, 10
Fixed effects:
Estimate Std. Error t value
(Intercept) 28.823909 5.517252 5.22
I(SHTS * HT) -0.211216 0.004152 -50.88
SIZE1-ha -1.832286 3.866764 -0.47
Correlation of Fixed Effects:
(Intr) I(SH*H
I(SHTS*HT) -0.075
SIZE1-ha -0.279 0.012
> attr(VarCorr(sm1.3)$Tag_PLOT,"correlation")
(Intercept) I(SHTS * HT)
(Intercept) 1 NaN
I(SHTS * HT) NaN 1
More information about the R-sig-mixed-models
mailing list