Run length,average run length and false alarm rate of shewhart x-bar chart: exact derivations by conditioning

The effects of estimation of the control limits on the performance of the popular Shewhart X-bar chart are examined via the average run length and the probability of a false alarm, when one or both of the process mean and variance are unknown. Exact expressions for the run length, the average run length (ARL) and the false alarm rate are obtained, in each case, using expectation by conditioning. Applying Jensen's inequality, together with expectation by conditioning, a simple lower bound to the ARL is obtained. This could be useful in designing the charts. The expressions for the exact ARL and the exact probabilities of false alarm are evaluated, using simulations, for various numbers of subgroups and shift sizes. The calculations throw new light on the performance of the Shewhart X-bar chart. Some recommendations are given.     

A class of distribution-free control charts

Summary.  A class of Shewhart-type distribution-free control charts is considered. A key advantage of these charts is that the in-control run length distribution is the same for all continuous process distributions. Exact expressions for the run length distribution and the average run length (ARL) are derived and properties of the charts are studied via evaluations of the run length distribution probabilities and the ARL. Tables are provided for implementation for some typical ARL values and false alarm rates. The charts proposed are preferable from a robustness point of view, have attractive ARL properties and would be particularly useful in situations where one uses a classical Shewhart   X -chart. A numerical illustration is given.     

ARL-Design of S Charts with k-of-k Runs Rules

The standard S chart signals an out-of-control condition when one point exceeds a control limit. It can be augmented with runs rules to improve its performance in detecting assignable causes. A commonly used rule signals when k consecutive points exceed a control limit. This rule can be used alone or to supplement the standard chart. In this article we derive ARL expressions for charts with the k-of-k runs rule. We show how to design S charts with this runs rule, compare their ARL performance, and make a control chart recommendation when it is important to monitor for both increases and decreases in process dispersion.     

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     

Nonlinear recursive estimation of volatility via estimating functions

For certain volatility models, the conditional moments that depend on the parameter are of interest. Following Godambe and Heyde (1987), the combined estimating function method has been used to study inference when the conditional mean and conditional variance are functions of the parameter of interest (See Ghahramani and Thavaneswaran [Combining Estimating Functions for Volatility. Journal of Statistical Planning and Inference, 2009, 139, 1449-1461] for details). However, for application purposes, the resulting estimates are nonlinear functions of the observations and no closed form expressions of the estimates are available. As an alternative, in this paper, a recursive estimation approach based on the combined estimating function is proposed and applied to various classes of time series models, including certain volatility models.     

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.     

Shewhart x-charts with estimated process variance

Properties of the Shewhart X-chart for controlling the mean of a process with a normal distribution are investigated for the situation where the process variance Ó²must be estimated from initial sample data. The control limits of the X-chart depend on the estimate of Ó²and thus, unlike the case when Ó²is known, the X-chart is not equivalent to a sequence of independent tests. When Ó²is estimated the distribution of the run length is not geometric and cannot be characterized simply in terms of the probability of a signal at a given point. The average run length (ARL) for the X-chart is expressed in terms of an integral involving the normal cdf, and it is shown that the chart signals with

probability one, but the ARL may not be finite if the size of the 2 sample used to estimate Ó²is sufficiently small. In addition, certain bounds for the ARL are also derived. Numerical integration is use to show that the effect of using small sample sizes in estimating Ó²is to increase the ARL and the variance of the run length distribution     

Zone control charts with estimated parameters for detecting prespecified changes in linear profiles

ABSTRACT

In profile monitoring, control charts are proposed to detect unanticipated changes, and it is usually assumed that the in-control parameters are known. However, due to the characteristics of a system or process, the prespecified changes would appear in the process. Moreover, in most applications, the in-control parameters are usually unknown. To overcome these issues, we develop the zone control charts with estimated parameters to detect small shifts of these prespecified changes. The effects of estimation error have been investigated on the performance of the proposed charts. To account for the practitioner-to-practitioner variability, the expected average run length (ARL) and the standard deviation of the average run length (SDARL) is used as the performance metrics. Our results show that the estimation error results in the significant variation in the ARL distribution. Furthermore, in order to adequately reduce the variability, more phase I samples are required in terms of the SDARL metric than that in terms of the expected ARL metric. In addition, more observations on each sampled profile are suggested to improve the charts' performance, especially for small phase I sample sizes. Finally, an illustrative example is given to show the performance of the proposed zone control charts.     

Transition probabilities and ARL performance of Six Sigma zone control charts

In this article, Six Sigma zone control charts (SSZCCs) are proposed for world class organizations. The transition probabilities are obtained using the Markov chain approach. The Average Run Length (ARL) values are then presented. The ARL performance of the proposed SSZCCs and the standard Six Sigma control chart (SSCC) without zones or run rules is studied. The ARL performance of these charts is then compared with those of the other standard zone control charts (ZCCs), the modified ZCC and the traditional Shewhart control chart (SCC) with common run rules. As expected, it is shown that the proposed SSZCC outperforms the standard SSCC without zones or run rules for process shifts of any magnitude. When compared to the other standard ZCCs and the Shewhart chart with common run rules, it is observed that the proposed SSZCCs have much higher false alarm rates for smaller shifts and hence they prevent unwanted process disturbances. The application of the proposed SSZCC is illustrated using a real time example.     

Modified inference about the mean of the exponential distribution using moving extreme ranked set sampling

The maximum likelihood estimator (MLE) and the likelihood ratio test (LRT) will be considered for making inference about the scale parameter of the exponential distribution in case of moving extreme ranked set sampling (MERSS). The MLE and LRT can not be written in closed form. Therefore, a modification of the MLE using the technique suggested by Maharota and Nanda (Biometrika 61:601–606, 1974) will be considered and this modified estimator will be used to modify the LRT to get a test in closed form for testing a simple hypothesis against one sided alternatives. The same idea will be used to modify the most powerful test (MPT) for testing a simple hypothesis versus a simple hypothesis to get a test in closed form for testing a simple hypothesis against one sided alternatives. Then it appears that the modified estimator is a good competitor of the MLE and the modified tests are good competitors of the LRT using MERSS and simple random sampling (SRS).     

On a problem in two-stage sampling

Beattie 1962 proposed an alternate procedure for monitoring continuous processes which, unlike existing continuous samplingplans under Mil-Std-1235B, does not require periods of 100% inspection. The procedure, which is a composite of two Cumulative (Cusum) procedures, will be analyzed in terms of Cusum average run length (ARL) performance. An algorithm is proposed for generating sample plans based upon specifications on ARL performance in each of the Cusum sampling zones and the overall sampling risks of the procedure. The algorithm is fully described, and implemented using exact expressions for the Cusum ARL. An illustrated example is provided, describing how a continuous sampling maybe selected which matches performances about the AQL for a related plan in ISO 3951.     

Sample size determination for 2k-r experiments with a binomial response

This paper provides closed form expressions for the sample size for two-level factorial experiments when the response is the number of defectives. The sample sizes are obtained by approximating the two-sided test for no effect through tests for the mean of a normal distribution, and borrowing the classical sample size solution for that problem. The proposals are appraised relative to the exact sample sizes computed numerically, without appealing to any approximation to the binomial distribution, and the use of the sample size tables provided is illustrated through an example.     

Linear estimation for the extended exponential power distribution

Tiao and Lund [The use of OLUMV estimators in inference robustness studies of the location parameter of a class of symmetric distributions. J Amer Statist Assoc. 1970;65(329):370–386] tabulated the coefficients of the best linear unbiased estimators (BLUEs) of location and scale for a particular family of symmetric distributions. This family was a reparameterization of the extended exponential power distribution (EEPD) with the shape parameter restricted to be greater than or equal to one. In this work, we consider the BLU estimation of the location and scale parameters of the EEPD when the shape parameter is one-third and one-half. We obtain closed-form expressions for the single and product moments of the order statistics when the shape parameter is in general in the form of a reciprocal of an integer. These expressions are then used to determine the BLUEs and the corresponding variances for complete samples of size 20 and less. We consider some other linear estimators of the location and scale parameters and then compare them with the BLUEs. Finally, we present a numerical example to illustrate the developed results.     

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.     

Alternative tests for heteroscedasticity of disturbances:a comparative study

This paper presents three small sample tests for testing the heteroscedasticity among regression disturbances. The power of these tests are compared with two of the leading tests for this hypothesis, one by Goldfeld and Quandt [5] and the other by Theil [17]. We also provide a heuristic method of selecting the number of middle observations to be deleted for the Goldfeld-Quandt type of tests.     

Implementing optimal allocation for sequential continuous responses with multiple treatments

In practice, it is important to find optimal allocation strategies for continuous response with multiple treatments under some optimization criteria. In this article, we focus on exponential responses. For a multivariate test of homogeneity, we obtain the optimal allocation strategies to maximize power while (1) fixing sample size and (2) fixing expected total responses. Then the doubly adaptive biased coin design [Hu, F., Zhang, L.-X., 2004. Asymptotic properties of doubly adaptive biased coin designs for multi-treatment clinical trials. The Annals of Statistics 21, 268–301] is used to implement the optimal allocation strategies. Simulation results show that the proposed procedures have advantages over complete randomization with respect to both inferential (power) and ethical standpoints on average. It is important to note that one can usually implement optimal allocation strategies numerically for other continuous responses, though it is usually not easy to get the closed form of the optimal allocation theoretically.     

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.     

On the Characteristic Function for Asymmetric Exponential Power Distributions

The econometric literature has seen a surge of developments in the theory and applications of asymmetric exponential power distributions (AEPDs). Here, for the first time, we derive explicit closed form expressions for the characteristic function of AEPDs. The expressions involve the complex parameter Wright generalized hypergeometric function.     

Tests for randomness of observations that involve two factors

Tests for randomness of observations that involve one factor have been considered by many authors among them are Mosteller [4], Bateman [21], Barton and David [1] and Shaughnessy [7]. However, on many occasions, data involve two different factors such as time and location, temperature and pressure, or levels of doses and responses of patients and so on. In this paper, we consider tests for randomness of observations that involve two factors for which data are given in a matrix form. Some new definitions of runs of a matrix of data are defined and discussed. A special kind of run is proposed for the test for randomness. Distributions and properties of this run are studied. Some critical regions are tabulated.     

Likelihood ratio tests for variance components in linear mixed models

Although the asymptotic distributions of the likelihood ratio for testing hypotheses of null variance components in linear mixed models derived by Stram and Lee [1994. Variance components testing in longitudinal mixed effects model. Biometrics 50, 1171–1177] are valid, their proof is based on the work of Self and Liang [1987. Asymptotic properties of maximum likelihood estimators and likelihood tests under nonstandard conditions. J. Amer. Statist. Assoc. 82, 605–610] which requires identically distributed random variables, an assumption not always valid in longitudinal data problems. We use the less restrictive results of Vu and Zhou [1997. Generalization of likelihood ratio tests under nonstandard conditions. Ann. Statist. 25, 897–916] to prove that the proposed mixture of chi-squared distributions is the actual asymptotic distribution of such likelihood ratios used as test statistics for null variance components in models with one or two random effects. We also consider a limited simulation study to evaluate the appropriateness of the asymptotic distribution of such likelihood ratios in moderately sized samples.