如何计算与二项式比例增加相关的置信区间?

机器算法验证 r 置信区间 二项分布 部分
2022-03-21 11:01:07

不知道你是否可以帮助我,但这是问题所在。

我有两个二项式比例 A 和 B (95% CI) -

A=2%±0.2%,B=3%±0.2%.

换言之,B 的比例比 A 高 50%。

商界人士倾向于以百分比增长来谈论这种差异。换言之,B 的比例比 A 高 50%。我如何计算 B 的比例比 A 增加 50% 的置信区间。

我对 Fieller 的公式进行了一些研究,据说该公式可以提供精确的百分比效果。但是,除了不理解 Fieller 定理之外,我不知道如何在 R 中计算这个。

任何人都可以帮忙吗?

编辑:例如,如果这些间隔的 CI 为 95%,则 50% 增加的 CI = 98.52%。

1个回答

按照 whuber 的链接到维基百科,你有

假设是联合正态分布的,并且不太接近零(即更具体地说,与相比,的标准误差很小)abbbb

Var(ab)=(ab)2(Var(a)a2+Var(b)b2).

虽然实际上你想要Var(BA)

如果你的 95% CI 是的方差,所以取平方根并乘以 1.96 得到±0.002AB(0.002/1.96)20.00000104Var(BA)0.00846

BA1.5±0.18

如果您必须将其转换为百分比(我认为它令人困惑而不是启发),那么它就变成了

B 的比例比 A 高 50%,正负 18%,即高出 32% 到 68%。

在 R 中,你可以通过类似的东西来模拟它

> n <- 1000000
> A <- 0.02 + (0.002 / qnorm(0.975)) * rnorm(n)
> B <- 0.03 + (0.002 / qnorm(0.975)) * rnorm(n)
> C <- B / A
> quantile(C,  probs = c(0.025, 0.5, 0.975))
    2.5%      50%    97.5% 
1.333514 1.499955 1.697418 

这是相当接近的。