Re: [R] Is using glht with "Tukey" for lme post-hoc comparisons an appropriate substitute to TukeyHSD?

riate interactions found -- default contrast might be inappropriate ## ## From ?glht, read about the interaction_average argument. ## With the glht argument interaction_average=TRUE, they give the same answer ## With your example, summary(glht(model3, linfct=mcp(Treatment="Tukey", interaction_average=

Re: [R] Is using glht with "Tukey" for lme post-hoc comparisons an appropriate substitute to TukeyHSD?

. Both TukeyHSD and glht give essesntially identical confidence intervals for the example in ?glht. What aren't you satisfied with? amod <- aov(breaks ~ tension, data = warpbreaks) confint(glht(amod, linfct = mcp(tension = "Tukey"))) TukeyHSD(amod) On Mon, Jan 2, 2012 at 6:19 PM,

[R] Is using glht with "Tukey" for lme post-hoc comparisons an appropriate substitute to TukeyHSD?

Hello, I am trying to determine the most appropriate way to run post-hoc comparisons on my lme model. I had originally planned to use Tukey HSD method as I am interested in all possible comparisons between my treatment levels. TukeyHSD, however, does not work with lme. The only other c