[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

>m0 <- gls(Reaction ~ Days, sleepstudy)
>m1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)

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


Frank Lawrence

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