常用的GRR分析方法有极差法、均值极差法、方差分析法,下面以一个具体的例子来说明基于方差分析法的GRR指标计算过程。
案例背景
· 被测物料: 金属棒端面直径(单位:mm)
· 测量系统: 直径测量软件
· 操作者: 3人(A, B, C)
· 物料样本: 10个
· 试验流程: 每个操作者测量每片物料各3次
第一步:收集数据
操作者A的测量数据
物料 | 试验1 | 试验2 | 试验3 | 求和 |
1 | 225.00 | 225.22 | 224.89 | 675.11 |
2 | 234.01 | 234.02 | 234.20 | 702.23 |
3 | 236.63 | 236.42 | 236.78 | 709.83 |
4 | 226.49 | 226.50 | 226.32 | 679.31 |
5 | 239.45 | 239.41 | 239.43 | 718.29 |
6 | 225.12 | 225.22 | 224.98 | 675.32 |
7 | 234.45 | 234.38 | 234.29 | 703.12 |
8 | 236.48 | 236.33 | 236.41 | 709.22 |
9 | 226.53 | 226.31 | 226.67 | 679.51 |
10 | 239.43 | 239.18 | 239.35 | 717.96 |
操作者B的测量数据
物料 | 试验1 | 试验2 | 试验3 | 求和 |
1 | 225.03 | 225.17 | 224.87 | 675.07 |
2 | 234.48 | 234.33 | 234.25 | 703.06 |
3 | 236.45 | 236.34 | 236.34 | 709.13 |
4 | 226.55 | 226.16 | 226.67 | 679.38 |
5 | 239.42 | 239.18 | 239.45 | 718.05 |
6 | 225.01 | 225.18 | 224.93 | 675.12 |
7 | 234.12 | 234.07 | 234.20 | 702.39 |
8 | 236.55 | 236.44 | 236.78 | 709.78 |
9 | 226.44 | 226.45 | 226.32 | 679.21 |
10 | 239.45 | 239.40 | 239.34 | 718.19 |
操作者C的测量数据
物料 | 试验1 | 试验2 | 试验3 | 求和 |
1 | 225.27 | 224.83 | 225.23 | 675.33 |
2 | 234.11 | 234.00 | 235.25 | 702.36 |
3 | 236.11 | 236.39 | 236.32 | 708.82 |
4 | 226.32 | 226.40 | 226.56 | 679.28 |
5 | 239.34 | 239.56 | 239.31 | 718.21 |
6 | 225.17 | 225.22 | 224.89 | 675.28 |
7 | 234.13 | 234.11 | 234.20 | 702.44 |
8 | 236.56 | 236.43 | 236.78 | 709.77 |
9 | 226.47 | 226.45 | 226.32 | 679.24 |
10 | 239.45 | 239.44 | 239.40 | 718.29 |
第二步:计算关键统计量
1. 计算每个操作者30次测量值的和XA、XB、XC
· XA = 675.11 + 702.23 + … + 679.51 + 717.96 = 6969.9
· XB = 675.07 + 703.06 + … + 679.21 + 718.19 = 6969.38
· XC = 675.33 + 702.36 + … + 679.24 + 718.29 = 6969.02
2. 计算每片物料9次测量值的和Xp1到Xp10
· Xp1 = 675.11 + 675.07 + 675.33 = 2025.51
· Xp2 = 702.23 + 703.06 + 702.36 = 2107.65
· Xp3 = 709.83 + 709.13 + 708.82= 2127.78
· Xp4 = 679.31 + 679.38 + 679.28 = 2037.97
· Xp5 = 718.29 + 718.05 + 718.21 = 2154.55
· Xp6 = 675.32 + 675.12 + 675.28 = 2025.72
· Xp7 = 703.12 + 702.39 + 702.44 = 2107.95
· Xp8 = 709.22 + 709.78 + 709.77 = 2128.77
· Xp9 = 679.51 + 679.21 + 679.24 = 2037.96
· Xp10 = 717.96 + 718.19 + 718.29 = 2154.44
3. 计算所有操作者所有物料90次测量值的总和X
· X = 6969.9 + 6969.38 + 6969.02 = 20908.3
4. 计算基准变差Rb,物料变差Rp,操作者变差Ro,重复性变差Re,总变差Ra
· Rb = X * X / 10 / 3 / 3 = 4857300.099
· Rp = SUMSQ(Xp1, Xp2, …, Xp10) / 3 / 3= 4860142.616
· Ro = SUMSQ(XA,XB,XC) / 10 / 3 = 4857300.112
· Re = SUMSQ(XA1, XA2, …, XC10) / 3 = 4860143.176
· Ra = SUMSQ(X1,X1,X3, …, X90) = 4860144.282
5. 计算实际P值XP
· XP = F.DIST.RT(((Ra - Rb) – (Rp - Rb) – (Ro - Rb) – (Ra – Re)) / 18 / ((Ra – Re) / 60), 18, 60) = 0.076
由于XP > 0.05,因此不存在交互作用
6. 方差分析
· 物料离差平方和SSp = Rp - Rb = 2842.52
· 操作者离差平方和SSo = Ro - Rb = 0.01
· 重复性离差平方和SSe = Ra - Rb - SSp - SSo= 1.65
· 物料均方MSp = SSp / 9 = 315.835
· 操作者均方MSo = SSo / 2 = 0.007
· 重复性均方MSe = SSe / 78= 0.021
第三步:计算方差分量
计算各变差
1. 重复性EV
· EV = MSe = 0.0212
2. 再现性AV
· 由于MSo < MSe
所以AV = 0
3. 测量系统GRR
· GRR = EV + AV = 0.0212
4. 物料间变差PV
· PV = (MSp - MSe) / 3 / 3 = 35.09
5. 总变差TV
· TV = GRR + PV = 35.11
第四步:计算百分比贡献
1. 基于总变差的百分比

2. 数据分类数

第五步:结论判定
判定标准对照
指标 | 结果 | 标准 | 判定 |
%GRR(TV) | 2.46% | <10%:可接受 | 可接受 |
ndc | 57 | ≥5:可接受 | 可接受 |