[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