Here we'll demonstrate the use of anova() to compare two models fit by lme() - note that the models must be nested and the both must be fit by ML rather than REML.