################################################################## #Cohen,J. (1988). #Statistical power analysis for the behavioral sciences (2nd ed.). #Hillsdale, NJ:Erlbaum ################################################################## source("http://blue.zero.jp/yokumura/R/power/Power.R") #----------------------------------------# # Chapter2: The t test for means (pp.19--74) #----------------------------------------# ##ex 2.1 power.t(sig.level=.05, delta=0.5, n1=30, n2=30) power.t(sig.level=.05, delta=0.70, n1=30, n2=30) power.t(sig.level=.05, delta=0.71, n1=30, n2=30) power.t(sig.level=.05*2, delta=0.50, n1=30, n2=30) power.t(sig.level=.05*2, delta=0.71, n1=30, n2=30) power.t(sig.level=.05, delta=0.8, n1=30, n2=30) power.t(sig.level=.05, delta=0.2, n1=30, n2=30) ##ex 2.2 power.t(sig.level=.01, delta=0.2, n1=500, n2=500) power.t(sig.level=.05, delta=0.2, n1=500, n2=500) ##ex 2.3 power.t(sig.level=.05*2, delta=0.6, n1=72, n2=72) power.t(sig.level=.05*2, delta=0.6, n1=90, n2=60) power.t(sig.level=.05, delta=0.6, n1=72, n2=72) power.t(sig.level=.05, delta=0.6, n1=90, n2=60) ##ex 2.4 power.t(sig.level=.01, delta=0.5, n1=100, n2=100) power.t(sig.level=.05, delta=0.5, n1=100, n2=100) power.t(sig.level=.01*2, delta=0.5, n1=100, n2=100) ##ex 2.5 power.t(sig.level=.10, delta=0.20, n1=60, n2=60) power.t(sig.level=.10, delta=0.30, n1=60, n2=60) power.t(sig.level=.10, delta=0.28, n1=60, n2=60) ##ex 2.6 power.t(sig.level=.05, delta=0.60, n1=50, n2=50) power.t(sig.level=.05, delta=0.40, n1=50, n2=50) ##ex 2.7 power.t(sig.level=.01*2, delta=0.70, n1=60, n2=60) power.t(sig.level=.01*2, delta=0.80, n1=60, n2=60) power.t(sig.level=.01*2, delta=0.74, n1=60, n2=60) #----------------------------------------# # Chapter3: Significance of a product moment r_s (pp.75--107) #----------------------------------------# ##ex3.1 power.cor2(sig.level=0.05, delta=0.3, n=50) power.cor2(sig.level=0.05*2, delta=0.3, n=50) #directional ##ex3.2 power.cor2(sig.level=0.01, delta=0.1, n=500) power.cor2(sig.level=0.05, delta=0.1, n=500) power.cor2(sig.level=0.10, delta=0.1, n=500) power.cor2(sig.level=0.01, delta=0.1, n=1000) power.cor2(sig.level=0.05, delta=0.1, n=1000) power.cor2(sig.level=0.01, delta=0.2, n=500) power.cor2(sig.level=0.05, delta=0.2, n=500) power.cor2(sig.level=0.10, delta=0.2, n=500) power.cor2(sig.level=0.01, delta=0.2, n=1000) power.cor2(sig.level=0.05, delta=0.2, n=1000) power.cor2(sig.level=0.10, delta=0.2, n=1000) ##ex3.3 power.cor2(sig.level=0.01*2, delta=0.4, n=120) #----------------------------------------# # Chapter6: Differences between proportions (pp.179--213) #----------------------------------------# ##ex 6.1 power.prop(h=.3, n=c(80,80), sig.level=0.05, one.sided=T) ##ex 6.2 power.prop(h=.4, n=c(100,100), sig.level=0.05, one.sided=F) power.prop(h=.4, n=c(120,120), sig.level=0.05, one.sided=F) power.prop(h=.4, n=c(140,140), sig.level=0.05, one.sided=F) power.prop(h=.4, n=c(125,125), sig.level=0.05, one.sided=F) power.prop(h=.3, n=c(120,120), sig.level=0.05, one.sided=F) power.prop(h=.3, n=c(140,140), sig.level=0.05, one.sided=F) power.prop(h=.3, n=c(125,125), sig.level=0.05, one.sided=F) ##ex 6.3 power.prop(h=.3, n=c(120,120), sig.level=0.05, one.sided=T) power.prop(h=.3, n=c(80,245), sig.level=0.05, one.sided=T) ##ex 6.4 power.prop(h=.5, n=c(85,450), sig.level=0.01, one.sided=T) power.prop(h=.5, n=c(85,450), sig.level=0.01, one.sided=F) ##ex 6.5 power.prop(h=.20, n=60, sig.level=0.05, one.sided=F) #----------------------------------------# # Chapter7: Chi-square tests for goodness of fit and contingency tables (pp.215--271) #----------------------------------------# ##ex 7.1 power.chi(sig.level=.05,u=3,n=100,delta=.2) power.chi(sig.level=.05,u=3,n=100,delta=.3) power.chi(sig.level=.05,u=3,n=100,delta=.289) ##ex 7.2 power.chi(sig.level=.1,u=6,n=200,delta=.2) ##ex 7.3 power.chi(sig.level=.01,u=2,n=140,delta=.3) power.chi(sig.level=.01,u=2,n=140,delta=.4) power.chi(sig.level=.01,u=2,n=140,delta=.346) ##ex 7.4 power.chi(sig.level=.05,u=4,n=80,delta=.2) power.chi(sig.level=.05,u=4,n=80,delta=.3) power.chi(sig.level=.05,u=4,n=80,delta=.283) ##ex 7.5 power.chi(sig.level=.01,u=20,n=400,delta=.1) power.chi(sig.level=.10,u=20,n=400,delta=.1) power.chi(sig.level=.01,u=20,n=400,delta=.3) power.chi(sig.level=.05,u=20,n=400,delta=.3) #----------------------------------------# # Chapter8: The analysis variance (pp.273--406) #----------------------------------------# ##ex 8.1 power.f(sig.level=.05,u=3,delta=.20, n=20) power.f(sig.level=.05,u=3,delta=.25, n=20) power.f(sig.level=.05,u=3,delta=.28, n=20) power.f(sig.level=.10,u=3,delta=.20, n=20) power.f(sig.level=.10,u=3,delta=.25, n=20) power.f(sig.level=.10,u=3,delta=.28, n=20) ##ex 8.2 power.f(sig.level=.05,u=2,delta=.20, n=200) power.f(sig.level=.05,u=2,delta=.25, n=200) power.f(sig.level=.05,u=2,delta=.23, n=200) power.f(sig.level=.05,u=3,delta=.33, n=200) power.f(sig.level=.01,u=3,delta=.23, n=140) power.f(sig.level=.01,u=3,delta=.33, n=140) power.f(sig.level=.01,u=2,delta=.15, n=200) power.f(sig.level=.01,u=3,delta=.15, n=150) ##ex 8.3 power.f(sig.level=.05,u=5,delta=.15, n=50) power.f(sig.level=.10,u=5,delta=.15, n=50) power.f(sig.level=.05,u=2,delta=.15, n=100) power.f(sig.level=.05,u=1,delta=.15, n=150) ##ex 8.4 power.f(sig.level=.05,u=10,delta=.25, n=27) power.f(sig.level=.05,u=12,delta=.25, n=27) power.f(sig.level=.05,u=11,delta=.25, n=27) ##ex 8.5 power.f2(df1=2,df2=60,sig.level=0.05,delta=.25) power.f2(df1=3,df2=60,sig.level=0.05,delta=.40) ##ex 8.6 power.f2(df1=1,df2=96,sig.level=0.05,delta=.10) power.f2(df1=2,df2=96,sig.level=0.05,delta=.25) power.f2(df1=3,df2=96,sig.level=0.05,delta=.40) power.f2(df1=1,df2=120,sig.level=0.05,delta=.10) power.f2(df1=2,df2=120,sig.level=0.05,delta=.25) power.f2(df1=3,df2=120,sig.level=0.05,delta=.40) ##ex 8.7 power.f2(df1=1,df2=96,sig.level=0.05,delta=.10) #S power.f2(df1=2,df2=96,sig.level=0.05,delta=.25) #A power.f2(df1=3,df2=96,sig.level=0.05,delta=.40) #C power.f2(df1=2,df2=96,sig.level=0.05,delta=.25) #S * A power.f2(df1=3,df2=96,sig.level=0.05,delta=.25) #S * C power.f2(df1=6,df2=96,sig.level=0.05,delta=.25) #A * C power.f2(df1=6,df2=96,sig.level=0.05,delta=.25) #S * A * C ##ex 8.8 power.f2(df1=2,df2=120,sig.level=0.01,delta=.1) #A power.f2(df1=2,df2=120,sig.level=0.05,delta=.1) #A power.f2(df1=2,df2=120,sig.level=0.10,delta=.1) #A power.f2(df1=2,df2=120,sig.level=0.01,delta=.25) #A power.f2(df1=2,df2=120,sig.level=0.05,delta=.25) #A power.f2(df1=2,df2=120,sig.level=0.10,delta=.25) #A power.f2(df1=2,df2=120,sig.level=0.01,delta=.40) #A power.f2(df1=2,df2=120,sig.level=0.05,delta=.40) #A power.f2(df1=2,df2=120,sig.level=0.10,delta=.40) #A power.f2(df1=3,df2=120,sig.level=0.01,delta=.1) #B power.f2(df1=3,df2=120,sig.level=0.05,delta=.1) #B power.f2(df1=3,df2=120,sig.level=0.10,delta=.1) #B power.f2(df1=3,df2=120,sig.level=0.01,delta=.25) #B power.f2(df1=3,df2=120,sig.level=0.05,delta=.25) #B power.f2(df1=3,df2=120,sig.level=0.10,delta=.25) #B power.f2(df1=3,df2=120,sig.level=0.01,delta=.40) #B power.f2(df1=3,df2=120,sig.level=0.05,delta=.40) #B power.f2(df1=3,df2=120,sig.level=0.10,delta=.40) #B power.f2(df1=4,df2=120,sig.level=0.01,delta=.1) #C power.f2(df1=4,df2=120,sig.level=0.05,delta=.1) #C power.f2(df1=4,df2=120,sig.level=0.10,delta=.1) #C power.f2(df1=4,df2=120,sig.level=0.01,delta=.25) #C power.f2(df1=4,df2=120,sig.level=0.05,delta=.25) #C power.f2(df1=4,df2=120,sig.level=0.10,delta=.25) #C power.f2(df1=4,df2=120,sig.level=0.01,delta=.40) #C power.f2(df1=4,df2=120,sig.level=0.05,delta=.40) #C power.f2(df1=4,df2=120,sig.level=0.10,delta=.40) #C power.f2(df1=6,df2=120,sig.level=0.01,delta=.1) #A * B power.f2(df1=6,df2=120,sig.level=0.05,delta=.1) #A * B power.f2(df1=6,df2=120,sig.level=0.10,delta=.1) #A * B power.f2(df1=6,df2=120,sig.level=0.01,delta=.25) #A * B power.f2(df1=6,df2=120,sig.level=0.05,delta=.25) #A * B power.f2(df1=6,df2=120,sig.level=0.10,delta=.25) #A * B power.f2(df1=6,df2=120,sig.level=0.01,delta=.40) #A * B power.f2(df1=6,df2=120,sig.level=0.05,delta=.40) #A * B power.f2(df1=6,df2=120,sig.level=0.10,delta=.40) #A * B power.f2(df1=8,df2=120,sig.level=0.01,delta=.1) #A * C power.f2(df1=8,df2=120,sig.level=0.05,delta=.1) #A * C power.f2(df1=8,df2=120,sig.level=0.10,delta=.1) #A * C power.f2(df1=8,df2=120,sig.level=0.01,delta=.25) #A * C power.f2(df1=8,df2=120,sig.level=0.05,delta=.25) #A * C power.f2(df1=8,df2=120,sig.level=0.10,delta=.25) #A * C power.f2(df1=8,df2=120,sig.level=0.01,delta=.40) #A * C power.f2(df1=8,df2=120,sig.level=0.05,delta=.40) #A * C power.f2(df1=8,df2=120,sig.level=0.10,delta=.40) #A * C power.f2(df1=12,df2=120,sig.level=0.01,delta=.1) #B * C power.f2(df1=12,df2=120,sig.level=0.05,delta=.1) #B * C power.f2(df1=12,df2=120,sig.level=0.10,delta=.1) #B * C power.f2(df1=12,df2=120,sig.level=0.01,delta=.25) #B * C power.f2(df1=12,df2=120,sig.level=0.05,delta=.25) #B * C power.f2(df1=12,df2=120,sig.level=0.10,delta=.25) #B * C power.f2(df1=12,df2=120,sig.level=0.01,delta=.40) #B * C power.f2(df1=12,df2=120,sig.level=0.05,delta=.40) #B * C power.f2(df1=12,df2=120,sig.level=0.10,delta=.40) #B * C power.f2(df1=24,df2=120,sig.level=0.01,delta=.1) #A * B * C power.f2(df1=24,df2=120,sig.level=0.05,delta=.1) #A * B * C power.f2(df1=24,df2=120,sig.level=0.10,delta=.1) #A * B * C power.f2(df1=24,df2=120,sig.level=0.01,delta=.25) #A * B * C power.f2(df1=24,df2=120,sig.level=0.05,delta=.25) #A * B * C power.f2(df1=24,df2=120,sig.level=0.10,delta=.25) #A * B * C power.f2(df1=24,df2=120,sig.level=0.01,delta=.40) #A * B * C power.f2(df1=24,df2=120,sig.level=0.05,delta=.40) #A * B * C power.f2(df1=24,df2=120,sig.level=0.10,delta=.40) #A * B * C ##ex 8.8 power.f2(df1=3,df2=88,sig.level=0.05,delta=.293) #R * C power.f2(df1=1,df2=88,sig.level=0.05,delta=.375) #R power.f2(df1=3,df2=88,sig.level=0.05,delta=.523) #C power.f2(df1=3,df2=128,sig.level=0.05,delta=.293) #R * C power.f2(df1=3,df2=88,sig.level=0.10,delta=.293) #R * C #----------------------------------------# # Chapter9: Multiple regression and correlation analysis (pp.407-465) #----------------------------------------# ##ex 9.1 power.multi(n=95,n.ind=5,sig.level=.05,delta=.1111) ##ex 9.2 power.multi(n=90,n.ind=3,sig.level=.01,delta=.15) ##ex 9.3 power.multi(n=90,n.ind=3,sig.level=.01,delta=.15) ##ex 9.4 power.f(sig.level=.05,u=11,delta=.25, n=27) power.multi(n=326,n.ind=11,sig.level=.05,delta=.0625) ##ex 9.5 power.f(sig.level=.05,u=3,delta=.28, n=20) power.multi(n=80,n.ind=3,sig.level=.05,delta=.0784) ##ex 9.6 power.cor2(sig.level=0.05, delta=0.3, n=50) power.multi(n=50,n.ind=1,sig.level=.05,delta=.0989) ##ex 9.7 power.t(sig.level=.05, delta=0.6, n1=90, n2=60) power.multi(n=150,n.ind=1,sig.level=.05,delta=.0864) ##ex 9.8 power.multi(n=90,n.ind=2,sig.level=.01,delta=.16) ##ex 9.9 power.multi(n=80,n.ind=3,sig.level=.05,delta=.0784) power.multi(n=80,n.ind=3,sig.level=.05,delta=.1228) ##ex 9.10 power.multi(n=90,n.ind=4,sig.level=.01,delta=.0714) power.multi(n=180,n.ind=4,sig.level=.01,delta=.0714) ##ex 9.11 power.multi(n=148, n.ind=2,sig.level=.01,delta=.1111) ##ex 9.12 power.multi(n=95, n.ind=2,sig.level=.05, delta=.0444) power.multi(n=95, n.ind=1,sig.level=.05, delta=.0333) ##ex 9.13 power.multi(n=200, n.ind=3,sig.level=.01, delta=.075) power.multi(n=200, n.ind=5,sig.level=.01, delta=.125) power.multi(n=200, n.ind=4,sig.level=.01, delta=.05) ##ex 9.14 power.multi(n=120, n.ind=6,sig.level=.05, delta=.0625) power.multi(n=120, n.ind=3,sig.level=.05, delta=.16) ##ex 9.15 power.multi(n=90, n.ind=1, sig.level=.01, delta=.092) power.multi(n=90, n.ind=1, sig.level=.05, delta=.092) power.multi(n=90, n.ind=2, sig.level=.05, delta=.0575) power.multi(n=90, n.ind=1, sig.level=.05, delta=.0575)