Multivariate CUSUM and EWMA Control Charts for Skewed Populations Using Weighted Standard Deviations

This article proposes a heuristic method of constructing multivariate cumulative sum and exponentially weighted moving average control charts for skewed populations based on the weighted standard deviation method which adjusts the variance–covariance matrix of quality characteristics and approximates the probability density function using several multivariate normal distributions. These control charts, however, reduce to the conventional control charts when the underlying distribution is symmetric. In-control and out-of-control average run lengths of the proposed control charts are compared with those of the conventional control charts for multivariate lognormal and Weibull distributions. Simulation results show that considerable improvements over the standard method can be achieved when the underlying distribution is skewed.     

Multivariate Control Charts Based on Hybrid Novelty Scores

We propose a new nonparametric multivariate control chart that integrates a novelty score. The proposed control chart uses as its monitoring statistic a hybrid novelty score, calculated based on the distance to local observations as well as on the distance to the convex hull constructed by its neighbors. The control limits of the proposed control chart were established based on a bootstrap method. A rigorous simulation study was conducted to examine the properties of the proposed control chart under various scenarios and compare it with existing multivariate control charts in terms of average run length (ARL) performance. The simulation results showed that the proposed control chart outperformed both the parametric and nonparametric Hotelling's T ² control charts, especially in nonnormal situations. Moreover, experimental results with real semiconductor data demonstrated the applicability and effectiveness of the proposed control chart. To increase the capability to detect small mean shift, we propose an exponentially weighted hybrid novelty score control chart. Simulation results indicated that exponentially weighted hybrid score charts outperformed the hybrid novelty score based control charts.     

The Design of the Multivariate Synthetic Exponentially Weighted Moving Average Control Chart

A multivariate synthetic exponentially weighted moving average (MSEWMA) control chart is presented in this study. The MSEWMA control chart consists of a multivariate exponentially weighted moving average (MEWMA) control chart and a conforming run length control chart. The average run length of the MSEWMA control chart is obtained using a Markov chain approach. From the numerical comparisons, it is shown that the MSEWMA control chart is more efficient than the multivariate synthetic T ² control chart and the MEWMA control chart for detecting shifts in the process mean vector.     

Bootstrap-Based T 2 Multivariate Control Charts

Control charts have been used effectively for years to monitor processes and detect abnormal behaviors. However, most control charts require a specific distribution to establish their control limits. The bootstrap method is a nonparametric technique that does not rely on the assumption of a parametric distribution of the observed data. Although the bootstrap technique has been used to develop univariate control charts to monitor a single process, no effort has been made to integrate the effectiveness of the bootstrap technique with multivariate control charts. In the present study, we propose a bootstrap-based multivariate T ² control chart that can efficiently monitor a process when the distribution of observed data is nonnormal or unknown. A simulation study was conducted to evaluate the performance of the proposed control chart and compare it with a traditional Hotelling's T ² control chart and the kernel density estimation (KDE)-based T ² control chart. The results showed that the proposed chart performed better than the traditional T ² control chart and performed comparably with the KDE-based T ² control chart. Furthermore, we present a case study to demonstrate the applicability of the proposed control chart to real situations.     

Profile Monitoring Using Nonparametric Bootstrap T 2 Control Chart

In certain statistical process control applications, performance of a product or process can be monitored effectively using a linear profile or a linear relationship between a response variable and one or more explanatory variables. In this article, we design a nonparametric bootstrap control chart for monitoring simple linear profiles based on T ² statistic. We evaluate the performance of the proposed method in phase II. The average and standard deviation of the run length under different shifts in the intercept, slope, and standard deviation are considered as the performance measures. Simulation results show that the performance of the proposed bootstrap control chart improves as the size of the available data increases.     

A Multivariate Synthetic Control Chart for Monitoring Process Mean Vector

In this article, a multivariate synthetic control chart is developed for monitoring the mean vector of a normally distributed process. The proposed chart is a combination of the Hotelling's T ² chart and Conforming Run Length chart. The operation, design, and performance of the chart are described. Average run length comparisons between some other existing control charts and the synthetic T ² chart are presented. They indicate that the synthetic T ² chart outperforms Hotelling's T ² chart and T ² chart with supplementary runs rules.     

The Design and Performance of the Multivariate Synthetic-T 2 Control Chart

One of the objectives of research in statistical process control is to obtain control charts that show few false alarms but, at the same time, are able to detect quickly the shifts in the distribution of the quality variables employed to monitor a productive process. In this article, the synthetic-T ² control chart is developed, which consists of the simultaneous use of a CRL chart and a Hotelling's T ² control chart. The ARL is calculated employing Markov chains for steady and zero-state scenarios. A procedure of optimization has been developed to obtain the optimum parameters of the synthetic-T ², for zero and steady cases, given the values of in-control ARL and magnitude of shift which needs to be detected rapidly. A comparison between (standard T ², MEWMA, T ² with variable sample size, and T ² with double sampling) charts reveals that the synthetic-T ² chart always performs better than the standard T ² chart. The comparison with the remaining charts demonstrate in which cases the performance of this new chart makes it interesting to employ in real applications.     

Monitoring Multi-Attribute Processes Based on NORTA Inverse Transformed Vectors

Although multivariate statistical process control has been receiving a well-deserved attention in the literature, little work has been done to deal with multi-attribute processes. While by the NORTA algorithm one can generate an arbitrary multi-dimensional random vector by transforming a multi-dimensional standard normal vector, in this article, using inverse transformation method, we initially transform a multi-attribute random vector so that the marginal probability distributions associated with the transformed random variables are approximately normal. Then, we estimate the covariance matrix of the transformed vector via simulation. Finally, we apply the well-known T ² control chart to the transformed vector. We use some simulation experiments to illustrate the proposed method and to compare its performance with that of the deleted-Y method. The results show that the proposed method works better than the deleted-Y method in terms of the out-of-control average run length criterion.     

Multivariate control chart based on the highest possibility region

The T ² control chart is widely adopted in multivariate statistical process control. However, when dealing with asymmetrical or multimodal distributions using the traditional T ² control chart, some points with relatively high occurrence possibility might be excluded, while some points with relatively low occurrence possibility might be accepted. Motived by the thought of the highest posterior density credible region, we develop a control chart based on the highest possibility region to solve this problem. It is shown that the proposed multivariate control chart will not only meet the false alarm requirement, but also ensure that all the in-control points are with relatively high occurrence possibility. The advantages and effectiveness of the proposed control chart are demonstrated by some numerical examples in the end.     

A Multivariate Adaptive Exponentially Weighted Moving Average Control Chart

A multivariate extension of the adaptive exponentially weighted moving average (AEWMA) control chart is proposed. The new multivariate scheme can detect small and large shifts in the process mean vector effectively. The proposed scheme can be viewed as a smooth combination of a multivariate exponentially weighted moving average (MEWMA) chart and a Shewhart χ²-chart. The optimal design of the proposed chart is given according to a pre-specified in-control average run length and two shift sizes; a small and large shift each measured in terms of the non centrality parameter. The signal resistance of the newly proposed multivariate chart is also given. Comparisons among the new chart, the MEWMA chart, and the combined Shewhart-MEWMA (S-MEWMA) chart in terms of the standard and worst-case average run length profiles are presented. In addition, the three charts are compared with respect to their worst-case signal resistance values. The proposed chart gives somewhat better worst-case ARL and signal resistance values than the competing charts. It also gives better standard ARL performance especially for moderate and large shifts. The effectiveness of our proposed chart is illustrated through an example with simulated data set.     

Economic-statistical design of synthetic double sampling T2 chart

Abstract

An economic-statistical design of the synthetic double sampling (synDS) T² chart is presented in this study. The cost function is minimized to obtain the optimal design parameters of the synDS T² chart by incorporating the statistical constraints or the constraints on the average number of samples. An example is provided and a sensitivity analysis is conducted to study the effect of model parameters on the optimal solution of the design. The numerical comparison shows that the synDS T² chart performs better than the synthetic T² chart and the multivariate exponentially weighted moving average chart, in terms of the cost.     

THE CONTROL CHART FOR INDIVIDUAL OBSERVATIONS FROM A MULTIVARIATE NON-NORMAL DISTRIBUTION

The Hotelling's T²statistic has been used in constructing a multivariate control chart for individual observations. In Phase II operations, the distribution of the T²statistic is related to the F distribution provided the underlying population is multivariate normal. Thus, the upper control limit (UCL) is proportional to a percentile of the F distribution. However, if the process data show sufficient evidence of a marked departure from multivariate normality, the UCL based on the F distribution may be very inaccurate. In such situations, it will usually be helpful to determine the UCL based on the percentile of the estimated distribution for T². In this paper, we use a kernel smoothing technique to estimate the distribution of the T²statistic as well as of the UCL of the T²chart, when the process data are taken from a multivariate non-normal distribution. Through simulations, we examine the sample size requirement and the in-control average run length of the T²control chart for sample observations taken from a multivariate exponential distribution. The paper focuses on the Phase II situation with individual observations.     

Mixed multivariate EWMA-CUSUM control charts for an improved process monitoring

Multivariate exponential weighted moving average and cumulative sum charts are the most common memory type multivariate control charts. They make use of the present and past information to detect small shifts in the process parameter(s). In this article, we propose two new multivariate control charts using a mixed version of their design setups. The plotting statistics of the proposed charts are based on the cumulative sum of the multivariate exponentially weighted moving averages. The performances of these schemes are evaluated in terms of average run length. The proposals are compared with their existing counterparts, including HotellingT², MCUSUM, MEWMA, and MC1 charts. An application example is also presented for practical considerations using a real dataset.     

Integration of support vector machines and control charts for multivariate process monitoring

Statistical process control tools have been used routinely to improve process capabilities through reliable on-line monitoring and diagnostic processes. In the present paper, we propose a novel multivariate control chart that integrates a support vector machine (SVM) algorithm, a bootstrap method, and a control chart technique to improve multivariate process monitoring. The proposed chart uses as the monitoring statistic the predicted probability of class (PoC) values from an SVM algorithm. The control limits of SVM-PoC charts are obtained by a bootstrap approach. A simulation study was conducted to evaluate the performance of the proposed SVM–PoC chart and to compare it with other data mining-based control charts and Hotelling's T ² control charts under various scenarios. The results showed that the proposed SVM–PoC charts outperformed other multivariate control charts in nonnormal situations. Further, we developed an exponential weighed moving average version of the SVM–PoC charts for increasing sensitivity to small shifts.     

A Comparison Between Two Multivariate EWMA Schemes

The signal issued by a control chart triggers the process professionals to investigate the special cause. Change point methods simplify the efforts to search for and identify the special cause. In this study, using maximum likelihood estimation, a multivariate joint change point estimation procedure for monitoring both location and dispersion simultaneously is proposed. After a signal is generated by the simultaneously used Hotelling's T ² and/or generalized variance control charts, the procedure starts detecting the time of the change. The performance of the proposed method for several structural changes for the mean vector and covariance matrix is discussed.     

Skewness Reduction Approach in Multi-Attribute Process Monitoring

Since the product quality of many industrial processes depends upon more than one dependent variable or attribute, they are either multivariate or multi-attribute in nature. Although multivariate statistical process control is receiving increased attention in the literature, little work has been done to deal with multi-attribute processes. In this article, we develop a new methodology to monitor multi-attribute processes. To do this, first we transform multi-attribute data in a way that their marginal probability distributions have almost zero skewness. Then, we estimate the transformed covariance matrix and apply the well-known T ² control chart. In order to illustrate the proposed method and evaluate its performance, we use two simulation experiments and compare the results with the ones from both MNP chart and the χ² control chart.     

Combined synthetic and |S| chart for monitoring process dispersion

This article proposes a multivariate control chart, the syn-|S| chart, which comprises a standard |S| subchart and a multivariate synthetic sample generalized variance |S| (synthetic |S|) subchart, for detecting shifts in the covariance matrix of a multivariate normally distributed process. A procedure for the optimal design of the syn-|S| chart by minimizing the average extra quadratic loss is provided. The syn-|S| chart has better overall performance compared to the synthetic |S| chart and the standard |S| chart, based on the zero-state and steady-state modes. An example is given to illustrate the operation of the synthetic |S| chart.     

Hotelling’s T2 control chart with double warning lines

Recent studies have shown that the T ² control chart with variable sampling intervals (VSI) and/or variable sample sizes (VSS) detects process shifts faster than the traditional T ² chart. This article extends these studies for processes that are monitored with VSI and VSS using double warning lines (T ² —DWL). It is assumed that the length of time the process remains in control has exponential distribution. The properties of T ² —DWL chart are obtained using Markov chains. The results show that the T ² —DWL chart is quicker than VSI and/or VSS charts in detecting almost all shifts in the process mean.     

Economic design based on a multivariate Bayesian VSI chart with a dual control limit

In this paper, a multivariate Bayesian variable sampling interval (VSI) control chart for the economic design and optimization of statistical parameters is designed. Based on the VSI sampling strategy of a multivariate Bayesian control chart with dual control limits, the optimal expected cost function is constructed. The proposed model allows the determination of the scheme parameters that minimize the expected cost per time of the process. The effectiveness of the Bayesian VSI chart is estimated through economic comparisons with the Bayesian fixed sampling interval and the Hotelling's T² chart. This study is an in-depth study on a Bayesian multivariate control chart with variable parameter. Furthermore, it is shown that significant cost improvement may be realized through the new model.