[R-sig-ME] lmer vs glmmPQL

Daniel Ezra Johnson danielezrajohnson at gmail.com
Wed Jun 24 15:47:04 CEST 2009


Regarding this thread, what about the method of fitting nested models
and using anova() to estimate a p-value.

For example:

mod2 = lmer(y ~ genotype + (1|block), family = binomial, data)
mod0 = lmer(y ~ (1|block), family = binomial, data)

anova(mod0,mod2)

How does that p-value (which is one value for the whole term
"genotype") relate to the individual coefficient p-values derived from
the Z-scores inside summary(mod2)?

Thanks,
Dan




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