Testing for lack of fit in blocked, split-plot and other multi-stratum designs

TY - JOUR

T1 - Testing for lack of fit in blocked, split-plot and other multi-stratum designs

AU - Goos, Peter

AU - Gilmour, Steven G.

PY - 2017/11

Y1 - 2017/11

N2 - Textbooks on response surface methodology emphasize the importance of lack-of-fit tests when fitting response surface models, and stress that, to be able to test for lack of fit, designed experiments should have replication and allow for pure-error estimation. In this paper, we show how to obtain pure-error estimates and how to carry out a lack-of-fit test when the experiment is not completely randomized, but a blocked experiment, a split-plot experiment, or any other multi-stratum experiment. Our approach to calculating pure-error estimates is based on residual maximum likelihood (REML) estimation of the variance components in a full treatment model (sometimes also referred to as a cell means model). It generalizes the one suggested by Vining et al. (2005) in the sense that it works for a broader set of designs and for replicates other than center point replicates. Our lack-of-fit test also generalizes the test proposed by Khuri (1992) for data from blocked experiments because it exploits replicates other than center point replicates and works for split-plot and other multi-stratum designs as well. We provide analytical expressions for the test statistic and the corresponding degrees of freedom, and demonstrate how to performthe lack-of-fit test in the SAS procedure MIXED. We re-analyze several publisheddata sets and discover a few instances in which the usual response surface model exhibits significant lack of fit.

AB - Textbooks on response surface methodology emphasize the importance of lack-of-fit tests when fitting response surface models, and stress that, to be able to test for lack of fit, designed experiments should have replication and allow for pure-error estimation. In this paper, we show how to obtain pure-error estimates and how to carry out a lack-of-fit test when the experiment is not completely randomized, but a blocked experiment, a split-plot experiment, or any other multi-stratum experiment. Our approach to calculating pure-error estimates is based on residual maximum likelihood (REML) estimation of the variance components in a full treatment model (sometimes also referred to as a cell means model). It generalizes the one suggested by Vining et al. (2005) in the sense that it works for a broader set of designs and for replicates other than center point replicates. Our lack-of-fit test also generalizes the test proposed by Khuri (1992) for data from blocked experiments because it exploits replicates other than center point replicates and works for split-plot and other multi-stratum designs as well. We provide analytical expressions for the test statistic and the corresponding degrees of freedom, and demonstrate how to performthe lack-of-fit test in the SAS procedure MIXED. We re-analyze several publisheddata sets and discover a few instances in which the usual response surface model exhibits significant lack of fit.

U2 - 10.1080/00224065.2017.11918000

DO - 10.1080/00224065.2017.11918000

M3 - Article

SN - 0022-4065

VL - 49

SP - 320

EP - 336

JO - JOURNAL OF QUALITY TECHNOLOGY

JF - JOURNAL OF QUALITY TECHNOLOGY

IS - 4

ER -