[R] why incomplete lmer output?

From: Maria Meyen <maria_meyen_at_yahoo.de>
Date: Tue, 15 Apr 2008 10:51:46 +0200 (CEST)


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 doesnt 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. Im 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. Received on Tue 15 Apr 2008 - 08:55:30 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 15 Apr 2008 - 09:30:30 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive