我正在尝试计算 Fisher 在 R 中对 2x2 列联表的精确检验的功效。但不确定我是否使用了正确的参数。谁能看看我下面的逻辑并告诉我这是否是正确的方法?
使用的 R 函数(在 中statmod
):
power.fisher.test(p1, p2, n1, n2, alpha=0.05, nsim=100, alternative="two.sided")
p1
- 第一个要比较的比例。
p2
- 要比较的第二个比例。
n1
- 第一个样本量。
n2
- 第二个样本量
我的示例列联表:
Yes|No
Team1 a | b
Team2 c | d
我计算上面列联表的比例和样本大小的方式。这些是正确的吗?还是我需要做任何修改?
p1=a/(a+c)
p2=b/(b+d)
n1=a+c
n2=b+d