《数据分析在新版GMP中的灵活应用》由会员分享，可在线阅读，更多相关《数据分析在新版GMP中的灵活应用（168页珍藏版）》请在金锄头文库上搜索。
1、 数据分析 在新版在新版GMPGMP中的中的 灵活应用 苏勤 2010.07 研讨提纲一.新版GMP实施箭在弦上,药企有备才能无患二.数据分析在企业内控标准制订中的应用三.数据分析在确认与验证中的应用四.数据分析在生产控制中的应用五. 数据分析在质量控制和质量保证中的应用 1. 在质量控制管理中的应用 2.在产品放行中的应用 3.在稳定性考察中的应用 4.在偏差管理中的应用 5.在对供应商的审核中的应用 6.在产品质量回顾分析中的应用新版GMP实施箭在弦上 药企有备才能无患 我国新版GMP实施在即，从法规层面上保障了万千百姓的安全用药，对广大中国药品生产企业来说，是一种挑战，更是一个难得的机遇
2、;是企业淘汰落后，实现升级换代，在更高平台上展翅高飞的起点。 建立与国际标准相当的GMP，有助于促进国家间药品GMP检查结果的互认。我国药企只有步入国际标准的轨道，才有无限延伸的市场空间。 我国新版GMP即将颁布。从内容看，新版GMP借鉴了国际先进标准和药品监管经验，更加注重科学性，其中着重细化了软件要求，在硬件上也提出了较高的技术要求。 新版GMP引入受权人 、质量风险管理、变更、偏差、产品的年度审核以及纠正预防措施等概念的同时，全篇贯穿了生产工艺同注册工艺相一致的要求 ， 这是“质量源于设计”理念的体现。 QbD(QualityQbD(Quality By By Design) :Desi
3、gn) :质量源于设计 - 理念的推行l质量是检验出来的 产生了标准l质量是生产出来的 产生了GMP和ISO标准l质量是设计出来的 催生了风险管理lFDA官方的推行标志着QbD从概念走上了法规要求 一切离不开数字，数字是美妙的。这对在数字和公式面前畏缩的人来说，这听起来也许是荒谬的。但看自然界，从最小的贝壳到最大的漩涡星系，无不和宇宙中大部分美丽的东西都呈现出数学特性.古代最伟大的哲学家之一亚里斯多德曾说过: 数学格外地展现了秩序、对称和极限，而这些是美的极致形式。 数美丽的天性， 也赋予我们在制药生产中的管理艺术. 善待数据 珍惜管理的无形资源 掌握现状， 工序调节， 工序管理， 检查和评价
4、， 分析和改进。不要轻信数据数据的来源：历史的，现在的虚假的数据 数据与事实不符 人为的虚假数据 经过修改的数据 因无知造成的错误数据错误的数据 抽取数 据 的 方法不好 数据的抄写、计算错误 异常值取不到的数据妥当安排提提有 科 学 依 据数据优良恰如其分 介 绍FDA检查要求法则2005年12月FDA CGMP对中国培训中提到，不可接受或错误的数据和记录. 不诚实数据的特征: 不真实，编造，虚假，无系统. 被认为是审查程序的材料之一的递呈书中重要数据缺失.应该让上报的数据未被上报. 不准确(例:第一次数据未达标，重测的数据达标准，但未经实验室调查，将此数据上报用于申请中.)要坚持诚信诚信是
5、立人之本诚信是交友之道诚信是为政之基诚信是经商之魂诚信是心灵良药诚信是中华民族的传统美德算:特征数比:统计值找:相关因素看:动态变化 数据分析与 新版GMP管理 数据分析在企业内控标准制订中的应用 数据分析在确认与验证中的应用 数据分析在生产控制中的应用 数据分析在质量控制管理中的应用 数据分析在放行中的应用 数据分析在稳定性考察中的应用 数据分析在偏差管理中的应用 数据分析在对供应商的审核中的应用 数据分析在产品质量回顾分析中的应用 数据分析在企业内控标准 -制订中的应用 C CPP值值值值数据特数据特征值征值可运用的数学模型控制图数据的数据的筛选筛选 数据的筛选 在实际测试中常常会发现有几
6、个测量值很接近，而有一个离得很远，表面看它已超出界限了。这个 值可能是原分布中的一个，是有效的;也可能是未发现的误差而引入的应弃去。 目前广泛采用Q法和格鲁布斯法,现再介绍一法,按美国EJ鲍尔推荐的方法进行处理。 步骤如下： 1.计算这群检测值的平均值 2.计算极差R 3.计算可疑值Xi与平均值 之差的绝对值，再用极差R除，得出ti ，与规定附表的临界值比较，若ti比表上的t值大，则应弃去此可疑值。 注：用公式注：用公式t=|Xt=|Xi i- |/R- |/R计算计算t t，如计算值超过表上的值时，则所调查的值是无效的。此概率约为0.950.95。例例1 1某分析者对一样品检测，得，得（1
7、1 1 1）93.393.393.393.3，（2 2 2 2）93.393.393.393.3，（，（，（，（3 3 3 3）93.493.493.493.4，（，（，（，（4 4 4 4）93.493.493.493.4，（，（，（，（5 5 5 5）93.393.393.393.3，（，（，（，（6 6 6 6）94.094.094.094.0。 问：第六个结果有效吗？n3456789t1.531.050.860.760.690.640.60n10111213141520t0.580.560.540.520.510.500.46附表如下：附表如下：抛弃无效测量的抛弃无效测量的临界值临界值
8、临界值临界值 解：1.计算方法六个结果的平均值93.45 2.计算极差R94.093.30.7 3.计算可疑值与平均值之差的绝 对值再用极差除：ti i|X |R （94.0-93.45）0.7079 4.与临界值t=0.76(n6)比较 ti i0.79 t=0.76 94.0是一个离群数据，应弃去。 数据特征值 位置特征量: 子样平均值 子样中位数 差异特征量: 极差 R=XmaxmaxmaxmaxX Xminminminmin 标准差 变异系数 （相对标准差（相对标准差RSDRSD）平均值的二种计算算术平均值加权平均值产品质量性能产品质量性能比较比较生产工艺生产工艺、环境环境状况比较状况
9、比较操作水平操作水平比较比较数据数据特性值特性值设备能力比较数据特性值的应用 用控制图的思路 制订企业内控标准 不论原辅料还是中间体及产品的质量的企业内控标准，从所使用的原辅料质量特性的检测数据、中间体质量特性的检测数据、生产的产品质量特性的检测数据，取连续20个数，弃去异常点，作控制图，计算平均值，标准差S与控制线UCL、CL、LCL的值，即可确定此标准质量特性值的范围为 3S。时间或组号3倍标准偏差3倍标准偏差质量特性值UCLCLLCL控制图特点样本的平均值 为 中心线 CL +3s 为上控制线 UCL -3s 为下控制线 LCL 区间概率为99.73% 用工序能力指数(Cp值) ) -
10、-验证所制订标准的可行性工序能力（B） 含义：在一定时间内工序处于稳定状态下的实际加工能力。 测定：B= 6S 1.是工序的一种可以量度的特性； 2. 反映与公差无关； 3.在正态分布情况下，6范围内的概率 为99.73%。工序能力指数Cp含义：产品公差范围（T）与工序能力（B） 之比。符号：样本均值；S样本标准差；TU规范上限；TL规范下限；T规范范围M 规范范围中心值中心偏移量k 偏移系数工序能力判断准则 工序能力状态工序能力状态工序能力指数工序能力指数严重不足严重不足C CP P(C(CPKPK、C CPUPU、C CPLPL) ) 0.670.67不足不足0.67 C0.67 CP P
11、(C(CPKPK、C CPUPU、C CPLPL) ) 1 1正常正常1 C1 CP P(C(CPKPK、C CPUPU、C CPL PL ) ) 1.331.33理想理想1.33 C1.33 CP P(C(CPKPK、C CPUPU、C CPLPL) 1.67) 1.67过剩过剩C CP P(C(CPKPK、C CPUPU、C CPLPL) )1.671.67备注：收集观测值个数备注：收集观测值个数n30 n30 工序能力指数对应的不合格品率 C CP P不合格品率不合格品率C CP P不合格品率不合格品率1.671.676/10006/1000万万1.11.11/10001/10001.5
12、1.57/1007/100万万1 13/10003/10001.331.336/106/10万万0.670.674.55/1004.55/1001.21.23/13/1万万0.330.3331.75/10031.75/100 制订生产过程的中间体标准 在成熟的工艺条件下，收集近20个特性观测值，统计其平均值、极差与标准偏差与产品法定标准规定对特性值比较，以极差及标准偏差来考虑工艺波动范围，综合比较来确定其质量特性值。 例例2 2制定抗生素BL溶配岗位中间体（乙醇溶液中抗生素B的效价）控制标准。收集某月正常生产的连续20批数据并进行数据筛选如下:序号序号检测值检测值序号序号检测值检测值11663
13、52111623992166944121694273159847131612114166623141647465167098151686666165963161607127160320171664878161123181709239161652191722541016676520172604 对序号3，18，19，20的数筛选结果舍去X19与X20(略) 计算:平均值 =164848 标准偏差S=3388 根据控制图原理计算 CL = 164848 UCL= +3S=175011 LCL= -3S=154684 因此，此工艺点的内控标准初步定为154684175011，考虑到测定方法的误差，故
14、将乙醇溶液中抗生素B的效价的内控标准定为155000175000 。 制订成品的企业内控标准 例例3 3 USP-29规定抗生素P中抗生素G含量（HPLC法）为51.0%-59.6%，现制定此产品的质量特性值的企业内控标准。 收集某年正常生产的连续12个月的数据（平均值）并进行筛选如下：序号序号检测值检测值序号序号检测值检测值155.855.8756.256.2255.955.9856.256.2356.156.1956.356.3456.456.41056.956.9556.156.11155.255.2655.455.41255.655.6 计算 Cpk= 2.41.67，工序能力过高。在
15、保证产品质量的前提下，为减少生产成本，将生产过程中的工序能力适当放宽，将计算出的控制范围进行调整为54.058.0%，作为此产品抗生素G含量的企业内控标准。质量标准管理的 关健在 取样样本量与标准差的关系从图可知，当测定次数n5时，则减小的速度变慢，并随着测定次数的继续增加曲线趋于平坦。在实际工作中，通常进行5次左右的测定就可以了，最多不超过10次左右。 从图可知，当测定次数n5时，则 减小的速度变慢，并随着测定次数的继续增加曲线趋于平坦。在实际工作中，通常进行5次左右的测定就可以了，最多不超过10次左右。 药品的检验是破坏性的检验，决定了样品的抽样形式，抽样过程中存在的弃真()的的第一类错误
16、和纳伪的第二类错误()，它们给企业带来了管理的风险。在管理中如何把握，看我们的水平发挥。数据分析数据分析在确认与验证中的应用在确认与验证中的应用数字的数字的数字的特性值特性值特性值区间估计区间估计区间估计假设假设假设检验检验检验 可运用的数学模型散点图 GMP的概念 确认 qualification ：证明任何厂房、设备及其配套系统正确安装之后，能够正常运行，并且得到预期结果的有记录的活动。 验证 validation: 能够高度确保某个工艺、方法或系统能够持续地产生符合规定标准的结果的有记录的活动。 “验证”的词义有时可以扩展，使之包括“确认” 的概念的概念。1.验证是一个证实或确认设计的过
17、程；2.验证是一个确立文件的过程；3.3.验证是一个提前发现问题的过程。 （验证时可安排最差条件试验、极限试验和挑战性试验） 美国FDA分析方法验证指南 要求:根据实际使用各种确 认药典分析方法的适用性。 ( ( 21CFR211 . 194 (a) (2 )(a) (2 ) ) )指出: 需要对分析方法验证资料进行统计分析(比如 : 线性回归分析，相对标准偏差)以说明方法的正确性。在开始分析方法验证之前，应当就要确定用于验证资料分析用的统计方法。 FDA 对颗粒含量均匀度的判定 至少10个位置点， 每个位置一次取3个样 品 ， 测定其中一个样品结果表示为标示量的 百分含量 ， 每批含量的相对
18、标准偏差满足RSD5%各批单点含量与平均值的绝对差10%则判断终混颗粒通过此判定。片剂、胶囊剂的含量均匀度判定 一.FDA Guidance 的标准 1.选取20个取样点，每点至少抽取7粒，分析其中3粒，共60粒；2.全部药粒含量（折重，每批共60粒）的RSD，RSD4.0%；3.每个取样点的含量平均值在90%-110% (表示为含量相当于标示量的）；4.每一个单个结果应在75%-125%之间。 二.Bergums 片剂、胶囊含量均匀度判定1.RSD上限标准表（90%置信度下，未来样品至少有95%的概率通过USP含量均匀度的2步判定）。2.计算样品的RSD，如其小于表中查得的RSD标准，则通过
19、。3.还必须满足药典、法规规定的超出限度的样品数量。 J. S. Bergum and M. L. Utter， “Statistical Methods for Uniformity and Dissolution Testing，” in Pharmaceutical Process Validation， 3rd ed.， Eds. R. Nash， R. and A. Wachter (Marcel Dekker， New York， 2003).FDAFDA色色谱方法验证审评指南指出:均数来估计样品分析准确度，RSD是样品分析的精密度。方法验证中，建议最少10次进样，RSD1%，对于
20、药物或药物制剂中的活性成分出厂和稳定性研究，其RSD1%，系统适用性试验的精密度RSD至少5次进样。对低浓度杂质，较高的变异系数可以接受。 假设检验假设检验 基本思想 根据是小概率原理，所谓小概率原理：即“在一次试验中，概率很小的事件认为是实际上不可能发生的事件，而概率接近于1的事件认为是实际上必然发生的事件。” 我们可以根据所获样本，运用统计分析方法，对总体X的某种假设H0 0作出接受或拒绝的判断。 基本步骤1.建立假设 H0：2.选择检验统计量，给出拒绝的形式3.给出显著性水平4.由样本值计算出统计量的值5.根据统计量的不同分布查相应的分布表6.作出判断 当则拒绝H0则接受H0 基本类型
21、1. 检验母体平均数 2. 检验二个母体平均数相等 3. 检验母体方差 4.检验二个母体方差相等 自由度dfn1+n22n1+n22t检验用于新方法可用性.例例3用含氮量为10.0 %的标准物来检验氮分析新方法可用性1.H0:0 H1 : 0 0 10.0%2.选用t检验 t(0 )(S )3.选0.05 ，df=3(n-1)， 查表，则t 0.052.354.用新方法测定4次，所得到的平均值为9.8%，标 准差是0.2% 5.按公式计算，得t(0 )(S ) (10.09.8)(0.2 )2.0 6.比较 2.02.35 接受H0 即新方法可用. 数据特征值发出可信的技术报告 中国药典中药品
22、质量标准分析方法验证指导原则对数据要求:在规定范围内，至少有9个测定结果进行评价.例如，设计3个不同浓度，每个浓度各分别制备3份供试品溶液，进行测定，应报告已知加入量的回收率()或测定结果平均值与真实值之差及其相对标准偏差或可信限。例例4 4含量测试方法验证报后要:a)记录标准溶液中被测物峰面积，计算平均值与RSD%。b)记算在每份溶液的含量及相应的RSD%。 区间估计验收器具例例5 5XXXX丁基胶塞首模质量的验收数据19.4419.5219.5219.5019.5019.4819.4619.4819.5219.5019.5219.4419.4419.4419.5419.4419.4619.
23、4619.4419.4419.5019.4819.4419.5019.4419.4619.4419.4819.4419.5219.5019.4819.4619.4819.4819.4619.4419.4619.4619.5219.4819.4819.4819.4819.4419.5219.4619.5019.5219.4819.5019.5419.5219.5019.4819.4619.5019.5219.4819.5019.4819.5219.5219.5419.4619.4819.4219.5219.5419.4819.5019.4819.4819.4619.5219.4419.5019.
24、4219.5019.5219.5219.4819.5019.4419.4819.4619.5219.5219.4419.4419.4619.5219.5219.46 计算19.48，S0.032RXmaxXmin19.5419.420.12冠部直径标准为19.50.2又3S19.4830.03219.480.096即区间为(19.3819.58)在标准区间(19.319.7)中此件模俱验收合格，可在生产中用。380392生长周期（生长周期（h）70006700摇摇瓶瓶效效价价u/g3866400例例6 6 XXXX斜斜面面培培养养周期由其外周期由其外 观观及及生生产产上上反反馈馈来来确确定定的
25、的, ,根根据据0808年年4 4月月-09-09年年1 1月月的的5656批批斜斜 面面的的生生长长周周期期和和摇摇瓶瓶效效价价，作作左左图图，结结果果两两者者不不相关。相关。 斜面生长周期不合适是要因 结论结论生长周期与摇瓶效价相关图生长周期与摇瓶效价相关图散布图确定斜面生长周期不合适t检验确认设备工作状态例例7 7某药片标准重为100毫克，抽25片分别称重，检查压片机工作是否正常？计算得样本均数为98.4毫克，标准偏差为4.5毫克。 = 98.4 S= 4.5 0 =100 n=251.1.检验假设检验假设：H0: = 0 H1 :02.2.计算统计量计算统计量： 3. 3. 确定概率确
26、定概率：=0.05， f = 25-1=24 查t值表，t0.05(24)=1.711 ， 1.778 1.711 ， 4.4.判断结果判断结果：样本均数98.4毫克与总体均数100毫克有显著性差别，压片机工作不正常。新版GMP230条要求，符合下列情况之一的，应对检验方法进行验证: : (1)采用新的检验方法; (2)检验方法需变更的; (3)采用中华人民共和国药典及其它法定标准 未收载消息的检验方法; (4)法规规定的其它需要验证的检验方法。对不需要进行验证的检验方法，企业应对检验方法进行确认，以确保检验数据准确、可靠。U检验确认二种方法的准确度例8 8用新方法与经典方法对比.测定一样品的
27、某特性，得测定数，试问新方法能否代经典方法.经典法经典法18.8919.2019.0019.7019.40-新方法新方法20.1020.5018.6519.2519.4019.99新方法能代替经典方法. 假设检验的联用 在两个总体均服从正态分布前提下，先用夏皮罗威尔克方法检验正态性，再用F检验，看两个正态总体方差是否齐性，在方差无显著性差异时，最后用T检验对两个总体均值进行比较，此法联用可广泛用于工艺改变前后产品质量的变化，验证二个检验方法对同一产品检验的特性值是否有差异。 假设检验确认二种方法的同等应用例9 9 XXX抗生素含量检验中用旋光法 代电位法的确认.批号E2E3E29E30E31E
28、34E43E452E528E529旋光法99.9199.5099.5510099.3799.9599.9199.6899.6499.70电位法99.8799.6999.7599.6299.5899.7699.6499.7099.8899.58批号E530E532E5807-E444Na37Na40Na42Na62Na64Na65旋光法99.6099.6099.699.6999.751001599.7599.75100.199.88电位法99.7099.5299.6499.94100.099.8799.9499.9499.9499.941.用S-W检验法检验二样本的正态性 经计算机统计后，数据为
29、=0.01， W(20， 0.01%)=0.868 W旋=9.716831282E-01 W电=9.121425475E-01 经比较数据后，二样本均值呈正态性 2.用F检验法检验两总体方差是否相等 0.75 0.76 =0.05 又F= 1.782.49 故两总体方差无显著性差异3.用T检验法检验两总体均值是否相等 查表 0.1772.024 两总体均值无显著性差异. 故旋光法可代替电位法检测.测定数据的收集和处理 统计验证样品B成分测试数据表样品编号第一批第二批第三批140.3740.4040.40240.1940.3940.49340.1840.5340.53440.2540.6240.
30、32540.2740.2740.37640.3640.2040.20740.4140.5340.53840.4240.5140.51940.4740.2740.271040.2940.5240.5240.3240.4240.41s0.09970.14020.1205u0.2850.1260.142 假设检验作产品均一性验证的分析例例1010产品的均一性分析1) 1) 抗生素P产品中，某化合物B含量的测定得到上述数据表。在长期生产过程中，可知化合物B的含量服从正态分布N（40.50，2.0），取=0.05，查正态分布函数的临界值为K0.05=1.96因为因为u u1 1=0.285 =0.285
31、 K K0.050.05=1.96=1.96 u u2 2=0.126 =0.126 K K0.050.05=1.96=1.96 u u3 3=0.142 =0.142 K K0.050.05=1.96=1.96 故接受原假设H成立，在给定的显著性水平 =0.05时，所取3批的混粉每批各是均匀的。 哈特利(Hartley)检验法 -多个总体方差的检验 假设 HO: =-= , m3的正态总体,独立取容量为K1=K2=Ki Ki3 , 计算 自由度df=K-1 查哈特利检验临界值H, 当HH时, 则HO成立。2 2) ) 又所取3批样品的样本为K1= K2= K3=10，假设H：12=22=32
32、进行统计检验，运用哈特利检验法进行计算H=1.977当m=3，f=k-1=10-1=9，=0.05查哈利特临界值表，得H0.05=5.34比较H=1.977 H0.05=5.34故在给定的显著性水平=0.05时，所取3批混粉后的产品的B含量的标准差没有显著性差异。又当m=3，f=k-1=10-1=9，=0.01查哈利特临界值表，得H0.01=8.5比较H=1.97 H0.01=8.5故在给定的显著性水平=0.01时，所取3批混粉后的产品的B含量的标准差没有显著性差异。根据对抗生素P的3批抽样测定B成分含量，可知抗生素P产品质量是均一的，可见抗生素P的混粉工艺是稳定可靠的。回归分析正交试验控制图
33、假设检验可运用的数学模型数据分析在生产控制中的应用正交试验设计 正交设计是用于安排多因素试验的一种科学方法.在多因素的情况下，选择最佳条件的比较方便方法.它适用于多因素的工艺参数的优化，我们将其应用于研究工作、生产试验的设计，极大地提高了试验的可靠性和效益。 常用正交表L4(23 )，L8(27 )，L9(34)，L8(41 24 ) 原理:利用数理统计的观点，运用正交性原理，从大量的试验点中，挑选适量的具有代表性典型性的试验点，用正交表来合理安排试验的一种方法. 用途:正交试验就是利用正交表的特点进行研究试验，找出最佳工艺条件为生产奠定基础.常用正交表常用正交表: : 正交试验设计步骤: 1
34、.1.明确试验目的，确定考核指标，明确试验目的，确定考核指标， 2.2.根据试验目的，挑选出试验的因素，根据试验目的，挑选出试验的因素， 3. 3. 选定每个因素变化的水平选定每个因素变化的水平 4. 4. 选择最合适的选择最合适的L L表表 5. 5. 根据因素和水平设计出试验方案表根据因素和水平设计出试验方案表 6. 6. 进行试验并记录好数据进行试验并记录好数据 7. 7. 设计计算表并计算、分析设计计算表并计算、分析 8. 8. 找到最佳搭配找到最佳搭配 用正交试验寻找最佳工艺点例例1111为了找到某产品中试生产中，影响产品收率的主要因素，从而为大生产确定最佳操作条件，我们又对结晶工序
35、的关键控制点进行了正交试验，确定了DBED溶液预加量（L）、结晶温度（）、结晶搅拌速度（r/min）这三因素，确定了各自的3个水平，如下表： 因素因素水平水平1 1水平水平2 2水平水平3 3A (DBEDA (DBED预加量，预加量，L)L)4.04.04.54.55.05.0B (B (结晶温度，结晶温度，) )5 515152525C(C(结晶搅拌速度结晶搅拌速度 r/min)r/min)100100150150200200选用选用3 3因素因素3 3水平所用的正交表水平所用的正交表L L9 9（3 34 4），安排试验，），安排试验，并得产品收率数据，分析如下表：并得产品收率数据，分析
36、如下表： 由上述图可见：DBED溶液最佳预加量为A2即4.5L；结晶温度对收率的影响不明显，为节约能源，选用B3即25比较好；最佳结晶搅拌速度为C2即150r/min，所以确定的最佳工艺条件为：A2 B3 C2 即DBED预加量为4.5L，结晶温度室温，最佳结晶搅拌速度为150r/min。 控制图控制图是对过程质量特性值的数据进行测定、记录、评估，从而监察过程是否处于控制状态的一种用统计方法设计的图.作用：对生产过程进行监控 发现异常，及时告警控制图的种类: 分析用控制图 控制用控制图 其结构:中心线 CL 上控制线 UCL 3 下控制线 LCL 3 控制图 的原理质量由人机料法环五个因素，它
37、又可分偶然因素(简称偶因)与异常因素(简称异因)。偶因是始终存在的，对质量的波动影响微小，但难以除去.偶波是不可避免的，但对质量波动影响微小，异波对质量波动的影响大，且采取措施不难消除，故在生产过程中造成异波的异因是需监控的对象，一旦发生，应尽快 采取措施加以消除，当生产过程中只存在偶波时，产品质量将形成典型分布，如果除了偶波还有异波，产品质量的分布必将偏离原来的典型分布。根据典型分布是否偏离就能判断异波即异因是否发生，而典型分布的偏离可由控制图检出，控制图的实质是区分偶因与异因，控制界限就是区分偶波与异波的科学界限。 作图步骤1.数据表 记入必要事项与数据。2.控制图 纵坐标为样本质量特性，
38、横坐标 为样本号；平均值为中心线，正负三倍标准偏差处为上下控制线；注明相关内容和数据。正常正常 控制图的判断 必须同时满足以下两条: 1. 所有点子在上下控制线内； 2. 线内点子随机排列无缺陷。 控制图的异常判定准则 点出界 界内点排列不随机 一点落在A区外 连续9点落在中心线同一侧 连续6点递增或递减 连续14点中相邻点上下交替 连续3点中有2点落在中心线同一侧的B区以外 连续5点中有4点落在中心线同一侧的C区以外 连续18点在中心线两侧，但无一在C区中. 注意问题1.确定控制对象（质量指标）应有数值。2.一张控制图分析数值只是质量特性值中的一项。3.控制图能提前起报警作用，仍要分析异常原
39、因。4.条件变化或长时间后应重作图分析。5.控制图中分析的质量标准宜采用内控标准。6.不能用规格界限来代替控制界限。7.控制图分析全过程的记录应归档保管。 控制图呈现生产的脉搏例例1212在在产品分装工序中，为保证成品装量处于产品分装工序中，为保证成品装量处于 控制状态，规定每班每台机随机抽样控制状态，规定每班每台机随机抽样4 4次，每次次，每次抽抽8 8支，测装量值，质量要求工艺控制范围为支，测装量值，质量要求工艺控制范围为10101010 914mg914mg，生产现场作出下图。，生产现场作出下图。 从图分析，从图分析，本机装量差本机装量差异异100%100%合格，合格，工序稳定。工序稳定
40、。 新版GMP提出的限度线 行动线： action lever 对于超标结果立即采取纠正措施（如：洁净区微生物和微粒监测超标。） 警戒线 ：Alert lever 比行动线低，超过后，进行调查，不一定采取纠正措施，有可能采取预防措施。(如：洁净区微生物和微粒监测；培养基罐装试验。) 新版GMP237条：超标(OOS)结果调查 用管理限报警关键的控制要素在企业质量控制中的警戒限和行动限，如:对工艺用水与洁净区的环境在控制中根据正常工序情况，制定出可行的指标，即警戒限2S，超出此界限意味着发出了报警。行动限：3S超出此界限时，质量管理上要立即采取行动。例例13-1洁净区环境的控制(1)-菌落数洁净
41、级别测试状状态表面菌落数表面菌落数CFU/24-30cm2标标标标准准准准行行行行动动限限限限警戒限警戒限警戒限警戒限ISOISOEGMPEGMPFDAFDA5A100级静态测试111动态测试7B与100级相连10000级静态测试542动态测试与100级不相连10000级静态测试554动态测试8C100000级静态测试252018动态测试9D-动态测试504035例例13-213-2洁净区环境的控制(2)-沉降菌沉降菌洁净级别测试状状态沉降菌沉降菌 个个/ /（90mm 4h90mm 4h）沉降菌沉降菌计数数方法方法标准标准标准标准行行行行动动限限限限报报警限警限警限警限ISOISOEGMPE
42、GMPFDAFDA5A100级静态测试111取最高菌落数7B与100级相连10000级静态测试542动态测试与100级不相连10000级静态测试553动态测试8C100000级静态测试504035动态测试9D- 动态测试1008070例例13-313-3洁净区环境的浮游菌控制(3)-浮游菌浮游菌洁净级别测试状状态浮游菌（浮游菌（CFU/mCFU/m3 3）取取样体体积标准标准标准标准行行行行动动限限限限报报警限警限警限警限ISOISOEGMPEGMPFDAFDA5A100级静态测试1111000L7B与100级相连10000级静态测试10821000L动态测试与100级不相连10000级静态测
43、试101071000L动态测试8C100000级静态测试1008070100L动态测试9D-动态测试200160140100L 例例13-413-4洁净区环境的控制(4)-尘埃尘埃洁净度度级别测试状状态尘埃粒子数埃粒子数 个个/m/m3 30.5m0.5m标标准准准准0.5m0.5m行行行行动动限限限限0.5m0.5m报报警限警限警限警限5m5m标标准准准准5m5m行行行行动动限限限限5m5m报报警限警限警限警限ISOISOFDAFDAEU EU GMPGMP5100级A静态352035002800202016动态352035002800202016710000级B静态352035002800
44、292924动态3520003500002800002900200016008100000级C静态352000350000280000290020001600动态3520000350000028000002900020000160009D静态352000035000002800000290002000016000动态-例例13-513-5工艺用水的水质控制水质控制纯化水质量控制标准纯化水质量控制标准项项 目目标标 准准报警限报警限行动限行动限电导率(25) s/cm1.31.0 1.1总有机碳(TOC) ppb500 250 400微生物总数 CFU/mL 80 15 60细菌内毒素 EU/m
45、L0.25 0.1250.125硝酸盐 ppm0.2 0.20.2 回归分析 研究变量间的定量关系式时称为回归分析. 直线回归:通过试验可以测定一个变数对另一个变数的影响 在数学上可写成Yf(X) ， 如果二者的关系是直线关系，则此式就呈YabX，这时叫直线回归. 回归分析的主要内容:从一组数据出发确定这些变量之间的定量关系式对这些定量关系式进行统计检验从影响着某一个量的许多变量中，判断哪些量的影响是显著的，哪些是不显著的利用所求得的关系式对生产过程进行予报和控制.判断两个变量有无相关关系的方法: 简易法 相关系数法 回归分析法找出影响产品质量的主因回归分析法找出影响产品质量的主因每种产品都有
46、几个质量特性值，它们间是互相影响并互相制约的、找出影响产品质量的主因，就可以抓住主要矛盾，在生产过程中加以控制，使产品质量保持稳定。运用电子计算机为工具，以多元线性回归法、求出因变量Y与自变量Xi(i=1、2n)的方程Y=b0+b1x1 + b2x2+ bnxn，用方差分析检验回归方程，当FF时，方程效果显著，求出F1、F2Fn比较，按F值大小，即可找出影响Y的因素顺序。回归分析找到影响质量的主因例例1414找出影响青霉素钠含量的主要因素 现收集青霉素钠35批的有关质量特性值数据(略)，根据经验，该产品含量与pH值、264nm吸收值、280 nm吸收值之间关系密切，为此进行验证，并分析各因素影
47、响的程度。经计算机计算，可得：y=70.03223-1.003950X1+43.112129X2-51.44648X3 上述数据经计算机处理后，可得： 复相关系数R=0.7761；剩余标准差S=0.3098607； 回归平方和U=4.507692 残差平方和Q=2.976423； 方差比 F=15.64948； F1=11.43947； F2=31.19659； F3=13.87310。 当=0.01时，F(0.01)=4.51；因F=15.64948F(0.01)=4.51； 故此回归方程效果显著。 各个自变量Xi显著性的检验 当=0.01时，F(0.01)=7.56； 因 F1=11.43
48、947F(0.01)=7.56； F2=31.19659F(0.01)=7.56； F3=13.87310F(0.01)=7.56； 故X1、X2、X3对Y均有显著影响。 又31.1965913.8731011.43947，故 F2F3 F1 可知青霉素钠内在质量中pH值、264nm吸收值、 280 nm吸收值对其含量均有显著影响，按影响程 度，第一位的是264nm吸收值，其次是280nm吸收 值，第三是pH值。生产工艺中要严格控制对应的有 关工艺操作点，以求产品的高含量与高纯度。 线性回归寻找工艺控制条件例例1515 寻找产品的工艺控制条件在产品的实验室试验中，DBED溶液流加速度对结晶粒度
49、影响较大，故用线性回归定量分析。测试数据见下表： 回归方程为：y=143.27-9.91x相关系数为：-0.94从上图可知：DBED溶液流加速度与结晶粒度呈线形关系，以此为依据，可以确定不同结晶粒度所需不同的DBED溶液流加速度。 假设检验确认配方改变的可能性例1616某橡胶配方，原用氧化锌5克，现减为1克，现对两种配方，测量橡胶的伸长力如下：原原配配方方540533525520545531541529534S=7.99新新配配方方565577580575556542560532570561S=15.39问两种配方效果是否相同解：假设H0：12=22， H1： F=S22/S12=236.8/
50、63.86=3.7，取=0.10，根据自由度9和8查F分布，得临界值F0.05=3.39。因为F=3.73.39，所以H0不成立，两总体的标准差不等，所以不用再比较均值。 数据分析在质量控制中的应用 数据筛选数据修约数据特征值假设检验线性回归控制图可运用的数学模型数字的修约4 4舍，舍，6 6入入 3.54243.5 2.36532.43.54243.5 2.36532.4遇5要考虑，5后非 0 应进1，5后皆0看奇偶，5前偶数应舍去，5前奇数则进1。2.15072.22.15072.2不得连续修约。15.456415正确15.456415.45615.4615.5 15 错 新版中国药典的检
51、验方法和限度 本版药典中规定的各种纯度和限度数值以及制剂的重（装）量差异，系包括上限和下限两个数值本身及中间数值。规定的这些数值不论是百分数还是绝对数字，其最后一位数字都是有效位。 试验结果在运算过程中，可比规定的有效数字多保留一位数，而后根据有效数字的修约规则进舍至规定有效位。计算所得的最后数值或测定读数值均可按修约规则进舍至规定的有效位，取此数值与标准中规定的限度数值比较，以判断是否符合规定的限度。 数值修约规则应按照现行版国家标准进行。 (GB/T8170-2008)化学药物质量控制分析方法验证技术指导原则【H】GPH 3-1要求: 杂质检测结果应提供具体的试验数据杂质含量小于1.0%，
52、则报告应精确到小数点后第二位，如杂质含量大于1.0%，则报告的数据可精确到小数点后第一位。 数据特征值确立分析结果 在一般的分析测试中，对一个试样平等地测34次，出报告时，先算出平均值 ，再算出结果的平均偏差d，相对平均偏差D，当D0.2时，可取 出数据，否则应重试.例例1717测一溶液浓度，三次测定结果为0.2041， 0.2039，0.2043 。计算 0.2041; d0.0001; D0.05 D0.050.2 可用0.2041，molL报告分析结果。 若测定在4次以上，则要对总体均值的可能存在区间作出估计，要进行统计处理。 数据特征值表示检验误差欧州药典对容量分析滴定要求至少平行 滴
53、定6份，对其分析的确定规定如下:容量滴定容量滴定重复性重复性(RSD)相相对对准确度准确度酸碱滴定酸碱滴定0.330.67非水滴定非水滴定0.330.67碱的共碱的共轭轭酸酸0.330.67氧化氧化还还原法原法0.51.0银银量法量法0.51.0络络合滴定法合滴定法0.671.33 我国药品检验的精密度要求分析方法分析方法标标定或原料定或原料测测定定制制剂剂分析分析重量法重量法0.10.2中和法中和法0.20.5氧化氧化还还原法原法0.20.5综综合量法合量法0.20.5银银量法量法0.20.5比色法比色法2.02.0紫外分光光度法紫外分光光度法2.03.0气相色气相色谱谱法法-5.0液相色液
54、相色谱谱法法-5.0旋光法旋光法23.0微生物法微生物法-5.0 用数据特征值来考评操作者的水平例例1818对气相色谱的操作人员进行技术考核， 进样10次，每次0.5l，得色谱峰高为： 142.1 147.0 146.2 145.2 143.8 146.2 147.3 150.3 149.9 151.8（mm） =146.98 S=3.00 CV V=2.04% % 有经验的色谱工作人员很容易将CV V控制在1% 以内，可认为该实验人员的技术还不够稳定，操作不够熟练。 t检验估计特性值区间例例1919试样测CI-的含量为:47.6447.64% % ， 47.9647.96% % ， 47.5
55、247.52% % ，47.5547.55% % ， 当0.05时，求平均值的区间. 解: 47.67% S0.20% 当n10时， 且未知，用t统计量来求 0.05 ， dfn-1413 查表得t0.0253.18 47.670.318 有95%的把握，认为CI-的含量范围为 ( ( 47.3547.3547.99 47.99 % % ) ) F检验两种方法的精密度例例2020要确定两种方法是否具有同样的精密度.分析人员用每种方法对9份样品检测，得到数据，进行统计.1. H1. H0 0:2 21 12 22 2 H H1 1:2 21 12 22 2 2. 2. 选用F检验3. 3. 选0
56、.05 0.05 查表，则F F 0.05 0.053.443.444. 4. 计算 S S2 21 1 16 16 ， S S2 22 2 1010 按公式计算 F FS S2 2大大S S2 2小小16161010 1.6 5.比较比较 1.6 1.63.443.44 接受H H0 0 线性回归确立 标准比色液的 工作曲线例例21铂钴标准色度溶液标准 工作曲线 测定铂钴标准色度溶液的吸光度，用波长420nm，10cm比色皿，得以下数值:铂钴色度铂钴色度(Hazen)(Hazen)吸光度吸光度(A)(A)铂钴色度铂钴色度(Hazen)(Hazen)吸光度吸光度(A)(A)铂钴色度铂钴色度(H
57、azen)(Hazen)吸光度吸光度(A)(A)50.022300.0971000.301100.035350.1141500.452150.046400.1272500.738200.065500.1523000.901250.084600.1964001.180用回归分析法求直线方程式为Y(Hazen)ab(吸光度A) Y(Hazen)340(吸光度A)3绘制标准工作曲线，根据抗生素检测所需，仅需色度值060，即可.用线性回归作杂质含量的稳定性分析例例22关于BLT产品杂质C的分析根据生产车间2006年8月至2007年2月期间，BLT生产工序能力降低，导致出厂产品的杂质C含量较2005年平
58、均高出0.2%，根据18个月的留样产品中杂质C的分析，数据如右。留样时间留样时间（月）（月）杂质杂质C C含量含量00.0139330.0151460.0140990.01479120.01556180.01855R=0.863，Y=0.013+0.0002X，由于BLT的有效期为2年，当X=24时，Y=0.0183；S=0.00154，所以 =0.013+0.0002X3S，当X=24月时，Y的预报值为（0.02290.0137），即BLT杂质C含量为2.23%-1.37%，有99.7%的可能。新GMP第二百二十八条要求:3.宜采用便于趋势分析的方法保存某些数据(如检验数据、环境监察数据、制
59、药用水的微生物监测数据)。 控制图看趋势例例2323洁净环境趋势分析 例例2424工艺用水质量趋势分析数据分析在产品放行中的应用均均值置信置信区区间数数据据特性特性值线性性回回归正正态曲曲线 可运用的数学模型GMP的有关概念 可接受标准：Acceptance Criteria 指建立在相应的取样方法之上的药品质量检验标准和接受拒收的的标准（例如可接受不可接受的质量水平），是决定批准或拒收一批药品的必需因素。 超标结果：OOSout of specification 测试结果超过了预先设定的可接收标准。 超常结果：OOT out of trend测试结果超过了企业日常的水平。U检验确定产品的性能
60、值例例2525检验某产品的含氯量，根据经验已知含氯量服从正态分布N(60.20 ， 0.12)现抽查新产品六次，数据如下: 60.04 ， 60.63 ， 60.67 ， 60.66 ， 60.70 ， 60.71 .试问这批产品含氯量有否变化.解:已知 060.20 0.12 n6 假设H0:0 H 1: 0 =60.57=7.55取0.05， 查正态分布表临界值K0.05 1.96 | u |7.55K0.05 1.96 这批产品含氯量有显著性变化. 用 检验法检查物料改换的可行性 例2626车间生产某种药品的某成分的含量平均值为52.8%，标准差为1.6%，为降成本，改换 了原材料，生产
61、后抽取9 批样品，测得某成分的含量如下： 51.9，53.0，52.7，54.1，53.2，52.3，52.5，51.1，54.7 问 原材料能否改换? 解：假设H0：2=02 ，H1:2 02 选用2统计量，2(n1)S22， 当=0.10，f=n-1=8 查2分布表， 得12 1-/2 =2 0.95 =2.733， 22 /2 =2 0.05 =15.52 舍弃域为22.733 或 215.51 计算2 =(n1)S22=3.73。2.7332=3.7315.51 ， 2 值未落入舍弃区域 接受H0，方差没有显著变化。 又H0：=0=52.8 请朋友们自己完成。数据分析在稳定性考察中的应
62、用 线性回归控制图散点图因果图 可运用的数学模型 例例2727 2009年度进厂的 盐酸普鲁卡因含量趋势分析 用控制图考察进厂原料质量的稳定性例例2828对XX产品5年的留样某质量特性值(略)进行检 测,数值Y，(数据略) 按回归方程Y=aX + b， 其中X为留样时间（Y检测值省略）， 分别得到： Y普=0.0000756X + 0.3938 Y青=0.0000144X + 0.5568 Y水=0.00000865X +0.0322 YPH=0.0021X + 5.2543 Y效=0.1043X + 984.153用线性回归法考察产品的稳定性线性回归法确定药品的有效期例例29某药品在温度25
63、2，相对湿度605%条件下进行长期实验，得各时间的含量如下表，以时间为自变量X，含量（%）为因变量进行回归。时间（月）时间（月）0 03 36 69 912121818含量（含量（% %）99.399.397.697.697.397.398.498.496.096.094.094.0得回归方程：y=99.18-0.26x， r=0.8970， 因果图排列图控制图正交试验线性回归 可运用的数学模型数据分析在偏差管理中的应用新版GMP关于偏差评估第第2 25 57 7条 任何偏差都应评估对产品质量的潜在影响。企业可以根据偏差的性质、范围、对产品质量潜在影响的程度将偏差分类（如重大、次要偏差）。对重
64、大偏差的评估还应考虑是否需要对产品进行额外的检验以及对产品有效期的影响，必要时，应对生产或质量控制有重大偏差的产品进行稳定性考察。新版GMP: 260条 企业应建立纠正和预防措施系统，对投诉、产品缺陷、召回、偏差、自检或外部检查结果、工艺性能和产品质量监测趋势等进行调查并采取纠正和预防措施。调查的深度和形式应与风险的级别相适应。纠正和预防措施系统采用的方法应能改进产品和工艺，增进对产品和工艺的理解。你如何知道出了问题-识别问题Diagram 2Diagram3 获取数据分析数据追溯数据? ? ? ?坚 持 纠 正调查原因因因因果果果图图图控控控制制制图图图采取有效措施正正交交试试验验线线性性回
65、回归归 纠 正 问 题工工工序序序分分分析析析系系统统图图排排列列图图分分层层法法对对策策表表因果图因果图 让枯木让枯木逢春逢春抓住生产抓住生产抓住生产抓住生产现场的六现场的六现场的六现场的六大因素大因素大因素大因素人人料料法法环环测测机机注意事项：1.确定原因时应集思广益，以免疏漏。2.确定原因，尽可能具体。3.有多少质量特性，就要绘制多少张因果图。4.要把影响的主因在图上表出来。5.主因的确定要通过验证或确认。6.对确定的主因要采取对策整改环环人人料料法法测测责任心不强操作不熟练培训不到位考核制度不完善机机检测结果不准确空气湿度超标苄苄星星青青霉霉素素水水份份控控制制能能力力低低检测有偏差
66、洗涤乙脂水份高高循环水更换周期长真空偏低真空系统不畅通真空滤芯堵塞干燥温度控制不准确检测样品不具有代表性溶媒小罐消理不及时干燥时间控制不准确干燥工艺不完善 用因果图寻找产生偏差的原因 例例3030例例3131经过对11个末端因素进行分析(略)，得出影响水份控制能力的要因有3个。用系统图总结主因主主要要原原因因循环水更换周期长循环水更换周期长溶媒小罐清理不及时溶媒小罐清理不及时干燥工艺参数不准确干燥工艺参数不准确使用对策表 实现军中无戏言序号 目标问题点目标值对策与措施负责人 时间 备注用对策表进行纠正例例32制定对策序序号号要因要因对策对策目标目标措施措施地点地点负责负责人人完成完成时间时间1
67、循环水循环水更换周更换周期长期长统计每天真统计每天真空泵出口的空泵出口的真空度变化真空度变化使真空泵出使真空泵出口真空小于口真空小于-0.092Mpa缩短循环水更换周期缩短循环水更换周期真空泵真空泵房间房间张哲全张哲全2007.72真空小真空小罐清理罐清理不及时不及时对小罐进行对小罐进行改造，便于改造，便于操作操作小罐内残余小罐内残余溶媒不超过溶媒不超过2L1、在小罐底部焊接阀门、在小罐底部焊接阀门2、规定每批产品干燥完、规定每批产品干燥完毕清理小罐毕清理小罐二车间二车间何东召何东召2007.85干燥工干燥工艺参数艺参数不准确不准确采用正交试采用正交试验表，筛选验表，筛选最佳工艺参最佳工艺参数
68、数优化工艺优化工艺参数参数调节干燥时间、温度调节干燥时间、温度干燥间干燥间何华何华2007.8因素因素位级位级AB干燥温度干燥温度干燥时间干燥时间1706h2726.5h3747.0h 用正交试验 纠正工艺参数例例3333 降低苄星青霉素水分含量 所制定的水平因素分析表：因素因素试验号试验号AB试验结果试验结果添加折线图(略)温度温度时间时间水份（）水份（）1117.92217.63317.84127.85227.56327.77137.48237.29336.6K17.7%7.8%通过试验最佳组合为：通过试验最佳组合为：A2B3K27.2%7.7%K37.5%7.0%干燥温度为干燥温度为72
69、，干燥时间为干燥时间为7h为最佳工艺为最佳工艺 排列图 使管理对症下药 排列图 是为寻找主要质量问题或影响质量主要问题所使用的图. 排列图的原理巴雷特分析法:“关键的少数”和“次要的多数”. 帕累托分析:关键的少数;次要的多数排列图的 用途1.分析现象用，与不良结果有关:质量， 成本， 交货期，安全。2.分析原因用，与过程控制有关:操作者， 机器， 原材料， 作业方法。排列图的 分析通常占总频 数 80% 以 上的 项目是主要问题，占总频数 10% 以上的项目是次要 问 题 ，余下的占总 频 数 10% 左右的项目是更次要的一般问题。排列图 的作图步骤1. 收集数据.2. 选题.3. 设计数据
70、记录表.4. 填表并计算，按数量从大到小顺序填入表中，其他项的数据由许多数据很小的项目合并在一起，将其列在最后.5. 画两根纵轴和一根横轴，左纵轴标上频数的刻度，最大刻度主总频数;右纵轴，标上频率的刻度，最大5.画两根纵轴和一根横轴，左纵轴标上频数的刻度，最大刻度主总频数;右纵轴，标上频率的刻度，最大刻度为100%.左右刻度高度相等.6. 在横轴上按频数大小画上矩形，其高代表各项目频数的大小.7. 在每个直方柱右侧上方，标上累计值，描点，用实线连接，成巴雷特曲线.8. 在图上记全有关事项.排列图 在制药生产企业中的在制药生产企业中的应用应用1.寻找产品存在的质量问题2.寻找影响产品质量问题的主
71、因3.进行用户信息分析4.进行产品市场分析5.进行产品成本分析6.不合格品分析7.生产现场存在问题分析8.各类事故运用排列图进行投诉信息分析排列图适用于找出主要的质量问题，以便确定质量改进项目。我们对用户反馈的信息分别按品种、质量特性项及用户地区加以分层排列，找出主要产品，及该产品主要质量问题和存在问题多的所在地区加以针对性整改，努力降低质量问题，力求用户满意。 用排列图分析偏差的主因例例3434对分公司年用户来信的品种作出分层排列图如图所示。 从图可知A产品和B产品是用户反馈质量问题的主要品种。 产品产品A产品B产品C产品D产品E产品反馈批数反馈批数7315652累计批数累计批数738894
72、99101比率比率(%)72.314.86.04.92.0批数累积批数累积百分比率百分比率(%)72.387.193.198.0100.0xxxx年全年从图中可知，澄明度及吸潮是用户从图中可知，澄明度及吸潮是用户反映反映A A产品的主要质量特性。产品的主要质量特性。特性特性澄明度吸潮过敏包装松盖其它批数批数232110658累计批数累计批数234454606573比率比率(%)31.528.813.78.26.811.0批数累积批数累积百分比率百分比率(%)31.560.374.082.289.0100xxxx年全年从图看，湖北、四川、湖南是由于气从图看，湖北、四川、湖南是由于气候潮影响而吸潮
73、，河北是由于当时发候潮影响而吸潮，河北是由于当时发大水而影响产品吸潮，为此生产工艺大水而影响产品吸潮，为此生产工艺上要进一步控制上要进一步控制A A产品粉针的严密度。产品粉针的严密度。地区地区湖北河北四川湖南海南其他省反馈批数反馈批数653223累计批数累计批数61114161821比率比率(%)28.623.814.39.59.514.3批数累积批数累积百分比率百分比率(%)28.652.466.776.285.7100xxxx年全年用线性回归-对重大偏差产品的稳定性分析例例35关于XX产品杂质C的分析根据生产车间2006年8月至2007年2月期间，XX生产工序能力降低，导致出厂产品的杂质C
74、含量较2005年平均高出0.2%，根据18个月的留样产品中杂质C的分析，数据如右。留样时间留样时间（月）（月）杂质杂质C C含量含量00.0139330.0151460.0140990.01479120.01556180.01855R=0.863，Y=0.013+0.0002X，由于BLT的有效期为2年，当X=24时，Y=0.0183；S=0.00154，所以Y=0.013+0.0002X3S，当X=24月时，Y的预报值为（0.02290.0137），即BLT杂质C含量为2.23%-1.37%，有99.7%的可能。CP值计算排列图区间估计数据分析在供应商审核中的应用 可运用的数学模型 审核原材
75、料供应方的 质量保证能力 例例3636华药生产用的工业盐酸由大成、沙岭庄、石市电化厂供应，收集9798年盐酸含量数据，比较它们工序能 力 指 数 ，其中标准规定盐酸含量 31.0%。N为批数 TL规范下限 为盐酸含量平均值S为含量标准偏差 CPL为规定下限时，工序能力指数计算。 盐酸含量二年数据统计如下表盐酸含量二年数据统计如下表 供方名称供方名称N NX Xi i平均值平均值S SC CPLPL大成大成121233.94633.9460.7380.7381.331.33沙岭庄沙岭庄343431.91231.9120.5830.5830.520.52石市石市电化厂电化厂434332.26832
76、.2680.7650.7650.550.55 可知，大成厂工序能力正常，其它二厂生产工序能力严重不足，需供应方进行整改。 用区间判定设备使用质量例例373720B2模板A区抽检数据表冠部直径T:19.50.2 塞径直径T:130.1总高T:8.70.3冠部厚度T:30.2519.5012.968.753.0319.4612.948.753.1219.4812.948.703.0719.4812.948.662.9319.4912.938.703.0219.4912.928.743.0019.5012.948.763.0419.4712.938.743.0619.5112.938.703.031
77、9.5012.968.662.9719.4612.928.672.9419.4912.948.662.9619.49S0.01612.94S0.0138.71S0.0393.01S0.057TuTL:19.319.7TuTL:12.913.1TuTL:8.49.0TuTL:2.753.2519.419.512.913.08.68.82.843.18新版GMP 产品质量回顾分析: 第273条 应每年对所有生产的药品按品种进行产品质量回顾分析，以确认工艺稳定可靠，以及原辅料、成品现行质量标准的适用性，及时发现不良趋势，确定产品及工艺改进的方向。应考虑以往回顾分析的历史数据，回顾分析应有文件记录。
78、回顾分析报告的内容：1.产品所用原辅料的所有变更，尤其是来自新供应商的原辅料；2.关键中间控制点及成品的检验结果；3.所有不符合质量标准的批次及其调查；4.所有重大偏差（包括不符合要求）及相关的调查、所采取的整改措施和预防措施的有效性；5.生产工艺或检验方法的所有变更；6.药品注册所有变更的申报、批准或退审；7.稳定性考察的结果及任何不良趋势；8.所有因质量原因造成的退货、投诉召回其当时的调查；9.其它以往产品工艺或设备的整改措施是否完善；10.新获得注册批准的药品和注册批准有变更的药品上市后的质量状况；11.相关设备和设施，如空调净化系统、水系统、压缩空气等的确认状态；12.对技术协议的回顾
79、分析，以确保内容更新。特征值因果图排列图控制图CP值正交试验假设检验线性回归原辅料变更中间体检验成品检验变更稳定性考察偏差返工批报废批投诉水，空气质量性能趋势分析验证与确认方法应用项目数据分析在回顾分析报告 的灵活应用 (仅供参考)注:应根据具体情况选用最适宜的检验工作数数据据特特征征值值的的应应用用生产管理质量风险 回顾总结(1)数据特征值的应用 确立分析结果设备使用质量判定可信的技术报告的关键同一事物的比较中间体质量标准制定质量特性置信区间的估计检验误差的表示COA数据的发出确定标准的分析方法质量标准分析方法验证要求工艺用水的水质控制洁净区环境的控制质量企内标准指标的确定回顾总结(2):
80、工序能力指数的应用监督工序质量的一种方法企业质量标准制订后的验证比较几个总体的工序质量某个产品质量问题的予防审核原材料供应商的质量C CP值 可运用的数学模型 寻找问题的因素寻找问题的因素纠正工艺参数纠正工艺参数 确定产品配方确定产品配方 正交试验的应用 回顾总结(3):检验标准品对照仪器使用比较专业水平比较.质量与标准值的比较实验室数据比对新老工艺新老工艺比较比较验证数据处理供应商产品质量对比回顾总结(4): 假设检验的应用回顾总结(5):回归分析的应用确立工作曲线找工艺条件质量特性关系趋势分析找主因找主因回顾总结(6): 因果图的应用的应用物料控制 工艺控制成品质量成品质量设备运行 确认与验证 因果图产品市场分析产品的质量问题产品产品成本成本不合格品分析 生产现场 问题投诉信息分析排列图排列图各类各类事故事故各类偏差分析回顾总结(7):(7): 排列图的的应用 生产工序现生产工序现 场控制场控制 企内标准制订洁净环境 趋势分析回顾总结(8):(8):应 用控制图
