The performance of multivariate CUSUM control charts with estimated parameters

In this paper, we study the effect of estimating the vector of means and the variance–covariance matrix on the performance of two of the most widely used multivariate cumulative sum (CUSUM) control charts, the MCUSUM chart proposed by Crosier [Multivariate generalizations of cumulative sum quality-control schemes, Technometrics 30 (1988), pp. 291–303] and the MC1 chart proposed by Pignatiello and Runger [Comparisons of multivariate CUSUM charts, J. Qual. Technol. 22 (1990), pp. 173–186]. Using simulation, we investigate and compare the in-control and out-of-control performances of the competing charts in terms of the average run length measure. The in-control and out-of-control performances of the competing charts deteriorate significantly if the estimated parameters are used with control limits intended for known parameters, especially when only a few Phase I samples are used to estimate the parameters. We recommend the use of the MC1 chart over that of the MCUSUM chart if the parameters are estimated from a small number of Phase I samples.     

Reflections on Changes in the Library and Serials Literature

Abstract

Johnson describes changes in library journal literature topics over the past one hundred years, with a focus on the period between her first Balance Point column in 2003 and this one, her last. She reminisces on her experience as a Balance Point column editor.     

Estimated Data Collection Organization Effect in the National Medical Care Expenditure Survey

The National Medical Care Expenditure Survey, which has a complex survey design, is further complicated by combining two independently drawn national samples of households from the Research Triangle Institute (RTI) and the National Opinion Research Center (NORC). It is assumed that since the structures of these national-area samples are similar, they are thereby compatible and allow for the derivation of statistically equivalent unbiased national estimates of relevant health parameters. In this paper, national parameter estimates of relevant demographic and health measures were separately produced for each of the two survey organizations and a determination was made regarding statistical equivalence across samples. In addition, the data-collection organization effect was estimated and tested within the framework of a multivariate analysis appropriate for complex survey data. The findings consistently revealed a statistically nonsignificant data-collection organization effect when testing for differences in domain estimates of the relevant health parameters under consideration.     

Weighted Searching Probability for Classes of Equivalent Search Designs Comparison

Search design is searching and estimating for a few non zero effects in a large set of effects along with estimation of elements in a set of unknown parameters. In presence of noise, the probability of discrimination between the true non zero effect from an alternative one depends on the design and an unknown parameter, say ρ. We develop a new criterion for design comparison which is independent of ρ and for a family density weight function show that it discriminates and ranks the designs precisely. This criterion is invariance to the variable noise which may be present between designs due to noise factors. This allows us to extend the design comparison to classes of equivalent designs.     

Akaike Information Criterion for Selecting Variables in the Nested Error Regression Model

The Akaike Information Criterion (AIC) is developed for selecting the variables of the nested error regression model where an unobservable random effect is present. Using the idea of decomposing the likelihood into two parts of “within” and “between” analysis of variance, we derive the AIC when the number of groups is large and the ratio of the variances of the random effects and the random errors is an unknown parameter. The proposed AIC is compared, using simulation, with Mallows' C _p , Akaike's AIC, and Sugiura's exact AIC. Based on the rates of selecting the true model, it is shown that the proposed AIC performs better.     

Capturing the Spillover Effect With Multiplicative Error Models

Recent statistical models for the analysis of volatility in financial markets serve the purpose of incorporating the effect of other markets in their structure, in order to study the spillover or the contagion phenomena. Extending the Multiplicative Error Model we are able to capture these characteristics, under the assumption that the conditional mean of the volatility can be decomposed into the sum of one component representing the proper volatility of the time series analyzed, and other components, each representing the volatility transmitted from one other market. Each component follows a proper dynamics with elements that can be usefully interpreted. This particular decomposition allows to establish, each time, the contribution brought by each individual market to the global volatility of the market object of the analysis. We experiment this model with four stock indices.     

Analyzing Binomial Data in a Split-Plot Design: Classical Approach or Modern Techniques?

Modeling data that are non-normally distributed with random effects is the major challenge in analyzing binomial data in split-plot designs. Seven methods for analyzing such data using mixed, generalized linear, or generalized linear mixed models are compared for the size and power of the tests. This study shows that analyzing random effects properly is more important than adjusting the analysis for non-normality. Methods based on mixed and generalized linear mixed models hold Type I error rates better than generalized linear models. Mixed model methods tend to have higher power than generalized linear mixed models when the sample size is small.     

The Effect of Methods for Handling Missing Values on the Performance of the MEWMA Control Chart

This article is concerned with the effect of the methods for handling missing values in multivariate control charts. We discuss the complete case, mean substitution, regression, stochastic regression, and the expectation–maximization algorithm methods for handling missing values. Estimates of mean vector and variance–covariance matrix from the treated data set are used to build the multivariate exponentially weighted moving average (MEWMA) control chart. Based on a Monte Carlo simulation study, the performance of each of the five methods is investigated in terms of its ability to obtain the nominal in-control and out-of-control average run length (ARL). We consider three sample sizes, five levels of the percentage of missing values, and three types of variable numbers. Our simulation results show that imputation methods produce better performance than case deletion methods. The regression-based imputation methods have the best overall performance among all the competing methods.     

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.     

Fractional replication in simulation studies

It is shown how to use fractional replication in simulation studies. Examples are given. Considerable savings in number of runs required can be achieved through the use of fractional replication ideas.