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

No.08071　Rによるサンプル数検定について　　【金木犀】　2008/10/24(Fri) 18:16

3つの測定機器で測定する実験を行なっています。
過去のデータを元にサンプル数を求めたいと考えています。

Rで

> power.t.test(delta=25,sd=30,power=0.8,sig.level=0.05/3)

    Two-sample t test power calculation

             n = 31.6212
         delta = 25
            sd = 30
     sig.level = 0.01666667

         power = 0.8
   alternative = two.sided

 NOTE: n is number in *each* group

このeach groupという意味がよく分からないのですが…
3つの測定機器で測定するので3群となるので
それぞれの群が32人だから31.6*3=94.8としたらいいという事なのでしょうか？？
ペアで比べているので31.6*2=63.2とした方がいいのでしょうか？？

No.08074　Re: Rによるサンプル数検定について　　【青木繁伸】　2008/10/24(Fri) 19:19

3つの機械で，同じサンプルを測っているのでしょうか？

別なら，each というのは，t.test が二群の検定だから各群が32個のデータ
同じものを測っているのなら，power.t.test にオプションを指定しないと。。。
> power.t.test(delta=25,sd=30,power=0.8,sig.level=0.05/3, type="paired")

     Paired t test power calculation 

              n = 18.03077
          delta = 25
             sd = 30
      sig.level = 0.01666667
          power = 0.8
    alternative = two.sided

 NOTE: n is number of *pairs*, sd is std.dev. of *differences* within pairs 

No.08075　Re: Rによるサンプル数検定について　　【金木犀】　2008/10/24(Fri) 22:30

青木先生いつも丁寧なご指導ありがとうございます。

＞3つの機械で，同じサンプルを測っているのでしょうか？

3つの機械で同じサンプルを測定しておりますので
n = 18.03077
n is number of *pairs*,
ということは各群18人のサンプルでよいという事でしょうか？？

No.08076　Re: Rによるサンプル数検定について　　【青木繁伸】　2008/10/24(Fri) 22:49

> ということは各群18人のサンプルでよいという事でしょうか？？

まあ，対応があるので，18サンプルを3つの方法で測ればよいと言うこと，だと思いますけど。。。

No.08086　Re: Rによるサンプル数検定について　　【金木犀】　2008/10/26(Sun) 19:30

青木先生ありがとうございます。
非常に勉強になりました。