「統計学関連なんでもあり」の過去ログ--- 042

No.07910　サンプルサイズの検出力分析　　【タロー】　2008/10/10(Fri) 22:25

論文の中で検出力を算出する必要が出てきました。トータル30人の被験者をランダムに10人づつに分けて三群作りました。果たして被験者が何人必要であったかを算出したいのですが，算出方法がわかりません。ご教授いただければありがたいです。いろいろ調べてみましたが理解できませんでした。どうかよろしくお願いします。困っています。

No.07911　Re: サンプルサイズの検出力分析　　【青木繁伸】　2008/10/10(Fri) 22:34

R だと，

power.anova.test {stats}	R Documentation

Power calculations for balanced one-way analysis of variance tests

Description

Compute power of test or determine parameters to obtain target power.

というのがあります。
SPSSとかにもあるかも知れません。そのほか，フリーソフトでは G*power でも出来たと思いますけど。

> example(power.anova.test)

pwr.n.> power.anova.test(groups=4, n=5, between.var=1, within.var=3)

     Balanced one-way analysis of variance power calculation 

         groups = 4
              n = 5
    between.var = 1
     within.var = 3
      sig.level = 0.05
          power = 0.3535594

 NOTE: n is number in each group 


pwr.n.> # Power = 0.3535594
pwr.n.> 
pwr.n.> power.anova.test(groups=4, between.var=1, within.var=3,
pwr.n.+                  power=.80)

     Balanced one-way analysis of variance power calculation 

         groups = 4
              n = 11.92613
    between.var = 1
     within.var = 3
      sig.level = 0.05
          power = 0.8

 NOTE: n is number in each group 


pwr.n.> # n = 11.92613
pwr.n.> 
pwr.n.> ## Assume we have prior knowledge of the group means:
pwr.n.> groupmeans <- c(120, 130, 140, 150)

pwr.n.> power.anova.test(groups = length(groupmeans),
pwr.n.+                  between.var=var(groupmeans),
pwr.n.+                  within.var=500, power=.90) # n = 15.18834

     Balanced one-way analysis of variance power calculation 

         groups = 4
              n = 15.18834
    between.var = 166.6667
     within.var = 500
      sig.level = 0.05
          power = 0.9

 NOTE: n is number in each group

No.07912　Re: サンプルサイズの検出力分析　　【タロー】　2008/10/10(Fri) 22:57

早々のご検診ありがとうございます。さらに質問させてください。こちらに示されているbetween.var = 1 within.var = 3は何の数字でしょうか？また，power =0.3535594とn=11.92613はどのように算出されたのでしょうか？

No.07913　Re: サンプルサイズの検出力分析　　【青木繁伸】　2008/10/10(Fri) 23:06

使用法は R のオンラインヘルプ参照。
計算方法は，R のプログラムを読む。
参考書は R のオンラインヘルプに以下のものが記載されている。と，書いてから，再度オンラインヘルプを参照したが，参考文献の記載はなかった。よほど一般的？そんなはずはないはずだが。R のプログラムを読むしかないかな。。というか，計算方法を知りたいのか，計算結果が得られればよいのか。。。

No.07914　Re: サンプルサイズの検出力分析　　【タロー】　2008/10/10(Fri) 23:12

ありがとうございます。プログラムを読むしかないですか...計算方法が知りたいのです。

No.07915　Re: サンプルサイズの検出力分析　　【青木繁伸】　2008/10/10(Fri) 23:17

G*power は，バージョン3が出たようだ。
http://www.psycho.uni-duesseldorf.de/abteilungen/aap/gpower3/
参考文献リストは以下にあるが，Cohen(1988) を見ろと書いてある。
http://www.psycho.uni-duesseldorf.de/aap/projects/gpower/literature.html

No.07916　Re: サンプルサイズの検出力分析　　【タロー】　2008/10/10(Fri) 23:23

ありがとうございました。