Christophe Pallier
2005-Oct-28 10:39 UTC
[R] aov() and lme() (repeated measures and ANOVA)
Hello,> Subject: [R] aov() and lme() > From: Jan Wiener <jan.wiener at 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