Re: [R] aov() and lme() (repeated measures and ANOVA)

From: Christophe Pallier <pallier_at_lscp.ehess.fr>
Date: Fri 28 Oct 2005 - 20:39:27 EST

Hello,

> Subject: [R] aov() and lme()
> From: Jan Wiener <jan.wiener@tuebingen.mpg.de>
> Date: Thu, 27 Oct 2005 13:14:59 +0200
>
> Sorry for reposting, but even after extensive search I still did not
> find any answers.
>
> using:
> summary(aov(pointErrorAbs~noOfSegments*turnAngle+Error(subj/(noOfSegments+turnAngle)),
> data=anovaAllData ))
>
> with subj being a random factor and noOfSegments and turnAngle being
> fixed factors, I get the following results:
>
> [...]
>
> No I trying to fit the same data with lme and using the following call:
>
> anova(lme(fixed=pointErrorAbs~noOfSegments*turnAngle, random=~1|subj,
> data=anovaAllData))
>
> Unfortunately the results are 'really' different from the aov()
> procedure (I guess I have the call wrong):
>

Maybe the following post concerning repeated measures and ANOVA can help:

(from http://tolstoy.newcastle.edu.au/~rking/R/help/03b/7663.html)

 > $ anova(lme(DV ~ GROUP*TRIAL,random= ~1|SUB, correlation=corCompSymm() ))
 >
 >(TRIAL and GROUP are fixed factors)
 >
 >Actually, you do not need lme for to run a repeated measure anova.
 >You could use the aov function:
 >
 > $ summary(aov(DV~GROUP*TRIAL+Error(SUB/TRIAL)))
 >
 > ... yields the same results (when data are balanced)


Christophe Pallier
www.pallier.org



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Fri Oct 28 20:43:56 2005

This archive was generated by hypermail 2.1.8 : Fri 28 Oct 2005 - 21:15:34 EST