> plot(count ~ spray, InsectSprays) > res.aov = aov(count ~ spray, InsectSprays) > summary(res.aov) > plot(TukeyHSD(res.aov))