(R) 找出接受概率为 0.1 的不良品总体比例

(R) Finding proportion of population defectives at probability 0.1 acceptance

我正在使用以下 R 代码:

library(AcceptanceSampling)

x <- OC2c(50, 2, type="hypergeom", N=4000)
plot(x, xlim=c(0,0.2))

生成图:

我想求出P(accept)(Y轴)为0.1时的比例。有没有办法在 R 中做到这一点?根据 https://cran.r-project.org/web/packages/AcceptanceSampling/AcceptanceSampling.pdf 上的文档,我感觉这个包不允许直接计算它。

据我所知,你不能准确地得到这个,但你可以接近。 x 是一个 S4 对象,您在这里感兴趣的插槽是 pacceptpd。我们可以得到 pd 值与 paccept 最接近 到所有 paccept 值的 0.1:

idx <- which.min(abs(x@paccept - 0.1))
round(cbind(paccept = x@paccept[idx], pd = x@pd[idx]), 3)
#      paccept  pd
# [1,]    0.11 0.1