Outlier detection by robust principal components analysis

The robust principal components analysis (RPCA) introduced by Campbell (Applied Statistics 1980, 29, 231–237) provides in addition to robust versions of the usual output of a principal components analysis, weights for the contribution of each point to the robust estimation of each component. Low weights may thus be used to indicate outliers. The present simulation study provides critical values for testing the kth smallest weight in the RPCA of a sample of n p-dimensional vectors, under the null hypothesis of a multivariate normal distribution. The cases p=2(2)10, 15, 20 for n=20, 30, 40, 50, 75, 100 subject to n≥p/2, are examined, with k≤√n.     

Estimation of gaussian mixtures with rotationally invariant covariance matrices

Homoscedastic and heteroscedastic Gaussian mixtures differ in the constraints placed on the covariance matrices of the mixture components. A new mixture, called herein a strophoscedastic mixture, is defined by a new constraint, This constraint requires the matrices to be identical under orthogonal trans¬formations, where different transformations are allowed for different matrices. It is shown that the M-step of the EM method for estimating the parameters of strophoscedastic mixtures from sample data is explicitly solvable using singular value decompositions. Consequently, the EM-based maximum likelihood estimation algorithm is as easily implemented for strophoscedastic mixtures as it is for homoscedastic and heteroscedastic mixtures. An example of a “noisy” Archimedian spiral is presented.     

An extension of bayesian measure of information to regression

This paper extends Lindley's measure of average information to the linear model, E(Y∣ß) = Xß. An expression which quantifies the average amount of information provided by the nxl vector of observations Y about the pxl vector of coefficient parameters ß will be derived. The effect of the structure of the regressor matrix, X, on the information measure is discussed. An information theoretic optimal design is characterized. Some applications are suggested.     

Partial Round-Robin Comparisons with Perfect Rankings

In teaching the development of uniformly most powerful unbiased (UMPU) tests, one rarely discusses the performance of alternative biased tests. It is shown, through the comparison of two independent Bernoulli proportions, that a biased test (the Z test) can be more powerful than the UMPU test (Fisher's exact test—randomized) in a large region of the alternative parameter space. A more general example is also given.     

A robust procedure for testing an assumed value of the population correlation coefficient

For testing that the population correlations coefficientp Q, Tiku and Balakrishnan1986) developed a robust test. This test is extended here to situcitions where one wants to test that p p , p being a specified non-zero value of p. o o     

Assessing convergence diagnostic tests for Bayesian Cox regression

The Markov chain Monte Carlo (MCMC) method generates samples from the posterior distribution and uses these samples to approximate expectations of quantities of interest. For the process, researchers have to decide whether the Markov chain has reached the desired posterior distribution. Using convergence diagnostic tests are very important to decide whether the Markov chain has reached the target distribution. Our interest in this study was to compare the performances of convergence diagnostic tests for all parameters of Bayesian Cox regression model with different number of iterations by using a simulation and a real lung cancer dataset.     

A class of distribution-free tests for two-way layout

In this article a class of distribution-free tests for the hypothesis of no row (treatment) effect in a two-way layout design, with several observations per cell, is proposed. The tests are based on U-statistics, constructed by considering minima of all possible subsamples of same size from each cell.The proposed class of tests is compared with the parametric test, Mack and Skillings test and Yate's test for two-way layout, in terms of Pitman ARE sense. It is seen that for the case of equal number of observations per cell, the proposed tests have better efficiency for exponential and uniform error distributions.     

The power functions of the likelihood ratio tests for a simply ordered trend in normal means

Likelihood ratio tests for the homogeneity of k normal means with the alternative restricted by an increasing trend are considered as well as the likelihood ratio tests of the null hypothesis that the means satisfy the trend. While the work is primarily a survey of results concerning the power functions of these tests, the extensions of some results to the case of not necessarily equal sample sizes are presented. For the case of known or unknown population variances, exact expressions are given for the power functions for k=3,4, and approximations are discussed for larger k. The topics of consistency, bias and monotonicity of the power functions are included. Also, Bartholomew's conjectures concerning minimal and maximal powers are investigated, with results of a new numerical study given.     

Elliptically contoured model and factorization of wllks1 A

In a series of papers, Kshirsagar (1964, 1971) and McHenry and Kshirsagar (1977), factorize Wilks' A into a number of factors and find the independent null multivariate beta densities of these factors. These factors are the likelihood ratio test criteria for testing the goodness of fit of certain assigned discriminant functions or canonical variables either in the space of independent or dependent variables. Essentially the factors of Wilks' A are the factors of certain multivariate beta distributed matrix or its determinant. The Bartlett decomposition of the underlying multivariate beta distribution into independent factors determines the distribution of these factors. The present paper generalizes Kshirsagar's (1971) normal theory to the elliptically contoured model, and shows that his results are null robust for the elliptically contoured model.