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.     

The economic design of multivariate binomial EWMA VSSI control charts

Since multi-attribute control charts have received little attention compared with multivariate variable control charts, this research is concerned with developing a new methodology to employ the multivariate exponentially weighted moving average (MEWMA) charts for m-attribute binomial processes; the attributes being the number of nonconforming items. Moreover, since the variable sample size and sampling interval (VSSI) MEWMA charts detect small process mean shifts faster than the traditional MEWMA, an economic design of the VSSI MEWMA chart is proposed to obtain the optimum design parameters of the chart. The sample size, the sampling interval, and the warning/action limit coefficients are obtained using a genetic algorithm such that the expected total cost per hour is minimized. At the end, a sensitivity analysis has been carried out to investigate the effects of the cost and the model parameters on the solution of the economic design of the VSSI MEWMA chart.     

A modified economic-statistical design of the T2 control chart with variable sample sizes and control limits

Recent studies have shown that using variable sampling size and control limits (VSSC) schemes result in charts with more statistical power than variable sampling size (VSS) when detecting small to moderate shifts in the process mean vector. This paper presents an economic-statistical design (ESD) of the VSSC T² control chart using the general model of Lorenzen and Vance [22]. The genetic algorithm approach is then employed to search for the optimal values of the six test parameters of the chart. We then compare the expected cost per unit of time of the optimally designed VSSC chart with optimally designed VSS and FRS (fixed ratio sampling) T² charts as well as MEWMA charts.     

Optimal Statistical Design of a Multivariate EWMA Chart Based on ARL and MRL

Statistical design is applied to a multivariate exponentially weighted moving average (MEWMA) control chart. The chart parameters are control limit H and smoothing constant r. The choices of the parameters depend on the number of variables p and the size of the process mean shift δ. The MEWMA statistic is modeled as a Markov chain and the Markov chain approach is used to determine the properties of the chart. Although average run length has become a traditional measure of the performance of control schemes, some authors have suggested other measures, such as median and other percentiles of the run length distribution to explain run length properties of a control scheme. This will allow a thorough study of the performance of the control scheme. Consequently, conclusions based on these measures would provide a better and comprehensive understanding of a scheme. In this article, we present the performance of the MEWMA control chart as measured by the average run length and median run length. Graphs are given so that the chart parameters of an optimal MEWMA chart can be determined easily.     

A Multivariate Group Runs Control Chart for Process Dispersion

In this article, control charts for bivariate as well as for multivariate normal data are proposed to detect a shift in the process variability. Methods of obtaining design parameters and procedures of implementing the proposed charts are discussed. Performance of the proposed charts is compared with some existing control charts. It is verified that the proposed charts significantly reduce the out of control “average run length” (ARL) as compared to other charts considered in the study. Also, when the process variability decreases (process improvement), it is verified that the ARL of the proposed multivariate control chart increases as compared to other charts considered in the study.     

Monitoring simple linear profiles using variable sample size schemes

ABSTRACT

Profile monitoring is one of the new research areas in statistical process control. Most of the control charts in this area are designed with fixed sampling rate which makes the control chart slow in detecting small to moderate shifts. In order to improve the performance of the conventional fixed control charts, adaptive features are proposed in which, one or more design parameters vary during the process. In this paper the variable sample size feature of EWMA₃ and MEWMA schemes are proposed for monitoring simple linear profiles. The EWMA₃ method is based on the combination of three exponentially weighted moving average (EWMA) charts for monitoring three parameters of a simple linear profile separately and the Multivariate EWMA (MEWMA) chart is based on the using a single chart to monitor the coefficients and variance of a general linear profile. Also a two-sided control chart is proposed for monitoring the standard deviation in the EWMA₃ method. The performance of the proposed charts is compared in terms of the average time to signal. Numerical examples show that using adaptive features increase the power of control charts in detecting the parameter shifts. Finally, the performance of the proposed variable sample size schemes is illustrated through a real case in the leather industry.     

Control charts for attributes with maxima nominated samples

We develop quality control charts for attributes using the maxima nomination sampling (MNS) method and compare them with the usual control charts based on simple random sampling (SRS) method, using average run length (ARL) performance, the required sample size in detecting quality improvement, and non-existence region for control limits. We study the effect of the sample size, the set size, and nonconformity proportion on the performance of MNS control charts using ARL curve. We show that MNS control chart can be used as a better benchmark for indicating quality improvement or quality deterioration relative to its SRS counterpart. We consider MNS charts from a cost perspective. We also develop MNS attribute control charts using randomized tests. A computer program is designed to determine the optimal control limits for an MNS p-chart such that, assuming known parameter values, the absolute deviation between the ARL and a specific nominal value is minimized. We provide good approximations for the optimal MNS control limits using regression analysis. Theoretical results are augmented with numerical evaluations. These show that MNS based control charts can yield substantial improvement over the usual control charts based on SRS.     

CUSUM control schemes for Gaussian processes

A CUSUM control scheme for detecting a change point in a Gaussian process is derived. An upper and a lower bound for the distribution of the run length and for its moments is given. In a Monte Carlo study the average run length (ARL) of this chart is compared with the ARL of two other CUSUM procedures which are based on approximations to the sequential probability ratio, and, moreover, with EWMA schemes for autocorrelated data. Results on the optimal choice of the reference value are presented. Furthermore it is investigated how these charts behave if the model parameters are estimated.     

Accounting for phase I sampling variability in the performance of the MEWMA control chart with estimated parameters

In this article, we assess the performance of the multivariate exponentially weighted moving average (MEWMA) control chart with estimated parameters while considering the practitioner-to-practitioner variability. We evaluate the chart performance in terms of the in-control average run length (ARL) distributional properties; mainly the average (AARL), the standard deviation (SDARL), and some percentiles. We show through simulations that using estimates in place of the in-control parameters may result in an in-control ARL distribution that almost completely lies below the desired value. We also show that even with the use of larger amounts of historical data, there is still a problem with the excessive false alarm rates. We recommend the use of a recently proposed bootstrap-based design technique for adjusting the control limits. The technique is quite effective in controlling the percentage of short in-control ARLs resulting from the estimation error.     

Multivariate EWMA Control Chart with Adaptive Sample Sizes

Standard multivariate control charts usually employ fixed sample sizes at equal sampling intervals to monitor a process. In this study, a multivariate exponential weighted moving average (MEWMA) chart with adaptive sample sizes is investigated. Performance measure of the adaptive-sample-size MEWMA chart is obtained through a Markov chain approach. The performance of the adaptive-sample-size MEWMA chart is compared with the fixed-sample-size control chart in terms of steady-state average run length for different magnitude of shifts in the process mean. It is shown that the adaptive-sample-size chart is more efficient than the fixed-sample-size MEWMA control chart in detecting shifts in the process mean.     

Optimal design of the adaptive exponentially weighted moving average control chart over a range of mean shifts

The adaptive exponentially weighted moving average (AEWMA) control chart is a smooth combination of the Shewhart and exponentially weighted moving average (EWMA) control charts. This chart was proposed by Cappizzi and Masarotto (2003) to achieve a reasonable performance for both small and large shifts. Cappizzi and Masarotto (2003) used a pair of shifts in designing their control chart. In this study, however, the process mean shift is considered as a random variable with a certain probability distribution and the AEWMA control chart is optimized for a wide range of mean shifts according to that probability distribution and not just for a pair of shifts. Using the Markov chain technique, the results show that the new optimization design can improve the performance of the AEWMA control chart from an overall point of view relative to the various designs presented by Cappizzi and Masarotto (2003). Optimal design parameters that achieve the desired in-control average run length (ARL) are computed in several cases and formulas used to find approximately their values are given. Using these formulas, the practitioner can compute the optimal design parameters corresponding to any desired in-control ARL without the need to apply the optimization procedure. The results obtained by these formulas are very promising and would particularly facilitate the design of the AEWMA control chart for any in-control ARL value.     

An ARL-unbiased design of time-between-events control charts with runs rules

In this paper, we consider incorporating the runs rules into the cumulative quantity control (CQC) chart for monitoring time-between-events data. We propose a simple and effective procedure to design a CQC chart coupled with runs rules that can yield average run length (ARL)-unbiased performance and meet the required in-control ARL. The proposed design involves determining a relation between the upper side and lower side false alarm probabilities. A Markov chain approach is used to evaluate the ARL performance of various control schemes studied in this paper. An extensive numerical comparison shows that the proposed design approach can result in a significant reduction in ARL for detecting increases in the occurrence rate of the event in comparison with the basic CQC charts.     

Synthetic exponential control charts with unknown parameter

The existing synthetic exponential control charts are based on the assumption of known in-control parameter. However, the in-control parameter has to be estimated from a Phase I dataset. In this article, we use the exact probability distribution, especially the percentiles, mean, and standard deviation of the conditional average run length (ARL) to evaluate the effect of parameter estimation on the performance of the Phase II synthetic exponential charts. This approach accounts for the variability in the conditional ARL values of the synthetic chart obtained by different practitioners. Since parameter estimation results in more false alarms than expected, we develop an exact method to design the adjusted synthetic charts with desired conditional in-control performance. Results of known and unknown in-control parameter cases show that the control limit of the conforming run length sub-chart of the synthetic chart should be as small as possible.     

A weighted likelihood ratio test-based chart for monitoring process mean and variability

In this paper, a new single exponentially weighted moving average (EWMA) control chart based on the weighted likelihood ratio test, referred to as the WLRT chart, is proposed for the problem of monitoring the mean and variance of a normally distributed process variable. It is easy to design, fast to compute, and quite effective for diverse cases including the detection of the decrease in variability and individual observation case. The optimal parameters that can be used as a design aid in selecting specific parameter values based on the average run length (ARL) and the sample size are provided. The in-control (IC) and out-of-control (OC) performance properties of the new chart are compared with some other existing EWMA-type charts. Our simulation results show that the IC run length distribution of the proposed chart is similar to that of a geometric distribution, and it provides quite a robust and satisfactory overall performance for detecting a wide range of shifts in the process mean and/or variability.     

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.     

The computation of average run length and average time to signal: an overview

Control charts are widely used in industries to monitor a process for quality improvement. Evaluation of the average run length (ARL) or average time to signal (ATS) plays an important role in the design of control charts and performance comparison. In this paper, we review several basic and popular procedures, including the Markov chain and integral equation methods for computing ARL, ATS and associated run length distributions for cumulative sum charts, exponentially weighted moving average charts and combined control charts, respectively. Some important references and key formulations are provided for practitioners.     

Some issues in the design of EWMA Charts

The traditional design procedure for selecting the parameters of EWMA charts is based on the average run length (ARL). It is shown that for some types of EWMA charts, such a procedure may lead to high probability of a false out-of-control signal. An alternative procedure based on both the ARL and the standard deviation of run length (SRL) is recommended. It is shown that, with the new procedure, the EWMA chart using its exact variance can detect moderate and large shifts of the process mean faster.     

An economic statistical design of the Poisson EWMA control charts for monitoring nonconformities

An economic statistical model of the exponentially weighted moving average (EWMA) control chart for the average number of nonconformities in the sample is proposed. The statistical and economic performance of proposed design are evaluated using the average run length (ARL) and the hourly expected cost, respectively. A Markov chain approach is applied to derive expressions for ARL. The cost model is established based on the general cost function given in Lorenzen and Vance [The economic design of control charts: a unified approach. Technometrics. 1986;28:3–11]. An example is provided to illustrate the application of the proposed model. A sensitivity analysis is also carried out to investigate the effects of model parameters on the solution of the economic statistical design by using the design of experiments (DOE) technique.     

A a comparison of the markov chain and the integral equation approaches for evaluating the run length distribution of quality control charts

Two methods that are often used to evaluate the run length distribution of quality control charts are the Markov chain and integral equation approaches. Both methods have been used to evaluate the cumulative sum (CUSUM) charts and the exponentially weighted moving average (EWMA) control charts. The Markov chain approach involves "discretiz-ing" the possible values which can be plotted. Using properties of finite Markov chains, expressions for the distribution of the run length, and for the average run length (ARL), can be obtained. For the CUSUM and EWMA charts there exist integral equations whose solution gives the ARL. Approximate methods can then be used to solve the integral equation. In this article we show that if the product midpoint rule is used to approximate the integral in the integral equation, then both approaches yield the same approximations for the ARL. In addition we show that the recursive expressions for the probability functions are the same for the two approaches. These results establish the integral equation approach as preferable whenever an integral equation can be found