### set up two-way ANOVA without interactions
amod <- aov(plates ~ block + detergent, data = detergent)
### set up all-pair comparisons
dht <- glht(amod, linfct = mcp(detergent = "Tukey"))
### see Westfall et al. (1999, p. 190)
confint(dht)
### see Westfall et al. (1999, p. 192)
summary(dht, test = univariate())
if (FALSE) {
summary(dht, test = adjusted("Shaffer"))
summary(dht, test = adjusted("Westfall"))
}
Run the code above in your browser using DataLab