Dear all,
I have problems with the lmer output (see
below)(2.6.2, WindowsXP), it is incomplete (no DFs, no
p-values)and also generating an anova-table doesn´t
work out. What is the fault? Does anybody know?
This example is about plant height within different
treatments (each with 5 replicates) measured several
times. I´m interested in the treatment effect.
> pap<-groupedData(heig~meas|box,outer=~treat,pap)
> mod<-lmer(heig~treat+(meas|box),pap)
> summary(mod)
Linear mixed-effects model fit by REML
Formula: heig ~ treat + (meas | box)
Data: pap
AIC BIC logLik MLdeviance REMLdeviance
134.1 156.6 -58.05 97.14 116.1
Random effects:
Groups Name Variance Std.Dev. Corr
box (Intercept) 0.621046 0.78806
meas 0.790570 0.88914 -1.000
Residual 0.037141 0.19272
number of obs: 90, groups: box, 30
Fixed effects:
Estimate Std. Error t value
(Intercept) 1.79473 0.08353 21.487
treatB1 0.23334 0.11813 1.975
treatB2 0.62677 0.11813 5.306
treatD1 0.07911 0.11813 0.670
treatD2 0.43942 0.11813 3.720
treatU -0.11933 0.11813 -1.010
Correlation of Fixed Effects:
(Intr) tretB1 tretB2 tretD1 tretD2
treatB1 -0.707
treatB2 -0.707 0.500
treatD1 -0.707 0.500 0.500
treatD2 -0.707 0.500 0.500 0.500
treatU -0.707 0.500 0.500 0.500 0.500
> anova(mod)
Analysis of Variance Table
Df Sum Sq Mean Sq
treat 5 2.11106 0.42221
Thanks a lot in advance
Maria Meyen
Department of Conservation Biology
Faculty of Biology
Philipps-University of Marburg
Karl-von-Frisch-Str.8
35043 Marburg
__________________________________________________________
Der Lieblings-Mailbox der Welt.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
