药品生产数据分析

蜗牛982015年04月08日1一、为什么要进行这次的培训二、数据的来源三、数据的描述性统计四、数据的图形化五、控制图六、警戒限与行动限七、生产统计应用21.法规符合性2.质量控制3.总结、报表4.工艺验证•回顾分析、培训需求•生产过程控制•季度、年中、年终总结，上交各种报表3·持续工艺验证1.物料与产品质量•物料来源、数量、质量•各阶段产品质量•各阶段产品生产•……2.环境和设备监控•环境监测•纯化水监测•设备参数•工艺参数•……3.其他•培训•偏差•变更•……•自己关注•领导需要4一、为什么要进行这次的培训二、数据的来源三、数据的描述性统计四、数据的图形化五、控制图六、警戒限与行动限七、生产统计应用5数据集中性数据离散性描述性统计67集中性众数中位数……均值离散性极差方差标准差相对标准偏差……8下面给出一组生产数据，来进行相应的计算。压片产量数据（kg）编号产量编号产量编号产量170.50971.221770.14271.601071.621871.06371.501171.701970.96470.751271.742070.96571.201370.682170.86670.841470.502270.98771.301570.26————871.301670.529均值（算术平均值）•Xbar=（X1+X2+……+Xn)/n=70.01众数•数据中出现次数最多的数据（数据可是一个，也可以是多个）：70.50、70.96、71.3010170.14270.26370.50470.50570.52670.68770.75870.84970.861070.961170.961270.981371.061471.201571.221671.301771.301871.501971.602071.622171.702271.7411中位数（70.97）判定计数数据排列升序/降序偶/奇数偶数中间两数均值奇数中间数（n+1）/2整数整数位数值小数两侧数据均值•数据排序•Xmax-Xmin=1.6极差•σ2=var（X）=Σ（Xj-Xbar）2/n•S2=Σ（Xj-Xbar）2/（n-1）=0.2159方差•σ•σ≈S=0.465标准差•RSD=σ（s）/Xbar×100%•RSD=0.66%相对标准偏差170.14270.26370.50470.50570.52670.68770.75870.84970.861070.961170.961270.981371.061471.201571.221671.301771.301871.501971.602071.622171.702271.7412通过上述计算，我们得出了素片产量数据的描述性统计结果项目平均值中位数众数最大值最小值数值70.0170.9770.50、70.96、71.3071.7470.14项目极差方差标准差相对标准偏差——数值1.60.21590.4650.66%——13为什么进行图形化？不直观难读懂浪费时间14一、为什么要进行这次的培训二、数据的来源三、数据的描述性统计四、数据的图形化五、控制图六、警戒限与行动限七、生产统计应用15图形散点图折线图柱形图饼图……164.1散点图70.0070.2070.4070.6070.8071.0071.2071.4071.6071.8072.000510152025产量174.2折线图70.0070.2070.4070.6070.8071.0071.2071.4071.6071.8072.000510152025产量18对于素片产量图表里面的数据，我们可以使用简单的柱形图。69.0069.5070.0070.5071.0071.5072.0012345678910111213141516171819202122产量194.4饼图20分类少，每个类别的数据差别较大数据大量重复出现其他合适的情况21224.4.2某一数据大量重复出现的饼图其他图形帕累托图（柏拉图）茎叶图箱线图直方图……23一、为什么要进行这次的培训二、数据的来源三、数据的描述性统计四、数据的图形化五、控制图六、警戒限与行动限七、生产统计应用24•区别引起关键质量特性（CTQ）或关键工序CTP）产生波动是正常原因和异常原因的一种工具。概念•是对关键质量特性或引起波动的正常原因的稳定性进行评估和监控。作用•从每个观察时段的子集或样本中采集数据采集2526产生原因波动类别处置方法正常波动异常波动过程产生过程外因素引起管理帮助采取措施1.防止将引起波动的正常原因错误的归于异常原因时，对过程进行过渡调整。2.防止将波动由异常原因归于是正常原因引起，275.3控制图条带28•1个点，距离中心线大于3个标准差1•连续9个点在中心线同一侧2•连续6个点，全部递增或全部递减3•连续14个点，上下交错4•3个点中有2个点，距离中心线（同侧）大于2个标准差5•5个点中有4个点，距离中心线（同侧）大于1个标准差6•连续15个点，距离中心线（任一侧）1个标准差以内7•连续6个点，距离中心线（任一侧）大于1个标准差8291•制定取样计划，收集数据2•判定数据类型3•根据数据选择合适的控制图4•制作控制图5•结果分析30数据计数型二项分布N图/NP图泊松分布C图/U图计量型正态分布控制图其他分布……3132计数数据N图/NP图样本数相同NP图P图样本数不同P图C图/U图取样单位相同C图U图取样单位不同U图33计量数据子组2≤n≤8Xbar-R图n≥9Xbar-S图单值I-MR图各图的名字R图、S图、MR图Xbar图、Xbar图、I图各图的作用检测过程稳定性检测数据集中趋势计量控制图过程的波动图过程集中趋势图3435363738控制图分析用控制图用于数据分析控制用控制图用于过程控制391.现在有一组上一年PVC硬片的使用数据。2.根据这组数据来制定一个合理的PVC使用限度来对之后PVC的用量进行控制。3.工具：Minitab4.方法：质量控制图40编号铝箔万片（kg）编号铝箔万片（kg）10.5280150.564920.5051160.512730.5649170.523740.5127180.507750.5237190.621960.5077200.539270.6219210.501180.5392220.478090.5011230.4919100.4780240.4877110.4919250.4659120.4877260.4732130.5280270.5518140.5051————4142•计数/计量1.数据类型•正太性检测、数据剔除2.数据处理•计量/计数控制图3.制作控制图•结果4.确定控制限度•计量数据数据类型•单值子组容量•I-MR图控制图43正太性检测数据剔除正太性再检测44使用正太概率图P=0.017＜0.05（α=0.05，显著性水平），数据不符合正太分布5.7.3数据的正太性检测45数据剔除的方法：四分位数法46编号铝箔万片（kg）编号铝箔万片（kg）10.5280150.564920.5051160.512730.5649170.523740.5127180.507750.5237190.621960.5077200.539270.6219210.501180.5392220.478090.5011230.4919100.4780240.4877110.4919250.4659120.4877260.4732130.5280270.5518140.5051————47P=0.627＞0.05，数据符合正太分布。485.7.4完成控制图49均值•0.51控制限•0.59•0.43实际控制上限•0.5550•对生产过程进行控制制作目的•防止异常数据影响使控制范围变宽剔除原因•质量改进是一个动态过程，控制图也是一个持续改进的过程持续改进511.目的：一直好奇压片的片重的稳定性，就弄了一个方案，准备对片重的稳定性和工序能力进行一下了解。2.数据收集计划：按照一个正常批次的压片时间，进行了25等分，每一个时间段取五片，分别称量片重，记录。3.方法：制作控制图，检查压片稳定性；过程能力分析，检查工序能力。4.工具：Minitab-控制图5253抽样简单随机抽样抽签随机号分层抽样分成类似的组，从每组中抽样，注意各组比例混合含量均匀度整群抽样包含各类情况的类似的群或组，然后抽取整个群成品取样系统抽样每k个单位抽取一次片重54Cp：过程能力指数,又称为潜在过程能力指数,为容差的宽度与过程波动范围之比。（短期）Cpk：过程能力指数,又称为实际过程能力指数,为过程中心μ与两个规范限最近的距离。Pp与Ppk：过程性能指数，计算方法与计算Cp和Cpk类似，所不同的是，它们是规范限与过程总波动的比值。（长期）子组编号片重检测结果（mg）S1S2S3S4S5均值1193192189192198192.802188193190189183188.603184196189187185188.204195186184184194188.605184184190189185186.406191187189186194189.407186190182186191187.008185184187190187186.609188194184189191189.2010190188182193191188.8011197192188192194192.6012184192189193196190.8013189188196188188189.8014193183185194191189.2015188184185190186186.6016193191191193194192.4017194192189195193192.6018191191190186192190.0019193189182196188189.6020185184194187196189.2021186192188192190189.6022188183186193194188.8023186188190190189188.6024189191187190188189.0025190189192184190189.0055•计量数据数据类型•5子组容量•Xbar-R图控制图565758组内标准差：3.65，Cp=1.21，Cpk=1.15整体标准差：3.87，Pp=1.19，Ppk=1.1759范围等级判断Cp1.67特级工序能力过高1.67Cp1.33一级工序能力充分1.33Cp1.0二级工序能力尚可1.0Cp0.67三级工序能力不充分0.67Cp四级工序能力不足606162一、为什么要进行这次的培训二、数据的来源三、数据的描述性统计四、数据的图形化五、控制图六、警戒限与行动限七、生产统计应用63控制图•Xbar±2σ•Xbar±3σ百分位数•95%•99%规格限百分比•60%•80%64合适的才是最好的65一、为什么要进行这次的培训二、数据的来源三、数据的描述性统计四、数据的图形化五、控制图六、警戒限与行动限七、生产统计应用66年终总结•发现A片成品收率连续三年降低疑问•降低是真实的？•什么引起的原因？需要的工作•确定真实性•查找可能原因67A片成品收率数据项目2011年2012年2013年收率93.22%92.68%92.28%较上年——-0.54-0.40较前年————-0.94686970包衣包衣的影响？包衣产量数据压片压片的影响？压片产量和损耗数据制粒制粒的影响？制粒产量数据相关性颗粒产量素片产量包衣片产量成品收率颗粒产量Pearson相关性1.767**.209.406**显著性（双侧）.000.087.001N68686868素片产量Pearson相关性.767**1.400**.524**显著性（双侧）.000.001.000N68686868包衣片产量P