2x2の因子計画に必要なサンプルサイズを計算しようとしています。2つの質問があります。
1）パッケージpwrと一元配置分散分析関数を使用して、次のコードを使用して必要なサンプルサイズを計算しています
pwr.anova.test(k = , n = , f = , sig.level = , power = )
ただし、一元配置分散分析よりもグループ平均の推定に効率的であるため、二元配置分散分析を検討したいと思います。私が見つけることができた双方向の分散分析関数はありません。これを行うためのパッケージまたはルーチンが[R]にありますか？
2）さらに、一元配置分散分析を使用しているため、サンプルサイズがより保守的（つまり大きく）になると想定しても安全ですか？
因子 A、因子 B、および AxB を含む 2 x 2 ANOVA では、これら 3 つの効果のそれぞれについて個別の統計的検出力推定値が得られます。
G Power 3 は、要因計画における効果の検出力を推定するための無料のソフトウェアといくつかの明確なチュートリアルを提供します: http://www.psycho.uni-duesseldorf.de/abteilungen/aap/gpower3/user-guide-by-design
検索した後、オンラインでこれに対する解決策が見つかりませんでした。
私がお勧めするのは (方法を知っている場合)、シミュレーションを使用してこれをプログラムすることです。その方法がわからない場合は、「検出力分析を実現するために、双方向分散分析のシミュレーションを作成するにはどうすればよいですか」という SO の質問を書いて、人々があなたを助けてくれるものを見てください :)
また、ここでコードを確認することから始めることもできます。
http://www.rforge.net/doc/packages/NCStats/power.sim.html
シミュレーションによる電力計算の開始用。
ジェロミーが書いたことに注目してください - この検出力分析は複数の結果に対するものです。
興味深いテーマです。フォローアップしたいと思います。
一番、
タル