include `plot(TukeyHSD(lm_output))`
include
plot(TukeyHSD(lm_output))