[R] why incomplete lmer output?

[DAVID ARTETA GARCIA]

Maria,

you might get an answer in this thread, from Douglas Bates

"lmer, p-values and all that", dated May 19 2006

HTH

David


Maria Meyen <maria_meyen at yahoo.de> ha escrito:

> 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 at 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.
>
>