2x2 分割表の R でのフィッシャーの正確確率検定の検出力を計算しようとしています。しかし、正しいパラメーターを使用しているかどうかはわかりません。以下の私のロジックを見て、これが正しいアプローチであるかどうか教えてもらえますか?
使用される R 関数 (でstatmod
):
power.fisher.test(p1, p2, n1, n2, alpha=0.05, nsim=100, alternative="two.sided")
p1
- 比較する最初の比率。
p2
- 比較する 2 番目の比率。
n1
- 最初のサンプル サイズ。
n2
- 2 番目のサンプル サイズ
私の分割表の例:
Yes|No
Team1 a | b
Team2 c | d
上記の分割表の比率とサンプルサイズを計算する方法。これらは正しいですか?または、何らかの修正を加える必要がありますか?
p1=a/(a+c)
p2=b/(b+d)
n1=a+c
n2=b+d