anov3 <- lm(logMaxAbund ~ Diet + Aquatic, data=bird) anova(anov3, anov2)