パワーアナリシス(一元配置分散分析)power.anova.test

*機能

一元配置分散分析のパワーアナリシスを行う。~
一元配置分散分析の検出力を求めること,また,所定の検出力を得るためのパラメータ(標本サイズなど)を求める。

*使用法

 power.anova.test(groups = NULL, n = NULL, between.var = NULL,
                  within.var = NULL, sig.level = 0.05, power = NULL)

*引数

 groups		群の数(水準の数)
 n	        群あたりのサンプルサイズ
 between.var    群間分散/n
 within.var     群内分散
 sig.level      有意水準(第1種の過誤の確率 α)デフォルトは 0.05
 power          検出力(1から第2種の過誤を引いた確率 1-β)

注:どういうわけか,between.var は群間分散そのものではないので注意!!!~
注:groups, n, between.var, power, within.var, 及び sig.level のうち,どれか一つだけ NULL にして関数を呼ぶ。そうすると,NULL であったパラメータの値が計算される。sig.level はデフォルトで 0.05 なので,sig.level を計算したいときにはそれを明示的に NULL にして関数を呼ばなければならない。

*戻り値

**関数への入力

 $ groups
 $ n
 $ between.var
 $ within.var
 $ sig.level
 $ power

**関数からの出力
 $ note      : 注
 $ method    : 方法

*例1 4 群,群あたりのサンプルサイズが10,群間分散が 304.2,群内分散が 100.87,有意水準 0.05 のときの検出力を求める
between.var は,304.2 ではなく,304.2/10=30.42
 > power.anova.test(groups=4, n=10, between.var=30.42, within.var=100.87, sig.level=0.05, power=NULL)
 
      Balanced one-way analysis of variance power calculation 
 
          groups = 4
               n = 10
     between.var = 30.42
      within.var = 100.87
       sig.level = 0.05
           power = 0.66122356442945618
 
  NOTE: n is number in each group 

*例2 4 群,群間分散が 304.2,群内分散が 100.87,有意水準 0.05,検出力=0.8 のときに必要なサンプルサイズを求める
between.var は,304.2 ではなく,304.2/10=30.42
 > power.anova.test(groups=4, n=NULL, between.var=30.42, within.var=100.87, sig.level=0.05, power=0.8)
 
      Balanced one-way analysis of variance power calculation 
 
          groups = 4
               n = 13.069996325957677
     between.var = 30.42
      within.var = 100.87
       sig.level = 0.05
           power = 0.8
 
  NOTE: n is number in each group


トップ   編集 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS