# NOT RUN {
library(ggplot2)
library(broom)
ggplot(salinity, aes(x = salinity_ppt)) +
geom_dotplot() +
facet_wrap(~site_number, ncol = 1)
tidy(aov(salinity_ppt ~ site_number, data = salinity))
# }
Run the code above in your browser using DataLab