data(depression)
head(depression)
# Two-way mixed ANOVA: treatment (between) x time (within)
depression_long <- reshape(depression, varying = c("t0", "t1", "t2", "t3"),
v.names = "score", timevar = "time",
direction = "long")
summary(aov(score ~ treatment * factor(time) + Error(factor(id)/factor(time)),
data = depression_long))
Run the code above in your browser using DataLab