Locally best rotation-invariant rank tests for modal location

For a general class of unipolar, rotationally symmetric distributions on the multi-dimensional unit spherical surface, a characterization of locally best rotation-invariant test statistics is exploited in the construction of locally best rotation-invariant rank tests for modal location. Allied statistical distributional problems are appraised, and in the light of these assessments, asymptotic relative efficiency of a class of rotation-invariant rank tests (with respect to some of their parametric counterparts) is studied. Finite sample permutational distributional perspectives are also appraised.     

A test of uniformity on shape spaces

A test of uniformity on the shape space Σ_m^k is presented, together with modifications of the test statistic which bring its null distribution close to the large-sample asymptotic distribution. The asymptotic distribution under suitable local alternatives to uniformity is given. A family of distributions on Σ_m^k is proposed, which is suitable for modelling shapes given by landmarks which are almost collinear.     

Affine-invariant aligned rank tests for the multivariate general linear model with VARMA errors

We develop optimal rank-based procedures for testing affine-invariant linear hypotheses on the parameters of a multivariate general linear model with elliptical VARMA errors. We propose a class of optimal procedures that are based either on residual (pseudo-)Mahalanobis signs and ranks, or on absolute interdirections and lift-interdirection ranks, i.e., on hyperplane-based signs and ranks. The Mahalanobis versions of these procedures are strictly affine-invariant, while the hyperplane-based ones are asymptotically affine-invariant. Both versions generalize the univariate signed rank procedures proposed by Hallin and Puri (J. Multivar. Anal. 50 (1994) 175), and are locally asymptotically most stringent under correctly specified radial densities. Their AREs with respect to Gaussian procedures are shown to be convex linear combinations of the AREs obtained in Hallin and Paindaveine (Ann. Statist. 30 (2002) 1103; Bernoulli 8 (2002) 787) for the pure location and purely serial models, respectively. The resulting test statistics are provided under closed form for several important particular cases, including multivariate Durbin-Watson tests, VARMA order identification tests, etc. The key technical result is a multivariate asymptotic linearity result proved in Hallin and Paindaveine (Asymptotic linearity of serial and nonserial multivariate signed rank statistics, submitted).     

Third-order efficiency of conditional tests in exponential models: The lattice case

As is well known, in full rank multivariate exponential families, tests of Neyman structure are uniformly most powerful unbiased for one-sided problems. For the case of lattice distributions, the power of these tests—evaluated at contiguous alternatives—is approximated by asymptotic expansions up to errors of order o(n^?1). Surprisingly the tests with Neyman structure are not third-order efficient in the class of all asymptotically similar tests unless the problem is univariate.     

Smoothed nonparametric tests and approximations of <Emphasis Type="Italic">p</Emphasis>-values

We propose new smoothed sign and Wilcoxon’s signed rank tests that are based on kernel estimators of the underlying distribution function of the data. We discuss the approximations of the p-values and asymptotic properties of these tests. The new smoothed tests are equivalent to the ordinary sign and Wilcoxon’s tests in the sense of Pitman’s asymptotic relative efficiency, and the differences between the ordinary and new tests converge to zero in probability. Under the null hypothesis, the main terms of the asymptotic expectations and variances of the tests do not depend on the underlying distribution. Although the smoothed tests are not distribution-free, making use of the specific kernel enables us to obtain the Edgeworth expansions, being free of the underlying distribution.     

An asymptotic expansion of the distribution of Rao's U-statistic under a general condition

In this paper we consider the problem of testing the hypothesis about the sub-mean vector. For this propose, the asymptotic expansion of the null distribution of Rao's U-statistic under a general condition is obtained up to order of n^-1. The same problem in the k-sample case is also investigated. We find that the asymptotic distribution of generalized U-statistic in the k-sample case is identical to that of the generalized Hotelling's T² distribution up to n^-1. A simulation experiment is carried out and its results are presented. It shows that the asymptotic distributions have significant improvement when comparing with the limiting distributions both in the small sample case and the large sample case. It also demonstrates the equivalence of two testing statistics mentioned above.     

Local efficiency of a Cramér-von Mises test of independence

Deheuvels proposed a rank test of independence based on a Cramér-von Mises functional of the empirical copula process. Using a general result on the asymptotic distribution of this process under sequences of contiguous alternatives, the local power curve of Deheuvels’ test is computed in the bivariate case and compared to that of competing procedures based on linear rank statistics. The Gil-Pelaez inversion formula is used to make additional comparisons in terms of a natural extension of Pitman's measure of asymptotic relative efficiency.     

Practical tests for randomized complete block designs

A class of affine-invariant test statistics, including a sign test and a related family of signed-rank tests, is proposed for randomized complete block designs with one observation per treatment. This class is obtained by using the transformation-retransformation approach of Chakraborty, Chaudhuri and Oja along with a directional transformation due to Tyler. Under the minimal assumption of directional symmetry of the underlying distribution, the null asymptotic distribution of the sign test statistic is shown to be chi-square with p-1 degrees of freedom. The same null distribution is also proved for the family of signed-rank statistics under the assumption of symmetry of the underlying distribution. The Pitman asymptotic relative efficiencies of the tests, relative to Hotelling-Hsu's T² are established. Several score functions are discussed including a simple linear score function and the optimal normal score function. The test based on the linear score function is compared to the other members of this family and other statistics in the literature through efficiency calculations and Monte Carlo simulations. This statistic has an excellent performance over a wide range of distributions and for small as well as large dimensions.     

Nonparametric methods for unbalanced multivariate data and many factor levels

We propose different nonparametric tests for multivariate data and derive their asymptotic distribution for unbalanced designs in which the number of factor levels tends to infinity (large a, small n_i case). Quasi gratis, some new parametric multivariate tests suitable for the large a asymptotic case are also obtained. Finite sample performances are investigated and compared in a simulation study. The nonparametric tests are based on separate rankings for the different variables. In the presence of outliers, the proposed nonparametric methods have better power than their parametric counterparts. Application of the new tests is demonstrated using data from plant pathology.     

On a new multivariate two-sample test

In this paper we propose a new test for the multivariate two-sample problem. The test statistic is the difference of the sum of all the Euclidean interpoint distances between the random variables from the two different samples and one-half of the two corresponding sums of distances of the variables within the same sample. The asymptotic null distribution of the test statistic is derived using the projection method and shown to be the limit of the bootstrap distribution. A simulation study includes the comparison of univariate and multivariate normal distributions for location and dispersion alternatives. For normal location alternatives the new test is shown to have power similar to that of the t- and T²-Test.     

Rank procedures for some two-population multivariate extended classification problems

Given independent samples from three multivariate populations with cumulative distribution functions F⁽¹⁾(x), F⁽²⁾(x), and F⁽⁰⁾(x) = θF⁽¹⁾(x) + (1 ? θ)F⁽²⁾(x), where 0 ≤ θ ≤ 1 is unknown, the three-action problem involving decision as to whether the value of θ is high, low, or intermediate, is considered. A class of consistent procedures based on the relative spacing of three sample averages of linearly compounded rank scores is formulated. The asymptotic operating characteristics of the procedures when F⁽¹⁾ and F⁽²⁾ come close together are studied and the best choice of the compounding coefficients in terms of these considered. The consequence of using estimates of the best coefficients on the asymptotic operating characteristics is also examined.     

Nonnull distribution of likelihood ratio criterion for reality of covariance matrix

In this paper the distribution of the likelihood ratio test for testing the reality of the covariance matrix of a complex multivariate normal distribution is investigated. Some simplifications in the noncentral distribution are made and the noncentral distribution is derived for the special case where the rank of the noncentrality matrix is two. In the null case exact expressions for the distribution are given up to p = 6, and percentage points are tabulated. These percentage points were compared with percentage points derived from an asymptotic expansion of the distribution, and the accuracy of the approximation was found to be sufficient for several practical situations.     

Dimension estimation in sufficient dimension reduction: A unifying approach

Sufficient Dimension Reduction (SDR) in regression comprises the estimation of the dimension of the smallest (central) dimension reduction subspace and its basis elements. For SDR methods based on a kernel matrix, such as SIR and SAVE, the dimension estimation is equivalent to the estimation of the rank of a random matrix which is the sample based estimate of the kernel. A test for the rank of a random matrix amounts to testing how many of its eigen or singular values are equal to zero. We propose two tests based on the smallest eigen or singular values of the estimated matrix: an asymptotic weighted chi-square test and a Wald-type asymptotic chi-square test. We also provide an asymptotic chi-square test for assessing whether elements of the left singular vectors of the random matrix are zero. These methods together constitute a unified approach for all SDR methods based on a kernel matrix that covers estimation of the central subspace and its dimension, as well as assessment of variable contribution to the lower-dimensional predictor projections with variable selection, a special case. A small power simulation study shows that the proposed and existing tests, specific to each SDR method, perform similarly with respect to power and achievement of the nominal level. Also, the importance of the choice of the number of slices as a tuning parameter is further exhibited.     

HAC estimation and strong linearity testing in weak ARMA models

In the framework of ARMA models, we consider testing the reliability of the standard asymptotic covariance matrix (ACM) of the least-squares estimator. The standard formula for this ACM is derived under the assumption that the errors are independent and identically distributed, and is in general invalid when the errors are only uncorrelated. The test statistic is based on the difference between a conventional estimator of the ACM of the least-squares estimator of the ARMA coefficients and its robust HAC-type version. The asymptotic distribution of the HAC estimator is established under the null hypothesis of independence, and under a large class of alternatives. The asymptotic distribution of the proposed statistic is shown to be a standard χ² under the null, and a noncentral χ² under the alternatives. The choice of the HAC estimator is discussed through asymptotic power comparisons. The finite sample properties of the test are analyzed via Monte Carlo simulation.     

Asymptotic distributions of robust shape matrices and scales

It has been frequently observed in the literature that many multivariate statistical methods require the covariance or dispersion matrix Σ of an elliptical distribution only up to some scaling constant. If the topic of interest is not the scale but only the shape of the elliptical distribution, it is not meaningful to focus on the asymptotic distribution of an estimator for Σ or another matrix Γ∝Σ. In the present work, robust estimators for the shape matrix and the associated scale are investigated. Explicit expressions for their joint asymptotic distributions are derived. It turns out that if the joint asymptotic distribution is normal, the estimators presented are asymptotically independent for one and only one specific choice of the scale function. If it is non-normal (this holds for example if the estimators for the shape matrix and scale are based on the minimum volume ellipsoid estimator) only the scale function presented leads to asymptotically uncorrelated estimators. This is a generalization of a result obtained by Paindaveine [D. Paindaveine, A canonical definition of shape, Statistics and Probability Letters 78 (2008) 2240-2247] in the context of local asymptotic normality theory.     

Local asymptotic normality in a stationary model for spatial extremes

De Haan and Pereira (2006) [6] provided models for spatial extremes in the case of stationarity, which depend on just one parameter β>0 measuring tail dependence, and they proposed different estimators for this parameter. We supplement this framework by establishing local asymptotic normality (LAN) of a corresponding point process of exceedances above a high multivariate threshold. Standard arguments from LAN theory then provide the asymptotic minimum variance within the class of regular estimators of β. It turns out that the relative frequency of exceedances is a regular estimator sequence with asymptotic minimum variance, if the underlying observations follow a multivariate extreme value distribution or a multivariate generalized Pareto distribution.     

Nonparametric tests of independence between random vectors

A nonparametric test of the mutual independence between many numerical random vectors is proposed. This test is based on a characterization of mutual independence defined from probabilities of half-spaces in a combinatorial formula of Möbius. As such, it is a natural generalization of tests of independence between univariate random variables using the empirical distribution function. If the number of vectors is p and there are n observations, the test is defined from a collection of processes R_n,A, where A is a subset of {1,…,p} of cardinality |A|>1, which are asymptotically independent and Gaussian. Without the assumption that each vector is one-dimensional with a continuous cumulative distribution function, any test of independence cannot be distribution free. The critical values of the proposed test are thus computed with the bootstrap which is shown to be consistent. Another similar test, with the same asymptotic properties, for the serial independence of a multivariate stationary sequence is also proposed. The proposed test works when some or all of the marginal distributions are singular with respect to Lebesgue measure. Moreover, in singular cases described in Section 4, the test inherits useful invariance properties from the general affine invariance property.     

Local powers of two-sample and multi-sample rank tests for Lehmann's contaminated alternative

Summary Local powers of two- andk-sample rank tests under alternatives of contaminated distributions are investigated. It is shown that the rank tests based on normal scores and Wilcoxon scores are superior to thet-test or theF-test for many choices of alternatives of contaminated distributions and that the values of the asymptotic relative efficiency of the rank test based on Wilcoxon scores with respect to the normal scores are about one in all the investigated cases.     

Optimal adaptive generalized Pólya urn design for multi-arm clinical trials

A class of optimal adaptive multi-arm clinical trial designs is proposed based on an extended generalized Pólya urn (GPU) model. The design is applicable to both the qualitative and quantitative responses and achieves, asymptotically, some pre-specified optimality criterion. Such criterion is specified by a functional of the response distributions and is implemented through the relationship between the design matrix and its first eigenvector. The asymptotic properties of the design are studied using the existing methods on GPU. Some examples for commonly used clinical designs are given as illustration.     

Some high-dimensional tests for a one-way MANOVA

A statistic is proposed for testing the equality of the mean vectors in a one-way multivariate analysis of variance. The asymptotic null distribution of this statistic, as both the sample size and the number of variables go to infinity, is shown to be normal. Thus, this test can be used when the number of variables is not small relative to the sample size. In particular, it can be used when the number of variables exceeds the degrees of freedom for error, a situation in which standard MANOVA tests are invalid. A related statistic, also having an asymptotic normal distribution, is developed for tests concerning the dimensionality of the hyperplane formed by the population mean vectors. The finite sample size performances of the normal approximations are evaluated in a simulation study.