我正在尝试计算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
答案 0 :(得分:1)
Fisher精确检验用于检验组间差异。因此,您必须计算每个团队中“是”的概率:
p1 = a/(a + b)
p2 = c/(c + d)
n1 = a + b
n2 = c + d