计算R中Fisher测试功率函数的正确参数

时间:2012-12-06 13:14:59

标签: r statistics

我正在尝试计算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

1 个答案:

答案 0 :(得分:1)

Fisher精确检验用于检验组间差异。因此,您必须计算每个团队中“是”的概率:

p1 = a/(a + b)

p2 = c/(c + d)

n1 = a + b
n2 = c + d