[R-sig-ME] Residuals plot for a linear mixed model
Manuel Spínola
mspinola10 at gmail.com
Mon Mar 21 14:14:51 CET 2011
Dear list members,
I got the attached residuals plot for linear mixed model after
logtransformationn of the response variable.
Is the plot telling me that I have a problem with the model fitted? If
logtransformation didn't help, what are thealternatives?
> mode = lmer(log(Swiftness.1 + 0.05) ~ Lure + Sex + Facility.Size +
(1|Subject), REML = F, data = otter)
> summary(mode)
Linear mixed model fit by maximum likelihood
Formula: log(Swiftness.1 + 0.05) ~ Lure + Sex + Facility.Size + (1 |
Subject)
Data: otter
AIC BIC logLik deviance REMLdev
546 574.9 -262 524 511.5
Random effects:
Groups Name Variance Std.Dev.
Subject (Intercept) 0.75907 0.87125
Residual 9.33305 3.05500
Number of obs: 102, groups: Subject, 17
Fixed effects:
Estimate Std. Error t value
(Intercept) 1.1968 1.0832 1.105
Lure[T.2] -0.4375 1.0479 -0.418
Lure[T.3] 0.9928 1.0479 0.947
Lure[T.4] 0.4381 1.0479 0.418
Lure[T.5] -0.3059 1.0479 -0.292
Lure[T.6] -1.1340 1.0479 -1.082
Sex[T.1] -0.2301 0.7411 -0.311
Facility.Size[T.2] 1.9357 0.9550 2.027
Facility.Size[T.3] 1.2889 0.9820 1.312
Correlation of Fixed Effects:
(Intr) L[T.2] L[T.3] L[T.4] L[T.5] L[T.6] S[T.1] F.S[T.2
Lure[T.2] -0.484
Lure[T.3] -0.484 0.500
Lure[T.4] -0.484 0.500 0.500
Lure[T.5] -0.484 0.500 0.500 0.500
Lure[T.6] -0.484 0.500 0.500 0.500 0.500
Sex[T.1] -0.342 0.000 0.000 0.000 0.000 0.000
Fclt.S[T.2] -0.540 0.000 0.000 0.000 0.000 0.000 -0.055
Fclt.S[T.3] -0.544 0.000 0.000 0.000 0.000 0.000 0.000 0.617
predicted<-fitted(mode)
raw.residuals<-residuals(mode)
plot(predicted,raw.residuals)
abline(h=c(-1.96*sd(raw.residuals),0,1.96*sd(raw.residuals)))
Best,
Manuel
--
*Manuel Spínola, Ph.D.*
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
mspinola at una.ac.cr
mspinola10 at gmail.com
Teléfono: (506) 2277-3598
Fax: (506) 2237-7036
Personal website: Lobito de río
<https://sites.google.com/site/lobitoderio/>
Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: residuals.pdf
Type: application/pdf
Size: 9555 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20110321/8f0aea10/attachment-0003.pdf>
More information about the R-sig-mixed-models
mailing list