Optimal designs of the multivariate synthetic chart for monitoring the process mean vector based on median run length

The average run length (ARL) is usually used as a sole measure of performance of a multivariate control chart. The Hotelling's T², multivariate exponentially weighted moving average (MEWMA) and multivariate cumulative sum (MCUSUM) charts are commonly optimally designed based on the ARL. Similar to the case of univariate quality control, in multivariate quality control, the shape of the run length distribution changes in accordance to the magnitude of the shift in the mean vector, from highly skewed when the process is in‐control to nearly symmetric for large shifts. Because the shape of the run length distribution changes with the magnitude of the shift in the mean vector, the median run length (MRL) provides additional and more meaningful information about the in‐control and out‐of‐control performances of multivariate charts, not given by the ARL. This paper provides a procedure for optimal designs of the multivariate synthetic T² chart for the process mean, based on MRL, for both the zero and steady‐state modes. Two Mathematica programs, each for the zero state and steady‐state modes are given for a quick computation of the optimal parameters of the synthetic T² chart, designed based on MRL. These optimal parameters are provided in the paper, for the bivariate case with sample sizes, nin{4, 7, 10}. The MRL performances of the synthetic T², MEWMA and Hotelling's T² charts are also compared. Copyright © 2011 John Wiley & Sons, Ltd.     

Economic statistical design of a T2 control chart with double warning lines

Recent studies have shown that enhancing the common T² control chart by using variable sample sizes (VSS) and variable sample intervals (VSI) sampling policies with a double warning line scheme (DWL) yields improvements in shift detection times over either pure VSI or VSS schemes in detecting almost all shifts in the process mean. In this paper, we look at this problem from an economical perspective, certainly at least as an important criterion as shift detection time if one considers what occurs in the industry today. Our method is to first construct a cost model to find the economic statistical design (ESD) of the DWL T² control chart using the general model of Lorenzen and Vance (Technometrics 1986; 28 :3–11). Subsequently, we find the values of the chart parameters which minimize the cost model using a genetic algorithm optimization method. Cost comparisons of Fixed ratio sampling, VSI, VSS, VSIVSS with DWL, and multivariate exponentially weighted moving average (MEWMA) charts are made, which indicate the economic efficacy of using either VSIVSS with DWL or MEWMA charts in practice if cost minimization is of interest to the control chart user. Copyright © 2010 John Wiley & Sons, Ltd.     

Model‐based Multivariate Monitoring Charts for Autocorrelated Processes

Autocorrelation or nonstationarity may seriously impact the performance of conventional Hotelling's T² charts. We suggest modeling processes with multivariate autoregressive integrated moving average time series models and propose two model‐based monitoring charts. One monitors the predicted value and provides information about the need for mean adjustments. The other is a Hotelling's T² control chart applied to the residuals. The average run length performance of the residual‐based Hotelling's T² chart is compared with the observed data‐based Hotelling's T² chart for a group of first‐order vector autoregressive models. We show that the new chart in most cases performs well. Copyright © 2013 John Wiley & Sons, Ltd.     

Multivariate EWMA control charts using individual observations for process mean and variance monitoring and diagnosis

Most multivariate control charts in the literature are designed to detect either mean or variation shifts rather than both. A simultaneous use of the Hotelling T ² and |S| control charts has been proposed but the Hotelling T ² reacts to mean shifts, dispersion changes, and changes of correlations among responses. The combination of two multivariate control charts into one chart sometimes loses the ability to provide detailed diagnostic information when a process is out-of-control. In this research a new multivariate control chart procedure based on exponentially weighted moving average (EWMA) statistics is proposed to monitor process mean and variance simultaneously to identify proper sources of variations. Two multivariate EWMA control charts using individual observations are proposed to achieve a quick detection of mean or variance shifts or both. Simulation studies show that the proposed charts are capable of identifying appropriate types of shifts in terms of correct detection percentages. A manufacturing example is used to demonstrate how the proposed charts can be properly set-up based on average run length values via simulations. In addition, correct detection rates of the proposed charts are explored.     

Optimum Multiple and Multivariate Poisson Statistical Control Charts

This paper deals with the simultaneous statistical process control of several Poisson variables. The practitioner of this type of monitoring may employ a multiple scheme, i.e. one chart for controlling each variable, or may use a multivariate scheme, based on monitoring all the variables with a single control chart. If the user employs the multivariate schemes, he or she can choose from, for example, three options: (i) a control chart based on the sum of the different Poisson variables; (ii) a control chart on the maximum value of the different Poisson variables; and (iii) in the case of only two variables, a chart that monitors the difference between them. In this paper, the previous control charts are studied when applied to the control of p = 2, 3 and 4 variables. In addition, the optimization of a set of univariate Poisson control charts (multiple scheme) is studied. The main purpose of this paper is to help the practitioner to select the most adequate scheme for her/his production process. Towards this goal, a friendly Windows© computer program has been developed. The program returns the best control limits for each control chart and makes a complete comparison of performance among all the previous schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

To shrink or not to shrink: Hotelling's T2 chart based on shrunken covariance estimates

Several authors have studied the effect of parameter estimation on the performance of Phase II control charts and shown that large in‐control reference samples are necessary for the Phase II control charts to perform as desired. For higher dimensional data, even larger reference samples are required to achieve stable estimation of the in‐control parameters. Shrinkage estimation has been widely studied as a method to achieve stable estimation of the covariance matrix for high‐dimensional data. We investigate the average run length (ARL) distribution of the Hotelling T² chart when using a shrunken covariance matrix. Specifically, we explore the following questions: (1) Does the use of a shrinkage estimator of the covariance matrix result in reduced variability in the ARL performance of the T² chart? (2) Does the use of a shrinkage estimator of the covariance matrix result in a reduced occurrence of “strictly multivariate” false alarms on the T²chart? (3) How does shrinkage of the covariance matrix affect the out‐of‐control performance of the T² chart? We use a simulation study to investigate the use of shrinkage estimation with the Hotelling T² chart in Phase II. Our results indicate that, while shrinkage estimation affects the ARL performance of the T² chart, the benefits are small and occur in fairly specific circumstances. The benefits of shrinking may not justify the use of more advanced techniques.     

Hotelling's T2 control chart with adaptive sample sizes

Some quality control schemes have been developed when several related quality characteristics are to be monitored: simultaneous X¯ charts, Hotelling's T² chart, multivariate CUSUM and multivariate EWMA. Hotelling's T² control chart has the advantage of its simplicity but it is slow in detecting small process shifts. The latest developments in variable sample sizes for univariate control charts are applied in this paper to define an adaptive sample sizes T² control chart. As occurs in the univariate case the ARL improvements are very important particularly for small process shifts. An example is given to illustrate the use of the proposed scheme.     

Multivariate double sampling |S| charts for controlling process variability

A new multivariate statistical process variability control scheme, the multivariate double sampling (MDS) |S| control chart scheme for controlling shifts in a variance–covariance matrix, is proposed. An MDS |S| chart is a double-stage extension of a single-stage control chart based on the sample generalized variance. The statistical efficiency in terms of average run length of the MDS |S| charts is compared with that of traditional and adaptive sample size control charts based on the sample generalized variance for a bivariate case. The results of the comparison show that the MDS |S| charts provide a significant improvement in efficiency for detecting shifts of any range over these schemes. The performance of the MDS |S| charts in terms of detecting power is also compared with that of the modified likelihood ratio test, the sum of standardized variances for principal components charts and a decomposition procedure. The results of the comparison show that for certain out-of-control situations, the performance of the MDS |S| chart is among those of the best schemes.     

Optimum variable-dimension EWMA chart for multivariate statistical process control

The variable-dimension T² control chart (VDT² chart) was recently proposed for monitoring the mean of multivariate processes in which some of the quality variables are easy and inexpensive to measure while other variables require substantially more effort or expense. The chart requires most of the times that only the inexpensive variables be sampled, switching to sampling all the variables only when a warning is triggered. It has good ARL performance compared with the standard T² chart, while significantly reducing the sampling cost. However, like the T² chart, it has limited sensitivity to small and moderate mean shifts. To detect such shifts faster, we developed an exponentially weighted moving average (EWMA) version of the VDT² chart, along with Markov chain models for ARL calculation, and software (made available) for optimizing the chart design. The optimization software, which is based on genetic algorithms, runs in Windows^© and has a friendly user interface. The performance analysis shows the great gain in performance achieved by the incorporation of the EWMA procedure.     

Adaptive multivariate double sampling and variable sampling interval Hotelling's T2 charts

In this article, two adaptive multivariate charts, which combine the double sampling (DS) and variable sampling interval (VSI) features, called the adaptive multivariate double sampling variable sampling interval T² (AMDSVSI T²) and the adaptive multivariate double sampling variable sampling interval combined T² (AMDSVSIC T²) charts, are proposed. The real purpose of using the proposed charts is to provide flexibility by enabling the sampling interval length of the DS T² chart to be varied so that the chart's sensitivity can be enhanced. The fundamental difference between the two proposed charts is that when a second sample is taken, the AMDSVSI T² chart uses the information of the combined sample mean vectors while the AMDSVSIC T² chart uses the information of the combined T² statistics, in deciding about the process status. This research is motivated by existing combined DS and VSI charts in the literature, which show convincing performance improvement over the standard DS chart. Consequently, it is believed that adopting this existing approach in the multivariate case will enable superior multivariate DS charts to be proposed. Numerical results show that the proposed charts outperform the existing standard T² and other adaptive multivariate charts, in detecting shifts in the mean vector, for the zero‐state and steady‐state cases. The performances of both charts when the shift sizes in the mean vector are unknown are also measured. The application of the AMDSVSI T² chart is illustrated with an example.     

A boosting approach for understanding out-of-control signals in multivariate control charts

The most widely used tools in statistical quality control are control charts. However, the main problem of multivariate control charts, including Hotelling's T ² control chart, lies in that they indicate that a change in the process has happened, but do not show which variable or variables are the source of this shift. Although a number of methods have been proposed in the literature for tackling this problem, the most usual approach consists of decomposing the T ² statistic. In this paper, we propose an alternative method interpreting this task as a classification problem and solving it through the application of boosting with classification trees. The classifier is then used to determine which variable or variables caused the change in the process. The results prove this method to be a powerful tool for interpreting multivariate control charts.     

The Effect of Autocorrelation on the Hotelling T2 Control Chart

One of the basic assumptions for traditional univariate and multivariate control charts is that the data are independent in time. For the latter, in many cases, the data are serially dependent (autocorrelated) and cross‐correlated because of, for example, frequent sampling and process dynamics. It is well known that the autocorrelation affects the false alarm rate and the shift‐detection ability of the traditional univariate control charts. However, how the false alarm rate and the shift‐detection ability of the Hotelling T² control chart are affected by various autocorrelation and cross‐correlation structures for different magnitudes of shifts in the process mean is not fully explored in the literature. In this article, the performance of the Hotelling T² control chart for different shift sizes and various autocorrelation and cross‐correlation structures are compared based on the average run length using simulated data. Three different approaches in constructing the Hotelling T² chart are studied for two different estimates of the covariance matrix: (i) ignoring the autocorrelation and using the raw data with theoretical upper control limits; (ii) ignoring the autocorrelation and using the raw data with adjusted control limits calculated through Monte Carlo simulations; and (iii) constructing the control chart for the residuals from a multivariate time series model fitted to the raw data. To limit the complexity, we use a first‐order vector autoregressive process and focus mainly on bivariate data. © 2014 The Authors. Quality and Reliability Engineering International Published by John Wiley & Sons Ltd.     

An assessment of the kernel‐distance‐based multivariate control chart through an industrial application

Traditional multivariate quality control charts assume that quality characteristics follow a multivariate normal distribution. However, in many industrial applications the process distribution is not known, implying the need to construct a flexible control chart appropriate for real applications. A promising approach is to use support vector machines in statistical process control. This paper focuses on the application of the ‘kernel‐distance‐based multivariate control chart’, also known as the ‘k‐chart’, to a real industrial process, and its assessment by comparing it to Hotelling's T² control chart, based on the number of out‐of‐control observations and on the Average Run Length. The industrial application showed that the k‐chart is sensitive to small shifts in mean vector and outperforms the T² control chart in terms of Average Run Length. Copyright © 2010 John Wiley & Sons, Ltd.     

M‐ATTRIVAR: An attribute‐variable chart to monitor multivariate process means

In this paper, an attribute‐variable control chart, namely, M‐ATTRIVAR, is introduced to monitor possible shifts in a vector of means. The monitoring starts using an attribute chart (classifying the units as approved or not using a gauge) and continues in such a way until a warning signal is given, shifting the control to a variable chart for the next sampling. If the variable chart does not confirm the warning, the monitoring returns to an attribute control. Otherwise, the monitoring remains with the variable chart. Whenever any of the charts (attribute or variable) signals an alarm, the control scheme triggers an alarm. The main advantage of this new proposal is the possibility of judging the state of the process only by the attribute chart most of the time (normally more economical and faster). The performance of the M‐ATTRIVAR control chart is compared versus the main competitor (T² control chart) in terms of performance detection (out‐of‐control average run length) but also economically (average sampling cost). The M‐ATTRIVAR is always cheaper than T², and in many scenarios, it detects quicker process shifts than the T² control chart. A numerical example illustrates a practical situation.     

Alternative Hotelling's T2 Charts using Winsorized Modified One‐Step M‐estimator

Hotelling's T² chart is a popular tool for monitoring statistical process control. However, this chart is sensitive in the presence of outliers. To alleviate the problem, this paper proposed alternative Hotelling's T² charts for individual observations using robust location and scale matrix instead of the usual mean vector and the covariance matrix, respectively. The usual mean vector in the Hotelling T² chart is replaced by the winsorized modified one‐step M‐estimator (MOM) whereas the usual covariance matrix is replaced by the winsorized covariance matrix. MOM empirically trims the data based on the shape of the data distribution. This study also investigated on the different trimming criteria used in MOM. Two robust scale estimators with highest breakdown point, namely S_n and T_n were selected to suit the criteria. The upper control limits for the proposed robust charts were calculated based on simulated data. The performance of each control chart is based on the false alarm and the probability of outlier's detection. In general, the performance of an alternative robust Hotelling's T² charts is better than the performance of the traditional Hotelling's T² chart. Copyright © 2012 John Wiley & Sons, Ltd.     

Multivariate statistical process control charts based on the approximate sequential χ2 test

Similar to the univariate CUSUM chart, a multivariate CUSUM (MCUSUM) chart can be designed to detect a particular size of the mean shift optimally based on the scheme of a sequential likelihood ratio test for the noncentrality parameter. However, in multivariate case, the probability ratio of a sequential test is intractable mathematically and the test statistic based on the ratio does not have a closed form expression which makes it impractical for real application. We drive an approximate log-likelihood ratio and propose a multivariate statistical process control chart based on a sequential χ² test to detect a change in the noncentrality parameter. The statistical properties of the proposed test statistic are explored. The average runs length (ARL) performance of the proposed charts is compared with other MCUSUM charts for process mean monitoring. The experimental results reveal that the proposed charts achieve superior, both zero-state and steady-state, ARL performance over a wide range of mean shifts, especially when the dimension of measurements is large.     

The variable sample size variable dimension T2 control chart

In this paper, we describe the development of the variable dimension and variable sample size T² control chart (VSSVDT²), which is an enhancement of the variable dimension T² chart (VDT²). In the VDT² control chart, the number of variables that are measured to compute the T² statistic is made variable. Some of the variables are easy or inexpensive to measure and are always monitored. The variables that are more difficult or expensive to measure are measured only when the T² value from the previous sample exceeds a specified threshold. The VDT² control chart performs well for moderate and large shifts in the mean vector. However, its performance for small shifts is poor. To improve the chart’s performance in detecting such shifts, we propose the application of the variable sample size technique to the VDT² control chart, resulting in the VSSVDT² control chart. To promote the use of the VSSVDT² chart, a user-friendly software has been developed, which the final user can use to find the best parameters of the VSSVDT² chart for a particular process.     

Hotelling's T2 control charts with fixed and variable sample sizes for monitoring short production runs

Short production runs are common in enterprises that require a high degree of flexibility and variety in manufacturing processes. To date, past research on short production runs has little focus on the multivariate control charts. In view of this, fixed sample size (FSS) and variable sample size (VSS) Hotelling's T² charts are designed to monitor the process mean when the production horizon is finite. Optimal parameters to minimize the out‐of‐control (1) truncated average run length (TARL) and (2) expected TARL (ETARL) are provided such that the in‐control TARL is equal to the number of inspections (say I). The numerical study considers the run length performances of the FSS and VSS T² short‐run charts for both known and unknown shift sizes. The VSS T² short‐run chart performs well in swiftly detecting various mean shifts in comparison with the FSS T² short‐run chart. Additionally, the VSS T² short‐run chart is superior to the FSS T² short‐run chart, in terms of the truncated standard deviation of the run length, expected truncated standard deviation of the run length, probability that the chart signals an alarm within the I inspections, ie, P(I) and expected P(I). A case study on the impurity profile of a crystalline drug substance illustrates the implementation of the VSS T² short‐run chart.     

A comparative study of phase II robust multivariate control charts for individual observations

Use of Hotelling's T² charts with high breakdown robust estimates to monitor multivariate individual observations are the recent trend in the control chart methodology. Vargas (J. Qual. Tech. 2003; 35: 367‐376) introduced Hotelling's T² charts based on the minimum volume ellipsoid (MVE) and the minimum covariance determinant (MCD) estimates to identify outliers in Phase I data. Studies carried out by Jensen et al. (Qual. Rel. Eng. Int. 2007; 23: 615‐629) indicated that the performance of these charts heavily depends on the sample size, amount of outliers and the dimensionality of the Phase I data. Chenouri et al. (J. Qual. Tech. 2009; 41: 259‐271) recently proposed robust Hotelling's T² control charts for monitoring Phase II data based on the reweighted MCD (RMCD) estimates of the mean vector and covariance matrix from Phase I. They showed that Phase II RMCD charts have better performance compared with Phase II standard Hotelling's T² charts based on outlier free Phase I data, where the outlier free Phase I data were obtained by applying MCD and MVE T² charts to historical data. Reweighted MVE (RMVE) and S‐estimators are two competitors of the RMCD estimators and it is a natural question whether the performance of Phase II Hotelling's T² charts with RMCD and RMVE estimates exhibits similar pattern observed by Jensen et al. (Qual. Rel. Eng. Int. 2007; 23: 615‐629) in the case of MCD and MVE‐based Phase I Hotelling's T² charts. In this paper, we conduct a comparative study to assess the performance of Hotelling's T² charts with RMCD, RMVE and S‐estimators using large number of Monte Carlo simulations by considering different data scenarios. Our results are generally in favor of the RMCD‐based charts irrespective of sample size, outliers and dimensionality of Phase I data. Copyright © 2010 John Wiley & Sons, Ltd.     

A clustering algorithm-based control chart for inhomogeneously distributed TFT-LCD processes

Statistical process control techniques have been widely used to improve processes by reducing variations and defects. In the present paper, we propose a multivariate control chart technique based on a clustering algorithm that can effectively handle a situation in which the distribution of in-control observations is inhomogeneous. A simulation study was conducted to examine the characteristics of the proposed control chart and to compare them with Hotelling’s T ² multivariate control charts that are widely used in real-world processes. Moreover, an experiment with real data from the thin film transistor liquid crystal display (TFT-LCD) manufacturing process demonstrated the effectiveness and accuracy of the proposed control chart.