New cumulative sum control charts for monitoring process variability

In this article, we propose new cumulative sum (CUSUM) control charts using the ordered ranked set sampling (RSS) and ordered double RSS schemes, with the perfect and imperfect rankings, for monitoring the variability of a normally distributed process. The run length characteristics of the proposed CUSUM charts are computed using the Monte Carlo simulations. The proposed CUSUM charts are compared in terms of the average and standard deviation of run lengths with their existing competitor CUSUM charts based on simple random sampling. It turns out that the proposed CUSUM charts with the perfect and imperfect rankings are more sensitive than the existing CUSUM charts based on the sample range and standard deviation. A similar trend is present when these CUSUM charts are compared with the fast initial response features. An example is also used to demonstrate the implementation and working of the proposed CUSUM charts.     

An EWMA chart for monitoring process mean

In the statistical process control literature, there exists several improved quality control charts based on cost-effective sampling schemes, including the ranked set sampling (RSS) and median RSS (MRSS). A generalized cost-effective RSS scheme has been recently introduced for efficiently estimating the population mean, namely varied L RSS (VLRSS). In this article, we propose a new exponentially weighted moving average (EWMA) control chart for monitoring the process mean using VLRSS, named the EWMA-VLRSS chart, under both perfect and imperfect rankings. The EWMA-VLRSS chart encompasses the existing EWMA charts based on RSS and MRSS (named the EWMA-RSS and EWMA-MRSS charts). We use extensive Monte Carlo simulations to compute the run length characteristics of the EWMA-VLRSS chart. The proposed chart is then compared with the existing EWMA charts. It is found that, with either perfect or imperfect rankings, the EWMA-VLRSS chart is more sensitive than the EWMA-RSS and EWMA-MRSS charts in detecting small to large shifts in the process mean. A real dataset is also used to explain the working of the EWMA-VLRSS chart.     

A new multivariate CUSUM chart using principal components with a revision of Crosier's chart

In this article, we introduce a new multivariate cumulative sum chart, where the target mean shift is assumed to be a weighted sum of principal directions of the population covariance matrix. This chart provides an attractive performance in terms of average run length (ARL) for large-dimensional data and it also compares favorably to existing multivariate charts including Crosier's benchmark chart with updated values of the upper control limit and the associated ARL function. In addition, Monte Carlo simulations are conducted to assess the accuracy of the well-known Siegmund's approximation of the average ARL function when observations are normal distributed. As a byproduct of the article, we provide updated values of upper control limits and the associated ARL function for Crosier's multivariate CUSUM chart.     

A Comparison of Spatiotemporal Surveillance Methods for Nonhomogeneous Change Size

Spatiotemporal surveillance, especially in detection of emerging outbreaks is of particular importance. When an outbreak spreads across some areas, the incidence rate at the center of the outbreak area might be expected to be much higher than the rate at its edge. However, to the best of our knowledge, all existing methods assume a uniformly increasing rate across the entire area of the outbreak. The purpose of this study is to compare the performance of the spatiotemporal surveillance methods such as multivariate cumulative sum (MCUSUM) or multivariate exponentially weighted moving average (MEWMA) when the changes in size are nonhomogeneous. Monte Carlo simulations were conducted to examine the properties of these spatiotemporal surveillance methods and compared them in terms of the detection speed and the identification rate under various scenarios. The results showed that when nonhomogeneous change sizes are involved, the MCUSUM method taking into account spatial nonhomogeneity of increase rates yields a better identification than the method ignoring such change size pattern although the detection speeds are similar. Further, a case study for the detection of male thyroid cancer data in New Mexico in the United States was performed to demonstrate the applicability of these methods.     

The probability distribution and the expected value of a stopping variable associated with one-sided cusum procedures for non-negative integer valued random variables

The structure of a stopping variable N based on one-sided CUSUM procedures is analyzed. Stopping occurs when a Markovian sequence of maxima of partial sums {M } crosses a certain boundary. On the basis of a recursive relationship between the M_n+1 and M_n a recursive equation is derived for the determination of the defective distributions K_n(x) = P{M ≤ x, N ≤n} . This recursive equation yields a recursive algorithm for the determination of P {N > n} . The paper studies the case when the basic random variables are non-negative integers-valued. In these cases the values of P{N > n} and E{N} can be determined by solving proper systems of linear equations.     

Hybrid regions for post-change mean in a sequence of normal variables

The hybrid bootstrap uses resampling ideas to extend the duality approach to the interval estimation for a parameter of interest when there are nuisance parameters. The confidence region constructed by the hybrid bootstrap may perform much better than the ordinary bootstrap region in a situation where the data provide substantial information about the nuisance parameter, but limited information about the parameter of interest. We apply this method to estimate the post-change mean after a change is detected by a stopping procedure in a sequence of independent normal variables. Since distribution theory in change point problems is generally a challenge, we use bootstrap simulation to find empirical distributions of test statistics and calculate critical thresholds. Both likelihood ratio and Bayesian test statistics are considered to set confidence regions for post-change means in the normal model. In the simulation studies, the performance of hybrid regions are compared with that of ordinary bootstrap regions in terms of the widths and coverage probabilities of confidence intervals.     

A synthetic control chart for monitoring the small shifts in a process mean based on an attribute inspection

Abstract

In this paper, a synthetic control chart is proposed by integrating the salient features of the np_x chart and the CRL chart. The synthetic chart achieves higher detection effectiveness on both small and large mean shifts while retaining the operational simplicity of the attribute charts owing to only using attribute inspection. Both statistical and economic design of the synthetic chart are considered and numerical tests have indicated that the synthetic chart has a higher power for detecting mean shifts than the np_x chart, MON chart and CUSUM chart. In addition, sensitivity analyses are also performed under both the statistical and economic design model.     

On the Run Length of a State-Space Control Chart for Multivariate Autocorrelated Data

The literature on statistical process control (SPC) describes the negative effects of autocorrelation in terms of the increase in false alarms. This has been treated by the individual modeling of each series or the application of VAR models. In the former case, the analysis of the cross correlation structure between the variables is altered. In the latter, if the cross correlation is not strong, the filtering process may modify the weakest relations. In order to improve these aspects, state-space models have been introduced in multivariate statistical process control (MSPC). This article presents a proposal for building a control chart for innovations, estimating its average run length to highlight its advantages over the VAR approach mentioned above.     

CUSUM charts for monitoring the mean of a multivariate Gaussian process

In this paper control charts for the mean of a multivariate Gaussian process are considered. Using the generalized likelihood ratio approach and the sequential probability ratio test under an additional constraint on the magnitude of the change various types of CUSUM control charts are derived. It is analyzed under which conditions these schemes are directionally invariant. These charts are compared with several other control schemes proposed in literature. The performance of the charts is studied based on the maximum average delay.     

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     

Two tests for sequential detection of a change-point in a nonlinear model

In this paper, two tests, based on weighted CUSUM of the least squares residuals, are studied to detect in real time a change-point in a nonlinear model. A first test statistic is proposed by extension of a method already used in the literature but for the linear models. It is tested under the null hypothesis, at each sequential observation, that there is no change in the model against a change presence. The asymptotic distribution of the test statistic under the null hypothesis is given and its convergence in probability to infinity is proved when a change occurs. These results will allow to build an asymptotic critical region. Next, in order to decrease the type I error probability, a bootstrapped critical value is proposed and a modified test is studied in a similar way. A generalization of the Hájek–Rényi inequality is established.     

Construction of a criterion for testing hypothesis about covariance function of a stationary Gaussian stochastic process with unknown mean

In this paper, a new criterion is constructed for testing hypothesis about covariance function of Gaussian stationary stochastic process with an unknown mean. This criterion is based on the fact that we can estimate the deviation of covariance function from its estimator with a given accuracy and reliability in L_p metric.     

Wald's approximations to the average run length in cusum procedures

Wald's approximation to the ARL(average run length in cusum) (cumulative sum) procedures are given for an exponential family of densities. From these approximations it is shown that Page's (1954) cusum procedure is (in a sense) identical with a cusum procedure defined in terms of likelihood ratios. Moreover, these approximations are improved by estimating the excess over the boundary and their closeness is examined by numerical comparisons with some exact results. Some examples are also given.     

Phase II monitoring of binary profiles in the presence of within-profile autocorrelation based on Markov Model

This paper introduces a Markov model in Phase II profile monitoring with autocorrelated binary response variable. In the proposed approach, a logistic regression model is extended to describe the within-profile autocorrelation. The likelihood function is constructed and then a particle swarm optimization algorithm (PSO) is tuned and utilized to estimate the model parameters. Furthermore, two control charts are extended in which the covariance matrix is derived based on the Fisher information matrix. Simulation studies are conducted to evaluate the detecting capability of the proposed control charts. A numerical example is also given to illustrate the application of the proposed method.     

On control charts for monitoring the variance of a time series

In this paper we derive control charts for the variance of a Gaussian process using the likelihood ratio approach, the generalized likelihood ratio approach, the sequential probability ratio method and a generalized sequential probability ratio procedure, the Shiryaev–Roberts procedure and a generalized modified Shiryaev–Roberts approach. Recursive presentations for the calculation of the control statistics are given for autoregressive processes of order 1. In an extensive simulation study these schemes are compared with existing control charts for the variance. In order to asses the performance of the schemes both the average run length and the average delay are used.     

Optimal Detection of a Change Point in a Poisson Process for Different Observation Schemes

Abstract.  Change point problems are considered where at some unobservable time the intensity of a point process ( T_n ), n ∈  N , has a jump. For a given reward functional we detect the change point optimally for different information schemes. These schemes differ in the available information. We consider three information levels, namely sequential observation of ( T_n ), ex post decision after observing the point process up to a fixed time t * and a combination of both observation schemes. In all of these cases the detection problem is viewed as an optimal stopping problem which can be solved by deriving a semimartingale representation of the gain process and applying tools from filtering theory.     

Asymptotic behaviour of a test statistic for detection of change in mean of vectors

The paper gives an asymptotic distribution of a test statistic for detecting a change in a mean of random vectors with dependent components. The studied test statistic has a form of a maximum of a square Euclidean norms of vectors with components being standardized partial cumulative sums of deviations from means. The limit distribution was obtained using a result of Piterbarg [1994. High deviations for multidimensional stationary Gaussian processes with independent components. In: Zolotarev, V.M. (Ed.), Stability Problems for Stochastic Models, pp. 197–210].     

Detection of a change-point in variance by a weighted sum of powers of variances test

This paper proposes a weighted sum of powers of variances test for detecting changes in variance of a data sequence. Asymptotic critical value formulas are derived for this test. The modified weighted sum of powers of variances test is also introduced so that the accuracy of change-point detection is highly improved for a sample of small size. Simulation studies and real data analysis are presented to assess the proposed tests.     

Identifying the time of change in the mean of a two-stage nested process

Statistical process control charts are used to distinguish between common cause and special cause sources of variability. Once a control chart signals, a search to find the special cause should be initiated. If process analysts had knowledge of the change point, the search to find the special cause could be easily facilitated. Relevant literature contains an array of solutions to the change-point problem; however, these solutions are most appropriate when the samples are assumed to be independent. Unfortunately, the assumption of independence is often violated in practice. This work considers one such case of non-independence that frequently occurs in practice as a result of multi-stage sampling. Due to its commonality in practice, we assume a two-stage nested random model as the underlying process model and derive and evaluate a maximum-likelihood estimator for the change point in the fixed-effects component of this model. The estimator is applied to electron microscopy data obtained following a genuine control chart signal and from a real machining process where the important quality characteristic is the size of the surface grains produced by the machining operation. We conduct a simulation study to compare relative performances between the proposed change-point estimator and a commonly used alternative developed under the assumption of independent observations. The results suggest that both estimators are approximately unbiased; however, the proposed estimator yields smaller variance. The implication is that the proposed estimator is more precise, and thus, the quality of the estimator is improved relative to the alternative.