### set up one-way ANOVA
amod <- aov(minutes ~ blanket, data = recovery)
### set up multiple comparisons: one-sided Dunnett contrasts
rht <- glht(amod, linfct = mcp(blanket = "Dunnett"),
alternative = "less")
### cf. Westfall et al. (1999, p. 80)
confint(rht, level = 0.9)
### the same
rht <- glht(amod, linfct = mcp(blanket = c("b1 - b0 >= 0",
"b2 - b0 >= 0",
"b3 - b0 >= 0")))
confint(rht, level = 0.9)
Run the code above in your browser using DataLab