set.seed(113)
power.dunnett.test(n = 9, groups = 5, delta = 30,
within.var = 333.7)
## compare with t-test, bonferroni corrected
power.t.test(n = 9, delta = 30, sd = sqrt(333.7),
sig.level = 0.05 / 4)
if (FALSE) {
## asymptotic Monte-Carlo power analysis
set.seed(113)
powerMCTests(mu = c(rep(0,4), 30), n = 9,
parms = list(mean = 0, sd = sqrt(333.7)),
test = "dunnettTest", alternative = "two.sided")
}
Run the code above in your browser using DataLab