[R] anova(lme object)

Liaw, Andy andy_liaw at merck.com
Thu Aug 21 19:16:03 CEST 2003


Yes.  One way is to use anova(res1, type="marginal").  Read the help page
and the book (or any decent linear models book).

Andy

> -----Original Message-----
> From: Mahbub Latif [mailto:ahmlatif at yahoo.com] 
> Sent: Thursday, August 21, 2003 12:51 PM
> To: r-help at stat.math.ethz.ch
> Subject: [R] anova(lme object)
> 
> 
> Hi,
> 
> 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
> formula)
> 
> 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?
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list 
> https://www.stat.math.ethz.ch/mailman/listinfo> /r-help
> 

------------------------------------------------------------------------------
Notice:  This e-mail message, together with any attachments, contains
information of Merck & Co., Inc. (Whitehouse Station, New Jersey, USA), and/or
its affiliates (which may be known outside the United States as Merck Frosst,
Merck Sharp & Dohme or MSD) that may be confidential, proprietary copyrighted
and/or legally privileged, and is intended solely for the use of the
individual or entity named on this message.  If you are not the intended
recipient, and have received this message in error, please immediately return
this by e-mail and then delete it.




More information about the R-help mailing list