waste {multcomp}R Documentation

Industrial Waste Data Set

Description

Industrial waste output in a manufactoring plant.

Usage

data("waste")

Format

This data frame contains the following variables

temp

temperature, a factor at three levels: low, medium, high.

envir

environment, a factor at five levels: env1 ... env5.

waste

response variable: waste output in a manufacturing plant.

Details

The data are from an experiment designed to study the effect of temperature (temp) and environment (envir) on waste output in a manufactoring plant. Two replicate measurements were taken at each temperature / environment combination. Data were obtained from Westfall, Tobias, Rom, Wolfinger, and Hochberg (1999), page 177.

References

Westfall P. H., Tobias R. D., Rom D., Wolfinger R. D., Hochberg Y. (1999). Multiple Comparisons and Multiple Tests Using the SAS System. SAS Institute Inc., Cary, NC.

Examples


  ### set up two-way ANOVA with interactions
  amod <- aov(waste ~ temp * envir, data=waste)

  ### comparisons of main effects only
  K <- glht(amod, linfct = mcp(temp = "Tukey"))$linfct
  K
  glht(amod, K)

  ### comparisons of means (by averaging interaction effects)
  low <- grep("low:envi", colnames(K))
  med <- grep("medium:envi", colnames(K))
  K[1, low] <- 1 / (length(low) + 1)
  K[2, med] <- 1 / (length(low) + 1)
  K[3, med] <- 1 / (length(low) + 1)
  K[3, low] <- - 1 / (length(low) + 1)
  K
  confint(glht(amod, K))

  ### same as TukeyHSD
  TukeyHSD(amod, "temp")

  ### set up linear hypotheses for all-pairs of both factors
  wht <- glht(amod, linfct = mcp(temp = "Tukey", envir = "Tukey"))

  ### cf. Westfall et al. (1999, page 181)
  summary(wht, test = adjusted("Shaffer"))


[Package multcomp version 1.4-30 Index]