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

No.14165　３群のサンプルサイズの求め方　　【キラキラ】　2011/01/16(Sun) 14:41

初めまして。いつもこのサイトを見て，勉強させていただいております。本日は，3群のサンプルサイズについて，教えていただきたく，質問をさせていただきます。
ある尺度の得点を，先行研究からA群66点，B群54点，C群54点，標準偏差13.5と想定し，有意水準5％（両側）で検出力を80％で計算した場合，サンプルサイズは各群21名としました。この計算は，2群の平均値の比較（66点と54点）として，フリーのソフトで計算したのですが，2群の式を当てはめる考え方は間違っていますでしょうか。間違っているとすれば，上記のサンプルサイズは，何名になるのでしょうか。計算式と一緒に教えていただけると，大変助かります。よろしくお願い申し上げます。

No.14166　Re: 3群のサンプルサイズの求め方　　【青木繁伸】　2011/01/16(Sun) 14:52

> 2群の式を当てはめる考え方は間違っていますでしょうか

たぶん，間違っていると思います。
R の power.anova.test だと以下のようになります。
使用方法は，example(power.anova.test) を参考に。

> groupmeans <- c(66, 54, 54)
> power.anova.test(groups=length(groupmeans),
+ between.var=var(groupmeans),
+ within.var=13.5^2, power=0.8)

     Balanced one-way analysis of variance power calculation 

         groups = 3
              n = 19.32532
    between.var = 48
     within.var = 182.25
      sig.level = 0.05
          power = 0.8

 NOTE: n is number in each group

ちなみに，2群比較の設計は

> power.t.test(delta=66-54, sd=13.5, power=0.8, strict=TRUE)

     Two-sample t test power calculation 

              n = 20.87274
          delta = 12
             sd = 13.5
      sig.level = 0.05
          power = 0.8
    alternative = two.sided

 NOTE: n is number in *each* group

No.14167　Re: 3群のサンプルサイズの求め方　　【キラキラ】　2011/01/16(Sun) 15:38

青木先生　お忙しいところ，早々に，ご返事をいただき，ありがとうございました。
3群のサンプルサイズが，きちんと計算式で求められることがわかりました。今後，example(power.anova.test) を参考にして，サンプルサイズを求めていきたいと思います。