パワーアナリシス(一元配置分散分析)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 群間分散 within.var 群内分散 sig.level 有意水準(第1種の過誤の確率 α)デフォルトは 0.05 power 検出力(1から第2種の過誤を引いた確率 1-β) 注: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 **関数からの出力 $ note : 注 $ method : 方法 *例1 4 群,群あたりのサンプルサイズが10,群間分散が 1.28,群間分散が 0.35,有意水準 0.05 のときの検出力を求める > power.anova.test(groups=4, n=10, between.var=0.37, within.var=0.35, sig.level=0.05, power=NULL) Balanced one-way analysis of variance power calculation groups = 4 n = 10 between.var = 0.37 within.var = 0.35 sig.level = 0.05 power = 0.9973492 NOTE: n is number in each group *例2 4 群,群間分散が 1.28,群間分散が 0.35,有意水準 0.05,検出力=0.8 のときに必要なサンプルサイズを求める > power.anova.test(groups=4, n=NULL, between.var=0.37, within.var=0.35, sig.level=0.05, power=0.8) Balanced one-way analysis of variance power calculation groups = 4 n = 4.543592 between.var = 0.37 within.var = 0.35 sig.level = 0.05 power = 0.8 NOTE: n is number in each group