Gauge repeatability and reproducibility studies are significant to quality improvement and quality control. The approaches are always applied to determine the capability of the measurement system. Much of the literature in this field mainly focuses on univariate and multivariate measurement systems. However, the state-of-the-art methods are not appropriate when the quality of a product is characterized by a profile. Therefore, this paper proposes a method for the measured values which can be characterized by a simple linear profile. In addition, the slopes and intercepts of these profiles often vary due to measurement error. Thus, the simple linear profile gauge studies can be considered as a two-response (slope and intercept) problem. X-values transformation is used to make the slope and intercept of each profile independent. ANOVA is utilized to estimate the variance component of measurement error and other sources of variation. Then, the criteria precision to tolerance ratio and percent R&R are introduced to assess the simple linear profile measurement system capability. Finally, the proposed approach is applied to the spring length and elasticity measurement which demonstrates how to implement the method.
1. Introduction
Gauge repeatability and reproducibility (GR&R) studies are usually adopted for determining the capability of the measurement systems. Burdick et al. [1] and Montgomery [2] proposed that one of the objectives of measurement systems capability analysis (MSCA) is to determine whether the gauge is capable. Evaluating the capability of measurement system is necessary for other researches, such as process capability analysis, control charts, and design of experiment. In practice, a reliable measurement system is significant to quality improvement and quality control.
Burdick et al. [1] reviewed the measurement system capability analysis. In their paper, the criteria, typical model, and GR&R experiment were all introduced. Attribute data are widely applied in industry; Lyu and Chen [3] evaluated the R&R of a measurement system for attribute data based on the generalized linear models (GLMs). Some papers considered confidence intervals for gauge capability studies. Borror et al. [4] presented two methods to construct confidence intervals for variance components: one is the restricted maximum likelihood (REML) method and the other is the modified large sample (MLS) method. Gong et al. [5] proposed an approach by applying unweighted sums of squares method when confidence intervals for unbalanced two-factor gauge study were considered.
In practice, many products need several critical characteristics to describe their quality. For these cases, the univariate MSCA method may not be appropriate especially when the measured characteristics have some correlations. Thus, numerous studies in this field considered multivariate MSCA. An approach to analyze the two-dimensional GR&R of imbalance measurement was introduced by Sweeney [6]. He concluded that the variation can be underestimated via treating such data as one-response. Majeske [7] proposed the method and the criteria for multivariate MSCA by regarding the measured values as a vector and assuming that the measured values follow a multivariate normal distribution. MANOA is applied to estimate the variance-component matrices.
Gauge repeatability and reproducibility studies are widely used in practice. Li and Al-Refaie [8] utilized the define-measure-analyze-improve-control (DMAIC) procedure to improve the quality system involving measurements. Erdmann et al. [9] gave an example of GR&R in a hospital, and the experiment is about an ear thermometer for temperature measurements.
The existing research proposed many methods for univariate and multivariate MSCA. However, the measured values in some situations can be characterized by profiles. The measurement for the elasticity and length of springs is a good example. In this case, the elasticity of spring is a simple linear profile of an independent variable: the length of spring. In this paper, an approach for the situation where the measured values can be expressed as a simple linear profile is proposed.
The remainder of this paper is organized as follows. A brief review of gauge R&R is presented in Section 2. Section 3 proposes the approach and criteria for simple linear profile gauge R&R. A case study is introduced in Section 4. The last section contains our concluding remarks.
2. The Model and Criteria for GR&R
Typically, gauge studies are always designed with two factors, i parts and j operators, and each operator measures each part k times. The measured values are represented as Mijk. These values are always treated as true values by the operators because they usually contain errors. The typical model with a two-factor design for variance analysis of Mijk is as follows, see Burdick et al. [1] and Montgomery [2]:
(1)Mijk=μ+pi+oj+(op)ij+εijk,
where i=1,2,…,p, j=1,2,…,o, and k=1,2,…,r. μ is the true value. pi, oj, (op)ij, and εijk are all independent random variables which represent the effect of part, operator, operator-part interaction, and random error. Assume that they are all normally distributed, pi~N(0,σp2), oj~N(0,σo2), (op)ij~N(0,σ(op)2), and εijk~N(0,σe2). The variance of Mijk is given as follows:
(2)σMijk2=σp2+σo2+σ(op)2+σe2,σG2=σo2+σ(op)2+σe2,
where σG2 is the variance component which represents the effect of gauge.
Two criteria which are always utilized to assess the gauge capability are introduced in AIAG [10]. One is the ratio between the measurement precision estimate and the tolerance of the parts, P/T%. The other is the ratio of the measurement precision estimate and the process width, R&R%. The criteria are formulated as follows:
(3)PT%=KσGUSL-LSL×100%,R&R%=σGσG2+σp2×100%,
where USL and LSL are upper specification limit and lower specification limit of the part. K is either 5.15 or 6. According to AIAG [10], the measurement system is unacceptable if the value exceeds 30%, and it is acceptable if the value is lower than 10%. If the value is between 10% and 30%, the measurement system needs to be improved. References considering these criteria include Wheeler and Lyday [11] and Woodall and Borror [12].
3. The GR&R Method for Simple Linear Profile
In practice, the quality of some products can be described by simple linear profiles when they have a collection of critical characteristics to be measured. In simple linear profile gauge study, a two-factor (i parts and j operators) design is considered. Each operator measures all of the parts k times. The measured values are characterized by a simple linear profile:
(4)Yijk=A0(ijk)+A1(ijk)Xijk+εijk,
where i=1,2,…,p, j=1,2,…,o, and k=1,2,…,r. A0(ijk) is the intercept, and A1(ijk) is the slope. The random variable εijk is independent and normally distributed, εijk~(0,σ2). Due to measurement error, the profiles always vary, and the variation is reflected by the slopes and the intercepts. As is shown in Figure 1, where one profile represents one time measurement, the slopes and intercepts of these profiles are different.
The differences of the profiles.
Thus, when considering measurement error, (4) is given as follows:
(5)Yijk=[A0+piA0+ojA0+(op)ijA0+εijkA0]+[A1+piA1+ojA1+(op)ijA1+εijkA1]Xijk+εijk,
where A0 and A1 are constants. piA0, ojA0, (op)ijA0, and εijkA0 are independent random variables which represent the effect of part, operator, operator-part interaction, and the random error on A0(ijk), respectively. piA1, ojA1, (op)ijA1, and εijkA1are independent random variables which represent the effect of part, operator, operator-part interaction, and the random error on A1(ijk), respectively.
3.1. The Transforming Model
For simple linear profile gauge study, A0(ijk) and A1(ijk) can be obtained easily when the measured results are profiles directly. However, it is necessary to estimate the two parameters in some cases. The least-square method, see Johnson and Wichern [13], is used to estimate parameters and fit profiles. The fitted profile is the following:
(6)Y^ijk=A^0(ijk)+A^1(ijk)Xijk,
where Y^ijk is the fitted value. A^0(ijk) and A^1(ijk) are the least-square estimates of A0(ijk) and A1(ijk)which are normally distributed. A^0(ijk) and A^1(ijk) are correlative, and the covariance between them is σ012.
Kim et al. [14] introduced a method which can make the intercept and the slope of each profile independent via transforming the X-values, so that the model can be simplified. An alternative form of the model in (6) is the following:
(7)Yijk=B0(ijk)+B1(ijk)Xijk*+εijk,
where Xijk*=Xijk-X-, B0(ijk)=A0(ijk)+A1(ijk)X-, and B1(ijk)=A1(ijk). In this situation, the least-squares estimator of B0(ijk) is B^0(ijk)=A^0(ijk)+A^1(ijk)X-, and the least-squares estimator of B1(ijk) is B^1(ijk)=A^1(ijk). Both of B^0(ijk) and B^1(ijk) are normally distributed, and the covariance between them is zero, so B^0(ijk) and B^1(ijk) for each profile are independent. A separate gauge R&R analysis can be applied to B^0(ijk) and B^1(ijk), respectively.
3.2. The Model for Variance Analysis
In the alternative form, the measurement error is reflected by B^0(ijk) and B^1(ijk). When a two-factor design experiment is considered, the model for variance analysis of measurement error on B^0(ijk) and B^1(ijk) is the following:
(8)B^0(ijk)=μ0+piB0+ojB0+(op)ijB0+εijkB0,B^1(ijk)=μ1+piB1+ojB1+(op)ijB1+εijkB1,
where μ0 and μ1 are constants. piB0 and piB1 are independent random variables which represent the part effect on B^0(ijk) and B^1(ijk), respectively. ojB0 and ojB1 are independent random variables which represent the operator effect on B^0(ijk) and B^1(ijk), respectively. (op)ijB0 and (op)ijB1 are independent random variables which represent the operator-part interaction effect on B^0(ijk) and B^1(ijk), respectively. εijkB0 and εijkB1 are independent random variables which represent the random error effect on B^0(ijk) and B^1(ijk), respectively. Assume that piB0, ojB0, (op)ijB0, εijkB0, piB1, ojB1, (op)ijB1, and εijkB1 are all normally distributed, where piB0~N(0,σp02), ojB0~N(0,σo02), (op)ijB0~N(0,σop02), εijkB0~N(0,σ02), piB1~N(0,σp12), ojB1~N(0,σo12), (op)ijB1~N(0,σop12), and εijkB1~N(0,σ12).
The variance of B^0(ijk) considering the measurement error is
(9)V[B^0(ijk)]=σp02+σo02+σop02+σ02.
The variance component of B^0(ijk) for the gauge is
(10)σG02=σo02+σop02+σ02.
The variance component of B^0(ijk) for the measured values is
(11)σM02=σp02+σG02.
The variance of B^1(ijk) considering the measurement error is
(12)V[B^1(ijk)]=σp12+σo12+σop12+σ12.
The variance component of B^1(ijk) for the gauge is
(13)σG12=σo12+σop12+σ12.
The variance component of B^1(ijk) for the measured values is
(14)σM12=σp12+σG12.
A two-factor ANOVA with an interaction term is applied to estimate these variance components.
3.3. The Assessment Criteria
Two responses (B^0(ijk) and B^1(ijk)) are obtained in the gauge study for simple linear profile. The univariate criteria P/T% and R&R% can be calculated for B^0(ijk) and B^1(ijk), respectively. Four equations are utilized to assess the gauge capability for simple linear profile:
(15)PT0%=KσG0USL0-LSL0×100%,R&R0%=σG0σG02+σp02×100%,PT1%=KσG1USL1-LSL1×100%,R&R1%=σG1σG12+σp12×100%.
According to AIAG [10], the measurement system is unacceptable if one of the above assessed values exceeds 30%. It is acceptable if the values are all lower than 10%. If the values are all between 10% and 30%, the measurement system needs to be improved.
4. The Introductory Case
An example of gauge study about spring measurement is introduced in this section. In the quality improvement, it is necessary to measure the elasticity and the length of spring. According to Hooke’s law, when the spring has reached a state of equilibrium, its elasticity is a simple linear profile of the amount by which the free end of the spring is displaced from its relaxed position (when it is not stretched). In this case, three operators made three measurements on each of ten springs, Xijk is the length of spring after compression or elongation, and Yijk is the different elasticity when the spring is of different length.
The least-square method is used to fit profile and estimate the parameters A^0(ijk) and A^1(ijk) after the measurement. Then, B^0(ijk) and B^1(ijk) in (7) can be obtained via transforming X-values. A two-factor ANOVA with an interaction term is utilized to estimate the variance components of B^0(ijk) and B^1(ijk), respectively. The variance components of B^0(ijk) are shown in Table 1. The variance components of B^1(ijk) are shown in Table 2.
Profile variance component of B^0(ijk).
σp02
0.007894
σo02
0.000009
σop02
0
σ02
0.000149
Profile variance component of B^1(ijk).
σp12
0.000554
σo12
0.000001
σop12
0.000002
σ12
0.000010
According to (10), (11), (13), and (14), the variance components for the gauge and the measured values of B^0(ijk) and B^1(ijk) are the following:
(16)σG02=0.000158,σM02=0.000090,σG12=0.000013,σM12=0.000567.
The upper and lower specifications on B^0(ijk) are USL0=1.3283 and LSL0=0.9717. Then, the criteria for B^0(ijk) are the following:
(17)PT0%=21.18%,R&R0%=14.03%.
The upper and lower specifications on B^1(ijk) are USL1=-0.2464 and LSL1=-0.3223. Then, the criteria for B^1(ijk) are estimated as follows:
(18)PT1%=28.30%,R&R1%=15.16%.
In this case, all of the criteria are between 10% and 30%. According to AIAG [10], the measurement system needs to be improved. The improvement will be made based on the gauge R&R experiment. In this case, the approach used to improve the capability of measurement system provides a clearer instruction. The instruction includes the precise location where the spring should be placed and the measurement operations that should be conducted. This improvement can decrease the differences among the operators.
5. Conclusion
Gauge repeatability and reproducibility studies are important to guarantee the validity of data, which is essential to other researches. Numerous existing studies in this area are about univariate and multivariate measurement systems, but these methods may not be suitable when the quality of a product should be characterized by a profile. This paper proposes an approach to assess the gauge capability when a simple linear profile is used to reflect product quality. Our proposed method can simplify the measurement problems effectively, especially for those with the multi-dimensional measured values. The example of spring measurement is presented in this paper which shows how to implement the proposed method. Further, more and more methods and criteria should be proposed to assess the capability of the measurement system.
Appendix
In this appendix, the original measured values of the introductory case are provided as Table 3.
The measured data of spring case.
Part
Operator
Xijk
17
16
15
13.5
12.5
11
Yijk
1
1
0.26
0.59
0.86
1.37
1.68
2.19
1
1
0.24
0.57
0.86
1.37
1.66
2.17
1
1
0.24
0.55
0.86
1.35
1.66
2.17
1
2
0.26
0.57
0.86
1.35
1.66
2.17
1
2
0.24
0.55
0.86
1.37
1.66
2.17
1
2
0.26
0.57
0.86
1.37
1.66
2.17
1
3
0.26
0.59
0.86
1.35
1.66
2.17
1
3
0.24
0.55
0.86
1.37
1.68
2.17
1
3
0.24
0.57
0.86
1.35
1.66
2.17
2
1
0.46
0.75
1.06
1.53
1.84
2.33
2
1
0.39
0.73
0.99
1.48
1.77
2.26
2
1
0.39
0.73
0.99
1.48
1.79
2.24
2
2
0.44
0.73
0.99
1.48
1.79
2.24
2
2
0.44
0.75
0.99
1.46
1.79
2.26
2
2
0.39
0.75
1.04
1.51
1.79
2.26
2
3
0.46
0.75
1.06
1.48
1.79
2.26
2
3
0.44
0.73
1.04
1.48
1.79
2.26
2
3
0.39
0.75
1.02
1.48
1.79
2.24
3
1
0.22
0.51
0.79
1.26
1.59
2.06
3
1
0.22
0.48
0.77
1.24
1.57
2.04
3
1
0.19
0.48
0.77
1.24
1.57
2.06
3
2
0.22
0.46
0.75
1.24
1.55
2.04
3
2
0.22
0.46
0.75
1.24
1.55
2.04
3
2
0.22
0.46
0.77
1.24
1.55
2.06
3
3
0.22
0.48
0.79
1.26
1.59
2.06
3
3
0.22
0.48
0.77
1.22
1.55
2.04
3
3
0.19
0.46
0.77
1.26
1.57
2.06
4
1
0.26
0.55
0.82
1.24
1.53
1.97
4
1
0.26
0.53
0.82
1.24
1.53
1.95
4
1
0.26
0.53
0.79
1.22
1.53
1.95
4
2
0.26
0.53
0.79
1.22
1.51
1.93
4
2
0.26
0.53
0.79
1.22
1.48
1.88
4
2
0.26
0.53
0.79
1.22
1.51
1.93
4
3
0.26
0.55
0.82
1.22
1.48
1.93
4
3
0.26
0.53
0.79
1.22
1.53
1.95
4
3
0.26
0.53
0.79
1.22
1.48
1.93
5
1
0.19
0.48
0.79
1.26
1.59
2.06
5
1
0.19
0.48
0.79
1.26
1.55
2.04
5
1
0.19
0.48
0.79
1.26
1.55
2.02
5
2
0.19
0.51
0.77
1.24
1.55
2.04
5
2
0.19
0.48
0.79
1.24
1.55
2.02
5
2
0.19
0.42
0.73
1.19
1.48
2.04
5
3
0.19
0.48
0.79
1.26
1.57
2.06
5
3
0.19
0.46
0.79
1.24
1.55
2.02
5
3
0.19
0.46
0.77
1.24
1.55
2.02
6
1
0.28
0.57
0.82
1.22
1.46
1.86
6
1
0.28
0.55
0.79
1.19
1.42
1.79
6
1
0.28
0.55
0.79
1.19
1.46
1.84
6
2
0.28
0.55
0.82
1.22
1.46
1.84
6
2
0.28
0.55
0.79
1.22
1.46
1.84
6
2
0.28
0.57
0.82
1.22
1.46
1.84
6
3
0.28
0.55
0.79
1.22
1.46
1.84
6
3
0.28
0.55
0.79
1.19
1.44
1.82
6
3
0.26
0.53
0.79
1.19
1.42
1.79
7
1
0.19
0.48
0.77
1.24
1.53
1.99
7
1
0.19
0.46
0.77
1.22
1.53
1.97
7
1
0.17
0.46
0.75
1.22
1.53
1.97
7
2
0.19
0.46
0.77
1.22
1.53
1.97
7
2
0.17
0.46
0.77
1.22
1.53
1.95
7
2
0.19
0.46
0.77
1.22
1.53
1.99
7
3
0.19
0.46
0.75
1.22
1.53
1.97
7
3
0.19
0.46
0.77
1.22
1.55
1.97
7
3
0.17
0.46
0.75
1.22
1.51
1.97
8
1
0.26
0.55
0.84
1.33
1.64
2.13
8
1
0.24
0.53
0.82
1.31
1.62
2.08
8
1
0.24
0.53
0.82
1.33
1.64
2.08
8
2
0.26
0.55
0.82
1.28
1.62
2.08
8
2
0.24
0.53
0.82
1.28
1.59
2.06
8
2
0.26
0.57
0.82
1.33
1.62
2.08
8
3
0.26
0.53
0.82
1.33
1.62
2.08
8
3
0.26
0.53
0.82
1.33
1.64
2.08
8
3
0.24
0.53
0.82
1.28
1.62
2.06
9
1
0.15
0.44
0.75
1.22
1.53
1.99
9
1
0.15
0.39
0.75
1.22
1.53
1.99
9
1
0.13
0.39
0.73
1.22
1.53
1.99
9
2
0.15
0.39
0.73
1.22
1.53
2.02
9
2
0.17
0.42
0.73
1.22
1.53
2.02
9
2
0.15
0.42
0.75
1.22
1.53
1.99
9
3
0.15
0.39
0.73
1.22
1.53
2.02
9
3
0.15
0.39
0.75
1.22
1.53
1.99
9
3
0.15
0.39
0.73
1.19
1.51
1.97
10
1
0.26
0.55
0.79
1.19
1.42
1.82
10
1
0.26
0.53
0.79
1.19
1.42
1.82
10
1
0.26
0.53
0.77
1.17
1.39
1.75
10
2
0.28
0.55
0.79
1.19
1.42
1.82
10
2
0.26
0.53
0.77
1.15
1.37
1.75
10
2
0.28
0.57
0.79
1.19
1.44
1.82
10
3
0.26
0.53
0.79
1.15
1.39
1.77
10
3
0.26
0.53
0.79
1.15
1.39
1.75
10
3
0.26
0.53
0.77
1.13
1.37
1.77
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors wish to thank Dr. Jinyu Yang for introducing the elementary idea of the paper. This research is supported by the National Natural Science Foundation of China (71102140) and the National Research Foundation for the Doctoral Program of Higher Education of China (20100032120030).
BurdickR. K.BorrorC. M.MontgomeryD. C.A review of methods for measurement systems capability analysis20033543423542-s2.0-1442324285MontgomeryD. C.2007John Wiley and SonsLyuJ.ChenM.-N.Gauge capability studies for attribute data200824171822-s2.0-4054913129910.1002/qre.868BorrorC. M.MontgomeryD. C.RungerG. C.Confidence intervals for variance components from gauge capability studies19971363613692-s2.0-0031270476GongL.BurdickR. K.QuirozJ.Confidence intervals for unbalanced two-factor gauge R&R studies20052187277412-s2.0-2884444380010.1002/qre.682SweeneyS.Analysis of two-dimensional gage repeatability and reproducibility200719129372-s2.0-3384579484110.1080/08982110601057641MajeskeK. D.Approval criteria for multivariate measurement systems20084021401532-s2.0-41549145252LiM.-H. C.Al-RefaieA.Improving wooden parts' quality by adopting DMAIC procedure20082433513602-s2.0-4214917916810.1002/qre.905ErdmannT. P.DoesR. J. M. M.BisgaardS.Quality quandaries: a gage R&R study in a hospital201022146532-s2.0-7795121395310.1080/08982110903412924Automotive Industry Action Group2010Automotive Industry Action GroupWheelerD. J.LydayR. W.1984Statistical Process ControlsWoodallW. H.BorrorC. M.Some relationships between gage R&R criteria2008241991062-s2.0-4054911950910.1002/qre.870JohnsonR. A.WichernD. W.2002Upper Saddle River, NJ, USAPrentice hallKimK.MahmoudM. A.WoodallW. H.On the monitoring of linear profiles20033533173282-s2.0-1642519894