library(binom)
x <- 10
n <- 20
binom.confint(x, n, methods = "prop.test", conf.level = 0.80)
# method x n mean lower upper
# 1 prop.test 10 20 0.5 0.299298 0.700702
binom.confint(x, n, methods = "prop.test", conf.level = 0.95)
# method x n mean lower upper
# 1 prop.test 10 20 0.5 0.299298 0.700702
但它确实适用于binom.test
binom.test(x, n, conf.level=0.80)$conf.int
# [1] 0.3381709 0.6618291
binom.test(x, n, conf.level=0.95)$conf.int
# [1] 0.2719578 0.7280422
1条答案
按热度按时间fd3cxomn1#
这似乎是软件包中的一个bug。根据
methods
参数的docs:prop.test
:相当于prop.test(x = x, n = n, conf.level = conf.level)$conf.int
。但是,如果我们看看来源:
您可以看到
conf.level
参数没有传递给prop.test()
,因此它始终使用默认值(0.95
)。它应该说:这将导致不同CI的输出不同:
我不完全确定如何解决这个问题。Github repo似乎是只读的,因此您无法打开问题或合并请求。也许您可以使用CRAN page上的详细信息联系软件包维护人员。在此期间,你可以尝试下面的热修复。
热修复
这不是构建健壮代码库的方法,但如果你需要一个短期解决方案,你可以像这样修补函数:
然后我们可以运行
binom.confin2
: