R標準の関数で、検定力を求めます

 n=10, 平均値差delta=0.5, 標準偏差=1,  つまり効果量 d=0.5, 有意水準 .05 とした場合

power.t.test(n=10,delta=0.5,sd=1,sig.level=0.05)

​

> power.t.test(n=10,delta=0.5,sd=1,sig.level=0.05)

     Two-sample t test power calculation 

              n = 10
          delta = 0.5
             sd = 1
      sig.level = 0.05
          power = 0.1838375
    alternative = two.sided

NOTE: n is number in *each* group

> 

​

　検定力　power  1-β=0.1838375

​

​

​

検定力からサンプルサイズを決めます

　n=NULL として、power=0.8 とした場合、

power.t.test(n=NULL,delta=0.5,sd=1,sig.level=0.05,power=0.8)

​

> power.t.test(n=NULL,delta=0.5,sd=1,sig.level=0.05,power=0.8)

     Two-sample t test power calculation 

              n = 63.76576
          delta = 0.5
             sd = 1
      sig.level = 0.05
          power = 0.8
    alternative = two.sided

NOTE: n is number in *each* group
 

サンプルサイズ　n=64  （各群）

​

​

​

参考）　山田剛史・杉澤武俊・村井潤一郎 (2008)　Rによるやさしい統計学　　オーム社　pp.364-367