Monitoring Variation in a Multivariate Process When the Dimension is Large Relative to the Sample Size

A control procedure is presented for monitoring changes in variation for a multivariate normal process in a Phase II operation where the subgroup size, m, is less than p, the number of variates. The methodology is based on a form of Wilk' statistic, which can be expressed as a function of the ratio of the determinants of two separate estimates of the covariance matrix. One estimate is based on the historical data set from Phase I and the other is based on an augmented data set including new data obtained in Phase II. The proposed statistic is shown to be distributed as the product of independent beta distributions that can be approximated using either a chi-square or F-distribution. An ARL study of the statistic is presented for a range of conditions for the population covariance matrix. Cases are considered where a p-variate process is being monitored using a sample of m observations per subgroup and m < p. Data from an industrial multivariate process is used to illustrate the proposed technique.     

The Effectiveness of Multistage Ranked Set Sampling in Stratifying the Population

The spectral analysis of Gaussian linear time-series processes is usually based on uni-frequential tools because the spectral density functions of degree 2 and higher are identically zero and there is no polyspectrum in this case. In finite samples, such an approach does not allow the resolution of closely adjacent spectral lines, except by using autoregressive models of excessively high order in the method of maximum entropy. In this article, multi-frequential periodograms designed for the analysis of discrete and mixed spectra are defined and studied for their properties in finite samples. For a given vector of frequencies ω, the sum of squares of the corresponding trigonometric regression model fitted to a time series by unweighted least squares defines the multi-frequential periodogram statistic IM(ω). When ω is unknown, it follows from the properties of nonlinear models whose parameters separate (i.e., the frequencies and the cosine and sine coefficients here) that the least-squares estimator of frequencies is obtained by maximizing I ^M(ω). The first-order, second-order and distribution properties of I ^M(ω) are established theoretically in finite samples, and are compared with those of Schuster's uni-frequential periodogram statistic. In the multi-frequential periodogram analysis, the least-squares estimator of frequencies is proved to be theoretically unbiased in finite samples if the number of periodic components of the time series is correctly estimated. Here, this number is estimated at the end of a stepwise procedure based on pseudo-Flikelihood ratio tests. Simulations are used to compare the stepwise procedure involving I ^M(ω) with a stepwise procedure using Schuster's periodogram, to study an approximation of the asymptotic theory for the frequency estimators in finite samples in relation to the proximity and signal-to-noise ratio of the periodic components, and to assess the robustness of I ^M(ω) against autocorrelation in the analysis of mixed spectra. Overall, the results show an improvement of the new method over the classical approach when spectral lines are adjacent. Finally, three examples with real data illustrate specific aspects of the method, and extensions (i.e., unequally spaced observations, trend modeling, replicated time series, periodogram matrices) are outlined.     

Data-Driven Tests for Trend

This article proposes a new nonparametric test for the ordered alternatives problem in the k-sample setting for null hypothesis of lack of trend. This article further elaborates upon and extends the results of Ledwina and Wy?upek (2012a Ledwina , T. , Wy?upek , G. ( 2012a ). Two-sample test against one-sided alternatives . Scand. J. Statist. 39 : 358 – 381 .[Crossref], [Web of Science ®] , [Google Scholar]) obtained for k = 2. Simulations show that the new test has high and stable power and is able to control the Type I error to satisfactory extent, thus solving the problem posed in Terpstra and Magel (2003 Terpstra , J. T. , Magel , R. C. ( 2003 ). A new nonparametric test for the ordered alternative problem . J. Nonparametr. Statist. 15 : 289 – 301 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Our theoretical results say that asymptotic errors of both kinds do not exceed significance level, thus implying that the test is asymptotically unbiased.     

Shrinkage Estimation Under Multivariate Elliptic Models

The estimation of the location vector of a p-variate elliptically contoured distribution (ECD) is considered using independent random samples from two multivariate elliptically contoured populations when it is apriori suspected that the location vectors of the two populations are equal. For the setting where the covariance structure of the populations is the same, we define the maximum likelihood, Stein-type shrinkage and positive-rule shrinkage estimators. The exact expressions for the bias and quadratic risk functions of the estimators are derived. The comparison of the quadratic risk functions reveals the dominance of the Stein-type estimators if p ≥ 3. A graphical illustration of the risk functions under a “typical” member of the elliptically contoured family of distributions is provided to confirm the analytical results.     

Mixed sampling plan based on exponentially weighted moving average statistic

ABSTRACT

We will design a new mixed acceptance sampling plan based on the exponentially weighted moving average statistic in this article. The plan parameters of the proposed plan are determined by an optimization problem. The efficiency of the proposed plan is compared with the existing attribute sampling plan. An industrial example is given for illustration purpose.     

Approximation of distributions by using the Anderson Darling statistic

ABSTRACT

In practice, it is often not possible to find an appropriate family of distributions which can be used for fitting the sample distribution with high precision. In these cases, it seems to be opportune to search for the best approximation by a family of distributions instead of an exact fit. In this paper, we consider the Anderson–Darling statistic with plugged-in minimum distance estimator for the parameter vector. We prove asymptotic normality of the Anderson–Darling statistic which is used for a test of goodness of approximation. Moreover, we introduce a measure of discrepancy between the sample distribution and the model class.     

Restricted estimation and testing of hypothesis in linear measurement errors models

ABSTRACT

In this article, the linear models with measurement error both in the response and in the covariates are considered. Following Shalabh et al. (2007 Shalabh, Garg, G., Misra, N. (2007). Restricted regression estimation in measurement error models. Comput. Stat. Data Anal. 52:1149–1166.[Crossref], [Web of Science ®] , [Google Scholar], 2009 Shalabh, Garg, G., Misra, N. (2009). Use of prior information in the consistent estimation of regression coefficients in measurement error models. J. Multivariate Anal. 100:1498–1520.[Crossref], [Web of Science ®] , [Google Scholar]), we propose several restricted estimators for the regression coefficients. The consistency and asymptotic normality of the restricted estimators are established. Furthermore, we also discuss the superiority of the restricted estimators to unrestricted estimators under Pitman closeness criterion. We also develop several variance estimators and establish their asymptotic distributions. Wald-type statistics are constructed for testing the linear restrictions. Finally, Monte Carlo simulations are conducted to illustrate the finite-sample properties of the proposed estimators.     

Some Results on Dynamic Generalized Survival Entropy

Recently, Zografos and Nadarajah (2005 Zografos, K., Nadarajah, S. (2005). Survival exponential entropies. IEEE Trans. Inform. Theor. 51:1239–1246.[Crossref], [Web of Science ®] , [Google Scholar]) proposed two measures of uncertainty based on the survival function, called the survival exponential entropy and the generalized survival exponential entropy. In this article, we explore properties of the generalized survival entropy and the dynamic version of it. We study conditions under which the generalized survival entropy of first order statistic can uniquely determines the parent distribution. The exponential, Pareto, and finite range distributions, which are commonly used in reliability, have been characterized using this generalized measure. Another measure of entropy is also introduced in analogy with cumulative entropy which has been proposed by Di Crescenzo and Longobardi (2009) and some properties of it are given.     

A changepoint statistic with uniform type I error probabilities

ABSTRACT

Likelihood ratio tests for a change in mean in a sequence of independent, normal random variables are based on the maximum two-sample t-statistic, where the maximum is taken over all possible changepoints. The maximum t-statistic has the undesirable characteristic that Type I errors are not uniformly distributed across possible changepoints. False positives occur more frequently near the ends of the sequence and occur less frequently near the middle of the sequence. In this paper we describe an alternative statistic that is based upon a minimum p-value, where the minimum is taken over all possible changepoints. The p-value at any particular changepoint is based upon both the two-sample t-statistic at that changepoint and the probability that the maximum two-sample t-statistic is achieved at that changepoint. The new statistic has a more uniform distribution of Type I errors across potential changepoints and it compares favorably with respect to statistical power, false discovery rates, and the mean square error of changepoint estimates.     

An extension of the gamma distribution

ABSTRACT

The gamma distribution has been widely used in many research areas such as engineering and survival analysis. We present an extension of this distribution, called the Kummer beta gamma distribution, having greater flexibility to model scenarios involving skewed data. We derive analytical expressions for some mathematical quantities. The estimation of parameters is approached by the maximum likelihood method and Bayesian analysis. The likelihood ratio and formal goodness-of-fit tests are used to compare the presented distribution with some of its sub-models and non nested models. A real data set is used to illustrate the importance of the distribution.