[英]Confidence interval from binomial distribution: different results using `qbinom` and `binom.test`
我希望使用qbinom
和binom.test
计算比例的 CI 的结果相同,但它们实际上略有不同:
success <- 360
n <- 1226
lci <- qbinom(0.025, n, success/n)/n
uci <- qbinom(0.975, n, success/n)/n
c(lci, uci)
[1] 0.2683524 0.3189233
binom.test(success, n, success/n)$conf.int
[1] 0.2682571 0.3200123
我在这里想念什么?
binom.test
function 的源代码使用qbeta
而不是qbinom
,因为这是公认的精确二项式置信区间公式。 其中x
是成功的次数,由 binom.test 给出的 lci 和binom.test
是:
p.L <- function(x, alpha) {
if (x == 0)
0
else qbeta(alpha, x, n - x + 1)
}
p.U <- function(x, alpha) {
if (x == n)
1
else qbeta(1 - alpha, x + 1, n - x)
}
alpha <- (1 - 0.95)/2
lci <- p.L(x, alpha)
uci <- p.U(x, alpha)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.