[R] anova(lme object)
Prof Brian Ripley
ripley at stats.ox.ac.uk
Thu Aug 21 19:19:32 CEST 2003
On Thu, 21 Aug 2003, Mahbub Latif wrote:
> I use lme to fit models like
> R> res1 <- lme(y~A+B, data=mydata, random=~1|subject)
> R> res2 <- lme(y~B+A, data=mydata, random=~1|subject)
> (only difference between these two models are the
> sequence in which the indep variables are written in
> where y is continuous and A, B, and subject are
> factors. To get ANOVA table I used
> R> anova(res1)
> R> anova(res2)
> and found ANOVA tables corresponding to these two
> models are different. Is there any way I can get
> similar ANOVA tables from lme objects of this type?
Those are *sequential* anova tables, so should be different.
Read the help page ?anova.lme to find the answer to your question.
Hint: argument `type' may help.
Brian D. Ripley, ripley at stats.ox.ac.uk
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UK Fax: +44 1865 272595
More information about the R-help