跳到正文
PharmaBuddy· PharmaBuddy·· 6 小时前AI 评分71

解析GMP附录8取样件数规则与验收抽样方案的差异:批量300与301之间件数由18降至10

批量300与301:取样件数从18件掉到10件

AI 导读

一篇取样方法解析指出,GMP附录8的√n+1只规定取样件数、不含接收数与风险,批量由300增至301时取样件数由18件降至10件。文中按c=0、批不合格率2%复算,接收概率在N=300与N=301之间由68.1%升至81.7%,相差13.6个百分点。

正文

检察官翻到一份原料的取样 SOP,指着其中一行问:

“这个物料,批量 1200 桶,为什么取 35 桶?”

“附录 8 的规定,√n+1,35 桶。”

“好。那这个方案,当这批货的不合格率已经高到不能接受的时候,还有多大概率被你们放行?”

停顿。然后是一句几乎所有药企 QA 都说过的话:

“……我们是按规定取的。”

这个回答没有错。但它答的不是同一个问题。

√n+1 从来不是抽样方案。它是一条取样件数的规则。它可以告诉你从 1200 桶里该舀出多少桶,但它无法告诉你,这 35 桶的数据能不能支撑“整批合格”这个结论——更无法告诉你,这个结论错的时候,错得有多离谱。

一、取样里藏着两个问题,大多数人只答了一个

把“取样”这两个字拆开,里面其实是两件事。

第一件:取几件。
从 1200 桶里,取多少桶,才算“代表”了这一整批?这是代表性问题。它的评价标准是:你取的这些,和整批像不像。

第二件:收不收。
拿到数据以后,凭这 35 桶的结果,敢不敢判定 1200 桶合格?这是风险问题。它的评价标准是:你判错的概率有多大。

代表性靠取样件数解决,风险靠抽样方案解决。前者是统计学里的 sampling for testing,后者是 acceptance sampling——两个不同的分支,两套不同的数学工具。

问题就出在这里:GMP 附录 8 给的,是前者。而几乎所有的审计追问,问的是后者。

附录 8 第二十条的原文是这样的:

“被抽检的物料与产品是均匀的,且来源可靠,应按批取样。若总件数为 n,则当 n≤3 时,每件取样;当 3<n≤300 时,按√n+1 件随机取样;当 n>300 时,按√n/2+1 件随机取样。”

注意这句话的前半段——“被抽检的物料与产品是均匀的,且来源可靠”。这是整条规则能够成立的前提,也是被引用时最常省略的部分。

它的意思是:这条规则管的是“均匀物料怎么取才有代表性”。至于取回来怎么判、判错的风险有多大,它没有管,也不该由它管。

附录 8 这句话里的三个关键词

前提均匀——物料本身是均匀的

物料不均匀时,取再多样也不能代表整批

方法随机——不是随意,要可追溯

取样点、取样人、取样时间需可追溯

目的代表——这一批的缩影

目标是“像整批”,不是“能判收放”

没有接收数、风险、概率

这三样恰恰是审计要问的

二、√n+1 长什么样

先把它算出来。取整方式不同,结果会差一件,下表按“向上取整”给:

附录 8 规则下的取样件数

N=10011 件

√100+1=11

N=28918 件

√289+1=18

N=30018 件

√300+1≈18.3,有进位习惯的取 19

N=30110 件

√301/2+1≈9.7——批量更大,件数砍半

N=100017 件

√1000/2+1≈16.8

N=1000051 件

批量涨 100 倍,件数只涨 4.6 倍

N=100000160 件

此时件数占比只有 0.16%

表中数据按附录 8 公式复算

把这条曲线画出来,三个性质一目了然。

图 1 批量 300 两侧,取样件数出现断裂

一是次线性。批量从 100 涨到 10000,涨了 100 倍,件数只从 11 涨到 51,涨了 4.6 倍。这是它的优点,也是它在工业界活了这么多年的原因——取样成本不会随批量线性膨胀。

√n这个形式,和“样本均值的标准误与√n成反比”是同源的。想让估计精度翻一倍,样本量要翻四倍。用平方根,逻辑上说得通。

二是不连续。300 是一条硬边界。批量从 300 变成 301,件数从 18 掉到 10,几乎砍半。批量更大、责任更重的那一批,反而取得更少。

这个跳变很难用统计模型解释。√n有它的道理,但“+1”这个修正项、“300”这个分界点、以及超过 300 之后把系数从 1 改成 1/2——更像是经验拟合的结果,而不是从某个分布推导出来的。

三是与质量要求无关。公式里只有 N。没有 AQL,没有 LQL,没有 α,没有 β,也没有接收数 c。

一个只认批量、不认质量要求的规则,注定回答不了“判错的概率有多大”。因为“判错”这件事,必须先有“判”的标准。

三、给它补上一个 c,问题就全露出来了

√n+1 给的只有件数。要评价它,必须先替它假设一个接收数 c。

取多少合适?没有依据。讲义里的做法是c=0,也就是“一个不合格都不许有”。这不是苛刻,这是唯一能算的选择——因为规则本身没有给出任何关于 c 的信息。

补上 c=0,画 OC 曲线。

图 2 批不合格率固定 2%、c=0 时的接收概率

横轴是批量,纵轴是接收概率,批不合格率固定在 2%:

c=0、p=2% 时的接收概率

N=2588.6%

N=10080.1%

N=30068.1%

N=30181.7%

批量只多一件,接受概率跳升 13.6 个百分点

N=100070.9%

N=1000035.7%

N=1000004.0%

生产者风险 α 高达 96%

按二项分布复算,与文中数据一致

三条结论,一条比一条严重。

第一条:批量越大,方案越严,而且严到不合理。

批量 10 万件、不合格率 2% 的批,接收概率只有 4.0%。换句话说,生产者风险 α 高达 96%。一个合格率 98% 的供应商,在这套规则下几乎必被拒收。这不是在筛选质量,这是在劝退供应商。

第二条:300 这个断点,是真实的、可观测的、无法辩护的。

同一家供应商、同一批质量、同样是 2% 的不合格率:报 300 件,接收概率 68.1%;报 301 件,81.7%。批量多一件,严度松一档,一步跨过 13.6 个百分点。

如果审计员拿两个批量的物料来对比,这个跳变是解释不通的。因为它背后没有统计学,只有分段函数。

第三条:也是最根本的一条——它不在 OC 体系里。

OC 曲线有一个基本性质:方案的接收概率只由批质量 p 和方案 (n, c) 决定,与批量 N 基本无关(除极小批量)。这正是它被称为“特征曲线”的原因——它是方案的特征,不是批量的特征。

而√n+1 的接收概率,随 N 一路漂移。漂移本身就是证据:它不是抽样方案,它是取样件数规则,被误当成了抽样方案。

把取样件数当成抽样方案,是一个范畴错误。错的不是公式,是用法。

四、GB/T 2828.1 强在哪里,偏在哪里

工业界的标准答案是 GB/T 2828.1(等同采用 ISO 2859-1,往上溯源是美军标 MIL-STD-105)。它的结构是三段式:

定 AQL → 查样本量字码 → 查表得 (n, c)

AQL 这三个字容易被误读。它不是“合格质量水平”,而是“接收质量限”(Acceptable Quality Limit)。它的含义是:不合格率达到这个水平时,仍然认为值得以高概率接收。它是一条底线,不是一个目标。

先看它的形状。

图 3 样本量与接收数如何改变曲线的陡峭度

同样守一条底线,曲线的陡峭度可以差很多。越陡,好批和坏批的分界越清楚;越平,判定的含糊区间越宽。而陡峭度由 (n, c) 共同决定——这就是为什么“选了 AQL”还远远不够,你还得选方案。

再说它强在哪。

很多人以为 GB/T 2828.1 的精华是那张主表。不是。它的精华是转移规则。

四个检验状态与切换条件

加严正常 → 加严:连续 5 批中 2 批不接收

怀疑的对象从“这一批”转向“这个过程”

回正常加严 → 正常:连续 5 批接收

恢复的路径写死,不靠人情判断

放宽正常 → 放宽:转移得分 ≥30 且生产稳定

省检验时间与费用,但不是放松标准

暂停加严 → 暂停:累计 5 批不接收

得分规则:Ac≥2 且加严后接收 +3 分;Ac=0 或 1 时接收得 2 分

这是一个闭环。加严检验的价值不在于“这一批更严格”,而在于它把怀疑的对象从“这一批”转移到了“这个过程”——连续出问题,说明过程失控了,此时该做的是停下来查原因,而不是继续一批一批地挑。

这套动态反馈,是 GB/T 2828.1 真正被低估的部分。很多企业只用了那张表,把“加严/放宽”当成了摆设。

再说它偏在哪。

AQL 体系的设计目标是:让质量达到 AQL 的批,以高概率被接收。也就是说,它在优化生产者风险 α(通常在 5% 量级)。

消费者风险 β 呢?是查表时顺带出来的结果,不是一个被预先锁定的设计目标。你选定 AQL 和检验水平,表给你一个 (n, c),β 是多少,事前并不知道——画完 OC 曲线才知道。

在工业来料检验里,这个取舍是合理的:供应商和采购方是长期博弈,冤枉好供应商的代价很高。但在药品里,这个取舍要翻过来。

五、药品的答案:把 β 摆到桌面上

药品的最终消费者是患者。

AQL 思路问的是:怎么才能不冤枉合格的工厂?
GMP 思路问的是:怎么才能不把坏批次放到患者手上?

两个问题的答案不一样。

两种风险,两种后果

α生产者风险批本身合格,抽样判了拒收

后果:工厂受损,成本上升,好批被退回

β消费者风险=患者风险批本身不合格,抽样判了接收

后果:坏批次流向患者——这是 GMP 的噩梦

后者才是监管盯的那一头。

所以药品关键取样的方案设计顺序,应该反过来:

β 优先的四步

1定 LQL:多差算不可接受

例如“批阳性率达到 5% 就属于不能接受的坏批”

2定 β:坏批最多允许以多大概率放行

常规 ≤10%,高风险无菌产品 ≤5%

3反推 n 与 c

解出同时满足 β 的方案,画出 OC 曲线

4回头核算 α

在 AQL 水平上,好批被误拒的概率工厂能否承受

这个顺序和 GB/T 2828.1 正好相反。它是 β 优先,α 事后核算。

举一个能算出来的例子:N=1000,AQL=2%,RQL=8%,要求 α 和 β 都不超过 10%。

解出来是(82, 3)——取 82 件,允许 3 件不合格。
复核:p=2% 时接收概率 91.8%,α=8.2%;p=8% 时接收概率 9.8%,β=9.8%。
两个风险同时被锁在 10% 以内。

这就是“双向锁定”。而单纯查 GB/T 2828.1 的主表,只能锁住一个点——另一个点是表给你的,不是你要的。

反过来,如果直接拿工业品的 AQL 表套在无菌、微生物、关键原料上,最容易被质疑的一句话是:这套方案的设计没有以患者风险为首要目标,也没有论证 LQL 与 β,取样科学性不足。

Minitab 里对应的是“按属性抽样验收 → 创建抽样计划”,输入 AQL、RQL、α、β,输出 (n, c) 和 OC 曲线。这个功能大部分质量部门都有,但用得不多。

顺带给一段可以直接背下来的答法。当被问到 “Why did you select this sampling plan?”:

“The sampling plan is designed primarily to control consumer risk (patient risk): to limit the probability of accepting a poor-quality lot at the defined LQL level; producer risk α is a secondary consideration. The OC curve has been established to demonstrate both risks.”

译:本取样方案首要管控消费者(患者)风险,限定达到 LQL 不良水平的批次被误接收的概率;生产者风险作为次要考量,OC 曲线同时论证两项风险。

这句话之所以有效,是因为它把“我按规定取的”换成了“我为患者风险设计过”。前者是合规,后者是科学。

六、还有一条更聪明的路:不要数不合格品

前面说的都是计数型:把每件判成合格或不合格,然后数不合格的个数。

但药品的检验数据,大部分是连续的。

含量 96.5%,水分 0.32%,熔点 216.8℃,有关物质 0.08%——这些不是“合格/不合格”,这些是数值。

对数值,有更高效的工具:计量型抽样。

图 4 计量型判定看的是整条分布,不是某一件

它的判定不数个数,而是算一个比值:

QL =(均值 − L)/ s   QU =(U − 均值)/ s

L 是规范下限,U 是规范上限,s 是样本标准差。算出来的 QL 与一个临界值 k 比较,QL ≥ k 就接收。

k 从哪来?查表。GB/T 6378.1、ISO 3951-1、ANSI/ASQ Z1.9,都有对应的表。给定批量、检验水平、AQL,查出一个 (n, k)。

讲义里有一个例子,值得反复看:

某原料含量要求 ≥96.0%,检验水平 II,AQL=1%,批量 100 桶。查表得样本量字码 F,方案(10, 1.71)。
抽 10 个样:96.95、96.04、96.68、96.63、96.65、96.52、96.59、95.86、96.57、96.91。
均值 96.540,标准差 0.343。
QL =(96.540 − 96.0)/ 0.343 = 1.573。
k = 1.71。QL < k,拒收。

注意这个结论的形状:10 个数据里,9 个看起来都很漂亮,只有一个 95.86 略微擦线。凭直觉,这批“差不多是合格的”。但方案判了拒收。

为什么?因为均值 96.54 离 96.0 只有 0.54,而批内标准差是 0.343——均值只比下限高 1.57 个标准差。按正态分布推算,这个批里大约有 5.8% 的单件会低于 96.0。方案不接受这个比例。

更能说明问题的是反过来:就算你把 10 个数全部做到 96.0 以上,只要离散度够大,照样拒收。

比如一组均值 96.97、标准差 1.171 的数据——每一个都在 96.0 以上,全部“合格”——QL 只有 0.83,离 1.71 差得很远,整批拒收。

这就是计量型和计数型的分野:

计数型问:有几个不合格?
计量型问:整个分布,离规范线有多远?

同样 10 个数据,计数型只能告诉你“有没有不合格”,计量型还能告诉你“离规范线几个标准差”。信息量不同,判别力自然不同。这也是为什么计量型抽样通常只需要计数型几分之一的样本量。

但代价也很实在,三件事必须做到:

一,正态性。计量型抽样默认数据服从正态分布。含量、水分、熔点通常成立;有关物质这种严重右偏的数据,要先做变换,或者直接用计数型。

二,σ 的维护。如果用“已知 σ”法,这个 σ 必须是长期数据累积出来的,不是这一批算出来的。它需要定期更新(常见做法是每 10 批重算一次),并且要和一个上限比较:MPSD(最大过程标准差)。若 σ 已经大于 MPSD,说明这个过程本身就不可能稳定合格,此时不必抽样,直接判供应商整改。

三,双规范限的处理。含量 98.0%~102.0%、熔点 216.0~217.5℃ 这类双侧限,要先明确是分立控制、联合控制还是综合控制——只控上侧、只控下侧、还是控总和。三种算法的判定结果可能不同,SOP 里必须写清楚。

多标准的情况更复杂。一个原料可能同时有水分、熔点、有关物质 C、有关物质 D 四个标准,分别属于 A 类(关键)和 B 类。这时不能逐个判完了事,要算总的不合格率:

P = 1 −(1−P₁)(1−P₂)(1−P₃)……

然后要求 P 不超过该类别的限值 M。两个类别都满足,才可接收。

七、回到那张 SOP:三层楼

讲到这里,可以把整件事重新排一下。

取样的三层结构

第一层件数√n+1 可以留

身份是“取样件数”,管代表性,不是判定规则

第二层方案必须有 (n, c) 或 (n, k)

计数型 GB/T 2828.1;计量型 GB/T 6378.1、Z1.9

第三层风险LQL、β、OC 曲线

写给监管看,也写给自己看

一句话:件数管代表性,方案管判定,风险管说服力。

还有一件事常被忽略,它横跨第一层和第二层:样品怎么处理。

同样是取 18 桶,你可以每桶各取 10g 混成一份混合样,也可以每桶单独测。混合样给你的是“平均”,它把批内变异抹平了——而计量型判定恰恰需要 s,需要批内变异。

如果你永远只测混合样,你就永远得不到 s,也就永远做不了计量型判定。

所以“取几件”和“怎么处理这几件”,必须一起设计。只改前半句,等于没改。

八、五个可以马上做的动作

一,把取样 SOP 拆成两张表。

一张写“取几件”(代表性),一张写“怎么判”(方案与风险)。现在很多 SOP 把两件事揉在一段话里,结果谁也说不清自己在执行哪一条。

二,给关键物料补一份 LQL/β 论证。

不用很长。一页纸:这个物料的 LQL 定在多少,理由是什么;β 定为多少;对应的 (n, c) 是多少;OC 曲线附上。审计的时候,这一页比十页“按规定执行”有用。

三,每个关键取样方案,附一张 OC 曲线。

Minitab 三分钟能画出来。图一贴,α 和 β 都在上面,不用解释。

四,包材继续用 YBB 加 GB/T 2828.1,不要乱升级。

药包材标准里的取样,是按检验项目分级给 AQL 的(外观 4.0、开启力 4.0、凸边 2.5……),这套体系是合理的:它的风险模型是外观与尺寸,消费者是制剂厂,AQL 的思路对得上。真正需要改的,是无菌、微生物、关键原料这类以患者为直接风险终点的取样。

五,计量型项目优先,建立自己的 σ 库。

含量、水分、熔点、装量——这些项目本来就有连续数据,用计数型判是浪费信息。先把历史数据汇总出 σ,再谈 MPSD,再谈计量型方案。这件事做起来慢,但一旦建成,样本量能降下来,判别力反而上去。

√n+1 不是错的。

它只是被问了一个它回答不了的问题。

这么多年来,它承担了“取样件数”的职责,做得不算差——次线性、成本可控、操作简单。但当它被当成“抽样方案”,被用来支撑“整批合格”这个结论时,它就开始说一些它没有资格说的话。

取几件,是成本问题。
收不收,是风险问题。

把这两个问题分开,取样才从一门手艺,变成一门科学。

你们厂的取样 SOP,写的是件数,还是方案?

来源:PharmaBuddy · mp.weixin.qq.com