Statistical Monitoring of Dynamic Multivariate Processes: Part I: Modeling Auto- and Cross-Correlation

The work in this paper is of particular significance since it considers the problem of modelling cross- and auto-correlation in statistical process monitoring. The presence of both types of correlation can lead to fault insensitivity or false alarms, although in published literature to date, only autocorrelation has been broadly considered. The proposed method, which uses a Kalman innovation model, effectively removes both correlations. The paper (and Part 2 [2]) has emerged from work supported by EPSRC grant GR/S84354/01 and is of direct relevance to problems in several application areas including chemical, electrical, and mechanical process monitoring.

Statistical Monitoring of Dynamic Multivariate Processes: Part II: Identifying Fault Magnitude and Signature

This paper builds on work presented in the first paper, Part 1 [1] and is of equal significance. The paper proposes a novel compensation method to preserve the integrity of step-fault signatures prevalent in various processes that can be masked during the removal of both auto- and cross correlation. Using industrial data, the paper demonstrates the benefit of the proposed method, which is applicable to chemical, electrical, and mechanical process monitoring. This paper, (and Part 1 [1]), has led to further work supported by EPSRC grant GR/S84354/01 involving kernel PCA methods.

Estimating capability index in multivariate processes using bootstrap sequential sampling procedures

Capability indices in both univariate and multivariate processes are extensively employed in quality control to assess the quality status of production batches before their release for operational use. It is traditionally a measure of the ratio of the allowable process spread and the actual spread. In this paper, we will adopt a bootstrap and sequential sampling procedures to determine the optimal sample size for estimating a multivariate capability index introduced by Pearns et. al. [12]. Bootstrap techniques have the distinct advantage of placing very minimum requirement on the distributions of the underlying quality characteristics, thereby rendering them more relevant under a wide variety of situations. Finally, we provide several numerical examples where the sequential sampling procedures are evaluated and compared.

A control chart based on sample ranges for monitoring the covariance matrix of the multivariate processes

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A new chart based on sample variances for monitoring the covariance matrix of multivariate processes

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Monitoring the mean vector and the covariance ­matrix of multivariate processes with sample means and sample ranges

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The use of principal components and univariate charts to control multivariate processes

In this article, we evaluate the performance of the T2 chart based on the principal components (PC chart) and the simultaneous univariate control charts based on the original variables (SU X̄ charts) or based on the principal components (SUPC charts). The main reason to consider the PC chart lies on the dimensionality reduction. However, depending on the disturbance and on the way the original variables are related, the chart is very slow in signaling, except when all variables are negatively correlated and the principal component is wisely selected. Comparing the SU X̄, the SUPC and the T 2 charts we conclude that the SU X̄ charts (SUPC charts) have a better overall performance when the variables are positively (negatively) correlated. We also develop the expression to obtain the power of two S 2 charts designed for monitoring the covariance matrix. These joint S2 charts are, in the majority of the cases, more efficient than the generalized variance |S| chart.

A single chart with supplementary runs rules for monitoring the mean vector and the covariance matrix of multivariate processes

The MRMAX chart is a single chart based on the standardized sample means and sample ranges for monitoring the mean vector and the covariance matrix of multivariate processes. User's familiarity with the computation of these statistics is a point in favor of the MRMAX chart. As a single chart, the recently proposed MRMAX chart is very appropriate for supplementary runs rules. In this article, we compare the supplemented MRMAX chart and the synthetic MRMAX chart with the standard MRMAX chart. The supplementary and the synthetic runs rules enhance the performance of the MRMAX chart. © 2013 Elsevier Ltd.

Reliability Nonparametric Bayesian Estimation in Parallel Systems

Relevant results for (sub-)distribution functions related to parallel systems are discussed. The reverse hazard rate is defined using the product integral. Consequently, the restriction of absolute continuity for the involved distributions can be relaxed. The only restriction is that the sets of discontinuity points of the parallel distributions have to be disjointed. Nonparametric Bayesian estimators of all survival (sub-)distribution functions are derived. Dual to the series systems that use minimum life times as observations, the parallel systems record the maximum life times. Dirichlet multivariate processes forming a class of prior distributions are considered for the nonparametric Bayesian estimation of the component distribution functions, and the system reliability. For illustration, two striking numerical examples are presented.

Monitoring bivariate processes

The T2 chart and the generalized variance |S| chart are the usual tools for monitoring the mean vector and the covariance matrix of multivariate processes. The main drawback of these charts is the difficulty to obtain and to interpret the values of their monitoring statistics. In this paper, we study control charts for monitoring bivariate processes that only requires the computation of sample means (the ZMAX chart) for monitoring the mean vector, sample variances (the VMAX chart) for monitoring the covariance matrix, or both sample means and sample variances (the MCMAX chart) in the case of the joint control of the mean vector and the covariance matrix.

Multivariate monitoring of a biological wastewater treatment process: a case study at Melbourne Water's Western Treatment Plant

Biological wastewater treatment is a complex, multivariate process, in which a number of physical and biological processes occur simultaneously. In this study, principal component analysis (PCA) and parallel factor analysis (PARAFAC) were used to profile and characterise Lagoon 115E, a multistage biological lagoon treatment system at Melbourne Water's Western Treatment Plant (WTP) in Melbourne, Australia. In this study, the objective was to increase our understanding of the multivariate processes taking place in the lagoon. The data used in the study span a 7-year period during which samples were collected as often as weekly from the ponds of Lagoon 115E and subjected to analysis. The resulting database, involving 19 chemical and physical variables, was studied using the multivariate data analysis methods PCA and PARAFAC. With these methods, alterations in the state of the wastewater due to intrinsic and extrinsic factors could be discerned. The methods were effective in illustrating and visually representing the complex purification stages and cyclic changes occurring along the lagoon system. The two methods proved complementary, with each having its own beneficial features. (C) 2003 Elsevier B.V. All rights reserved.

A transformation-based multivariate chart to monitor process dispersion

Multivariate monitoring techniques such as multivariate control charts are used to control the processes that contain more than one correlated characteristic. Although the majority of previous researches are focused on controlling only the mean vector of multivariate processes, little work has been performed to monitor the covariance matrix. In this research, a new method is presented to detect possible shifts in the covariance matrix of multivariate processes. The basis of the proposed method is to eliminate the correlation structure between the quality characteristics by transformation technique and then use an S chart for each variable. The performance of the proposed method is then compared to the ones from other existing methods and a real case is presented.

O efeito da autocorrelação no desempenho do gráfico T2 de hotelling: caso bivariado

The chart of control of Hotelling T2 has been the main statistical device used in monitoring multivariate processes. Currently the technological development of control systems and automation enabled a high rate of collection of information of the production systems in very short time intervals, causing a dependency between the results of observations. This phenomenon known as auto correlation causes in the statistical control of the multivariate processes a high rate of false alarms, prejudicing in the chart performance. This entails the violation of the assumption of independence and normality of the distribution. In this thesis we considered not only the correlation between two variables, but also the dependence between observations of the same variable, that is, auto correlation. It was studied by simulation, the bi variate case and the effect of auto correlation on the performance of the T2 chart of Hotelling.

Gráfico de controle MCMAX para o monitoramento simultâneo do vetor de médias e da matriz de covariâncias

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Joint control charts for monitoring location and dispersion

The T-2 and the generalized variance vertical bar S vertical bar charts are used for monitoring the mean vector and the covariance matrix of multivariate processes. In this article, we propose for bivariate processes the use of the T-2 and the VMAX charts. The points plotted on the VMAX chart correspond to the maximum of the sample variances of the two quality characteristics. The reason to consider the VMAX statistic instead of the generalized variance vertical bar S vertical bar is the user's familiarity with the computation of simple sample variances; we can't say the same with regard to the computation of the generalized variance vertical bar S vertical bar.