1308 04异常值.pdf

异常值判断 格拉布斯检验法 Grubbs test method GB 4883-1985 概述 • 一组重复测量数据中，若个别数据与大多 数数据相差很大，暂定这个(些)数据为“可 疑值” 。 • 不能凭主观意愿剔除可疑值，只能根据客 观的统计检验方法，判断它们是否超过方 法规定的判据(临界值) 。 • 若超过判据，则判可疑值为“异常值”，将 其剔除，不参与任何计算。 • 若不超出判据，则不能判为异常值，不可 剔除，仍为正常值。 实测数据举例 • 重复测量某个量10次，按测量顺序获得以 下一组10个数据：8.2、5.4、14.0、7.3、4.7 、 9.0、6.5、10.1、7.7、6.0 。(每数前加一常数亦同) •计算两个统计量：①算术平均值 x =⎯ 7.89； ②单个测得值x 的实验标准偏差s(x )=2.704 。 i i • 原数据序列不能看出哪个数大，哪个数小。 将它们按照大小重新排列，得：4.7、5.4、 6.0、6.5、7.3、7.7、8.2、9.0、10.1、14.0。 计算可疑值的G值 • 先找出测量数据中“两个”最可能的可疑值： 最小值4.7、最大值14.0。 • 分别计算它们与算术平均值 x⎯之差，选择 其中差值较大者为可疑值。 • x⎯与4.7之差：7.89−4.7=3.19 ；14.0与 x⎯之差： 14.0−7.89=6.11 。因为6.11＞3.19，所以选择 序号为10的14.0为可疑值。 • 计算可疑值的G 值。公式：G =(x −⎯x)/s(x ) 。 i i i i • G =(x −⎯x)/s(x )=(14.0−7.89)/2.704=2.260 。 10 10 i 格拉布斯表临界值G (n) p • G (n)与测量次数n及置信概率p 有关，其中p p 和检出水平α互补：p +α=1 。α值自定。 • 若α取0.01 ，则p =0.99 ；若α取0.10 ，则p = 0.90；通常α取0.05，那么p =0.95 。 • 格拉布斯表横行为p ，竖列为n 。本例横竖 p =0.95 、n=10相交得临界值：G95(10)=2.176 。 • 判断。若G ＞G (n) ，异常值，剔除；否则 i p 非异常值，不剔除。本例Gi=G10=2.260 ＞ G95(10)=2.176，所以14.0为异常值，剔除。 格拉布斯表(部分) n p=0.95 p =0.99 n p=0.95 p =0.99 3 1.135 1.155 11 2.234 2.485 4 1.463 1.492 12 2.285 2.550 5 1.672 1.749 13 2.331 2.607 6 1.822 1.944 14 2.371 2.659 7 1.938 2.097 15 2.409 2.705 8 2.032 2.231 16