data(selfesteem2)
head(selfesteem2)
# Two-way repeated measures ANOVA (treatment x time), both within-subjects
selfesteem2_long <- reshape(selfesteem2, varying = c("t1", "t2", "t3"),
v.names = "score", timevar = "time",
idvar = c("id", "treatment"), direction = "long")
summary(aov(score ~ treatment * factor(time) +
Error(factor(id)/(treatment * factor(time))),
data = selfesteem2_long))
Run the code above in your browser using DataLab