### 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