anova-method function

Compare nested models with likelihood-based statistics