pwr.r

tchsieh — Apr 25, 2013, 10:58 PM

library(pwr)
samplesize <- c(2:200)
power <- pwr.t.test(n = samplesize, d = 0.5, sig.level = 0.05, power = NULL,
                    type = c("two.sample"))
names(power)

[1] "n"           "d"           "sig.level"   "power"       "alternative"
[6] "note"        "method"     

sample.size <- power$n
power.size <- power$power
plot(sample.size, power.size, type="l")