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.9196 2 57 0.00461 0.172 0.1429
Run the code above in your browser using DataLab