data(selfesteem)
head(selfesteem)
# One-way repeated measures ANOVA of the self-esteem score over time
selfesteem_long <- reshape(selfesteem, varying = c("t1", "t2", "t3"),
v.names = "score", timevar = "time",
direction = "long")
summary(aov(score ~ factor(time) + Error(factor(id)/factor(time)),
data = selfesteem_long))
Run the code above in your browser using DataLab