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.     

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.     

Statistically constrained economic design of the multivariate exponentially weighted moving average control chart

This paper shows that economic statistical design can provide for better statistical properties without significantly increasing optimal total costs. Cost comparisons between optimal economic statistical designs and optimal economic designs show no significant cost increases. The average run length (ARL) constraints added by economic statistical design significantly improve the statistical properties of the control chart scheme. False alarm frequency is limited while keeping good shift detection characteristics. In addition, the Multivariate Exponentially Weighted Moving Average (MEWMA) control schemes performed better from the cost standpoint than the benchmark pure statistical design—the Hotelling T² control chart. This improvement held for unconstrained and constrained designs. Finally, cost comparisons at small values of n showed significant advantage for the MEWMA schemes. Copyright © 2001 John Wiley & Sons, Ltd.     

Statistical process control for monitoring non-linear profiles using wavelet filtering and B-Spline approximation

A statistical process control framework is proposed to monitor non-linear profiles. The proposed methodology aims at identifying mean shifts or ‘shape changes’ in a profile. Discrete wavelet transformation (DWT) is applied to separate variation or noise from profile contours. B-splines are adopted to generate critical points to define the shape of a profile. The proposed method is innovative in that users can divide a profile into multiple segments to be monitored simultaneously. The high dimensionality problem that hinders the implementation of multivariate control charts can be solved by this framework. The distance difference statistic for each segment provides diagnostic information when the process of interest is out of control. These proposed statistics form a vector to be fed into any multivariate control chart such as the Hotelling T ² control chart. A decomposition method can also be applied on the T ² statistics when an out-of-control profile is detected. A simulation study applied to a forging process is conducted to demonstrate the property of the proposed method. The proposed method is capable of detecting profile shifts and identifying the exact location of problematic segments.     

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.     

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.     

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.     

Economic design of the integrated multivariate EPC and multivariate SPC charts

The goal of engineering process control (EPC) is to minimize variability by adjusting some manipulative process variables. The goal of statistical process control (SPC) is to reduce variability by monitoring and eliminating assignable causes of variation. As suggested by Box and Kramer and others, it is possible to reduce both special cause and common cause variations by integrating EPC and SPC. In the integrated multivariate EPC (MEPC) and multivariate SPC (MSPC) charts, we propose some statistical and economic criteria, such as the average Euclidean distance from the target vector and the average quality cost (AQC) to evaluate the performance of the MEPC/MSPC charts. The traditional average run length (ARL), average Euclidean distance and AQC of three MSPC charts are investigated and compared. The results of the simulations show that the MEPC/MGWMA chart is more effective and more economical than both the MEPC/MEWMA chart and the MEPC/Hotelling multivariate chart in detecting small shifts of the mean vector. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

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.     

Robust T2 control chart using median-based estimators

One of the most widely used multivariate control charts is the Hotelling T². In order to construct a Hotelling T² control chart, the mean vector (μ) and the variance–covariance matrix (Σ) must be first estimated. The classical estimators of μ and Σ are usually used to design Hotelling T² control chart. The classical estimators are sensitive to the presence of outliers. One way to deal with outliers is to use robust estimators. In this study, a robust T² control chart is proposed. The mean vector is obtained from the sample median. The median absolute deviation and the comedian are used as the estimates of the elements of the variance–covariance matrix. The proposed robust estimators of the mean vector and the variance–covariance matrix are compared with the sample mean vector and the sample variance–covariance matrix, and the M estimator of these parameters, through efficiency and robustness measures. The performances of the proposed robust T² control chart and the classical and the M estimators are also compared by means of average run length. Simulation results reveal that the proposed robust T² control chart has much better performance than the traditional Hotelling T² and similar performance to the M estimator in detecting shifts in process mean vector. Use of other robust estimators to estimate the process parameters is an area for further research.     

Performance of the MEWMA-CoDa control chart in the presence of measurement errors

This paper can be considered as an extension of the work of Tran et al (for monitoring compositional data using a multivariate exponentially weighted moving average MEWMA-compositional data [CoDa] chart) by taking into account potential measurement errors that are known to highly affect production processes. A linearly covariate error model with a constant error variance is used to study the impact of measurement errors on the MEWMA-CoDa control chart. In particular, the influence of the device parameters (σ_M,b), the number of independent observations m, and the the number of variables p are investigated in terms of the MEWMA optimal couples (r,H) as well as in terms of their corresponding ARLs. A comparison between the Hotelling-CoDa T² and the proposed chart is made in order to show that the MEWMA-CoDa chart is more efficient in detecting shifts in the presence of measurement errors. A real-life example of muesli production, using multiple measurements for each composition, is used to estimate the parameters and also to demonstrate how the MEWMA-CoDa can handle measurement errors to detect shifts in the process.     

MEWMA Control Chart and Process Capability Indices for Simple Linear Profiles with Within‐profile Autocorrelation

A multivariate exponentially weighted moving average (MEWMA) control chart is proposed for detecting process shifts during the phase II monitoring of simple linear profiles (SLPs) in the presence of within‐profile autocorrelation. The proposed control chart is called MEWMA‐SLP. Furthermore, two process capability indices are proposed for evaluating the capability of in‐control SLP processes, and their utilization is demonstrated through examples. Intensive simulations reveal that the MEWMA‐SLP chart is more sensitive than existing control charts in detecting profile shifts. Copyright © 2016 John Wiley & Sons, Ltd.     

Optimal T2 Control Chart with a Double Sampling Scheme – An Alternative to the MEWMA Chart

Recent studies have shown that a double sampling (DS) scheme yields improvements in detection times of process shifts over variable ratio sampling (VRS) methods that have been extensively studied in the literature. Additionally, a DS scheme is more practical than some of the VRS methods since the sampling interval is fixed. In this paper, we investigate the effect of double sampling on cost, a criterion as important as detection rate. We study economic statistical design of the DS T² chart (ESD DS T²) so that designs are found that are economically optimal but yet meet desired statistical properties such as having low probabilities of false searches and high probabilities of rapid detection of process shifts. Through an illustrative example, we show that relatively large benefits can be achieved in a comparison with the classical T² chart and the statistical DS T² charts with our ESD DS T² approach. Furthermore, the economic performance of the ESD DS T² charts is favorably compared to the MEWMA and other VRS T² control charts in the literature. Copyright © 2011 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.     

A new adaptive EWMA control chart for monitoring the process dispersion

The exponentially weighted moving average (EWMA), cumulative sum (CUSUM), and adaptive EWMA (AEWMA) control charts have had wide popularity because of their excellent speed in tracking infrequent process shifts, which are expected to lie within certain ranges. In this paper, we propose a new AEWMA dispersion chart that may achieve better performance over a range of dispersion shifts. The idea is to first consider an unbiased estimator of the dispersion shift using the EWMA statistic, and then based on the magnitude of this shift, select an appropriate value of the smoothing parameter to design an EWMA chart, named the AEWMA chart. The run length characteristics of the AEWMA chart are computed with the help of extensive Monte Carlo simulations. The AEWMA chart is compared with some of the existing powerful competitor control charts. It turns out that the AEWMA chart performs substantially and uniformly better than the EWMA‐S², CUSUM‐S², existing AEWMA, and HHW‐EWMA charts when detecting different kinds of shifts in the process dispersion. Moreover, an example is also used to explain the working and implementation of the proposed AEWMA chart.     

Multivariate multinomial T2 control chart using fuzzy approach

Quality of a product is often measured through various quality characteristics generally correlated. Multivariate control charts are a response to the need for quality control in such situations. If quality characteristics are qualitative, it sometimes happens that the product quality is defined by linguistic variables – where quality levels are represented by some specific words – and product units are classified into several linguistic forms categories, depending on the degree of fulfilment of expectations, creating a situation of fuzzy classifications. This study first reviews the concepts found in the literature on the development of fuzzy multivariate control charts. We propose a method to control these fuzzy quality evaluations, with correlated multiple attributes quality characteristics, through the use of a Hotelling T² control chart.     

The T2 chart with mixed samples to control bivariate autocorrelated processes

In this paper, we propose the use of the T² chart with the mixed sampling strategy (MS) to monitor the mean vector of bivariate processes with observations that fit to a first-order vector autoregressive model. With the MS, rational subgroups of size n are taken from the process and the selected units are regrouped to form the mixed samples. The units of the mixed samples are units selected from the last two rational subgroups. The aim of the proposed sampling strategy is to reduce the negative effect of the autocorrelation on the performance of the T² chart. When the two variables are autocorrelated, the MS always enhances the T² chart performance, however, the mixed samples are not recommended for bivariate processes with only one autocorrelated variable which is rarely affected by the assignable cause.     

Multivariate multiple sampling charts

A new multivariate statistical process control scheme, the Multivariate Multiple Sampling (MMS) control chart scheme, is proposed in this paper. A MMS chart is a multivariate extension of a double sampling X-bar control chart with at least two sampling stages. In the paper, a statistical design optimization procedure to design the MMS chart is presented and the performance of the MMS chart is investigated. The statistical efficiency in terms of average run length of the MMS chart is compared with that of the Hotelling chart both with and without variable sampling schemes, a multivariate CUMulative SUM (CUSUM) chart, and a multivariate Exponentially Weighted Moving Average (EWMA) chart. The ability of the MMS chart to handle the worst-case scenario is also investigated and compared with that of the multivariate EWMA and CUSUM charts. The results of the investigation show that even with only two sampling stages, the MMS chart provides an improvement in efficiency in detecting small shifts over the Hotelling chart without variable sampling schemes. When the number of sampling stages is equal to two, the MMS chart is better in detecting large shifts and the multivariate EWMA and CUSUM charts are better in detecting relatively small shifts. As the number of sampling stages is increased beyond two, the improvement in sensitivity of the MMS chart in detecting the small shifts increases. When the number of sampling stages ≥3, the MMS chart begins to give a better performance than a Hotelling chart with a variable sampling scheme for small shifts and is also better than a multivariate EWMA chart for both small and large shifts. As the number of sampling stages ≥4, the MMS chart begins to give a better performance than a multivariate CUSUM chart for both small and large shifts. The results of the investigation also show that the MMS chart outperforms the multivariate EWMA and CUSUM charts in the worst-case scenario situation.