Need a way to compare competing models other than AIC/BIC. Extending the anova() function is one way to do this.