[R] why incomplete lmer output?

[Maria Meyen]

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.