# Using aov() aov1 <- aov(logMaxAbund ~ Diet, data=bird) summary(aov1) # Using lm() anov1 <- lm(logMaxAbund ~ Diet, data=bird) anova(anov1)