#From the work of John Verzani
#One-way ANOVA

rater1 <- c(4,3,4,5,2,3,4,5)
rater2 <- c(4,4,5,5,4,5,4,4)
rater3 <- c(3,4,2,4,5,5,4,4)

scores <- data.frame(rater1,rater2,rater3)
boxplot(scores)

scores = stack(scores)
names(scores)

oneway.test(values ~ ind, data=scores, var.equal=T)