Journal Search Engine

ISSN : 1598-7248 (Print)
ISSN : 2234-6473 (Online)

Industrial Engineering & Management Systems Vol.17 No.2 pp.184-192
DOI : https://doi.org/10.7232/iems.2018.17.2.184

Interval Estimation for Mean Response in a Simple Nested Error Regression Model in Gauge R & R Study

Dong Joon Park*

Department of Statistics Pukyong National University, Busan, Republic of Korea

Corresponding Author, djpark@pknu.ac.kr

Received: January 31, 2017 Review June 21, 2017 Accepted: July 11, 2017

ABSTRACT

In a measurement system monitoring manufacturing processes, one often uses gauge repeatability and reproducibility (R & R) experiment to determine the amount of variability due to parts, final products, operators, or gauge. A measurement system that employs operators randomly chosen to conduct measurements on randomly selected parts can be statistically modelled. When a measurement is assumed to be linearly related to a predictor variable and a simple nested error regression model is applied for the measurement system, one might be interested in making inferences concerning the variability for the mean response in the model. In general, confidence intervals are uniformly more informative than hypothesis tests in making inferences to parametric values in statistical models. Confidence intervals for the mean response in a simple nested error regression model can be useful tools to determine whether the variability is appropriately managed in a manufacturing process. Several confidence intervals for the mean response in the model are proposed. The confidence intervals proposed in this article are based on a moderate large sample method. Computer simulation is performed to see if the proposed intervals maintain the stated confidence level and a numerical example is provided to explain the proposed intervals by calculating confidence intervals for the mean response in the model.

Keywords : Nested Error Regression Model , Confidence Intervals , Mean Response , Repeatability , Reproducibility

초록

키워드 :

This article has been cited by 0 article in crossref

Cited-By

Funding:

1. INTRODUCTION

A gauge is employed to collect replicated measurements on units by several different operators, setups, or time periods in measurement system. The measurements in manufacturing processes include variability due to parts, operators, and gauges. The data collected during measuring processes are used to monitor to decrease the variability generated in production. When a measurement is linearly related to a predictor variable and a simple nested error regression model is applied for the measurement system, one might be interested in making inferences concerning the variability for the mean response in the model. The simple nested regression model considered in this article is a simple regression model with nested error structure where the response variable is linearly related to a predictor variable in manufacturing processes for gauge R & R study. Burdick et al. (2005) define that repeatability represents the gauge variability when it is used to measure the same unit (with the same operator or the same setup or in the same time period). Reproducibility refers to the variability arising from different operators, setups, or time periods. Measurement system capability studies are often referred to as gauge repeatability and reproducibility (R & R) study. The sources of variability in gauge R & R study are expressed as variance components so that confidence intervals for the variance components are employed to determine the adequacy of manufacturing for gauge R & R study.

This article focuses on confidence intervals for the mean response in a simple nested regression model which models a manufacturing process. This simple nested regression model includes two error terms. One is associated with the first stage sampling unit which is an operator effect and the other with second stage sampling unit which is a measurement error. The two error terms are assumed to be independent and normally distributed with zero means and constant variances. This article derives several confidence intervals for the mean response in the simple nested regression model.

2. LITERATURE REVIEW

A simple nested error regression model applying for a manufacturing process in this article can be viewed as a statistical model for clustered observations where subjects are nested within groups in the field of educational or biological sciences. Nested data structures are easily found in two types of structures: repeated measures within individuals or students within schools. The data structures are commonly appeared when multi-stage sampling procedures are taken from a target population. The data with this structure are written as nested error structure models, hierarchical linear models, or multilevel models. Observations within clusters tend to be more homogeneous than those between clusters. Misuse of standard ordinary least squares regression when clustered data exist violates the assumptions of data structures.

A great deal of research on statistical models with hierarchically structured data has been conducting by researchers (Goldstein, 2011; Raudenbush, 1995; Longford, 1995). They started with the simplest multilevel model with each group’s intercept, $Y_{i j} = α_{j} + β_{j} X_{i j} + E_{i j}$ , which is at level 1. They then introduce level 2 random variables into the simplest multilevel model by replacing α_j by $β_{0 j} = β_{0} + O_{0 j}$ and β_j by $β_{1 j} = β_{1} + O_{1 j}$ . The level 1 variables are observations within clusters and level 2 variables observations between clusters. Goldstein (2011) referred to this model as a 2-level model with hierarchically structured data. This model becomes the simple nested error regression model in this article by letting β₀, O_0j and O_1j be α , O_i and 0 respectively. We assumed that a predictor variable X_ij and the regression coefficients α and β are fixed factors and O_i and E_ij are random factors where $O_{i} ~ N (0, σ_{O}^{2})$ and $E_{i j} ~ N (0, σ_{E}^{2})$ in the model.

Parameters, α and , β and variance components, $σ_{O}^{2}$ and $σ_{E}^{2}$ in multilevel models have been estimated by applying for variance component models, random slopes models, ordinary regression models, etc. (Huang, 2016; McNeish and Stapleton, 2016; Mass and Hox, 2004; Aitkin and Longford, 1986). They have also compared the results of estimation methods such as maximum likelihood estimation, restricted maximum likelihood estimation, ordinary least squares estimation, or robust estimation, depending on various situations and assumptions of multilevel models. Although maximum likelihood estimation and ordinary least squares estimation are not generalized to all situations, using ordinary least estimation is worth considering due to its parsimonious nature (Huang, 2017). The ordinary least squares regression that ignores the nested error structure can lead to substantially biased estimates of regression estimators and their standard errors (Goldstein et al., 1993; Longford, 1986). Their research on multilevel models has mostly concentrated on Point estimation of the parameters, bias of standard errors, and statistical significance.

Park and Burdick (1993, 1994) provided confidence intervals for regression coefficients and variance components in the regression model with one-fold nested regression model. Park (2013) and Park and Yoon (2016) conducted further study to propose confidence intervals for regression coefficients and variance components in a regression model with two-fold nested error structure using various estimation methods. To our knowledge, the inference concerning expected mean response in the nested error regression model for a gauge R & R study has not actively researched up until now. We make, therefore, an attempt to provide confidence intervals for mean response in the model by extending the researches of Park and Burdick (1993, 1994) and Park and Hwang (2002). We specifically aim to present alternative confidence intervals for mean response by relaxing the condition for the mean value of predictor variable, i.e., when an individual value of predictor variable X_ij is given rather than calculating the mean value of predictor variable X_i.

3. MODEL STATEMENT

Consider that a manufacturing process where parts are measured by operators and the response variable as a measurement is linearly related to a predictor variable. The simple nested error regression model is then written as

Y_{i j} = α + β X_{i j} + O_{i} + E_{i j}; i = 1, \dots, I; j = 1, \dots, J

(1)

where Y_ij is the j th measurement of a part measured by the ith operator as a response variable, α and β are regression parameters, X_ij is a fixed predictor variable related to the response variable Y_ij ,O_i is the ith randomly chosen operator effect, and E_ij is the jth measurement error of a part measured by the ith operator. Operator effect O_i is associated with the first stage sampling unit and measurement error effect E_ij with second stage sampling unit. Two error terms O_i and E_ij are jointly independent normal random variables with zero means and variances $σ_{O}^{2}$ and $σ_{P}^{2}$ , respectively. That is, $O_{i} ~ N (0, σ_{O}^{2})$ and $E_{i j} ~ N (0, σ_{E}^{2})$ . Since α ,β , and X_ij are fixed factors and O_i and E_ij are random factors, equation (1) is a mixed model.

The model (1) is written in matrix notation(2)

y = X β + Z o + e

(2)

where y is an IJ ×1 vector of measurements, X is an IJ × 2 matrix of predictor variables with a column of 1’s in the first column and a column of X_ij′s in the second column, β is a 2×1 vector of parameters with α and β as elements, Z is an IJ × I design matrix with 0’s and 1’s, i.e. $Z = \oplus_{i = 1}^{I} 1_{J}, 1_{J}$ is a vector of 1’s, o is an I × 1 vector of operator effects, e is an IJ × 1 vector of measurement errors. By the assumptions in model (1) the response variables have a multivariate normal distribution as follows:

y ~ N (X β, V)

(3)

where $V = σ_{O}^{2} Z Z^{'} + σ_{E}^{2} D_{I J}$ and D_IJ is an IJ × IJ identity matrix.

A possible partitioning for source of variability of model (1) is written in Table 1. The sums of squares in Table 1 is defined as follows: $R_{A} = S_{y y a} - {\hat{β}}_{A}^{2} S_{x x a}$ and $R_{L} = {\hat{β}}_{A}^{2} S_{x x a} + {\hat{β}}_{W}^{2} S_{x x w} - R_{1}$ where $S_{x x a} = J Σ_{i} {({\bar{X}}_{i .} - {\bar{X}}_{..})}^{2},$ $S_{y y a} = J Σ_{i} {({\bar{Y}}_{i .} - {\bar{Y}}_{..})}^{2},$ $S_{x y a} = J Σ_{i} ({\bar{X}}_{i .} - {\bar{X}}_{..}) ({\bar{Y}}_{i .} - {\bar{Y}}_{..}),$ $S_{x x w} {=Σ}_{i} Σ_{j} {(X_{i j} - {\bar{X}}_{i .})}^{2},$ $S_{y y w} = Σ_{i} Σ_{j} {(Y_{i j} - {\bar{Y}}_{i .})}^{2},$ $S_{x y w} = Σ_{i} Σ_{j} (X_{i j} - {\bar{X}}_{i .}) (Y_{i j} - {\bar{Y}}_{i .}),$ ${\bar{X}}_{i .} = Σ_{j} X_{i j} / J,$ ${\bar{Y}}_{i .} = Σ_{j} Y_{i j} / J,$ ${\bar{X}}_{..} = Σ_{i} Σ_{j} X_{i j} / I J,$ and ${\bar{Y}}_{..} = Σ_{i} Σ_{j} Y_{i j} / I J .$ Three estimators of regression coefficient β used in Table 1 are as follows: ${\hat{β}}_{c} = (S_{x y a} + S_{x y w}) / (S_{x x a} + S_{x x w}),$ ${\hat{β}}_{A} = S_{x y a} / S_{x x a},$ and ${\hat{β}}_{W} = S_{x y w} / S_{x x w}$ .

4. DISTRIBUTIONAL RESULTS OF ESTIMATORS OF MEAN RESPONSE AND SUMS OF SQUARES

In order to construct confidence intervals for the mean response E(Y_ij ) , the estimators of mean response and their distributional properties are derived. Three estimators of the mean response are presented. The independence between estimators of mean response and sums of squares in Table 1 is examined.

4.1. Within Group Estimator of Mean Response

The mean response of the simple nested error regression model (1) is defined as $E (Y_{i j}) = α + β X_{i j}$ where E(Y_ij ) means the expected value of the jth measurement of a part measured by the ith operator. Park and Burdick (1993) showed that within group ordinary least squares estimators (OLSE) ${\hat{α}}_{W}$ and ${\hat{β}}_{W}$ of parameters in model (1) are obtained from regression of Y_ij on the X_ij and the grouping variables. In matrix notation, the estimators are the first two elements of the vector ${(X^{*^{'}} X^{*})}^{G} X^{*^{'}} y$ where $X^{*} = [X Z]$ and ${(X^{*^{'}} X^{*})}^{G}$ is a generalized inverse of $X^{*^{'}} X^{*} .$ . The within group slope OLSE is shown as ${\hat{β}}_{W} = S_{x y w} / S_{x x w}$ and it is normally distributed, i.e. ${\hat{β}}_{W} ~ N (β, σ_{E}^{2} / S_{x x w}) .$ It can be shown by elementary algebra that the within group intercept OLSE is shown as ${\tilde{α}}_{W} = J {\hat{α}}_{W} / I = {\bar{Y}}_{..} - {\hat{β}}_{W} {\bar{X}}_{..}$ Since an unbiased estimator of E(Y_ij ) is ${\hat{Y}}_{i j}$ , the within group OLSE of the mean response E(Y_ij) is obtained by substituting within group OLSEs ${\tilde{α}}_{W}$ and ${\hat{β}}_{W}$ for parameters in model (1), i.e. ${\tilde{Y}}_{i j W} = {\tilde{α}}_{W} + {\hat{β}}_{W} X_{i j}$ . It can be shown by the assumptions of model (1) that within group OLSE of E(Y_ij ) is normally distributed as ${\tilde{Y}}_{i j w} ~ N (E (Y_{i j}), (J σ_{O}^{2} + σ_{E}^{2}) [1 / I J + {(X_{i j} - {\bar{X}}_{..})}^{2} ϕ / S_{x x w}])$ where $ϕ = σ_{E}^{2} / (J σ_{O}^{2} + σ_{E}^{2})$ .

4.2. Total Estimator of Mean Response

Total OLSEs ${\hat{α}}_{C}$ and ˆβC of parameters in model (1) are obtained from regression of Y_ij on X_ij . Park and Burdick (1993) showed that the total OLSEs are obtained by the vector ${(X^{'} X)}^{- 1} X^{'} y$ . The total slope OLSE is written as ${\hat{β}}_{C} = (S_{x y a} + S_{x y w}) / (S_{x x a} + S_{x x w})$ and it is normally distributed, i.e. ${\hat{β}}_{C} = (S_{x y a} + S_{x y w}) / (S_{x x a} + S_{x x w})$ where $r^{2} = S_{x x a} / (S_{x x a} + S_{x x w})$ . The total intercept OLSE is shown to be ${\hat{α}}_{C} = {\bar{Y}}_{..} - {\hat{β}}_{C} {\bar{X}}_{..} .$ The total OLSE of the mean response E(Y_ij) is therefore obtained by substituting total OLSEs ${\hat{α}}_{C}$ and ${\hat{β}}_{C}$ for parameters in model (1), i.e. ${\hat{Y}}_{i j C} = {\hat{α}}_{C} + {\hat{β}}_{C} X_{i j}$ . It can be shown by the assumptions of model (1) that total OLSE of E(Y_ij ) is normally distributed as ${\hat{Y}}_{i j C} ~ N (E (Y_{i j}), (J σ_{O}^{2} + σ_{E}^{2}) [1 / I J + (r^{2} J σ_{O}^{2} + σ_{E}^{2}) {(X_{i j} - {\bar{X}}_{..})}^{2} / {(J σ_{O}^{2} + σ_{E}^{2}) (S_{x x a} + S_{x x w})}]) .$

4.3. Best Linear Estimator of Mean Response

Park and Burdick (1994) derived among groups OLSEs ${\hat{α}}_{A}$ and ${\hat{β}}_{A}$ of parameters in model (1) by regression of ${\bar{Y}}_{i .}$ on ${\bar{X}}_{i .}$ . and the among groups OLSEs are obtained by the vector ${({M^{'}}_{X} M_{X})}^{- 1} {M^{'}}_{X} y$ where $M_{X} = M X$ and $M = (1 / J) Z .$ They showed that the best linear unbiased estimator (BLUE) ${\hat{β}}_{B}$ of slope parameter in model (1) is a linear combination of ${\hat{β}}_{A}$ and ${\hat{β}}_{W}$ . It can be shown by matrix algebra that the BLUEs ${\hat{α}}_{B}$ and ${\hat{β}}_{B}$ of two parameters in model (1) are calculated by the vector ${(X^{'} V^{- 1} X)}^{- 1} X^{'} V^{- 1} y$ where $V = V a r (y)$ in equation (3). The slope BLUE is written as ${\hat{β}}_{B} = (\hat{ϕ} S_{x y a} + S_{x y w}) / (\hat{ϕ} S_{x x a} + S_{x x w})$ where $\hat{ϕ} = {\hat{σ}}_{E}^{2} / (J {\hat{σ}}_{O}^{2} + {\hat{σ}}_{E}^{2})$ and it is normally distributed, i.e. ${\hat{β}}_{B} ~ N (β, σ_{E}^{2} / (ϕ S_{x x a} + S_{x x w})) .$ . The intercept BLUE is shown as ${\hat{α}}_{B} = {\bar{Y}}_{..} - {\hat{β}}_{B} {\bar{X}}_{..} .$ . The BLUE of the mean response E(Y_ij ) is therefore obtained by substituting the BLUEs ${\hat{α}}_{B}$ and ${\hat{β}}_{B}$ for two parameters in model (1), i.e. ${\hat{Y}}_{i j B} = {\hat{α}}_{B} + {\hat{β}}_{B} X_{i j} .$ It can be shown by the assumptions of model (1) that BLUE of E(Y_ij ) is normally distributed as ${\hat{Y}}_{i j B} ~ N (E (Y_{i j}), (J σ_{O}^{2} + σ_{E}^{2}) [1 / I J + {(X_{i j} - {\bar{X}}_{..})}^{2} ϕ / (ϕ S_{x x a} + S_{x x w})]) .$

4.4. Independence of Point Estimators and Sums of Squares

Park and Burdick (1994) showed that the sums of squares R_A and R_W are jointly independent chi-squared random variables with I − 2 and IJ − I −1 degrees of freedom, i.e. $R_{A} / (J σ_{O}^{2} + σ_{E}^{2}) ~ χ_{I - 2}^{2}$ and $R_{W} / σ_{E}^{2} ~ χ_{I J - I - 1}^{2}$ . Thus the sum of squares $R_{B} / (J σ_{O}^{2} + σ_{E}^{2})$ is a chisquared random variable with IJ − 3 degrees of freedom, $R_{B} / (J σ_{O}^{2} + σ_{E}^{2}) ~ χ_{I J - 3}^{2}$ where $R_{B} = R_{A} + R_{W} / ϕ$ since $R_{A} \sim (J σ_{O}^{2} + σ_{E}^{2}) χ_{I - 2}^{2}$ and $R_{W} / ϕ ~ (J σ_{O}^{2} + σ_{E}^{2}) χ_{I J - I - 1}^{2}$ and R_B is a linear combination of two independent chi-squared random variables. In order to construct confidence intervals for mean response E(Y_ij ) we need to show the independence between estimators, ${\tilde{Y}}_{i j W}, {\hat{Y}}_{i j C},$ and ${\hat{Y}}_{i j B}$ and the sums of squares, R_A, R_W, and R_B.

Theorem 1.The within group estimator ${\tilde{Y}}_{i j W}$ and the sum of squares R_A are independent.

Proof: In matrix notation the within group OLSE of mean response is written ${\tilde{Y}}_{i j w} = x_{1} {(X^{*^{'}} X^{*})}^{G} X^{*^{'}} y$ where $x_{1} = [J / I, X_{i j}, 0, 0, \dots, 0]$ , i.e., x₁ is a (2 + I)×1 vector of J / I, X_ij , and the number of 0’s with I as elements. The sum of squares R_A in Table 1 is defined as $B = [D_{I} - M_{X} {({M^{'}}_{X} M_{X})}^{- 1} {M^{'}}_{X}]$ in matrix notation. The matrix definitions of Y_ijw and R_A are now utilized to show the independence between them. It can be shown by matrix manipulation that $x_{1} {(X^{*^{'}} X^{*})}^{G} X^{*^{'}} \times (σ_{O}^{2} Z Z^{'} + σ_{E}^{2} D_{I J}) \times J A = 0$ since $Z Z^{'} A = J A$ and $x_{1} {(X^{*^{'}} X^{*})}^{G} X^{*^{'}} A = 0 .$ Thus Y_ijW and R_A are independent by Theorem 7.5 of Searle (1987).

Theorem 2. The within group estimator ${\tilde{Y}}_{i j w}$ and the sum of squares R_W are independent.

Proof: The sum of squares R_W is defined R_W = y′Wy where $W = D_{I J} - X^{*} {(X^{*^{'}} X^{*})}^{G} X^{*^{'}} .$ In order to show independence between ${\tilde{Y}}_{i j w}$ and R_W matrix definitions are used. It can be shown by matrix manipulation that $x_{1} {(X^{*^{'}} X^{*})}^{G} X^{*^{'}} \times (σ_{O}^{2} Z Z^{'} + σ_{E}^{2} D_{I J}) \times W = 0$ since $Z Z^{'} W = 0$ and ${(X^{*^{'}} X^{*})}^{G} X^{*^{'}} W = 0 .$ Thus ${\tilde{Y}}_{i j w}$ and R_W are independent by Theorem 7.5 of Searle (1987).

Theorem 3. The within group estimator ${\tilde{Y}}_{i j w}$ and the sum of squares R_Bare independent.

Proof: The sum of squares R_B was defined as $R_{B} = J y^{'} A y + y^{'} W y = y^{'} [J A + 1 / ϕ W] y$ . It can be written that R_B = J y′A y + y′Wy = y′[J A+1/ϕ W] y in matrix notation. The within group estimator ${\tilde{Y}}_{i j w}$ is therefore independent of the sum of squares R_B by Theorems 1 and 2.

Theorem 4.The total estimator ${\tilde{Y}}_{i j C}$ and the sum of squares R_A are independent.

Proof: Total OLSE of mean response is written ${\hat{Y}}_{i j C} = x_{2} {(X^{'} X)}^{- 1} X^{'} y$ where $x_{2} = [1 X_{i j}]$ in matrix notation. It can be shown by using the matrix definitions of ${\hat{Y}}_{i j C}$ and R_A that $x_{2} {(X^{'} X)}^{- 1} X^{'} \times (σ_{O}^{2} Z Z^{'} + σ_{E}^{2} D_{I J}) \times J A = 0$ since $R_{A} = J y^{'} A y, Z Z^{'} A = J A, and X^{'} A = 0 .$ . Thus ${\hat{Y}}_{i j C}$ and R_A are independent by Theorem 7.5 of Searle (1987).

LTheorem 5.The total estimator ${\hat{Y}}_{i j C}$ and the sum of squaresR_Ware independent.

Proof: It can be shown by using the matrix definitions of ${\hat{Y}}_{i j C}$ and R_W that $x_{2} {(X^{'} X)}^{- 1} X^{'} \times (σ_{O}^{2} Z Z^{'} + σ_{E}^{2} D_{I J}) \times W = 0$ using $Z Z^{'} W = 0$ and $X^{'} W = 0 .$ . Thus ${\hat{Y}}_{i j C}$ and R_W are independent by Theorem 7.5 of Searle (1987).

Theorem 6.The total estimator ${\hat{Y}}_{i j C}$ and the sum of squares R_Bare independent.

Proof: It follows that the total estimator ${\hat{Y}}_{i j C}$ and the sum of squares R_B are independent by Theorems 4 and 5.

Theorem 7. The best linear estimator ${\hat{Y}}_{i j B}$ and the sum of squares R_A are independent.

Proof: The matrix definitions of ${\hat{Y}}_{i j B}$ and R_A are utilized to show the independence between them. In matrix notation the BLUE of E(Y_ij ) is written as ${\hat{Y}}_{i j B} = x_{2} {(X^{'} V^{- 1} X)}^{- 1} X^{'} V^{- 1} y .$ It can be shown that $x_{2} {(X^{'} V^{- 1} X)}^{- 1} X^{'} V^{- 1} \times V \times J A = 0$ since $V = (σ_{O}^{2} Z Z^{'} + σ_{E}^{2} D_{I J})$ and $X^{'} A = 0$ . Thus ${\hat{Y}}_{i j B}$ and R_A are independent by Theorem 7.5 of Searle (1987).

Theorem 8.The best linear estimator ${\hat{Y}}_{i j B}$ and the sum of squares R_Ware independent.

Proof : It can be shown by using the matrix definitions of ${\hat{Y}}_{i j B}$ and R_W that $x_{2} {(X^{'} V^{- 1} X)}^{- 1} X^{'} V^{- 1} \times V \times W = 0$ since $V = (σ_{O}^{2} Z Z^{'} + σ_{E}^{2} D_{I J})$ and $X^{'} A = 0$ . Thus ${\hat{Y}}_{i j B}$ and R_W are independent by Theorem 7.5 of Searle (1987).

Theorem 9.The best linear estimator ${\hat{Y}}_{i j B}$ and the sum of squares R_B are independent.

Proof: It follows that the best linear estimator ${\hat{Y}}_{i j B}$ and the sum of squares R_B are independent by Theorems 7 and 8.

We showed the independence between three estimators of mean response and three sums of squares. In summary theorems say that ${\tilde{Y}}_{i j w}$ , R_A , R_W, and R_B are jointly independent, YˆijC , R_A, R_W, and R_B are jointly independent, and YˆijB, R_A, R_W, and R_B are jointly independent.

5. CONFIDENCE INTERVALS FOR MEAN RESPONSE

The confidence intervals for mean response are constructed using distributional properties of OLSEs and BLUE of E(Y_ij ) and independence between the estimators and the sums of squares described in Theorems in Section 4. Since ${\tilde{Y}}_{i j w}$ is normally distributed and ${\tilde{Y}}_{i j w}$ and R_A are independent by Theorem 1, it follows that an exact 100(1−α)% confidence interval for E(Y_ij ) is

{\tilde{Y}}_{i j w} \pm t_{(I - 2; α / 2)} \sqrt{S_{A}^{2} [\frac{1}{I J} + \frac{{(X_{i j} - {\bar{X}}_{..})}^{2} ϕ}{S_{x x w}}]}

(4)

where $S_{A}^{2}$ is the error mean square among groups and defined as $S_{A}^{2} = R_{A} / (I - 2)$ and $t_{(υ; δ)}$ is the t-value for ν degrees of freedom with δ area to the right. From two independent chi-squared random variables $R_{A} / (J σ_{O}^{2} + σ_{E}^{2}) ~ χ_{I - 2}^{2}$ and $R_{W} / σ_{E}^{2} ~ χ_{I J - I - 1}^{2}$ in Section 4.4, the error mean squares and their expected mean squares are obtained as follows:

E (S_{A}^{2}) = J σ_{O}^{2} + σ_{E}^{2} and

(5)

E (S_{W}^{2}) = σ_{E}^{2}

(6)

where $S_{W}^{2}$ is the error mean square within group and defined as $S_{W}^{2} = R_{W} / (I J - I - 1) .$ The unbiased estimators of the variances are respectively

{\hat{σ}}_{O}^{2} = (S_{A}^{2} - S_{W}^{2}) / J

and

{\hat{σ}}_{E}^{2} = S_{W}^{2}

from (5) and (6) and an estimator

\hat{ϕ} = S_{W}^{2} / S_{A}^{2}

is thus used to construct a confidence interval for E(Y_ij). It follows by substituting the unbiased estimator

\hat{ϕ}

for ϕ in interval (4) that an exact 100(1−α )% confidence interval for E(Y_ij) is

{\tilde{Y}}_{i j w} \pm t_{(I - 2; α / 2)} \sqrt{S_{A}^{2} [\frac{1}{I J} + \frac{{(X_{i j} - {\bar{X}}_{..})}^{2} S_{W}^{2}}{S_{x x w} S_{A}^{2}}]} .

(7)

In a similar manner, using the distributional property of ${\hat{Y}}_{i j C}$ and independence of ${\hat{Y}}_{i j C}$ and R_A in Theorem 4, it follows that an exact 100(1−α )% confidence interval for E(Y_ij ) is

{\hat{Y}}_{i j C} \pm t_{(I - 2; α / 2)} \sqrt{\frac{S_{A}^{2}}{I J} + \frac{{r^{2} S_{A}^{2} + (1 - r^{2}) S_{W}^{2}} {(X_{i j} - {\bar{X}}_{..})}^{2}}{S_{x x a} + S_{x x w}}} .

(8)

Using the distributional property of ${\hat{Y}}_{i j B}$ and independence of ${\hat{Y}}_{i j B}$ and R_A in Theorem 7, it follows that an exact 100(1−α )% confidence interval for E(Y_ij ) is

{\hat{Y}}_{i j B} \pm t_{(I - 2; α / 2)} \sqrt{\frac{S_{A}^{2}}{I J} + \frac{S_{A}^{2} S_{W}^{2} {(X_{i j} - {\bar{X}}_{..})}^{2}}{S_{W}^{2} S_{x x a} + S_{A}^{2} S_{x x w}}} .

(9)

On the other hand, one can construct confidence intervals for mean response by using independence between OLSEs and BLUE of E(Y_ij ) and the sum of squares within group R_W described in Theorems in Section 4. Since ${\tilde{Y}}_{i j w}$ is normally distributed and ${\tilde{Y}}_{i j w}$ and R_W are independent by Theorem 2, it follows that an exact 100(1−α )% confidence interval for E(Y_ij ) is

{\tilde{Y}}_{i j w} \pm t_{(I J - I - 1; α / 2)} \sqrt{\frac{S_{A}^{2}}{I J} + \frac{{(X_{i j} - {\bar{X}}_{..})}^{2} S_{W}^{2}}{S_{x x w}}} .

(10)

This is referred to as WG1 method.

Similarly, using the distributional property of ${\hat{Y}}_{i j C}$ and independence of ${\hat{Y}}_{i j C}$ and R_W in Theorem 5, it follows that an exact 100(1−α )% confidence interval for E(Y_ij ) is

{\hat{Y}}_{i j C} \pm t_{(I J - I - 1; α / 2)} \sqrt{\frac{S_{A}^{2}}{I J} + \frac{{r^{2} S_{A}^{2} + (1 - r^{2}) S_{W}^{2}} {(X_{i j} - {\bar{X}}_{..})}^{2}}{S_{W}^{2} (S_{x x a} + S_{x x w})}} .

(11)

This is referred to as T1 method.

Using the distributional property of ${\hat{Y}}_{i j B}$ and independence of ${\hat{Y}}_{i j B}$ and R_W in Theorem 8, it follows that an exact 100(1−α )% confidence interval for E(Y_ij ) is

{\hat{Y}}_{i j B} \pm t_{(I J - I - 1; α / 2)} \sqrt{\frac{S_{A}^{2}}{I J} + \frac{S_{A}^{2} S_{W}^{2} {(X_{i j} - {\bar{X}}_{..})}^{2}}{S_{W}^{2} S_{x x a} + S_{A}^{2} S_{x x w}}}

(12)

This is referred to as BL1 method.

One can also construct confidence intervals for mean response by using independence between OLSEs and BLUE of E(Y_ij ) and the sum of squares R_B described in Theorems in Section 4. Using the normality of ${\tilde{Y}}_{i j w}$ and independence of ${\tilde{Y}}_{i j w}$ and R_B in Theorem 3, it follows that an exact 100(1−α )% confidence interval for E(Y_ij ) is

{\tilde{Y}}_{i j w} \pm t_{(I J - 3; α / 2)} \sqrt{S_{B}^{2} [\frac{1}{I J} + \frac{{(X_{i j} - {\bar{X}}_{..})}^{2} S_{W}^{2}}{S_{x x w} S_{A}^{2}}]}

(13)

where $S_{B}^{2}$ is the error mean square and defined as $S_{B}^{2} = R_{B} / (I J - 3) .$ . This is referred to as WG2 method.

Using the normality of ${\hat{Y}}_{i j C}$ and independence of ${\hat{Y}}_{i j C}$ and R_B in Theorem 6, it follows that an exact 100(1−α )% confidence interval for E(Y_ij ) is

{\hat{Y}}_{i j C} \pm t_{(I J - 3; α / 2)} \sqrt{S_{B}^{2} [\frac{1}{I J} + \frac{{r^{2} + (1 - r^{2}) (S_{W}^{2} / S_{A}^{2})} {(X_{i j} - {\bar{X}}_{..})}^{2}}{(S_{x x a} + S_{x x w})}]} .

(14)

This is referred to as T2 method.

Similarly, using the normality of ${\hat{Y}}_{i j B}$ and independence of ${\hat{Y}}_{i j B}$ and R_B in Theorem 9, it follows that an exact 100(1−α )% confidence interval for E(Y_ij ) is

{\hat{Y}}_{i j B} \pm t_{(I J - 3; α / 2)} \sqrt{S_{B}^{2} [\frac{1}{I J} + \frac{S_{W}^{2} {(X_{i j} - {\bar{X}}_{..})}^{2}}{S_{W}^{2} S_{x x a} + S_{A}^{2} S_{x x w}}]} .

(15)

This is referred to as BL2 method.

We constructed nine exact confidence intervals for E(Y_ij ) in this Section. Confidence intervals (7), (8), and (9) use t-values with I − 2 degrees of freedom only whereas confidence intervals (10) to (15) use more degrees of freedom. Thus, confidence intervals (10) to (15) produce narrower interval lengths which are preferred in interval estimation than (7), (8), and (9) do.

6. SIMULATION STUDY

The performance of the confidence intervals proposed in Section 5 is examined by computer simulation. Twenty five designs are formed by taking all combinations of I = 3, 5, 10, 15, 20 and J = 3, 5, 10, 15, 20. The values of $σ_{O}^{2}$ are selected from the set of values (0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9) and the values of $σ_{E}^{2}$ are determined from $σ_{O}^{2} + σ_{E}^{2} = 1$ without loss of generality. Recall that the mean squares in Section 4 are chisquared random variables. In particular, $S_{A}^{2} ~ [(J σ_{O}^{2} + σ_{E}^{2}) / (I - 2)] χ_{(I - 2)}^{2},$ $S_{W}^{2} ~ [σ_{E}^{2} / (I J - I - 1)] χ_{(I j - I - 1)}^{2},$ and $S_{B}^{2} ~ [(J σ_{O}^{2} + σ_{E}^{2}) / (I J - 3)] χ_{(I J - 3)}^{2} .$ These mean squares are generated by the R_ANGAM function in SAS (Statistical Analysis System) by substituting the specific values of $σ_{O}^{2}$ and $σ_{E}^{2}$ that are selected respectively.

The values of S_xxa are chosen from the set of values (0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9) and the values of S_xxw are determined from $S_{x x a} + S_{x x w} = 1$ without loss of generality. The OLSEs and BLUE of E(Y_ij ) are generated by the distributional properties in Section 4. In particular, ${\tilde{Y}}_{i j w} ~ N (E (Y_{i j}), (J σ_{O}^{2} + σ_{E}^{2}) [1 / I J + {(X_{i j} - {\bar{X}}_{..})}^{2} ϕ / S_{x x w}]),$ and ${\hat{Y}}_{i j C} ~ N (E (Y_{i j}), (J σ_{O}^{2} + σ_{E}^{2}) [1 / I J + (r^{2} J σ_{O}^{2} + σ_{E}^{2}) {(X_{i j} - {\bar{X}}_{..})}^{2} / {(J σ_{O}^{2} + σ_{E}^{2}) (S_{x x a} + S_{x x w})}]),$ ${\hat{Y}}_{i j B} ~ N (E (J σ_{O}^{2} + σ_{E}^{2}) [1 / I J + {(X_{i j} - {\bar{X}}_{..})}^{2} ϕ / (ϕ S_{x x a} + S_{x x w})]) .$ 2 2 2 ${\tilde{Y}}_{i j w}$ . These estimators are generated by using R_ANNOR function SAS by substituting the specific values of S_xxa and S_xxw. The Simulated values of $S_{A}^{2}, S_{W}^{2}, S_{B}^{2}, {\tilde{Y}}_{i j w}, {\hat{Y}}_{i j C}, and {\hat{Y}}_{i j B}$ are substituted into confidence intervals (10) to (15). For each of 25 designs with different combinations of I and J, 2000 iterations are simulated and two sided confidence intervals for E(Y_ij ) are computed. Confidence coefficients are determined by counting the number of the intervals that contain E(Y_ij ). The average lengths of the two sided confidence intervals for E(Y_ij ) are computed.

Using the normal approximation to the binomial, if the true confidence coefficient is 0.90, there is less than a 2.5% chance that a simulated confidence coefficient based on 2000 replications will be less than 0.8868. The comparison criteria are: (i) the ability to maintain the stated confidence coefficient and (ii) the average length of two sided confidence intervals. Although narrower average interval lengths are preferable, it is necessary that an interval first maintain the stated confidence level.

Table 2 reports the range of simulated confidence coefficients for stated 90% two sided confidence intervals on E(Y_ij) with I = 3, 10, 20 and J = 3, 10, 20 as $σ_{O}^{2}$ and S_xxa range from 0.1 to 0.9, respectively. The WG1, T1, BL1, WG2, T2, and BL2 methods refer to the intervals (10), (11), (12), (13), (14), and (15), respectively. The WG2, T2, and BL2 methods generally maintain the stated confidence level 0.9 across all combinations of I and J. However, the WG1, T1, and BL1 methods give confidence coefficients that are much below 0.8868 when I = 3. The WG1 method slightly maintains the stated confidence level when I becomes greater than or equal to 10 whereas T1 method is too conservative since the simulated confidence coefficients are very close to 1.0. Thus WG1, T1, and BL1 methods are not generally recommended except WG1 method with I ≥ 10. The simulated confidence coefficients that are less than 0.8868 or abnormally greater than the stated confidence level are shown in boldface in Table 2.

Table 3 reports the range of the average interval lengths for WG1, WG2, T2, and BL2 methods with I = 3, 10, 20 and J = 3, 10, 20. The T1 and BL1 methods are eliminated because they do not generally maintain the simulated confidence level across all combinations of I and J in Table 2. The WG1 method is also eliminated for the same reason when I = 3. The three methods generally yield shorter average interval lengths as I and J increase since the degrees of freedom become large and the standard errors of OLSEs and BLUE of E(Y_ij ) become small. Although the average interval lengths of three methods change largely depending on the values of $σ_{O}^{2}, σ_{E}^{2}, S_{x x a} and S_{x x a}$ , BL2 method generally produces the shortest interval length than T2 and WG2 methods from Table 3.

7. EXAMPLE APPLICATION

One of the manufacturing processes of integrated circuits in semiconductor technology is to connect individual transistors, capacitors, and other circuit elements with a conducting metallic material. These connections are typically built by first depositing a thin blanket layer of metal over the entire silicon wafer and then etching away unnecessary portion. Czitrom and Spagon (1997) presented a data set of a designed experiment in seven factors and two blocks in 38 runs with 14 responses to perform this connection process. Reflectivity (%) as a response variable Y_ij and resistivity (U-ohm cm) as a predictor variable X_ij are selected from the data set and they are shown in Table 4. Three operators (I = 3) are chosen and five measurements (J = 5) are repeatedly conducted assuming a simple nested error regression model (1).

The data set in Table 4 is used to calculate confidence intervals for E(Y_ij ) . The overall means of reflectivity Y.. and resistivity X.. are respectively computed 12.047% and 87.333 U-ohm cm. Assume that a resistivity value X_ij = 85 U-ohm is given to estimate the value of reflectivity Y_ij. In that case practitioners often consider the variabilities of E(Y_ij) to see if there is any change in the connection manufacturing process. We suggest confidence intervals for mean response 3. The results E(Y_ij ) with certain degree of confidence by allowing the bound of errors rather than an estimated value ˆY_ij when a resistivity value is given. The WG2, T2, and BL2 methods are suggested when resistivity value is 85 because they keep the stated confidence level when I = 3. The results are presented in Table 5. In order to choose an appropriate confidence intervals, one generally prefers the shortest confidence interval that keeps the stated confidence level. The BL2 method yields the shortest interval length among three methods. This result is consistent with simulation study because Table 3 presents the same pattern of average interval lengths when I = 3.

8. CONCLUSIONS

This article presents statistical properties the OLSEs and BLUE of mean response E(Y_ij ) and independence of estimators and sums of squares appeared in a simple nested error regression model. A numerical example was illustrated to compute confidence intervals for E(Y_ij ) when a value of predictor variable X_ij is given. We first recommend to apply WG2, T2, and BL2 methods to compute a confidence interval for E(Y_ij ) when I is less then 10. We then suggest to choose the shortest confidence interval.

This article extends the researches of Park and Burdick (1993, 1994) and Park and Hwang (2002). This article proposes several confidence intervals for the mean response in a simple nested error regression model when a measurement is linearly related to a predictor variable in manufacturing processes for gauge R & R study. We specifically present alternative confidence intervals for mean response when an individual value of predictor variable X_ij is given rather than calculating the mean value of predictor variable X_i. Future research includes inference concerning expected mean response in a regression model with two-fold nested error structure using various estimation methods.

ACKNOWLEDGMENT

The author would like to thank the anonymous referees for their valuable time and comments on an earlier version of this article. This work was supported by a Research Grant of Pukyong National University (2017 year).

Figure

Table

Table 1.

A partition for source of variability of model (1)

Source of Variation	Degrees of Freedom	Sums of Squares
Predictor Variable after Mean	1	R1=β^C2 (Sxxa+Sxxw)
Operator Effect	I-1	RA+RL=Syya+β^W2 Sxxw −β^C2 (Sxxa+Sxxw)
Measurement Error Effect	IJ-I-1	RW= Syyw−β^W2 Sxxw
Corrected Total	IJ-1	RT=Syya+Syyw

Table 2.

Range of simulated confidence coefficients for 90% two sided intervals for E(Y_ij)

				Range	WG1	T1	BL1	WG2	T2	BL2
3	3	0.1	0.9	Max	0.9245	0.9685	0.8110	0.9010	0.9350	0.9070
		0.1	0.9	Min	0.8965	0.7990	0.6995	0.8790	0.9105	0.8855
		0.5	0.5	Max	0.9055	0.9840	0.8530	0.9170	0.9370	0.9225
		0.5	0.5	Min	0.8850	0.8255	0.7085	0.8905	0.9185	0.8770
		0.9	0.1	Max	0.9045	0.9970	0.8355	0.9425	0.9320	0.9340
		0.9	0.1	Min	0.8405	0.9085	0.7745	0.9200	0.9070	0.9125
3	10	0.1	0.9	Max	0.9070	0.9200	0.8025	0.9095	0.9435	0.9150
		0.1	0.9	Min	0.8870	0.7265	0.6875	0.8965	0.9145	0.8925
		0.5	0.5	Max	0.9015	0.9350	0.8360	0.9295	0.9400	0.9270
		0.5	0.5	Min	0.8845	0.7565	0.6955	0.9015	0.9100	0.9070
		0.9	0.1	Max	0.8875	0.9735	0.8045	0.9430	0.9400	0.9430
		0.9	0.1	Min	0.8010	0.8895	0.7780	0.9145	0.9010	0.9075
3	20	0.1	0.9	Max	0.9095	0.9075	0.8105	0.9195	0.9400	0.9145
		0.1	0.9	Min	0.8925	0.7210	0.6790	0.9025	0.9245	0.8990
		0.5	0.5	Max	0.9050	0.9010	0.8515	0.9345	0.9390	0.9265
		0.5	0.5	Min	0.8890	0.7475	0.7085	0.9045	0.9025	0.9100
		0.9	0.1	Max	0.8890	0.9485	0.8170	0.9395	0.9360	0.9480
		0.9	0.1	Min	0.8005	0.8835	0.7845	0.9245	0.8930	0.9110
10	3	0.1	0.9	Max	0.9195	0.9330	0.9035	0.9055	0.9135	0.9030
		0.1	0.9	Min	0.8815	0.8930	0.8760	0.8810	0.8930	0.8835
		0.5	0.5	Max	0.9170	0.9870	0.9025	0.9045	0.9100	0.9040
		0.5	0.5	Min	0.8940	0.9595	0.8730	0.8750	0.8935	0.8869
		0.9	0.1	Max	0.9185	1.0000	0.9105	0.9105	0.9060	0.9160
		0.9	0.1	Min	0.8940	0.9985	0.8910	0.8915	0.8990	0.8870
10	10	0.1	0.9	Max	0.9095	0.9140	0.8990	0.9030	0.9130	0.9025
		0.1	0.9	Min	0.8900	0.8885	0.8655	0.8890	0.8975	0.8820
		0.5	0.5	Max	0.9075	0.9770	0.9095	0.9125	0.9095	0.9120
		0.5	0.5	Min	0.8915	0.9445	0.8770	0.8895	0.8925	0.8905
		0.9	0.1	Max	0.9085	1.0000	0.9025	0.9195	0.9145	0.9140
		0.9	0.1	Min	0.8875	0.9985	0.8920	0.8965	0.8890	0.8975
10	20	0.1	0.9	Max	0.9125	0.9085	0.9055	0.9065	0.9180	0.9090
		0.1	0.9	Min	0.8870	0.8800	0.8690	0.8945	0.8915	0.8870
		0.5	0.5	Max	0.9100	0.9720	0.9060	0.9150	0.9090	0.9110
		0.5	0.5	Min	0.8895	0.9480	0.8795	0.8890	0.8890	0.8920
		0.9	0.1	Max	0.9045	1.0000	0.9065	0.9160	0.9125	0.9130
		0.9	0.1	Min	0.8855	0.9985	0.8850	0.8950	0.8935	0.9030
20	3	0.1	0.9	Max	0.9065	0.9230	0.9035	0.9075	0.9100	0.8995
		0.1	0.9	Min	0.8910	0.9065	0.8870	0.8885	0.8920	0.8890
		0.5	0.5	Max	0.9155	0.9810	0.9125	0.9070	0.9125	0.9055
		0.5	0.5	Min	0.8905	0.9635	0.8865	0.8810	0.8869	0.8850
		0.9	0.1	Max	0.9120	1.0000	0.9075	0.9110	0.9085	0.9055
		0.9	0.1	Min	0.8875	0.9990	0.8940	0.8875	0.8860	0.8875
20	10	0.1	0.9	Max	0.9130	0.9175	0.9005	0.9155	0.9095	0.9050
		0.1	0.9	Min	0.8955	0.9000	0.8765	0.8945	0.8930	0.8880
		0.5	0.5	Max	0.9145	0.9780	0.9055	0.9120	0.9060	0.9080
		0.5	0.5	Min	0.8910	0.9655	0.8815	0.8885	0.8885	0.8890
		0.9	0.1	Max	0.9110	1.0000	0.9150	0.9150	0.9095	0.9090
		0.9	0.1	Min	0.8895	0.9995	0.8885	0.8915	0.8880	0.8905
20	20	0.1	0.9	Max	0.9145	0.9160	0.9040	0.9205	0.9175	0.9035
		0.1	0.9	Min	0.8945	0.8995	0.8785	0.8950	0.8935	0.8880
		0.5	0.5	Max	0.9095	0.9770	0.9050	0.9110	0.9075	0.9070
		0.5	0.5	Min	0.8940	0.9640	0.8850	0.8900	0.8855	0.8910
		0.9	0.1	Max	0.9080	1.0000	0.9160	0.9120	0.9040	0.9120
		0.9	0.1	Min	0.8880	0.9995	0.8865	0.8905	0.8925	0.8920

Table 3.

Range of simulated average interval lengths for 90% two sided intervals for E(Y_ij)

I	J	σo2	σE2	Range	WG1	WG2	T2	BL2
3	3	0.1	0.9	Max	–	270.17	163.01	5.90
Min	–	20.95	11.16	4.16
0.5	0.5	Max	–	400.58	360.94	5.72
Min	–	34.63	7.91	4.71
0.9	0.1	Max	–	27.32	14.36	4.67
Min	–	9.64	6.99	3.90
3	10	0.1	0.9	Max	–	112.36	35.70	6.09
Min	–	18.92	12.24	4.49
0.5	0.5	Max	–	81.54	44.57	6.70
Min	–	14.17	11.06	5.50
0.9	0.1	Max	–	74.68	29.54	5.65
Min	–	8.65	9.32	4.12
3	20	0.1	0.9	Max	–	141.91	45.15	6.73
Min	–	18.34	12.83	5.13
0.5	0.5	Max	–	120.61	31.58	7.97
Min	–	18.05	13.18	6.28
0.9	0.1	Max	–	72.27	32.03	6.50
Min	–	9.52	11.53	4.46
10	3	0.1	0.9	Max	10.30	11.16	3.77	3.70
Min	3.50	3.77	3.62	3.52
0.5	0.5	Max	7.67	8.27	4.70	4.26
Min	2.70	2.91	3.07	2.85
0.9	0.1	Max	3.56	3.83	5.46	3.21
Min	1.54	1.61	2.35	1.61
10	10	0.1	0.9	Max	9.97	11.00	4.51	4.37
Min	3.36	3.71	3.66	3.57
0.5	0.5	Max	7.45	8.15	7.46	5.64
Min	2.59	2.83	3.59	2.81
0.9	0.1	Max	3.46	3.79	9.52	3.55
Min	1.48	1.58	3.49	1.58
10	20	0.1	0.9	Max	9.92	10.99	5.47	5.13
Min	3.34	3.69	3.80	3.60
0.5	0.5	Max	7.41	8.08	10.21	6.42
Min	2.57	2.82	4.27	2.80
0.9	0.1	Max	3.45	3.80	13.36	3.67
Min	1.47	1.58	4.68	1.58
20	3	0.1	0.9	Max	10.08	10.40	3.65	3.62
Min	3.38	3.48	3.35	3.32
0.5	0.5	Max	7.53	7.73	4.57	4.16
Min	2.58	2.66	2.85	2.62
0.9	0.1	Max	3.43	3.53	5.35	3.03
Min	1.33	1.35	2.17	1.35
20	10	0.1	0.9	Max	9.92	10.34	4.45	4.33
Min	3.32	3.45	3.43	3.35
0.5	0.5	Max	7.40	7.66	7.40	5.54
Min	2.52	2.61	2.61	2.59
0.9	0.1	Max	3.37	3.51	9.49	3.33
Min	1.31	1.35	3.39	1.35
20	20	0.1	0.9	Max	9.90	10.30	5.43	5.09
Min	3.31	3.43	3.58	3.37
0.5	0.5	Max	7.39	7.57	10.18	6.28
Min	2.51	2.60	4.13	2.59
0.9	0.1	Max	3.37	3.50	13.33	3.40
Min	1.30	1.35	4.61	1.35

Table 4.

A data set of reflectivity (Yij) and resistivity (Xij)

Measurement	Operator 1	Operator 2	Operator 3
1	11.8	85	12.4	90	12.1	88

2	12.2	89	11.8	85	12.5	90

3	11.8	85	12.0	87	12.2	87

4	12.0	88	11.8	85	11.8	86

5	11.9	85	12.1	88	12.3	92

Table 5.

90% confidence intervals for mean response E(Yij) when Xij=85

Confidence Interval	Estimated values	Confidence Interval	Interval Length
WG2	4.026	0.092	11.832	[11.766 ; 11.898]	0.133
T2	3.775	0.095	11.826	[11.771 ; 11.880]	0.109
BL2	2.914	0.105	11.803	[11.773 ; 11.832]	0.059

REFERENCES

M. Aitkin , N.T. Longford (1986) Statistical modeling issues in school effectiveness studies (with Discussion)., J. R. Stat. Soc. [Ser A], Vol.149 (1) ; pp.1-43
R.K. Burdick , C.M. Borror , D.C. Montgomery (2005) Design and Analysis of Gauge R & R Studies., SIAM,
V. Czitrom , P.D. Spagon (1997) Case Studies for Industrial Process Improvement, SIAM,Philadelphia, V.A.,
H. Goldstein , J. Rasbash , M. Yang , G. Woodhouse , H. Pan , D. Nuttall , S. Thomas (1993) A multilevelanalysis of school examination results., Oxf. Rev. Educ., Vol.19 (4) ; pp.425-433
H. Goldstein (2011) Multilevel Statistical Models., John Wiley & Sons,
F.L. Huang (2016) Alternatives to multilevel modeling for the analysis of clustered data., J. Exp. Educ., Vol.84 (1) ; pp.175-196
F.L. Huang (2017) Multilevel modeling and ordinary least squares regression: How comparable are they? The Journal of Experimental Education., Publishedonline, Vol.22 (Feb) ; pp.265-281
N.T. Longford (1986) A fast scoring algorithm for maximum likelihood estimation in unbalanced mixed models with nested random effects., Biometrika, Vol.74 (4) ; pp.817-827
N.T. Longford (1995) Random Coefficient Models, edited in Handbook of Statistical modeling for the Social and Behavioral Sciences., Plenum Press, ; pp.519-578
C.J.M. Mass , J.J. Hox (2004) The influence of violations of assumptions on multilevel parameter estimates and their standard errors., Comput. Stat. Data Anal., Vol.46 (3) ; pp.427-440
D.M. McNeish , L.M. Stapleton (2016) The effect of small sample size on two-level model estimates: A review and illustration., Educ. Psychol. Rev., Vol.28 (2) ; pp.295-314
D.J. Park , R.K. Burdick (1993) Confidence intervals on the among group variance component in a simple linear regression model with a balanced onefold nested error structure., Commun. Stat. Theory Methods, Vol.22 (12) ; pp.3435-3452
D.J. Park , R.K. Burdick (1994) Confidence intervals on the regression coefficient in a simple linear regression model with a balanced one-fold nested error structure., Commun. Stat. Simul. Comput., Vol.23 (1) ; pp.43-58
D.J. Park , H.M. Hwang (2002) Confidence intervals for the mean response in the simple linear regression model with balanced nested error structure., Commun. Stat. Theory Methods, Vol.31 (1) ; pp.107-117
D.J. Park (2013) Alternative confidence intervals on variance components in a simple regression model with a balanced two-fold nested error structure., Commun. Stat. Theory Methods, Vol.42 (13) ; pp.2281-2291
D.J. Park , M. Yoon (2016) Confidence intervals for the regression coefficient in a simple regression model with a balanced two-fold nested error structure., Commun. Stat. Theory Methods, Vol.45 (17) ; pp.5053-5065
S.W. Raudenbush (1995) Maximum likelihood estimation for unbalanced multilevel covariance structure models via the EM algorithm., Br. J. Math. Stat. Psychol., Vol.48 (2) ; pp.359-370
S.R. Searle (1987) Linear Models for Unbalanced Data., John Wiley & Sons, Inc.,

Confidence Interval	Estimated values			Confidence Interval	Interval Length
Confidence Interval	Intercept	Slope	Mean Response	Confidence Interval	Interval Length
WG2	4.026	0.092	11.832	[11.766 ; 11.898]	0.133
T2	3.775	0.095	11.826	[11.771 ; 11.880]	0.109
BL2	2.914	0.105	11.803	[11.773 ; 11.832]	0.059