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.     

EWMA control charts for monitoring correlated counts with finite range

In this work, we develop and study upper and lower one-sided EWMA control charts for monitoring correlated counts with finite range. Often in practice, data of that kind can be adequately described by a first-order binomial or beta-binomial autoregressive model. Especially, when there is evidence that data demonstrate extra-binomial variation, the latter model is preferable than the former. The proposed charts can be used for detecting upward or downward shifts in process mean level. Practical guidelines concerning the statistical design of the proposed charts are given, while the effect of the extra-binomial variation is investigated as well. Comparisons with existing control charting procedures are also provided. Finally, an illustrative real-data example is also given.     

Phase II Monitoring of Geometric Profiles

Geometric profiles can be modeled effectively by large and small scale components. In several articles, a regression model with spatial autoregressive error term is combined with control charts to monitor geometric profiles. However, once a signal occurs, control charts would not be able to determine whether the shift is due to the large or small scale component. In this article, a combination of a multivariate and an omnibus control charts is used to monitor the large scale and small scale components to determine whether the shift is due to the large scale or small scale components.     

Weighted Dickey–Fuller processes for detecting stationarity

Aiming at monitoring a time series to detect stationarity as soon as possible, we introduce monitoring procedures based on kernel-weighted sequential Dickey–Fuller (DF) processes, and related stopping times, which may be called weighted DF control charts. Under rather weak assumptions, (functional) central limit theorems are established under the unit root null hypothesis and local-to-unity alternatives. For general dependent and heterogeneous innovation sequences the limit processes depend on a nuisance parameter. In this case of practical interest, one can use estimated control limits obtained from the estimated asymptotic law. Another easy-to-use approach is to transform the DF processes to obtain limit laws which are invariant with respect to the nuisance parameter. We provide asymptotic theory for both approaches and compare their statistical behavior in finite samples by simulation.     

Monitoring Variance Change in Infinite Order Moving Average Processes and Nonstationary Autoregressive Processes

This article considers the sequential monitoring problem of variance change in stationary and non stationary time series. We suggest a CUSUM of squares procedure to detect variance change in infinite order moving average processes, and a residual CUSUM of squares procedure to detect variance change in non stationary autoregressive processes. Moreover, we introduce a bandwidth parameter to improve the monitoring power when change point does not occur at the early stage of monitoring. It is shown that both procedures have the same null distribution. The procedures are illustrated via a simulation study and an investigation of daily Mexico/US exchange rates.     

Monitoring Variability and Analyzing Multivariate Autocorrelated Processes

Traditional multivariate quality control charts are based on independent observations. In this paper, we explain how to extend univariate residual charts to multivariate cases and how to combine the traditional statistical process control (SPC) approaches to monitor changes in process variability in a dynamic environment. We propose using Alt's (1984) W chart on vector autoregressive (VAR) residuals to monitor the variability for multivariate processes in the presence of autocorrelation. We study examples jointly using the Hotelling T² chart on VAR residuals, the W chart, and the Portmanteau test to diagnose the types of shift in process parameters.     

Comparison of predictions by kriging and spatial autoregressive models

In geostatistics, the prediction of unknown quantities at given locations is commonly made by the kriging technique. In addition to the kriging technique for modeling regular lattice spatial data, the spatial autoregressive models can also be used. In this article, the spatial autoregressive model and the kriging technique are introduced. We extend prediction method proposed by Basu and Reinsel for SAR(2,1) model. Then, using a simulation study and real data, we compare prediction accuracy of the spatial autoregressive models with that of the kriging prediction. The results of simulation study show that predictions made by the autoregressive models are good competitor for the kriging method.     

A Bayesian hierarchical spatio-temporal rainfall model

A Bayesian hierarchical spatio-temporal rainfall model is presented and analysed. The model has the ability to deal with extensive missing or null values, uses a sophisticated variance stabilising rainfall pre-transformation, incorporates a new elevation model and can provide sub-catchment rainfall estimation and interpolation using a sequential kriging scheme. The model uses a vector autoregressive stochastic process to represent the time dependence of the rainfall field and an exponential covariogram to model the spatial correlation of the rainfall field. The model can be readily generalised to other types of stochastic processes. In this paper, some results of applying the model to a particular rainfall catchment are presented.     

Bootstrap-based unit root tests for higher order autoregressive models with GARCH(1, 1) errors

ABSTRACT

Bootstrap-based unit root tests are a viable alternative to asymptotic distribution-based procedures and, in some cases, are preferable because of the serious size distortions associated with the latter tests under certain situations. While several bootstrap-based unit root tests exist for autoregressive moving average processes with homoskedastic errors, only one such test is available when the innovations are conditionally heteroskedastic. The details for the exact implementation of this procedure are currently available only for the first order autoregressive processes. Monte-Carlo results are also published only for this limited case. In this paper we demonstrate how this procedure can be extended to higher order autoregressive processes through a transformed series used in augmented Dickey–Fuller unit root tests. We also investigate the finite sample properties for higher order processes through a Monte-Carlo study. Results show that the proposed tests have reasonable power and size properties.     

Control charts based on dependent count data with deflation or inflation of zeros

The process of serially dependent counts with deflation or inflation of zeros is commonly observed in many applications. This paper investigates the monitoring of such a process, the first-order zero-modified geometric integer-valued autoregressive process (ZMGINAR(1)). In particular, two control charts, the upper-sided and lower-sided CUSUM charts, are developed to detect the shifts in the mean process of the ZMGINAR(1). Both the average run length performance and the standard deviation of the run length performance of these two charts are investigated by using Markov chain approaches. Also, an extensive simulation is conducted to assess the effectiveness or performance of the charts, and the presented methods are applied to two sets of real data arising from a study on the drug use.     

New characteristics for portfolio surveillance

This paper elaborates the tools for the surveillance of the global minimum variance portfolio weights. Golosnoy and Schmid [V. Golosnoy and W. Schmid, EWMA control charts for optimal portfolio weights, Sequential Anal. 26 (2007), pp. 195–224] introduced exponentially weighted moving average-type control charts for this task based on the processes of the estimated weights as well as of their first differences. This paper proposes the new approximations to these processes exhibiting better stochastic properties for sequential monitoring purposes. The control schemes for the new processes are compared for different types of economically relevant changes using Monte Carlo simulations. The suggested procedures appear to be superior for the considered performance measures.     

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.     

Statistical Analysis Of Mixture Vector Autoregressive Models

In this paper, we reconsider the mixture vector autoregressive model, which was proposed in the literature for modelling non‐linear time series. We complete and extend the stationarity conditions, derive a matrix formula in closed form for the autocovariance function of the process and prove a result on stable vector autoregressive moving‐average representations of mixture vector autoregressive models. For these results, we apply techniques related to a Markovian representation of vector autoregressive moving‐average processes. Furthermore, we analyse maximum likelihood estimation of model parameters by using the expectation–maximization algorithm and propose a new iterative algorithm for getting the maximum likelihood estimates. Finally, we study the model selection problem and testing procedures. Several examples, simulation experiments and an empirical application based on monthly financial returns illustrate the proposed procedures.     

A Rejection Principle for Sequential Tests of Multiple Hypotheses Controlling Familywise Error Rates

We present a unifying approach to multiple testing procedures for sequential (or streaming) data by giving sufficient conditions for a sequential multiple testing procedure to control the familywise error rate (FWER). Together, we call these conditions a ‘rejection principle for sequential tests’, which we then apply to some existing sequential multiple testing procedures to give simplified understanding of their FWER control. Next, the principle is applied to derive two new sequential multiple testing procedures with provable FWER control, one for testing hypotheses in order and another for closed testing. Examples of these new procedures are given by applying them to a chromosome aberration data set and finding the maximum safe dose of a treatment.     

Monitoring the mean of autocorrelated observations with one generally weighted moving average control chart

Traditionally, using a control chart to monitor a process assumes that process observations are normally and independently distributed. In fact, for many processes, products are either connected or autocorrelated and, consequently, obtained observations are autocorrelative rather than independent. In this scenario, applying an independence assumption instead of autocorrelation for process monitoring is unsuitable. This study examines a generally weighted moving average (GWMA) with a time-varying control chart for monitoring the mean of a process based on autocorrelated observations from a first-order autoregressive process (AR(1)) with random error. Simulation is utilized to evaluate the average run length (ARL) of exponentially weighted moving average (EWMA) and GWMA control charts. Numerous comparisons of ARLs indicate that the GWMA control chart requires less time to detect various shifts at low levels of autocorrelation than those at high levels of autocorrelation. The GWMA control chart is more sensitive than the EWMA control chart for detecting small shifts in a process mean.     

On the run length of a Shewhart chart for correlated data

We consider an extension of the classical Shewhart control chart to correlated data which was introduced by Vasilopoulos/Stamboulis (1978). Inequalities for the moments of the run length are given under weak conditions. It is proved analytically that the average run length (ARL) in the in-control state of the correlated process is larger than that in the case of independent variables. The exact ARL is calculated for exchangeable normal variables and autoregressive processes (AR). Moreover, we compare this chart with residual charts. Especially, in the case of an AR(1)—process with positive coefficient, it turns out that the out-of-control ARL of the modified Shewhart chart is smaller than that of the Shewhart chart for the residuals.     

Parametric bootstrap process capability index control charts for both mean and dispersion

Abstract

Non-normal processes are common in practice. In this paper, we propose a novel approach to defining bootstrap process capability index (PCI) control charts to monitor the performance of in-control skew normal processes. We use a bootstrap method to calculate phase I control limits of the corresponding PCI control charts. The β-risk curves of the associated PCI control charts will be used to assess the performance of the PCI control charts. We use Monte-Carlo simulation to evaluate the performance of the proposed PCI control charts. A numerical example to illustrate the implementation of the proposed control charts.     

On the convergence of some adaptive estimation procedures

In the paper we consider some adaptive procedures for estimating the unknown parameters of autoregressive and moving average processes. In case of AK(p) and MA(1) processes sequences of estimators converging with probability one and In mean square are given     

Generalized Control Charts for Non-Normal Data Using g-and-k Distributions

Statistical control charts are often used in industry to monitor processes in the interests of quality improvement. Such charts assume independence and normality of the control statistic, but these assumptions are often violated in practice. To better capture the true shape of the underlying distribution of the control statistic, we utilize the g-and-k distributions to estimate probability limits, the true ARL, and the error in confidence that arises from incorrectly assuming normality. A sensitivity assessment reveals that the extent of error in confidence associated with control chart decision-making procedures increases more rapidly as the distribution becomes more skewed or as the tails of the distribution become longer than those of the normal distribution. These methods are illustrated using both a frequentist and computational Bayesian approach to estimate the g-and-k parameters in two different practical applications. The Bayesian approach is appealing because it can account for prior knowledge in the estimation procedure and yields posterior distributions of parameters of interest such as control limits.     

Self-Starting control charts for monitoring logistic regression profiles

In some applications, quality engineers cannot monitor the processes at the beginning of the production process. Because the process parameters are unknown and there are not enough initial samples to estimate the process parameters. Self-starting control charts are applied to monitor processes at the start-up stages with no enough initial samples. In this paper, we propose three self-starting control charts to monitor a logistic regression profile which models the relationship between a binomial response variable and explanatory variables. Also, we compare the proposed control charts with each other through simulation studies in terms of average run length (ARL) criterion.