r - 麦克尼马尔精确检验

Question

我正在 R 中对以下数据进行 McNemar 测试：

获得以下结果：

我理解结果，但是，有人可以向我解释如何计算置信区间？

Answer 1

您可以在这个小插图中阅读更多内容，也可以查看代码。使用此 wiki 图片进行说明：

优势比是 b / c，在你的情况下是150/86 = 1.744186。您可以围绕成功的比例构建一个二项式置信区间，将 b 视为成功，将试验次数视为 b + c。

在代码中，他们使用这段代码来计算：

library(exactci)
CI = binom.exact(150,86+150,tsmethod = "central")
CI

data:  150 and 86 + 150
number of successes = 150, number of trials = 236, p-value = 3.716e-05
alternative hypothesis: true probability of success is not equal to 0.5
95 percent confidence interval:
 0.5706732 0.6970596
sample estimates:
probability of success 
             0.6355932 

你有 b 的上限和下限，那么优势比是 p / 1- p ：

CI$conf.int/(1-CI$conf.int)
[1] 1.329228 2.300979

binom.exact 的小插曲指出：

'central' 方法给出了 Clopper-Pearson 区间，而 'minlike' 方法给出了 Stern (1954) 提出的置信区间（参见 Blaker, 2000）。

所以这是估计二项式置信区间的众多方法之一。