m <- c(12.4, 8.6, 10.5)
sd <- c(3.84, 3.12, 3.48)
n <- c(20, 20, 20)
test.anova.bs(m, sd, n)
# Should return:
# F dfA dfE p Eta-squared adj Eta-squared
# 5.919585 2 57 0.004614428 0.1719831 0.1429298
Run the code above in your browser using DataLab