glht & lmer (was: Post hoc tests with lme)
This also makes me wonder, is it possible to use glht to incorporate random effects as well? E.g., if a fixed effect differs by, say, block, can one use glht to test whether the effect differs within each block?
On Apr 16, 2008, at 9:16 PM, Michael Kubovy wrote:
I wonder if there's a problem with glht and the current version of lmer From http://epub.ub.uni-muenchen.de/2120/1/tr019.pdf : require(multcomp) require(lme4) data(trees513) mmod <- lmer(damage ~ species - 1 + (1 | plot) + (1 | lattice), data = trees513, family = binomial()) K <- diag(length(fixef(mmod))) (ci <- confint(glht(mmod, linfct = K))) Simultaneous Confidence Intervals Fit: Error in x$model$call : $ operator not defined for this S4 class
_____________________________
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400 Charlottesville, VA 22904-4400
Parcels: Room 102 Gilmer Hall
McCormick Road Charlottesville, VA 22903
Office: B011 +1-434-982-4729
Lab: B019 +1-434-982-4751
Fax: +1-434-982-4766
WWW: http://www.people.virginia.edu/~mk9y/
[[alternative HTML version deleted]]
_______________________________________________
R-sig-mixed-models at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models