m <- c(56.1, 51.2, 60.3, 68.2, 48.9, 70.5)
sd <- c(9.45, 8.79, 9.71, 8.90, 8.31, 9.75)
ci.random.anova(.05, m, sd, 20)
# Should return:
# Estimate LL UL
# Grand mean 59.200000 49.9363896 68.4636104
# Within SD: 9.166782 8.0509046 10.4373219
# Between SD: 8.585948 8.3239359 8.8562078
# Omega-squared: 0.467317 0.2284142 0.8480383
Run the code above in your browser using DataLab