# Plot for diagnostics
opar <- par(mfrow=c(2,2))
plot(anov1)
par(opar)

# Test assumption of normality of residuals
shapiro.test(resid(anov1))

# Test assumption of homogeneity of variance
bartlett.test(logMaxAbund ~ Diet, data=bird)