Displaying 1 result from an estimated 1 matches for "t3_c2".
Did you mean:
t1_c2
2009 Jan 12
0
Two-way repeated measures anova with lme
...s',1:9,sep=''),each=4)
),y=rnorm(108))
summary(aov(y ~ (Time*Cond) + Error(Subj/(Time*Cond)), data = df))
library(nlme)
anova(lme(y~Time*Cond,random=~1|Subj,data=df))
df2<-cast(df,Subj~Time+Cond)[-1]
mod<-lm(cbind(t1_c1,t1_c2,t1_c3,t1_c4,t2_c1,t2_c2,t2_c3,t2_c4,t3_c1,t3_c2,t3
_c3,t3_c4)~1,df2)
idata=data.frame(Time=gl(3,4),Cond=gl(4,1))
summary(Anova(mod,idata=idata,idesign=~Time*Cond),multivariate=F)
I would highly appreciate if anyone could explain why the results are
slightly different. I probably spent hours to find an answer on the
internet. I would lik...