Detection of Multiple Outbreaks Using Spatio-Temporal EWMA-Ordered Statistics

Spatio-temporal surveillance methods for detecting outbreaks of disease are fairly common in the literature with the scan statistic setting the benchmark. If the shape and size of the outbreaks are known in advance, then the scan approach can be designed to efficiently detect these, however, this is seldom true. Therefore we want to devise plans that are efficient at detecting a number of outbreaks that vary in size and shape. This paper examines plans which use the exponential weighted moving average statistic to build temporal memory into plans and tries to develop robust plans for detecting outbreaks of unknown shapes and sizes.     

Optimal exponentially weighted moving average (EWMA) plans for detecting seasonal epidemics when faced with non-homogeneous negative binomial counts

Exponentially weighted moving average (EWMA) plans for non-homogeneous negative binomial counts are developed for detecting the onset of seasonal disease outbreaks in public health surveillance. These plans are robust to changes in the in-control mean and over-dispersion parameter of the negative binomial distribution, and therefore are referred to as adaptive plans. They differ from the traditional approach of using standardized forecast errors based on the normality assumption. Plans are investigated in terms of early signal properties for seasonal epidemics. The paper demonstrates that the proposed EWMA plan has efficient early detection properties that can be useful to epidemiologists for communicable and other disease control and is compared with the CUSUM plan.     

A Modification on the Friedman Test Statistic

A modification on the well-known, nonparametric Friedman test statistic is suggested in this article. Probability distributions of the suggested test statistic under the null hypothesis are tabulated for some small sample cases. In addition to an example, simulation results for various sample sizes are presented. The simulation indicates that the modified test statistic performs better than the Friedman test in detecting treatment effects of small differences especially when the sample size is small.     

Recurrence relations for moments of different orders of extremes from gamma distribution

Certain recurrence relations for the moments of different orders of the largest order statistic from a gamma distribution with shape parameter p are obtained. By using this it is shown that for obtaining the moment of any order of each order statistic of a sample of size n from the gamma distribution, one has to evaluate at most n-2 single integrals.     

Time truncated acceptance sampling plans for generalized exponential distribution

Acceptance sampling plans for generalized exponential distribution when the lifetime experiment is truncated at a pre-determined time are provided in this article. The tables are provided for the minimum sample size required to ensure a certain median life of the experimental unit when the shape parameter is two. The operating characteristic function values of the sampling plans and the associated producer's risks are also presented. It is shown that the tables presented here can be used if instead of median life, other percentile life is chosen as the criterion or if the shape parameter is not two. Examples are provided for illustrative purposes.     

Improved attribute acceptance sampling plans based on maxima nomination sampling

This paper demonstrates the use of maxima nomination sampling (MNS) technique in design and evaluation of single AQL, LTPD, and EQL acceptance sampling plans for attributes. We exploit the effect of sample size and acceptance number on the performance of our proposed MNS plans using operating characteristic (OC) curve. Among other results, we show that MNS acceptance sampling plans with smaller sample size and bigger acceptance number perform better than commonly used acceptance sampling plans for attributes based on simple random sampling (SRS) technique. Indeed, MNS acceptance sampling plans result in OC curves which, compared to their SRS counterparts, are much closer to the ideal OC curve. A computer program is designed which can be used to specify the optimum MNS acceptance sampling plan and to show, visually, how the shape of the OC curve changes when parameters of the acceptance sampling plan vary. Theoretical results and numerical evaluations are given.     

Two‐Sample Bootstrap Hypothesis Tests for Three‐Dimensional Labelled Landmark Data

Abstract. We investigate resampling methodologies for testing the null hypothesis that two samples of labelled landmark data in three dimensions come from populations with a common mean reflection shape or mean reflection size‐and‐shape. The investigation includes comparisons between (i) two different test statistics that are functions of the projection onto tangent space of the data, namely the James statistic and an empirical likelihood statistic; (ii) bootstrap and permutation procedures; and (iii) three methods for resampling under the null hypothesis, namely translating in tangent space, resampling using weights determined by empirical likelihood and using a novel method to transform the original sample entirely within refection shape space. We present results of extensive numerical simulations, on which basis we recommend a bootstrap test procedure that we expect will work well in practise. We demonstrate the procedure using a data set of human faces, to test whether humans in different age groups have a common mean face shape.     

Nonparametric quality control charts based on the sign statistic

Nonparametric control chart are presented for the problem of detecting changes in the process median (or mean), or changes in the process variability when samples are taken at regular time intervals. The proposed procedures are based on sign-test statistics computed for each sample, and are used in Shewhart and cumulative sum control charts. When the process is in control the run length distributions for the proposed nonparametric control charts do not depend on the distribution of the observations. An additional advantage of the non-parametric control charts is that the variance of the process does not need to be established in order to set up a control chart for the mean. Comparisons with the corresponding parametric control charts are presented. It is also shown that curtailed sampling plans can considerably reduce the expected number of observations used in the Shewhart control schemes based on the sign statistic.     

An economic reliability test plan: Log-logistic distribution

Sampling plans in which items that are put to test, to collect the life of the items in order to decide upon accepting or rejecting a submitted lot, are called reliability test plans. The basic probability model of the life of the product is specified as the well-known log-logistic distribution with a known shape parameter. For a given producer's risk, sample size, termination number, and waiting time to terminate the test plan are computed. The preferability of the test plan over similar plans existing in the literature is established with respect to cost and time of the experiment.     

Modified stationarity tests with improved power in small samples

In a recent paper Kwiatkowski et al. (1992) propose the so-called KPSS statistic for testing the null hypothesis of stationarity against the alternative of a unit root process. The statistic employs a spectral estimator which can be shown to diverge with increasing sample size, given the alternative is true. Here, we suggest a modified spectral estimator which is shown to stabilize for moving average models. It is shown that this test statistic uniformly outperforms the KPSS statistic in an MA(1) model. Furthermore, a two-step nonparametric correction procedure is suggested, giving a test statistic with similar asymptotic properties as the original KPSS statistic. However, in small samples this correction performs better especially in detecting large random walk components. This paper was written while the author was a post-doctoral fellow at the University of Amsterdam. The author likes to thank Peter Boswijk, Inge van den Doel, Noud van Giersbergen and Jan F.Kiviet for their help during that time. Moreover, I would like to thank an anonymous referee for a number of helpful comments.     

The limiting distribution of a test for multivariate structure

We define a chi-squared statistic for p-dimensional data as follows. First, we transform the data to remove the correlations between the p variables. Then, we discretize each variable into groups of equal size and compute the cell counts in the resulting p-way contingency table. Our statistic is just the usual chi-squared statistic for testing independence in a contingency table. Because the cells have been chosen in a data-dependent manner, this statistic does not have the usual limiting distribution. We derive the limiting joint distribution of the cell counts and the limiting distribution of the chi-squared statistic when the data is sampled from a multivariate normal distribution. The chi-squared statistic is useful in detecting hidden structure in raw data or residuals. It can also be used as a test for multivariate normality.     

An unconditional space–time scan statistic for ZIP‐distributed data

A scan statistic is proposed for the prospective monitoring of spatiotemporal count data with an excess of zeros. The method that is based on an outbreak model for the zero‐inflated Poisson distribution is shown to be superior to traditional scan statistics based on the Poisson distribution in the presence of structural zeros. The spatial accuracy and the detection timeliness of the proposed scan statistic are investigated by means of simulation, and an application on the weekly cases of Campylobacteriosis in Germany illustrates how the scan statistic could be used to detect emerging disease outbreaks. An implementation of the method is provided in the open‐source R package scanstatistics available on the Comprehensive R Archive Network.     

A two-sample distribution-free scale test of the smirnov type

The two-sample, distribution-free statistics of Smirnov (1939) are used to define a new statistic. While the Smirnov statistics are used as a general goodness-of-fit test, a distribution-free scale test based on this new statistic is developed. It is shown that this new test has higher power than the two-sided Smirnov statistic in detecting differences in scale for some symmetric distributions with equal means/medians. The critical values of the proposed test statistic and its limiting distribution are given     

Multi - sample test of equal gamma distribution scale parameters in presence of unknown common shape parameter

Shiue and Bain proposed an approximate F statistic for testing equality of two gamma distribution scale parameters in presence of a common and unknown shape parameter. By generalizing Shiue and Bain's statistic we develop a new statistic for testing equality of L >= 2 gamma distribution scale parameters. We derive the distribution of the new statistic ESP for L = 2 and equal sample size situation. For other situations distribution of ESP is not known and test based on the ESP statistic has to be performed by using simulated critical values. We also derive a C(α) statistic CML and develop a likelihood ratio statistic, LR, two modified likelihood ratio statistics M and MLB and a quadratic statistic Q. The distribution of each of the statistics CML, LR, M, MLB and Q is asymptotically chi-square with L - 1 degrees of freedom. We then conducted a monte-carlo simulation study to compare the perfor- mance of the statistics ESP, LR, M, MLB, CML and Q in terms of size and power. The statistics LR, M, MLB and Q are in general liberal and do not show power advantage over other statistics. The statistic CML, based on its asymptotic chi-square distribution, in general, holds nominal level well. It is most powerful or nearly most powerful in most situations and is simple to use. Hence, we recommend the statistic CML for use in general. For better power the statistic ESP, based on its empirical distribution, is recommended for the special situation for which there is evidence in the data that λ₁ < … < λ_L and n₁ < … < n_L, where λ₁ …, λ_L are the scale parameters and n₁,…, n_L are the sample sizes.     

Application of skew-normal distribution for detecting differential expression to microRNA data

Traditional statistical modeling of continuous outcome variables relies heavily on the assumption of a normal distribution. However, in some applications, such as analysis of microRNA (miRNA) data, normality may not hold. Skewed distributions play an important role in such studies and might lead to robust results in the presence of extreme outliers. We apply a skew-normal (SN) distribution, which is indexed by three parameters (location, scale and shape), in the context of miRNA studies. We developed a test statistic for comparing means of two conditions replacing the normal assumption with SN distribution. We compared the performance of the statistic with other Wald-type statistics through simulations. Two real miRNA datasets are analyzed to illustrate the methods. Our simulation findings showed that the use of a SN distribution can result in improved identification of differentially expressed miRNAs, especially with markedly skewed data and when the two groups have different variances. It also appeared that the statistic with SN assumption performs comparably with other Wald-type statistics irrespective of the sample size or distribution. Moreover, the real dataset analyses suggest that the statistic with SN assumption can be used effectively for identification of important miRNAs. Overall, the statistic with SN distribution is useful when data are asymmetric and when the samples have different variances for the two groups.     

Multiple window discrete scan statistic for higher-order Markovian sequences

Accurate and efficient methods to detect unusual clusters of abnormal activity are needed in many fields such as medicine and business. Often the size of clusters is unknown; hence, multiple (variable) window scan statistics are used to identify clusters using a set of different potential cluster sizes. We give an efficient method to compute the exact distribution of multiple window discrete scan statistics for higher-order, multi-state Markovian sequences. We define a Markov chain to efficiently keep track of probabilities needed to compute p-values for the statistic. The state space of the Markov chain is set up by a criterion developed to identify strings that are associated with observing the specified values of the statistic. Using our algorithm, we identify cases where the available approximations do not perform well. We demonstrate our methods by detecting unusual clusters of made free throw shots by National Basketball Association players during the 2009–2010 regular season.     

Development of life-test sampling plans for exponential distributions based on accelerated life testing

This paper presents a procedure for developing life-test sampling plans for exponential distributions based upon accelerated life testing(ALT). Type II censoring is assumed at each overstress level. The derived test statistic is shown to be a quotient of two independent random variables, each of which is a rational power of a Chi-square random variable. The distribution of the test statistic is characterized by the H-function, which can be numerically evaluated to obtain desired sampling plans.     

A New Adaptive Variable Sample Size Approach in EWMA Control Chart

Adaptive control charts have been developed for improving the capability of control charts in detecting small shifts. In this article, we propose a new exponential weighted moving average control chart with variable sample size, in which the sample size is determined as an integer linear function by EWMA statistic value. The performance of the proposed VSS EWMA control chart is compared with FSS EWMA as well as traditional VSS EWMA control charts. The results show the better performance of the proposed VSS strategy respect to the traditional one and fixed sample size.     

Multivariate outbreak detection

Online monitoring is needed to detect outbreaks of diseases such as influenza. Surveillance is also needed for other kinds of outbreaks, in the sense of an increasing expected value after a constant period. Information on spatial location or other variables might be available and may be utilized. We adapted a robust method for outbreak detection to a multivariate case. The relation between the times of the onsets of the outbreaks at different locations (or some other variable) was used to determine the sufficient statistic for surveillance. The derived maximum-likelihood estimator of the outbreak regression was semi-parametric in the sense that the baseline and the slope were non-parametric while the distribution belonged to the one-parameter exponential family. The estimator was used in a generalized-likelihood ratio surveillance method. The method was evaluated with respect to robustness and efficiency in a simulation study and applied to spatial data for detection of influenza outbreaks in Sweden.     

Monitoring Process Mean and Variability with One Non-central Chi-square Chart

Traditionally, an X-chart is used to control the process mean and an R-chart to control the process variance. However, these charts are not sensitive to small changes in process parameters. A good alternative to these charts is the exponentially weighted moving average (EWMA) control chart for controlling the process mean and variability, which is very effective in detecting small process disturbances. In this paper, we propose a single chart that is based on the non-central chi-square statistic, which is more effective than the joint X and R charts in detecting assignable cause(s) that change the process mean and/or increase variability. It is also shown that the EWMA control chart based on a non-central chi-square statistic is more effective in detecting both increases and decreases in mean and/or variability.