# NOT RUN {
require(multcomp)
warp.lm <- lm(breaks ~ wool*tension, data = warpbreaks)
# Using 'emmGrid'
summary(glht(warp.lm, emmGrid(pairwise ~ tension | wool)))
# Same, but using an existing 'emmeans' result
warp.emmGrid <- emmeans(warp.lm, ~ tension | wool)
summary(as.glht(pairs(warp.emmGrid)))
# Same contrasts, but treat as one family
summary(as.glht(pairs(warp.emmGrid), by = NULL))
# }
Run the code above in your browser using DataLab