A multivariate dispersion ordering based on quantiles more widely separated

A multivariate dispersion ordering based on quantiles more widely separated is defined. This new multivariate dispersion ordering is a generalization of the classic univariate version. If we vary the ordering of the components in the multivariate random variable then the comparison could not be possible. We provide a characterization using a multivariate expansion function. The relationship among various multivariate orderings is also considered. Finally, several examples illustrate the method of this paper.     

The simplex dispersion ordering and its application to the evaluation of human corneal endothelia

A multivariate dispersion ordering based on random simplices is proposed in this paper. Given a R^d-valued random vector, we consider two random simplices determined by the convex hulls of two independent random samples of sizes d+1 of the vector. By means of the stochastic comparison of the Hausdorff distances between such simplices, a multivariate dispersion ordering is introduced. Main properties of the new ordering are studied. Relationships with other dispersion orderings are considered, placing emphasis on the univariate version. Some statistical tests for the new order are proposed. An application of such ordering to the clinical evaluation of human corneal endothelia is provided. Different analyses are included using an image database of human corneal endothelia.     

Preserving multivariate dispersion: An application to the Wishart distribution

Univariate dispersive ordering has been extensively characterized by many authors over the last two decades. However, the multivariate version lacks extensive analysis. In this paper, sufficient and necessary conditions are given to preserve the strong multivariate dispersion order through properties of the corresponding transformation. Finally, these results are applied to the Wishart distribution which can be viewed as “the spread of the dispersion”.     

Total positivity order and the normal distribution

Unlike the usual stochastic order, total positivity order is closed under conditioning. Here we provide a general formulation of the preservation properties of the order under conditioning; we study certain properties of the order including translation properties and the implications of having equality in the inequality defining the order. Specializing to the multivariate normal distribution, the study of total positivity order leads to new cones defined in terms of covariance M-matrices related to positive dependence, whose properties we study.     

Dependence structure of conditional Archimedean copulas

In this article, copulas associated to multivariate conditional distributions in an Archimedean model are characterized. It is shown that this popular class of dependence structures is closed under the operation of conditioning, but that the associated conditional copula has a different analytical form in general. It is also demonstrated that the extremal copula for conditional Archimedean distributions is no longer the Fréchet upper bound, but rather a member of the Clayton family. Properties of these conditional distributions as well as conditional versions of tail dependence indices are also considered.     

A new dependence ordering with applications

In this paper, we introduce a new copula-based dependence order to compare the relative degree of dependence between two pairs of random variables. Relationship of the new order to the existing dependence orders is investigated. In particular, the new ordering is stronger than the partial ordering, more monotone regression dependence as developed by Avérous et al. [J. Avérous, C. Genest, S.C. Kochar, On dependence structure of order statistics, Journal of Multivariate Analysis 94 (2005) 159-171]. Applications of this partial order to order statistics, k-record values and frailty models are given.     

Regression models for multivariate ordered responses via the Plackett distribution

We investigate the properties of a class of discrete multivariate distributions whose univariate marginals have ordered categories, all the bivariate marginals, like in the Plackett distribution, have log-odds ratios which do not depend on cut points and all higher-order interactions are constrained to 0. We show that this class of distributions may be interpreted as a discretized version of a multivariate continuous distribution having univariate logistic marginals. Convenient features of this class relative to the class of ordered probit models (the discretized version of the multivariate normal) are highlighted. Relevant properties of this distribution like quadratic log-linear expansion, invariance to collapsing of adjacent categories, properties related to positive dependence, marginalization and conditioning are discussed briefly. When continuous explanatory variables are available, regression models may be fitted to relate the univariate logits (as in a proportional odds model) and the log-odds ratios to covariates.     

Choosing joint distributions so that the variance of the sum is small

The paper considers how to choose the joint distribution of several random variables each with a given marginal distribution so that their sum has a variance as small as possible. A theorem is given that allows the solution of this and of related problems for normal random variables. Several specific applications are given. Additional results are provided for radially symmetric joint distributions of three random variables when the sum is identically zero.     

Some new results on multivariate dispersive ordering of generalized order statistics

Let be generalized order statistics based on a continuous distribution function F with parameters k and (m₁,…,m_n−1). Chen and Hu (2007) [8] investigated the sufficient conditions on F and on the parameters k and m_i’s such that , where , and is the Shaked-Shanthikumar multivariate dispersive order. Since the order does not possess the closure property under marginalization, one may naturally wonder whether the corresponding multivariate margins of the above random vectors are also ordered in the order . This is answered affirmatively in this paper. Some comparison results for generalized order statistics from two samples are presented. Potential applications are also mentioned.     

On the asymptotic power of some k-sample statistics based on the multivariate empirical process

Some k-sample Kolmogorov-Smirnov and Cramér-von Mises-type statistics, based on the multivariate empirical process, are studied. Expressions for their asymptotic power are obtained against various classes of alternative distribution functions.     

Dependence structures of multivariate Bernoulli random vectors

In some situations, it is difficult and tedious to check notions of dependence properties and dependence orders for multivariate distributions supported on a finite lattice. The purpose of this paper is to utilize a newly developed tool, majorization with respect to weighted trees, to lay out some general results that can be used to identify some dependence properties and dependence orders for multivariate Bernoulli random vectors. Such a study gives us some new insight into the relations between the concepts of dependence.     

A trivariate chi-squared distribution derived from the complex Wishart distribution

The joint density for a particular trivariate chi-squared distribution given by the diagonal elements of a complex Wishart matrix is derived. This distribution has applications in the processing of multilook synthetic aperture radar data. The expression for the density is in the form of an infinite series that converges rapidly and is simple and fast to compute. The expression is shown to reduce to known forms for a number of special cases and is validated by simulation. The characteristic function is also derived and used to relate joint moments of the trivariate distribution to the parameters of the density function.     

Bounds for a multivariate extension of range over standard deviation based on the Mahalanobis distance

The range over standard deviation of a set of univariate data points is given a natural multivariate extension through the Mahalanobis distance. The problem of finding extrema of this multivariate extension of “range over standard deviation” is investigated. The supremum (maximum) is found using Lagrangian methods and an interval is given for the infinimum. The independence of optimizing the Mahalanobis distance and the multivariate extension of range is demonstrated and connections are explored in several examples using an analogue of the “hat” matrix of linear regression.     

Probability and moment inequalities for sums of weakly dependent random variables,with applications

Doukhan and Louhichi [P. Doukhan, S. Louhichi, A new weak dependence condition and application to moment inequalities, Stochastic Process. Appl. 84 (1999) 313–342] introduced a new concept of weak dependence which is more general than mixing. Such conditions are particularly well suited for deriving estimates for the cumulants of sums of random variables. We employ such cumulant estimates to derive inequalities of Bernstein and Rosenthal type which both improve on previous results. Furthermore, we consider several classes of processes and show that they fulfill appropriate weak dependence conditions. We also sketch applications of our inequalities in probability and statistics.     

On the singularity of multivariate skew-symmetric models

In recent years, the skew-normal models introduced by Azzalini (1985) [1]-and their multivariate generalizations from Azzalini and Dalla Valle (1996) [4]-have enjoyed an amazing success, although an important literature has reported that they exhibit, in the vicinity of symmetry, singular Fisher information matrices and stationary points in the profile log-likelihood function for skewness, with the usual unpleasant consequences for inference. It has been shown (DiCiccio and Monti (2004) [23], DiCiccio and Monti (2009) [24] and Gómez et al. (2007) [25]) that these singularities, in some specific parametric extensions of skew-normal models (such as the classes of skew-t or skew-exponential power distributions), appear at skew-normal distributions only. Yet, an important question remains open: in broader semiparametric models of skewed distributions (such as the general skew-symmetric and skew-elliptical ones), which symmetric kernels lead to such singularities? The present paper provides an answer to this question. In very general (possibly multivariate) skew-symmetric models, we characterize, for each possible value of the rank of Fisher information matrices, the class of symmetric kernels achieving the corresponding rank. Our results show that, for strictly multivariate skew-symmetric models, not only Gaussian kernels yield singular Fisher information matrices. In contrast, we prove that systematic stationary points in the profile log-likelihood functions are obtained for (multi)normal kernels only. Finally, we also discuss the implications of such singularities on inference.     

The Stein phenomenon for monotone incomplete multivariate normal data

We establish the Stein phenomenon in the context of two-step, monotone incomplete data drawn from , a (p+q)-dimensional multivariate normal population with mean and covariance matrix . On the basis of data consisting of n observations on all p+q characteristics and an additional N−n observations on the last q characteristics, where all observations are mutually independent, denote by the maximum likelihood estimator of . We establish criteria which imply that shrinkage estimators of James-Stein type have lower risk than under Euclidean quadratic loss. Further, we show that the corresponding positive-part estimators have lower risk than their unrestricted counterparts, thereby rendering the latter estimators inadmissible. We derive results for the case in which is block-diagonal, the loss function is quadratic and non-spherical, and the shrinkage estimator is constructed by means of a nondecreasing, differentiable function of a quadratic form in . For the problem of shrinking to a vector whose components have a common value constructed from the data, we derive improved shrinkage estimators and again determine conditions under which the positive-part analogs have lower risk than their unrestricted counterparts.     

A novel trace test for the mean parameters in a multivariate growth curve model

A trace test for the mean parameters of the growth curve model is proposed. It is constructed using the restricted maximum likelihood followed by an estimated likelihood ratio approach. The statistic reduces to the Lawley-Hotelling trace test for the Multivariate Analysis of Variance (MANOVA) models. Our test statistic is, therefore, a natural extension of the classical trace test to GMANOVA models. We show that the distribution of the test under the null hypothesis does not depend on the unknown covariance matrix Σ. We also show that the distributions under the null and alternative hypotheses can be represented as sums of weighted central and non-central chi-square random variables, respectively. Under the null hypothesis, the Satterthwaite approximation is used to get an approximate critical point. A novel Satterthwaite type approximation is proposed to obtain an approximate power. A simulation study is performed to evaluate the performance of our proposed test and numerical examples are provided as illustrations.     

The distribution of the residual from a general elliptical multivariate linear model

Given a general multivariate linear model of full or less than full rank, we find the distributions of internally and externally studentised residuals, assuming normal and elliptical distributions.