基于HL-MAD组合估计量的常规控制图稳健改进.pdf

收稿日期:基金项目:国家社会科学基金项目(B G L );广东省普通高校特色创新类项目(K T S C X )作者简介:李熠玲(),女,湖北黄冈人,硕士研究生,主要从事稳健时间序列分析的研究.通信作者:王志坚,E m a i l:q q c o m.基于HLMAD组合估计量的常规控制图稳健改进李熠玲,任凯亮,王志坚李熠玲,任凯亮,王志坚(广东财经大学 统计与数学学院,广东 广州 ;康涅狄格大学 文理学院,美国 康涅狄格州 )摘要:常规控制图的控制中心及控制限是不稳健的,易受离群值影响而导致监控失效.采用HLMA D组合估计量对控制中心、控制限进行稳健改进,构建出稳健控制图.数值模拟与金融实证分析均表明,相比于传统控制图,稳健控制图能更好地对受控过程的异常现象进行识别及克服“漏报警”现象,且在大样本、高污染率下效果更优.另外,通过估计量的相对效比较研究发现:在控制图的稳健估计量选择问题上,组合估计量比单个估计量具有显著的优良性.关键词:统计过程控制;HLMA D组合估计量;稳健估计;相对效中图分类号:O 文献标识码:A文章编号:()D O I:/j i s s n 引言 引言统计过程控制是一项建立在数理统计学原理基础之上的过程质量管理技术.从诞生至今,经过近 年的发展与完善,统计过程控制技术已经被广泛应用于社会经济发展的各个领域,包括物流、数字经济、服务业等,取得了良好的社会和经济效益,其监控的过程也从传统单变量过程拓展到各种多变量场景.而在实际工作中,控制图则是常被用来对受控过程进行质量管理的工具之一,通过使用控制图,能够对过程绩效及其波动原因进行直观而实时的监控和分析,帮助人们及早地发现和识别异常现象并找出异常原因,从而提高产品或服务质量.质量控制图的基本原理主要包含以下四个方面:受控过程服从正态分布的假定;控制限的准则;在一次试验中小概率事件不发生的原理;数理统计学中的统计推断基本思想.受控过程的正态性假定是指受控过程在正常状态下,其相关的特性指标数据服从正态分布假定.控制限的准则是指当某个过程特性数据服从正态分布N(,)时,可以利用正态分布的准则得出控制图的设计原理即控制上限U C L,控制中心线C L,控制下限L C L.小概率事件不发生原理是指在一次观测试验中,正常情况下小概率事件通常认为不发生.统计推断的基本思想认为,若在过程控制中一旦出现了“小概率事件发生”的现象,则可以认为受控过程出现异常.同时,应注意到,一方面,统计过程控制技术的理论基础之一,控制限的准则即为常规控制图的控制设计原理,其中控制上限U C L、控制中心线C L、控制下限L C L所采用的统计量,均值与标准差对离群值敏感,缺乏稳健性;另一方面,在离群值存在时,过程质量特性数据很难满足正态分第 卷第期 年 月 上 饶 师 范 学 院 学 报J OUR NA LO FS HAN G R AONO RMA LUN I V E R S I T YV o l ,N o J u n 布假定.因此,常规控制图在对过程异常情况识别时,易受离群值干扰,常出现“遮蔽效应”或“漏报警”现象,从而极大地降低了控制图监测性能,如何解决此类问题即是本研究的出发点.文文献献综综述述近 年,稳健统计技术在各领域的过程控制中应用越来越多,包括经济、金融、工业制造、服务业等.刚开始多是对单变量控制图采用不同的稳健方法进行设计,如均值X b a r控制图.维贾雅(V i j a y a)和墨蒂(M u r t y)讨论了采用基于风险的方法来找到X b a r控制图的最佳参数对X b a r控制图进行稳健经济设计.周纯光等提出了一种基于小波的预分析稳健控制图,用于监测第一阶段过程控制中的均值漂移.吴纯杰和王兆军分析了休哈特X b a r控制图的不稳健性,并对其进行了稳健修正.阿卜秋勒(A b d u l)等研究了利用考虑不确定参数区间估计的稳健优化方法进行X b a r控制图的经济统计设计(E S D),开发了一种启发式算法来获得控制图的稳健方案,效果优于传统E S D.萨利赫(S a l i h)等对因素选择控制图的最优稳健设计进行了研究.随后,越来越多的学者开始关注多变量控制图的稳健性.哈布沙(H a b s h a h)和阿什坎(A s h k a n)提出了稳健多元C U S UM图和多元EWMA图,以解决散点异常值变化小的问题.范(F a n)等基于分层聚类树原理开发了一种新的稳健多变量控制图,该图可以有效地检测多维数据中的潜在异常值,同时控制遮蔽和淹没效应.阿索坎(A s o k a n)和贾亚尚卡尔(J a y a s a n k a r)提出了监测第一阶段多变量个体观测过程均值的稳健控制图.张(Z h a n g)等提出了一种无分布的多变量统计过程控制图(M S P C),以检测多变量过程变量的一般分布变化.安杰尼丝(A n g e l l y s)等 将截尾均值应用于稳健多变量控制图的异常值诊断.拉吉(R a j i)等 基于S t a h e lD o n o h o稳健估计量(S D R E)构建了一种稳健多变量控制图,同时从第一阶段估计过程参数.卡丽纳(C a b a n a)和利洛(L i l l o)基于稳健重加权收缩估计,提出了一种用于个体观测的稳健多变量质量控制技术.萨巴诺(S a b a h n o)和塞拉诺(C e l a n o)用可变参数控制图监测存在自相关的多变量变异系数.近年来,也有学者研究自相关过程残差控制图的稳健性.王志坚 通过权重函数对A RMA模型与GA R CH模型进行稳健建模,最后构建稳健残差控制图.王志坚等 通过构建稳健A R模型,为自相关过程残差控制图的设计提供了理论依据.萨里阿提(S h a r i a t i)提出了一种适用于自相关序列的稳健控制图新方法,该图对污染数据的影响具有稳健性.萨拉赫(S a l a h)等 研究了伽马回归模型下残差控制图的剖面监测.通过梳理文献发现,不少学者从控制图控制限的位置参数、尺度参数角度采用稳健估计量来构建稳健控制图.哈菲兹(H a f i z)等 建立了稳健S h e w h a r t位置参数控制特征的逐步筛选方法.纳迪娅(N a d i a)和沙希德(S h a h i d)比较了文献中的六种不同稳健尺度估计下的EWMA控制图性能,模拟研究结果表明,基于估计量Q_n的控制图在非正常过程中表现相对较好.吴纯杰等、庄芳等 分别分析了稳健似然比累积和控制图及EWMA方差控制图的不稳健性,并给出了稳健改进的方法.王志坚、苏拥英等 对常规过程控制图的敏感性进行了分析并给出了稳健化方法.卡奥(K a o)认为,当存在污染数据时,标准偏差的估计会由于其高偏差性而降低控制图的检测能力,于是提出了基于不同筛选的平方A估计量,该估计量在抗干扰方面表现相对最好.与前人不同的是,本研究尝试采用估计量H o d g e s L e h m a n n(简写HL)与中位绝对离差(MA D)相结合的方式对常规控制图控制中心及控制限进行稳健改进,并通过改进前后对比的研究方法验证稳健控制图的可行性和有效性.相相关关理理论论与与方方法法常规控制图的不稳健性研究常规休哈特控制图设计原理,主要包括控制中心、上下控制限,相应表达式如()式所示:第期李熠玲,等:基于HLMA D组合估计量的常规控制图稳健改进常规控制图U C L:k C L:L C L:k()()式中的k通常取,显然U C L、C L、L C L所采用的统计量缺乏稳健性,导致整个控制图对离群值过于敏感,当受控过程出现异常情况时会使得控制中心、控制限的值不能反映大多数样本数据的特征,因此控制图常会出现“漏报警”现象.接下来通过举例论证在对异常现象监控时常规控制图如何“失控”.先采用R软件生成 个服从N(,)分布的随机数作为过程的特性指标数据,再通过随机抽样技术在 个随机数里面随机抽取一个数用去代替之,我们将作为 个随机数里面的异常值.接下来基于不含异常值及含个异常值的序列构建两个常规控制图,第一个是不含异常值常规控制图,第二个是含个异常值常规控制图,两个控制图的监测结果如图所示.图不含异常值的常规控制图(左)与含个异常值的常规控制图(右)从图可以看出,两个控制图的中间实线为控制中心线,上下两条虚线为控制上下限.由于图左图中的受控过程不含异常值,所有的点均在上下控制限以内,没有出现“虚假报警”或者“漏报警”现象,监控成功.而图右图中的受控过程含有一个异常值,但从图右图中可以看出,异常值在上下控制限内,控制图并没有对异常值进行报警,而是将异常值误判正常值,监控失败.显然,监控失败的原因是由于控制限对异常值敏感,异常值的存在导致控制上下限间距拉大所致.因此,对常规控制图的控制中心及控制限进行稳健改进具有重要意义.稳健常规控制图构建的基本原理寻找均值、标准差的稳健组合估计量来改进控制中心及控制限,对构造稳健控制图具有重要作用.研究表明构建稳健控制图是一个较为复杂的系统工作,并不是仅仅寻找均值与标准差的稳健估计量.由于均值与标准差在控制图中是有机整体,若将两者分开研究有可能会导致控制图过于稳健或者缺乏稳健,其后果是:过于稳健会出现虚假报警,缺乏稳健会出现漏报警现象.常用的均值稳健估计量有:中位数(M e d i a n)、切尾均值(T r i mm e dM e a n)、缩尾均值(W i n s o r i z e dM e a n)、三均值(T r i p l eM e a n)以及H o d g e sL e h m a n n(HL、HL、HL)等.常用的标准差稳健估计量有:缩尾标准差(WS D)、中位绝对离差(M e d i a na b s o l u t ed e v i a t i o n)、平均绝对离差(M e a na b s o l u t ed e v i a t i o n)、四分位数间距(I Q R)以及S h a m o s估计量等.本研究经过反复模拟实验比较后,采用昌世凯(C h a n s e o k)等 提出的H o d g e s L e h m a n n估计量作为均值稳健估计量,选取中位绝对离差(MA D)作为标准差稳健估计量,该两个估计量组合在一起,称为HLMA D组合估计量.研究表明,采用该组合估计量来稳健改进常规控制图能得到相对最优的稳健监测效果.H o d g e s L e h m a n n估计量具有种形式,可分别简写为HL、HL 及HL,其表达式分别为:上 饶 师 范 学 院 学 报 (第 卷)HLm e d i a nij(XiXj),HLm e d i a nij(XiXj),HLm e d i a n(i,j)(XiXj),在此,经比较后选HL,其中位绝对离差(MA D)表达式为:MADME D|XME D(X)|,将该两个估计量作为控制图的稳健组合估计量来改进控制中心与控制限,可得到如下稳健控制图(图):常规控制图U C L:(x)k(x)(x)C L:(x)L C L:(x)k(x)(x)异常值干扰稳健控制图U C L:m e d i a nij(xixj)ME D xME D(x)C L:m e d i a nij(xixj)L C L:m e d i a nij(xixj)ME D xME D(x)图基于HLMA D组合估计量的稳健常规控制图构建原理模拟研究下面通过数值模拟仿真研究来说明稳健改进的可行性和有效性.分别模拟随机产生样本量为、的标准正态分布N(,)随机数,目的是从小样本、中样本和大样本三种情形来观测改进效果.三种情形的污染率分别为、,即分别覆盖轻污染、中污染、重污染.需要说明的是,在小样本 的情形下,当污染率为 时,异常值不足个,为保证各样本量、各污染率下至少有个异常值,本研究将小样本 的轻污染率由原来的 增加到,其他情形不变.污染分布选择均匀分布,污染分布构造过程如()式所示:P()N(n;,)U n i f(n;,)()其中为污染率,N(n;,)表示标准正态分布,U n i f(n;,)表示最小值为、最大值为的均匀分布.根据研究需要,有时要产生负异常值,这时可将污染分布()式变形为()式:P()N(n;,)U n i f(n;,)()基于()式和()式产生的随机数,得到各样本量、各污染率下的监控效果如表所示.从表可以看出,当样本量为 时,常规控制图在各污染率下的监控正确率全部为,稳健控制图为 .当样本量为 与 时,常规控制图与稳健控制图的报警正确率均随着污染率的增加而减少,但稳健控制图报警正确率一直显著高于常规控制图且数值相对稳定.表不同样本量不同污染率下的常规与稳健控制图监控效果对比表样样本本量量污污染染率率()异异常常点点个个数数常常规规控控制制图图正正确确率率()稳稳健健控控制制图图正正确确率率()实实证证检检验验为进一步检验本研究所提出的基于HL MA D组合估计量所改进的稳健控制图对含异常值序列的监第期李熠玲,等:基于HLMA D组合估计量的常规控制图稳健改进控效果,本文选取 年月 日 年月 日泰山石油(代码:S Z)收益率作为受控对象.数据来源于英为财情网站(h t t p s:/c n i n v e s t i n g c o m),有效样本量为 .泰山石油数据的探索性分析结果如图所示.图显示:第一,泰山石油的收盘价图波动幅度很大;第二,泰山石油的收益率图呈现异方差性且存在多个绝对值较大的数,这一特征与本研究目标相吻合.至于绝对值较大的数是否为异常值,则需要通过统计检验进行识别.图泰山石油股票收盘价图(左)与收益率图(右)采用王志坚、汪志红和王斌会等提出的时间序列异常值稳健检测法 对泰山石油样本数据进行异常值检测,结果为表所示.表泰山石油收益率异常值检测结果表序序号号日日期期类类型型序序号号日日期期类类型型 /I O /I O /AO /I O /AO /AO /I O /AO /I O /I O从表可以看出,稳健检测法检测到异常值共 个,其中,个为I O型,个为AO型.而泰山石油收益率的正态性检验结果显示为:W ,pv a l u e e ,即拒绝服从正态分布的原假设.可见,异常值的存在导致序列违背了受控过程满足正态分布假定.接下来,采用常规控制图与稳健改进控制图分别对收益率序列进行监控,试图通过监控将序列中异常值“报警”出来,监控结果见图.从图可以看出,常规控制图出现个异常值报警,序号分别是:、.稳健控制图出现 个异常值报警,序号分别是:、.常规控制图漏报警个异常值,报警正确率,稳健控制图所报警异常值个数与稳健检测法检测到异常点个数一致,报警正确率 .泰山石油收益率的实证结果进一步验证了稳健控制图的可行性与有效性.上 饶 师 范 学 院 学 报 (第 卷)图常规控制图(左)与稳健控制图(右)监控结果图相相对对效效比比较较一个估计量的有效性,通常用相对效率来进行评价,相对效率的概念在各种估计量进行评估比较时非常有用.关于的相对效定义为:R E(|)V a r()V a r()()()式中分子为估计量的最小方差,分母为实际方差,通常当该比值为或者接近时,才认为该估计量是有效的,或者称为相对最优.泰山石油样本量为 ,据此算出泰山石油收益率常用的位置参数与尺度参数估计量的相对效,结果如表所示.表各估计量的相对效比较估估计计量量相相对对效效估估计计量量相相对对效效位置估计m e a nm e d i a n HL HL HL 尺度估计s dr a n g e m a d s h a m o s 注:表中“”表示不存在从表可以看出,位置估计量相对效最小的是中位数m e d i a n,尺度估计量相对效最小的是极差r a n g e.显然中位数是位置估计量的一个很好的稳健估计,而极差并不是尺度估计量的一个好的估计量,但表显示,极差相对效最小.本文通过对比研究发现,由于控制图的构建涉及到样本均值与标准差两个统计量,仅仅考虑单个统计量的相对效难以构建一个监控效果满意的控制图.因此,需要将位置与尺度统计量两者结合起来构建组合估计量才能得到一个相对最优的稳健控制图.因此,在控制图的稳健估计量的选择问题上,研究者不能仅仅用单个估计量的相对效比较来作为稳健估计量的唯一选取标准,从某种意义上讲,此发现亦是本研究的学术贡献之一.第期李熠玲,等:基于HLMA D组合估计量的常规控制图稳健改进结结束束语语常规休哈特控制图的控制限表达式为,其中的样本均值、标准差是不稳健.离群值会导致控制中心、控制限发生改变,以此设置的控制图容易引起“漏报警”现象,从而严重影响人们对受控过程质量的判断.本研究通过比较并选取均值的稳健估计量HL 与标准差的稳健估计量MA D作为稳健组合估计量构建了稳健常规控制图,模拟与实证分析均表明本研究构建的稳健控制图能有效地对异常值进行监控.另外,本研究发现,在稳健控制图构建过程中,均值与标准差的稳健估计量不应该分开选取,而要作为一个整体来考虑才能达到更好的效果,这一点在相对效的比较研究中得到了进一步论证.参考文献:V I J AYABV,MUR T YSNS An e wa p p r o a c ht or o b u s te c o n o m i cd e s i g no fc o n t r o l c h a r t sJ A p p l i e ds o f tc o m p u t i n gj o u r n a l,():周纯光,邹长亮,王兆军基于小波的预分析稳健控制图(英文)J应用概率统计,():吴纯杰,王兆军基于C o r n i s h F i s h e r展开的稳健综合控制图设计J应用概率统计,():A B D U LSS,R E Z ABK,HE N G S OON G R o b u s te c o n o m i cs t a t i s t i c a ld e s i g no fX b a rc o n t r o lc h a r tJ I n t e r n a t i o n a lj o u r n a l o fp r o d u c t i o nr e s e a r c h,():S A L I H OD,AHME D M G,AHME D M A A m a t h e m a t i c a lm o d e l f o r t h eo p t i m a l r o b u s td e s i g no f c a u s es e l e c t i n gc o n t r o l c h a r t sJ E u r o p e a n j o u r n a l o f i n d u s t r i a l e n g i n e e r i n g,():HA B S HAH M,A S HKANS R o b u s tm u l t i v a r i a t e c o n t r o l c h a r t s t od e t e c t s m a l l s h i f t s i nm e a nJ M a t h e m a t i c a l p r o b l e m si ne n g i n e e r i n g,(p t):F ANSKS,HUANG H K,CHANGYJ R o b u s tm u l t i v a r i a t e c o n t r o l c h a r t f o ro u t l i e rd e t e c t i o nu s i n gh i e r a r c h i c a l c l u s t e rt r e e i nSW J Q u a l i t ya n dr e l i a b i l i t ye n g i n e e r i n g i n t e r n a t i o n a l,():A S OKAN M V,J AYA S ANKA RV R o b u s t c o n t r o l c h a r t s f o rm o n i t o r i n gp r o c e s sm e a no f p h a s e m u l t i v a r i a t e i n d i v i d u a l o b s e r v a t i o n sJ I n t e r n a t i o n a l j o u r n a l o fq u a l i t y,s t a t i s t i c s,&r e l i a b i l i t y,:Z HAN GC,CHE NN,Z OUCL R o b u s tm u l t i v a r i a t e c o n t r o l c h a r tb a s e do ng o o d n e s s o f f i t t e s tJ J o u r n a l o f q u a l i t yt e c h n o l o g y,():AN G E L L Y SPAG,R I S T E RBP,R O B E R T OJ,e t a l A p p l i c a t i o no f r o b u s tm u l t i v a r i a t e c o n t r o l c h a r tw i t hW i n s o r i z e dM e a n:ac a s es t u d yJ J o u r n a l o f i n d u s t r i a l e n g i n e e r i n g i n t e r n a t i o n a l,(s u p p l):R A J I IA,A B B A SN,A B U J I YA MR,e t a l R o b u s tm u l t i v a r i a t e s h e w h a r t c o n t r o l c h a r t b a s e do n t h e s t a h e l d o n o h o r o b u s t e s t i m a t o ra n dm a h a l a n o b i sd i s t a n c e f o rm u l t i v a r i a t eo u t l i e rd e t e c t i o nJ M a t h e m a t i c s,():C A B ANAE,L I L L ORE R o b u s tm u l t i v a r i a t e c o n t r o l c h a r tb a s e do ns h r i n k a g e f o r i n d i v i d u a l o b s e r v a t i o n sJ J o u r n a l o fq u a l i t yt e c h n o l o g y,():S A B AHNO H,C E L ANOG M o n i t o r i n g t h em u l t i v a r i a t e c o e f f i c i e n t o f v a r i a t i o n i np r e s e n c eo f a u t o c o r r e l a t i o nw i t hv a r i a b l ep a r a m e t e r sc o n t r o l c h a r t sJ Q u a l i t yt e c h n o l o g y&q u a n t i t a t i v em a n a g e m e n t,():王志坚稳健残差控制图的构建及在金融时序中的应用J数理统计与管理,():王志坚,王斌会稳健A R模型的构建及比较研究J数学的实践与认识,():王志坚稳健A R模型的构建及其在金融时序中的应用J统计与信息论坛,():S HA R I A T IN R o b u s t r e s i d u a l c o n t r o l c h a r t f o r c o n t a m i n a t e dt i m e s e r i e s:as o l u t i o nt ot h ee f f e c t so f o u t l i e r d r i v e np a r a m e t e rm i s e s t i m a t i o no nt h ec o n t r o l c h a r tp e r f o r m a n c eJ C o mm u n i c a t i o n s i ns t a t i s t i c st h e o r ya n dm e t h o d s,():上 饶 师 范 学 院 学 报 (第 卷)S A L AH M,E N G Y M,S HE R E E NAL P r o f i l em o n i t o r i n go f r e s i d u a l s c o n t r o l c h a r t su n d e r g a mm a r e g r e s s i o nm o d e lJE a s t e r n E u r o p e a n j o u r n a l o f e n t e r p r i s e t e c h n o l o g i e s,():HA F I ZZN,MA R I TS,MUHAMMA DR,e t a l Q u a l i t yq u a n d a r i e s:as t e p w i s ea p p r o a c hf o rs e t t i n gu par o b u s ts h e w h a r t l o c a t i o nc o n t r o l c h a r tJ Q u a l i t ye n g i n e e r i n g,():NA D I AS,S HAH I DK T h eEWMAc o n t r o l c h a r tb a s e do nr o b u s t s c a l ee s t i m a t o r sJ P a k i s t a n j o u r n a l o f s t a t i s t i c sa n do p e r a t i o nr e s e a r c h,():吴纯杰,魏一青,郁淼淼混合正态分布下的似然比累积和控制图及改进J中国科学(数学),():庄芳,郁淼淼,吴纯杰带参数估计的EWMA方差控制图稳健性分析J系统科学与数学,():王志坚,郑健松稳健S h e w h a r t控制图的构建及应用J统计与决策,():苏拥英,王志坚常规过程控制图的敏感性分析及稳健设计J统计与决策,():KAOSC Ar o b u s t s t a n d a r dd e v i a t i o nc o n t r o l c h a r tb a s e do ns q u a r eAe s t i m a t o rJ Q u a l i t ya n dr e l i a b i l i t ye n g i n e e r i n gi n t e r n a t i o n a l,():CHAN S E OKP,HA EWONK,M I N W I n v e s t i g a t i o no f f i n i t e s a m p l ep r o p e r t i e so f r o b u s t l o c a t i o na n ds c a l e e s t i m a t o r sJ C o mm u n i c a t i o n s i ns t a t i s t i c s s i m u l a t i o na n dc o m p u t a t i o n,():王志坚,王斌会稳健改进的AO型异常点检测法在金融时序中的应用J数理统计与管理,():汪志红,王志坚,王斌会一种新的时间序列T C型异常值稳健检测法J统计与决策,():汪志红,王志坚,王斌会时间序列新息异常值稳健诊断新方法J统计与决策,():王志坚,罗舒琪,王斌会基于稳健C o o k距离的时间序列异常值诊断J统计与决策,():R o b u s t I m p r o v e m e n t o fC o n v e n t i o n a lC o n t r o lC h a r t sB a s e do nH LMA DC o m b i n e dE s t i m a t o r sL IY i l i n g,R E NK a i l i a n g,WA N GZ h i j i a n(S c h o o l o fS t a t i s t i c s&M a t h e m a t i c s,G u a n g d o n gU n i v e r s i t yo fF i n a n c e&E c o n o m i c s,G u a n g z h o uG u a n g d o n g ,C h i n a;C o l l e g eo fL i b e r a lA r t sa n dS c i e n c e s,U n i v e r s i t yo fC o n n e c t i c u t,C o n n e c t i c u t ,U S A)A b s t r a c t:T h ec o n t r o l c e n t e r a n dc o n t r o l l i m i t o f c o n v e n t i o n a l c o n t r o l c h a r t s a r en o t r o b u s t a n da r e s u s c e p t i b l e t oo u t l i e r s,w h i c hm a y l e a dt om o n i t o r i n gf a i l u r e s I nt h i sp a p e r,t h eHLMA Dc o m b i n e de s t i m a t o r i su s e dt o i m p r o v e t h ec o n t r o l c e n t e ra n dc o n t r o l l i m i t i nar o b u s tw a y,a n dar o b u s t c o n t r o l c h a r t i sc o n s t r u c t e d B o t hn u m e r i c a l s i m u l a t i o na n df i n a n c i a l e m p i r i c a l a n a l y s i s s h o wt h a t c o m p a r e dw i t ht r a d i t i o n a l c o n t r o l c h a r t s,r o b u s t c o n t r o l c h a r t s c a nb e t t e r i d e n t i f ya b n o r m a l p h e n o m e n a i nc o n t r o l l e dp r o c e s s e sa n do v e r c o m et h ep h e n o m e n o no fm i s s i n ga l a r m s,a l s oh a v eb e t t e re f f e c t su n d e r l a r g e s a m p l e s a n dh i g hc o n t a m i n a t e r a t e s I na d d i t i o n,t h r o u g hc o m p a r a t i v e s t u d i e so n t h e r e l a t i v ee f f i c i e n c yo f e s t i m a t o r s,i t i s f o u n dt h a t i nt h es e l e c t i o no f r o b u s t e s t i m a t o r s f o r c o n t r o l c h a r t s,c o m b i n e de s t i m a t o r sh a ss i g n i f i c a n t s u p e r i o r i t yo v e r i n d i v i d u a l e s t i m a t o r s K e yw o r d s:s t a t i s t i c a l p r o c e s sc o n t r o l;HLMA Dc o m b i n e de s t i m a t o r;r o b u s t e s t i m a t i o n;r e l a t i v ee f f i c i e n c y第期李熠玲,等:基于HLMA D组合估计量的常规控制图稳健改进