Load required packages library(emmeans) library(multcomp) library(ggplot2) ## Load the data set data(waste) ## Sanity check ---- ?waste # brings up the help file, which obviously only works with...

Load required packages library(emmeans) library(multcomp) library(ggplot2) ## Load the data set data(waste) ## Sanity check ---- ?waste # brings up the help file, which obviously only works with built-in data sets str(waste) summary(waste) head(waste) ## Graphical data exploration ---- ##...using ggplot2 ggplot(data = waste, aes(x = temp, y = waste)) + geom_boxplot() ggplot(data = waste, aes(x = envir, y = waste)) + geom_boxplot() ggplot(data = waste, aes(x = temp, y = waste)) + geom_boxplot() + facet_wrap(facets = ~envir) ## ...same using traditional graphics boxplot(waste ~ temp, data = waste) boxplot(waste ~ envir, data = waste) boxplot(waste ~ temp * envir, data = waste) # Not as neat as the ggplot version ## Two-way ANOVA without interaction m1