[R-sig-ME] anova

Cougar Cougar_711 at msn.com
Wed Jun 6 22:15:42 CEST 2007


I wonder if it is possible to nest a generalized least squares model within
an lmer model for the purpose of testing the hypothesis that sigma.sq_beta_0
= 0.

Example: R 2.5.0, WinXP

>library(lme4)
>data(sleepstudy)
>library(nlme)
>m0 <- gls(Reaction ~ Days, sleepstudy)
>m1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)
>anova(m0,m1)

I do not get any result from the anova.  I assume it is because the model
objects are of different classes.  Is that correct?

Respectfully,

Frank Lawrence




More information about the R-sig-mixed-models mailing list