CNAS能力验证——稳健统计方法（算法A）计算稳健均值和标准差

1、案例介绍表E.1给出一轮能力验证计划的23个结果其中5个结果标示为“小于”某数量。由算法A得到的稳健均值x*)和标准差s*显示了三种不同的计算方法忽略“”作为定量数据进行分析删除带有“”的数据忽略“”数据的0.5倍作为定量结果的估计。在各种情形下超出可接受界限的结果标记为“#”。这假设定量部分超出x*±3s*的任何结果将被评定为“不可接受“行动信号。对于带有“”或“”号的结果的评定能力验证提供者也可使用其他规则。表E.1 包含截断结果的样本数据集以及调整结果的三个选项参加者结果忽略“”删除含“”的数据0.5×含“”的数据A1010——5B1010——5C12121212D19191919E2020——10F20202020G23232323H23232323J25252525K25252525L26262626M28282828N28282828P3030——15Q28282828R29292929S30303030T30303030U31313131V32323232W32323232Y454545#45Z5050#——25总结结果的数量23231823x*26.0126.8123.95s*7.235.298.602、算法A算法 A 来自 GB/T 6379.5。应用此算法计算得到数据平均值和标准差的稳健值。稳健性是估计算法的特点而不是其产生的估计值的特点因此严格来说称由此算法计算的平均值和标准差是稳健的是不确切的。然而为避免使用繁琐的术语“稳健均值”和“稳健标准差”应理解为利用稳健算法计算的总体均值和总体标准差的均值估计。从一个特定检测中得到的结果总数为p。按递增顺序排列p个检测数据表示为x1 ,x2 ,Lxi,L,xp。这些数据的稳健平均值和稳健标准差记为x* 和s*。计算x* 和s* 的初始值如下med 表示中位数根据以下步骤更新x* 和s* 的值。计算对每个xi(i 1,2,L,p) 计算再由下式计算x* 和s* 的新的取值其中求和符号对i求和。稳健估计值x* 和s* 可由迭代计算得出例如用已修改数据更新x* 和s* 直至过程收敛。当稳健标准差的第三位有效数字和稳健平均值相对应的数字在连续两次迭代中不再变化时即可认为过程是收敛的。3、示例中数据求取稳健均值和标准差的Matlab算法实现%% 代码中data为待求取稳健均值和标准差的一维向量数据可从excel中获取 x0median(data); s01.483*median(abs(data-x0)); x1zeros(length(data),1); while 1 delta1.5*s0; for i1:length(data) if data(i)x0-delta x1(i)x0-delta; elseif data(i)x0delta x1(i)x0delta; else x1(i)data(i); end end xnewmean(x1); snew1.134*std(x1); if abs(xnew-x0)1e-3abs(snew-s0)1e-3 break; else x0xnew; s0snew; end end fprintf(稳健均值 %.2f\n, x0); fprintf(稳健标准差 %.2f\n, s0);4、运行结果以第一列为例结果和表格中一致5、参考文献CNAS-GL0022018《能力验证结果的统计处理和能力评价指南》GB/T 28043-2019 利用实验室间比对进行能力验证的统计方法