[R] why incomplete lmer output?

[Andrew Robinson]

Hi Maria,

the output is not incomplete. See the FAQ here:

http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-are-p_002dvalues-not-displayed-when-using-lmer_0028_0029_003f

and the lengthy discussions about p-values and lmer (try Google).  

Andrew

On Tue, Apr 15, 2008 at 10:51:46AM +0200, Maria Meyen wrote:
> 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.

-- 
Andrew Robinson  
Department of Mathematics and Statistics            Tel: +61-3-8344-6410
University of Melbourne, VIC 3010 Australia         Fax: +61-3-8344-4599
http://www.ms.unimelb.edu.au/~andrewpr
http://blogs.mbs.edu/fishing-in-the-bay/