Determining Optimal Number of Samples for Constructing Multivariate Control Charts

Normally, an average run length (ARL) is used as a measure for evaluating the detecting performance of a multivariate control chart. This has a direct impact on the false alarm cost in Phase II. In this article, we first conduct a simulation study to calculate both in-control and out-of-control ARLs under various combinations of process shifts and number of samples. Then, a trade-off analysis between sampling inspection and false alarm costs is performed. Both the simulation results and trade-off analysis suggest that the optimal number of samples for constructing a multivariate control chart in Phase I can be determined.     

A Phase I nonparametric Shewhart-type control chart based on the median

A nonparametric Shewhart-type control chart is proposed for monitoring the location of a continuous variable in a Phase I process control setting. The chart is based on the pooled median of the available Phase I samples and the charting statistics are the counts (number of observations) in each sample that are less than the pooled median. An exact expression for the false alarm probability (FAP) is given in terms of the multivariate hypergeometric distribution and this is used to provide tables for the control limits for a specified nominal FAP value (of 0.01, 0.05 and 0.10, respectively) and for some values of the sample size (n) and the number of Phase I samples (m). Some approximations are discussed in terms of the univariate hypergeometric and the normal distributions. A simulation study shows that the proposed chart performs as well as, and in some cases better than, an existing Shewhart-type chart based on the normal distribution. Numerical examples are given to demonstrate the implementation of the new chart.     

Self-starting X-bar control chart based on Six Sigma quality and sometimes pooling procedure

Self-starting control charts have been proposed in the literature to allow process monitoring when only a small amount of relevant data is available. In fact, self-starting charts are useful in monitoring a process quickly, without having to collect a sizable Phase I sample for estimating the in-control process parameters. In this paper, a new self-starting control charting procedure is proposed in which first an effective initial sample is chosen from the perspective of Six Sigma quality, then the successive sample means are either pooled or not pooled (sometimes pooling procedure) for computing next Q-statistics depending upon its signal. It is observed that the sample statistics obtained so from this in-control Phase I situation can serve as more efficient estimators of unknown parameters for Phase II monitoring. An example is considered to illustrate the construction of the proposed chart and to compare its performance with the existing ones.     

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.     

Multivariate Control Chart Based on Multivariate Smirnov Test

Robust control charts are useful in statistical process control (SPC) when there is limited knowledge about the underlying process distribution, especially for multivariate observations. This article develops a new robust and self-starting multivariate procedure based on multivariate Smirnov test (MST), which integrates a multivariate two-sample goodness-of-fit (GOF) test based on multivariate empirical distribution function (MEDF) and the change-point model. As expected, simulation results show that our proposed control chart is robust to nonnormally distributed data, and moreover, it is efficient in detecting process shifts, especially large shifts, which is one of the main drawbacks of most robust control charts in the literature. As it avoids the need for a lengthy data-gathering step, the proposed chart is particularly useful in start-up or short-run situations. Comparison results and a real data example show that our proposed chart has great potential for application.     

The performance of multivariate CUSUM control charts with estimated parameters

In this paper, we study the effect of estimating the vector of means and the variance–covariance matrix on the performance of two of the most widely used multivariate cumulative sum (CUSUM) control charts, the MCUSUM chart proposed by Crosier [Multivariate generalizations of cumulative sum quality-control schemes, Technometrics 30 (1988), pp. 291–303] and the MC1 chart proposed by Pignatiello and Runger [Comparisons of multivariate CUSUM charts, J. Qual. Technol. 22 (1990), pp. 173–186]. Using simulation, we investigate and compare the in-control and out-of-control performances of the competing charts in terms of the average run length measure. The in-control and out-of-control performances of the competing charts deteriorate significantly if the estimated parameters are used with control limits intended for known parameters, especially when only a few Phase I samples are used to estimate the parameters. We recommend the use of the MC1 chart over that of the MCUSUM chart if the parameters are estimated from a small number of Phase I samples.     

Conditional design of the EWMA median chart with estimated parameters

The exponentially weighted moving average (EWMA) chart is often designed assuming the process parameters are known. In practice, the parameters are rarely known and need to be estimated from Phase I samples. Different Phase I samples are used when practitioners construct their own control chart's limits, which leads to the “Phase I between-practitioners” variability in the in-control average run length (ARL) of control charts. The standard deviation of the ARL (SDARL) is a good alternative to quantify this variability in control charts. Based on the SDARL metric, the performance of the EWMA median chart with estimated parameters is investigated in this paper. Some recommendations are given based on the SDARL metric. The results show that the EWMA median chart requires a much larger amount of Phase I data in order to reduce the variation in the in-control ARL up to a reasonable level. Due to the limitation of the amount of the Phase I data, the suggested EWMA median chart is designed with the bootstrap method which provides a good balance between the in-control and out-of-control ARL values.     

Multivariate Process Variability Monitoring

Understanding multivariate variability is a difficult task because there is no single measure that can be properly used. This article presents a new measure that features good properties. If this measure is simultaneously used with generalized variance, it will give a better understanding of multivariate variability. It can also efficiently be used for large data sets with high dimensions. Furthermore, when it is used for constructing a Shewhart-type chart to monitor multivariate variability, the resulting chart has a much better out-of-control ARL than the generalized variance chart. An example illustrates its advantage.     

A Semi‐empirical Bayesian Chart to Monitor Weibull Percentiles

This paper develops a Bayesian control chart for the percentiles of the Weibull distribution, when both its in‐control and out‐of‐control parameters are unknown. The Bayesian approach enhances parameter estimates for small sample sizes that occur when monitoring rare events such as in high‐reliability applications. The chart monitors the parameters of the Weibull distribution directly, instead of transforming the data as most Weibull‐based charts do in order to meet normality assumption. The chart uses accumulated knowledge resulting from the likelihood of the current sample combined with the information given by both the initial prior knowledge and all the past samples. The chart is adapting because its control limits change (e.g. narrow) during Phase I. An example is presented and good average run length properties are demonstrated.     

Control chart for monitoring multivariate COM-Poisson attributes

Statistical process control of multi-attribute count data has received much attention with modern data-acquisition equipment and online computers. The multivariate Poisson distribution is often used to monitor multivariate attributes count data. However, little work has been done so far on under- or over-dispersed multivariate count data, which is common in many industrial processes, with positive or negative correlation. In this study, a Shewhart-type multivariate control chart is constructed to monitor such kind of data, namely the multivariate COM-Poisson (MCP) chart, based on the MCP distribution. The performance of the MCP chart is evaluated by the average run length in simulation. The proposed chart generalizes some existing multivariate attribute charts as its special cases. A real-life bivariate process and a simulated trivariate Poisson process are used to illustrate the application of the MCP chart.     

A control chart for multivariate Poisson distribution using repetitive sampling

Control charts using repetitive group sampling have attracted a great deal of attention during the last few years. In the present article, we attempt to develop a control chart for the multivariate Poisson distribution using the repetitive group sampling scheme. In the proposed control chart, the monitoring statistic from the multivariate Poisson distribution has been used for the quick detection of the deteriorated process to avoid losses. The control coefficients have been estimated using the specified in-control average run lengths. The procedure of the proposed control chart has been explained by using the real-world example and a simulated data set. It has been observed that the proposed control chart is an efficient development for the quick detection of the nonrandom change in the manufacturing process.     

Addressing the effect of parameter estimation on phase II monitoring of multivariate multiple linear profiles via a new cluster-based approach

Abstract

Profile monitoring is applied when the quality of a product or a process can be determined by the relationship between a response variable and one or more independent variables. In most Phase II monitoring approaches, it is assumed that the process parameters are known. However, it is obvious that this assumption is not valid in many real-world applications. In fact, the process parameters should be estimated based on the in-control Phase I samples. In this study, the effect of parameter estimation on the performance of four Phase II control charts for monitoring multivariate multiple linear profiles is evaluated. In addition, since the accuracy of the parameter estimation has a significant impact on the performance of Phase II control charts, a new cluster-based approach is developed to address this effect. Moreover, we evaluate and compare the performance of the proposed approach with a previous approach in terms of two metrics, average of average run length and its standard deviation, which are used for considering practitioner-to-practitioner variability. In this approach, it is not necessary to know the distribution of the chart statistic. Therefore, in addition to ease of use, the proposed approach can be applied to other type of profiles. The superior performance of the proposed method compared to the competing one is shown in terms of all metrics. Based on the results obtained, our method yields less bias with small-variance Phase I estimates compared to the competing approach.     

Control charts for fraction nonconforming in a bivariate binomial process

Many multivariate quality control techniques are used for multivariate variable processes, but few work for multivariate attribute processes. To monitor multivariate attributes, controlling the false alarms (type I errors) and considering the correlation between attributes are two important issues. By taking into account these two issues, a new control chart is presented to monitor a bivariate binomial process. An example is illustrated for the proposed method. To evaluate the performance of the proposed method, a simulation study is conducted to compare the results with those using both the multivariate np chart and skewness reduction approaches. The results show that the correlation is taken into account in the designed chart and the overall false alarm is controlled at the nominal value. Moreover, the process shift can be quickly detected and the variable that is responsible for a signal can be determined.     

Monitoring Nonlinear Profiles with Random Effects by Nonparametric Regression

The monitoring of process/product profiles is presently a growing and promising area of research in statistical process control. This study is aimed at developing monitoring schemes for nonlinear profiles with random effects. We utilize the technique of principal components analysis to analyze the covariance structure of the profiles and propose monitoring schemes based on principal component (PC) scores. The number of the PC scores used in constructing control charts is crucial to the detecting power. In the Phase I analysis of historical data, due to the dependency of the PC-scores, we adopt the usual Hotelling T ² chart to check the stability. For Phase II monitoring, we study individual PC-score control charts, a combined chart scheme that combines all the PC-score charts, and a T ² chart. Although an individual PC-score chart may be perfect for monitoring a particular mode of variation, a chart that can detect general shifts, such as the T ² chart and the combined chart scheme, is more feasible in practice. The performances of the schemes under study are evaluated in terms of the average run length.     

A CUSUM control chart for monitoring the variance when parameters are estimated

CUSUM control chart has been widely used for monitoring the process variance. It is usually used assuming that the nominal process variance is known. However, several researchers have shown that the ability of control charts to signal when a process is out of control is seriously affected unless process parameters are estimated from a large in-control Phase I data set. In this paper we derive the run length properties of a CUSUM chart for monitoring dispersion with estimated process variance and we evaluate the performance of this chart by comparing it with the same chart but with assumed known process parameters.     

Conditional design of the CUSUM median chart for the process position when process parameters are unknown

In practice, different practitioners will use different Phase I samples to estimate the process parameters, which will lead to different Phase II control chart's performance. Researches refer to this variability as between-practitioners-variability of control charts. Since between-practitioners-variability is important in the design of the CUSUM median chart with estimated process parameters, the standard deviation of average run length (SDARL) will be used to study its properties. It is shown that the CUSUM median chart requires a larger amount of Phase I samples to sufficiently reduce the variation in the in-control ARL of the CUSUM median chart. Considering the limitation of the amount of the Phase I samples, a bootstrap approach is also used here to adjust the control limits of the CUSUM median chart. Comparisons are made for the CUSUM and Shewhart median charts with estimated parameters when using the adjusted- and unadjusted control limits and some conclusions are made.     

DDMA-charts: Nonparametric multivariate moving average control charts based on data depth

Summary: This paper studies the DDMA–chart, a data depth based moving–average control chart for monitoring multivariate data. This chart is nonparametric and it can detect simultaneously location and scale changes in the process. It improves upon the existing r– and Q–chart in the efficiency of detecting location changes. Both theoretical justifications and simulation studies are provided. Comparisons with some existing multivariate control charts via simulation results are also provided. Some applications of the DDMA–chart to the analysis of airline performance data (collected by the FAA) are demonstrated. The results indicate that the DDMA–chart is an effective nonparametric multivariate control chart.*Research supported in part by grants from the National Science Foundation, the National Security Agency, and the Federal Aviation Administration. The discussion on aviation safety in this paper reects the views of the authors, who are solely responsible for the accuracy of the analysis results presented herein, and does not necessarily reect the official view or policy of the FAA. The dataset used in this paper has been partially masked in order to protect confidentiality.     

Multivariate Synthetic |S| Control Chart with Variable Sampling Interval

A variable sampling interval (VSI) feature is introduced to the multivariate synthetic generalized sample variance |S| control chart. This multivariate synthetic control chart is a combination of the |S| sub-chart and the conforming run length sub-chart. The VSI feature enhances the performance of the multivariate synthetic control chart. The comparative results show that the VSI multivariate synthetic control chart performs better than other types of multivariate control charts for detecting shifts in the covariance matrix of a multivariate normally distributed process. An example is given to illustrate the operation of the VSI multivariate synthetic chart.     

Using the predictive distribution to determine control limits for the Bayesian MEWMA chart

Bayesian control charts have been proposed for monitoring multivariate processes with the multivariate exponentially weighted moving average (MEWMA) statistic. It has been suggested that we use limits based on the predictive distribution of the MEWMA statistic. This analysis, however is based on the erroneous result that the average run length (ARL) is a function of the exceedance probability, that is, the probability that the first point exceeds the control limit. We show how this result can be corrected and we discuss how the Bayesian MEWMA chart with limits based on the predictive distribution compares with other multivariate control chart procedures.     

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.