OC 曲线中实际缺陷率的估计我觉得这个问题有点奇怪。
无论如何，我的问题是；
-手数：10,000
-样本大小； 500
-AQL（可接受的失败）：0
- 样品的实际缺陷率：0% (0/500)
有没有办法估算剩余 9,000 个的缺陷率?
抽样检验中检验不良率为0时，如何估算不良率？
OC 曲线SPC回复：OC 曲线中实际缺陷率的估计嗨，迦勒，
感谢您的友好和频繁的评论。
现在似乎清楚了一些。
prob = 0;
K = 9000;
alpha = 0.05; // change to your preference
While( prob < alpha,
prob = Hypergeometric Probability( 10000, K, 1000, 85 );
K--;
);
Print( K + 1 );
//
回复：OC 曲线中实际缺陷率的估计在我的代码中发现了另一个错误。 对于上限，它应该从上往下看。 它适用于 0 和 1 缺陷案例的原因只是分布的侥幸。 我已经更新了初始帖子中的代码并添加了另一个注释。
回复：OC 曲线中实际缺陷率的估计嘿@IkjuShin ！
为长时间的延迟道歉。 一定是在某处错过了通知。
回顾我之前的回复，我可以确认错误不在代码中，而是在编码器中（这是典型的）。 在公式中，我应该有“prob>=alpha”而不是“prob>=1-alpha”。 我的意图是正确的，但我的公式是错误的。 我已经返回并进行了必要的更正，并为未来的读者添加了注释。
对于那个很抱歉！ 不过很好抓！
回复：OC 曲线中实际缺陷率的估计你好，卡金
感谢您的友好评论。 我知道你评论的意思。
但是，无论 X 的值如何（0、1 或 2），结果始终为 0。
我不知道脚本中有什么样的错误。
回复：OC 曲线中实际缺陷率的估计如果通过您的问题，您的意思是要计算剩余 9500 个单元中可能存在的缺陷数量的上限，因为我在 500 个样本中仅发现 1 个，那么这是我冗长的说法“是的！”。
回复：OC 曲线中实际缺陷率的估计你好，
这是非常有帮助的评论。 我想就此提出进一步的问题。
N=10,000 n = 500，n = 1 (1/500) 中的实际缺陷数，
在这种情况下，估计剩余 9500 个缺陷的脚本是否正确？
谢谢
prob=1;
K=0;
alpha=0.05; // change to your preference
While(prob>=1-alpha,
Print(K);
K++;
prob=Hypergeometric Probability(10000,K,500,1)
);
回复：OC 曲线中实际缺陷率的估计嘿@IkjuShin ！
这实际上是一个非常好的问题！ 一个比你想象的更经常被问到的问题。
我将把我的答案分成两部分。 首先，我将说明 OC 曲线实际上告诉您的内容。 接下来，我会直接回答你的问题。
首先，OC 曲线告诉您在给定一定的接受水平（在这种情况下，样本中的缺陷为 0）和假定的真实缺陷水平（假定的批次中存在的缺陷百分比）的情况下，批次被接受的机会。 换句话说，您正在玩曲线的假设游戏。 例如，根据您那里的输出，如果真正的缺陷百分比约为 0.58%（即 10,000 个批次中有 58 个缺陷），那么您将在大约 5% 的时间内接受该批次。 这似乎是合理的，因为您有相当大的样本量。 基本上，OC 曲线可以帮助您了解对于给定的缺陷率，您应该多久接受一次批次。 它有助于确定您可能有多少缺陷，但不是明确的。
所以这给我们带来了你的问题：如果我看到 0 个缺陷，我应该如何估计可能有多少缺陷？ 显然，公然说“必须有 0 个缺陷”似乎太冒险了。 解决这个问题的一种方法是计算真实缺陷率的置信区间。 假设一个二项式模型和 1-alpha 的置信水平，你会得到一个看起来像 [0,1-alpha^(1/N)] 的区间，其中 N 是样本数（我附上了一个简短的这篇文章供参考）。 对于您的示例，如果 1-alpha=0.95（95% 置信度），那么您的上限将为 0.00597（大约 0.6%）。 当您的点估计实际上为 0 时，这只是尝试和估计真实比率的众多方法之一，但它肯定是最简单的。
现在您还有一个问题，即在剩余的 9,000 个未抽样的单位中可能有多少缺陷。 这需要不同的答案。 为此，我们需要在超几何模型下进行操作（这里是维基百科链接了解更多信息）。 从技术上讲，这就是我们应该一直在努力的方向，但是对于我们正在使用的这些非常大的尺寸，之前的结果不会有太大的不同。 重新表述您的问题，对于 N=10,000 个总单位的超几何模型，我们对其中的 n=500 进行抽样，我们要问的是，考虑到我们有 k=0 个缺陷，K（批次中的缺陷数量）的可能值是多少我们的样品？ 一个快速回答这个问题的方法是对 K 做一些类似于我们在假设二项式模型时对 p 所做的事情：计算给定置信度 1-alpha 的上限。 下面是一段 JSL 代码，可以为我们做到这一点：
prob=0;
K=9500;
alpha=0.05 // change to your preference
While(prob<ALPHA></ALPHA>Print(K+1);
基本上，它寻找 K 的最大值，使得在 500 个单位的样本中观察到 0 个缺陷的概率不会低于 alpha 1-阿尔法.因此，对于您的示例，如果 1-alpha=0.95，则 K 的上限为 1。因此，对于您的示例，如果 1-alpha=0.95，则 K 的上限为 58。
注意：更正了公式中的一个错误；应该寻找 K 其中 Pr(0 个缺陷给定 K 总)=alpha。
注2：又一个修复；公式应该从上到下寻找 K 的上限
希望有帮助！
