library(TeachingDemos) power.examp(diff = 4, stdev = 1, alpha = .05, n = 10) power.examp(diff = 4, stdev = 1, alpha = .05, n = 10) power.examp(diff = 4, stdev = 1, alpha = .05, n = 10) power.examp(diff = 4, stdev = 1, alpha = .05, n = 10) power.examp(diff = .1, stdev = 1, alpha = .05, n = 10) power.examp(diff = .5, stdev = 1, alpha = .05, n = 10) 1.64 - sqrt(20)/30 1-pnorm(1.49) 1.64 - 10*sqrt(20)/30 1.64 - 10*sqrt(20)/30 1-pnorm(.15) 1.64 - 20*sqrt(20)/30 1-pnorm(-1.34) 1.64*30/sqrt(20) 1.64*30/sqrt(200) qnorm(.2) 30*(1.64 + .84)/sqrt(20) 30*(1.64 + .84)/sqrt(30) 30*(1.64 + .84)/sqrt(100) 30*(1.64 + .84)/sqrt(200) 30*(1.64 + .84)/sqrt(300) 30*(1.64 + .84)/sqrt(10000) 30*(1.64 + .84)/sqrt(100) 30^2*(-.84-1.64)^2/2^2 30^2*(-.84-1.64)^2/1^2 30^2*(-.84-1.64)^2/.5^2 30^2*(-.84-1.64)^2/10^2 30^2*(-.84-1.64)^2/10^2 30^2*(-.84-1.64)^2/8^2 30^2*(-.84-1.64)^2/.1^2 30^2*(-.84-1.64)^2/10^2 install.packages("TeachingDemos") library(TeachingDemos) power.examp() power.examp(diff=.1) power.examp(diff=.1,n =1)