第十九章 质量管理的统计方法

参考书目为安德森的《商务与经济统计》,以下为个人的学习总结,如果有错误欢迎指正。有需要本书pdf的,链接在本文末尾。(仅限个人学习使用,请勿牟利)

第十九章 质量管理的统计方法

德克萨斯州的陶氏化学公司,在生产过程中每隔一个周期抽取产品组成样本,得到样本均值,并记录在\bar x控制图上。如果有异常的均值可以追溯到生产操作员,并重新培训操作员。让产品质量大幅改进。

美国质量协会(American Society for Quality,ASQ)对质量的定义:质量时产品或服务的所有性质和特征,这些性质和特征使得产品或服务能够满足特定的需要。如今的竞争中,产品质量尤为重要。
从简单排除生产线上的缺陷产品到建立运用广泛的公司战略。扩大质量的范围自然导致全面质量(total quality,TQ)的概念。
全面质量在不同组织中实现方式不尽相同,主要依据三个基本原理:

  • 以顾客和股东为中心
  • 将参与精神与团队合作贯穿整个组织
  • 以不断改进和学习为中心

19.1 理念和框架

20世纪初,质量管理实践仅限定于检验已完成的产品并剔除有缺陷的,此后有沃尔特A.休哈特(Walter A. Shewhart)、W.爱德华·戴明(W. Edwards Deming)和约瑟夫·朱兰(Joseph Juran)做出了巨大贡献。

朱兰给质量下了定义:质量必须与实用性相匹配。
朱兰的质量方法关注三个质量过程:

  • 质量计划
  • 质量控制
  • 质量改进

朱兰与戴明的需要在组织中开展重大文化变革的理念相比,朱兰的方案旨在当前的组织系统内提高质量。

在质量改进运动中,许多学者都发挥了作用,其中两个最重要的项目是美国的马尔科姆·鲍德里奇国家质量奖和ISO 9000国际认证过程,近年来,六西格玛的使用也在增加。

19.1.1 马尔科姆·鲍德里奇国家质量奖

由美国总统颁发的马尔科姆·鲍德里奇国家质量奖(BNQP)授予在七个方面表现杰出的组织:

  1. 领导
  2. 战略策划
  3. 客户和市场
  4. 测量、分析和信息管理
  5. 人力资源
  6. 流程管理
  7. 经营成果

1998年首次颁奖,在连续8年中,假设股票指数(获得马尔科姆·鲍德里奇国家质量奖的美国上市公司)优于标准普尔500,在某一年里这个比是4.4:1。

19.1.2 ISO 9000

ISO 9000由5个国际标准组成的一个系列标准,由设在瑞士日内瓦的国际标准化组织(ISO)于1987年颁布。公司可以利用这些标准帮助确定,维持一个有效的质量一致性体系需要有哪些方面。
例如:标准描述一个有效质量体系的需要;对确保定期测量和检测设备使其符合标准的需要,以及对维持一个令人满意地记录保持系统地需要。
总之,ISO 9000认证覆盖地鲍德里奇质量奖标准不到10%。

19.1.3 六西格玛

二十世纪80年代后期,摩托罗拉他们的目标是达到非常高的质量水平(每百万产品出现缺陷的机会不超过3.4),这个质量水平被称为六西格玛质量水平,达到这种质量目标的方法被称为六西格玛(six sigma)
一个组织可以承担两类六西格玛方案:

  • DMAIC(Define, Measure, Analyze, Improve and Control,界定、测量、分析、改进和控制)
  • DFSS(六西格玛设计)用于设计新产品、过程或服务。

六西格玛线和每百万个机会中的缺陷数,每百万个机会中的缺陷数(dpmo),要求dpmo不超过3.4个。

举例:KJW包装公司经营一条填装谷物箱得生产线,填装过程中\mu=16.05盎司,\sigma=0.10盎司。并假设填装总量服从正态分布。

image

假设管理人员认为15.45~16.65盎司是可接受得质量限。此时只有0.0000002%的谷物不在这个区间。
如果我们调整均值到16.20,标准差不变,此刻当重量超过16.65盎司的概率为0.0000034,即百万箱中有3.4箱。
我们可以得到结论,只要均值在16.05的1.5个标准差之内,就可以控制质量是满足六西格玛质量水平。即均值满足15.9~16.2这个区间时,KJW填装过程将被视作六西格玛过程,则每百万填装箱中最多只有3.4箱存在缺陷。

image

质量控制包括一系列的检验和测量以确定是否满足质量标准。下面两节将介绍质量控制中使用的两种统计方法。

19.1.4 服务业中的质量

相比制造业,服务业关注确保消费者满意度和提升消费者体验。通常留住老客户比拉新成本要低,因此设计提升消费者服务的质量控制过程对服务业至关重要。
但是服务业的消费者满意度非常主观,测量极具困难。我们可以通过测量诸如提供服务的及时性和满意度进行调查。

19.2 统计过程控制

在一个连续不制造产品的生产过程中考虑质量控制程序,以产品生产量的抽样和检验为依据,可以做出是继续生产过程还是调整生产过程的决定。

生产的产品在质量上的变异来源可以分为可指出的原因一般原因,前者指机器磨损、错误操作、原材料有问题等;后者指温度、湿度等差异,这些一般无法控制,也没必要控制。

统计质量控制的主要目标是可指出原因时我们应该调整生产过程,达到可以接受的质量水平。其中过程控制统计程序的依据是假设检验:


image

19.2.1 控制图

控制图(control chart)对确定产品中的质量变异是来源于一般原因(在控)还是来源于可指出的原因(失控),提供了一个做出决策的基础。
控制图根据所包含的数据类型进行分类,这里介绍\bar x控制图

image

当过程处于在控状态时,\bar x值大概率位于上控制限(UCL,Upper Control Limit)下控制限(LCL,Lower Control Limit)之间(补充中心线(CL, Central Line))。当\bar x值位于控制限之外时有强有力的统计依据,过程处于失控状态。横轴是时间,每过一段时间新增加一个点时,我们都相当于做了一次假设检验来确定过程是否处于在控状态。

19.2.2 \bar x控制图:过程的均值和标准差已知

例子:KJW包装公司填装谷物箱时,平均重量\mu=16.05盎司,标准差\sigma=0.10盎司;其中抽样分布的标准差\sigma_{\bar x}=\frac{\sigma}{\sqrt{n}}

image

回顾正态概率分布,大约99.7%的随机变量的数值都位于均值\pm 3个标准差之间。因此\bar x控制图的控制限如下:
\bar x控制图的控制险:过程的均值和标准差已知:

  • UCL=\mu-3\sigma_{\bar x}
  • LCL=\mu+3\sigma_{\bar x}

回到KJW,我们通过计算得到UCL和LCL的值。在样本序号5时处于失控状态,经过调整后回到在控状态。


image

19.2.3 \bar x控制图:过程的均值和标准差未知

例子:Jensen计算机用品公司生产直径3.5英寸的固态硬盘,每小时抽取5个组成一个样本,直到抽取了20个样本为止,得到样本均值\bar x_j和极差R_j,数据如下:

image

总样本均值:\bar {\bar x}=\frac{\sum_{k=1}^{k}\bar x_k}{k}
平均极差:\bar R=\frac{\sum_{k=1}^{k} R_k}{k}

前面讲到\bar x控制图的控制上下限为:\bar x \pm 3\frac{\sigma}{\sqrt{n}},其中\mu的点估计为\bar {\bar x}\sigma的点估计可以使用极差。
统计学家已经证明:过程标准差\sigma的一个估计量为平均极差除以d_2d_2是一个依赖于样本容量n的常数,通过查控制图系数表得到),即\sigma的估计量=\frac{\bar R}{d_2}

\bar x控制图的控制限:\bar {\bar x}\pm 3\frac{\bar R/d_2}{\sqrt{n}}=\bar {\bar x}\pm\frac{3}{d_2\sqrt{n}}\bar R=\bar {\bar x}\pm A_2\bar R
其中A_2=3/(d_2\sqrt{n})是一个仅依赖于样本容量的常数。
最终计算得到UCL=3.514,LCL=3.485。我们将数据制作成控制图,发现抽样期间过程的均值总是在控状态。

image

19.2.4 R控制图

为了构造极差控制图(R控制图),我们把样本的极差看作一个随机变量,那么这个随机变量的估计值为\hat \sigma_R=d_3\frac{\bar R}{d_2},接下来尝试计算UCL和LCL。

  • \hat \sigma_R=d_3\frac{\bar R}{d_2}代入
    • UCL=\bar R+3\hat \sigma_R=\bar R(1+3\frac{d_3}{d_2})
    • LCL=\bar R-3\hat \sigma_R=\bar R(1-3\frac{d_3}{d_2})
  • 我们令
    • D_4=1+3\frac{d_3}{d_2}
    • D_3=1-3\frac{d_3}{d_2}
  • 简化后
    • UCL=\bar R D_4
    • LCL=\bar R D_3

回到Jensen计算机用品问题的R控制图,计算得到UCL为0.053,LCL为0,\bar R为0.0253。得到下述控制图:

image

由于20个样本的极差都在控制限以内,我们确定抽样期间过程的变异性处于在控状态。

19.2.5 P控制图

p控制图是使用有缺陷项目比例的数据构造的控制图。

例子:邮局自动分拣机器有3%的概率误投。同样\bar p的抽样分布可以用来确定\bar p值所期望的变异性。
在第七章我们知道,\sigma_{\bar p}=\sqrt{\frac{p(1-p)}{n}},且要求np和n(1-p)都大于等于5时抽样分布近似服从正态分布。

P控制图的控制限

  • UCL=p+3\sigma_{\bar p}
  • LCL=p-3\sigma_{\bar p}

回到邮局的例子,其n=200,计算得到\sigma_{\bar p}=0.0121,则UCL=0.0663,LCL=-0.0063(概率不能为负数,这里手动调整为0)

image

上图表明,所有样本都是在控状态时抽取的。当一个处于在控状态的过程,若缺陷比例未知,我们可以通过抽取k个容量为n的不同样本,来计算p以及标准误差。

19.2.6 np控制图

np控制图是针对样本中有缺陷项目的个数而构造的控制图。此时n为样本容量,p为当过程处于在控状态时所观测到的有缺陷项目的概率。
当np和n(1-p)都大于等于5时,有缺陷项目书的分布近似服从均值为np,标准差为\sqrt{np(1-p)}的正态分布。

np控制图的控制限

  • UCL=np+3\sqrt{np(1-p)}
  • LCL=np-3\sqrt{np(1-p)}

回到邮件分拣的例子中,我们计算得到UCL=13.2375,LCL=-1.2375(手动调整为0)。因此当超高13件被分错时我们认为过程处于失控状态。

np控制图和p控制图所提供的信息相同(判断标准一样),区别在于前者是缺陷项目数量,后者是缺陷项目比例。

19.2.7 控制图的解释

控制图帮确定一个过程是否处于统计控制状态,错判概率很小。如果数据点落在控制限之外就应该调整生产过程。
除了数据点位于控制限之外,

  • 当数据点集中分布在中心线的一侧时也是预警信号,我们应该对生产过程进行仔细检查。
  • 当控制图出现一个趋势时(如当6~7个点呈现单调上升/下降时),也应该引起警示,可能出现设备老化或其他原因。

总结:单纯的\bar x控制图依赖于平均极差的数值,所以控制限没有太多的意义,因此实践中通常在构造\bar x控制图之前构造R控制图,如果R控制图表明过程的变异性处于在控状态,然后再构造\bar x控制图

19.3 接受抽样

在采购或进货的项目中,我们希望以指定产品的质量特性为依据来决定是否接受一组产品项目。在质量控制的术语中,一组项目称为一批(lot)接受抽样是一种统计方法,让我们在检测一批项目样本的基础上做出接受或拒绝的决定。

下图为接受抽样的一般步骤


image

接受抽样的程序:

  • 假设
    • H_0:高质量批
    • H_a:低质量批
  • 得到结论


    image

其中针对第一和第二类错误的意义:

  • 生产者风险:错误的拒绝了一个高质量批
  • 消费者风险:错误的接受了一个低质量批

19.3.1 KALI有限公司:接受抽样的实例

KALI是一家生产家电的公司,其中空调产品的超载保护器需要采购。为了保障质量我们可以对收到的每一个部件进行检验,这种方法被称为100%检验法,不过这个部件的检验费时且费用高,不可能实现100%检验。
替代方法:接受抽样

  • 一个接受抽样方案包含样本容量n和接受准则c。
  • 接受准则是在样本中发现并能接受有缺陷项目的最大数量。
  • 举例:n=15,c=0。则当样本中完全没有缺陷项目才能选择接受该批,否则拒绝该批。
    n和c的选择要权衡生产者风险和消费者风险在一个合理水平。

19.3.2 计算接受一批的概率

还是KALI的例子,我们收到大批量超载保护器,有5%的缺陷项目,抽样方案取n=15,c=0。

接受抽样的二项概率函数f(x)=\frac{n!}{x!(n-x)!}p^x(1-p)^{(n-x)}
其中n为样本容量,p为批中缺陷项目比例,x为样本中有缺陷项目数量,f(x)为样本中有x个缺陷项目的概率。
我们通过计算可以得到f(0)=0.4633,则拒绝的概率为0.5367(EXCEL可以用BINOM.DIST(0,15,0.05,1)来计算)
当缺陷项目比例开始变动时,接受该批的概率也随之改变:

image

根据表格内的数据绘制抽样特性曲线(operating characteristic,OC)
image

我们使用不同的接受抽样方案可以得到不同的抽样特征曲线。
image

19.3.3 选择接受抽样方案

作为管理人员,我们需要对批中有缺陷项目的比例指定两个数值:

  • p_0用于控制生产者风险,生产者风险\alpha:拒绝有缺陷比例为p_0的一批的概率
  • p_1用于控制消费者风险,消费者风险\beta:接受有缺陷比例为p_1的一批的概率

假设KALI公司的管理人员规定p_0=0.03,p_1=0.15

  • 当p=0.03时,我们选择接受的概率为0.63,拒收的概率为1-0.63=0.37。即生产者风险为0.37。(风险过高)
  • 当p=0.15时,我们选择接受的概率为0.087,举手的概率为1-0.087=0.913。即消费者风险为0.087。
image

若管理人员要求生产者风险\alpha=0.1,消费者风险\beta=0.19
利用p_0=0.03,\alpha=0.10,p_1=0.15,\beta=0.20,在图19-13表明n=20,c=1的接受抽样方案接近于同时满足生产者风险合消费者风险的需求。

解和逾期的生产风险和消费者风险,我们可以使用一些抽样方案表(如美国军用标准表(MIL-STD-105D))。

19.3.4 多重抽样方案

本节介绍的KALI公司的接受抽样是单样本方案,是因为仅仅用了一个样本或一个抽样阶段。
多重抽样方案涉及两个或多个抽样阶段,每个阶段有三种可能的结论:

  • 停止抽样接受该批,
  • 停止抽样拒绝该批,
  • 继续抽样。

流程如下图所示:

image

双重抽样方案的建立更加困难,因为样本容量n_1,n_2,c_1,c_2,c_3必须同时满足生产者和消费者所预期的双重风险。

补充
对于接受抽样,使用二项分布的依据是假设批的容量较大,如果较小应使用超几何分布。


链接: https://pan.baidu.com/s/1fc0q-Q4kj3g-7Fr4MHZaqw 提取码: 333c 复制这段内容后打开百度网盘手机App,操作更方便哦

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 199,830评论 5 468
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 83,992评论 2 376
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 146,875评论 0 331
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 53,837评论 1 271
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 62,734评论 5 360
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,091评论 1 277
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,550评论 3 390
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,217评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,368评论 1 294
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,298评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,350评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,027评论 3 315
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,623评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,706评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,940评论 1 255
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,349评论 2 346
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 41,936评论 2 341

推荐阅读更多精彩内容