anov5 <- lm(logMaxAbund ~ Diet + Aquatic, data=bird) anova(anov5, anov4)