[R-sig-ME] glht & lmer (was: Post hoc tests with lme)
Jarrett Byrnes
jebyrnes at ucdavis.edu
Thu Apr 17 06:24:07 CEST 2008
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
More information about the R-sig-mixed-models
mailing list