# NOT RUN {
## These should all give the same results
if (require(mosaicData)) {
model <- lm(age ~ substance, data=HELPrct)
TukeyHSD(model)
TukeyHSD( age ~ substance, data=HELPrct)
TukeyHSD(aov(age ~ substance, data=HELPrct))
}
# }
Run the code above in your browser using DataCamp Workspace