# Where does the Diet difference lie? TukeyHSD(aov(anov1),ordered=T) # or equivalently TukeyHSD(aov1,ordered=T)