Flexible bivariate beta distributions

Bivariate beta distributions which can be used to model data sets exhibiting positive or negative correlation are introduced. Properties of these bivariate beta distributions and their applications in Bayesian analysis are discussed. Three methods for parameter estimation are presented. The performance of these estimators is evaluated based on Monte Carlo simulations. Examples are provided to illustrate how additional parameters can be introduced to gain even more modeling flexibility. A possible extension of the proposed bivariate beta model and a multivariate generalization are also discussed.     

Characterizations of Arnold and Strauss’ and related bivariate exponential models

Characterizations of probability distributions is a topic of great popularity in applied probability and reliability literature for over last 30 years. Beside the intrinsic mathematical interest (often related to functional equations) the results in this area are helpful for probabilistic and statistical modelling, especially in engineering and biostatistical problems. A substantial number of characterizations has been devoted to a legion of variants of exponential distributions. The main reliability measures associated with a random vector X are the conditional moment function defined by m_φ(x)=E(φ(X)|X?x) (which is equivalent to the mean residual life function e(x)=m_φ(x)-x when φ(x)=x) and the hazard gradient function h(x)=-∇logR(x), where R(x) is the reliability (survival) function, R(x)=Pr(X?x), and ∇ is the operator . In this paper we study the consequences of a linear relationship between the hazard gradient and the conditional moment functions for continuous bivariate and multivariate distributions. We obtain a general characterization result which is the applied to characterize Arnold and Strauss’ bivariate exponential distribution and some related models.     

Smooth depth contours characterize the underlying distribution

The Tukey depth is an innovative concept in multivariate data analysis. It can be utilized to extend the univariate order concept and advantages to a multivariate setting. While it is still an open question as to whether the depth contours uniquely determine the underlying distribution, some positive answers have been provided. We extend these results to distributions with smooth depth contours, with elliptically symmetric distributions as special cases. The key ingredient of our proofs is the well-known Cramér-Wold theorem.     

Moment properties of multivariate infinitely divisible laws and criteria for multivariate self-decomposability

Ramachandran (1969) [9, Theorem 8] has shown that for any univariate infinitely divisible distribution and any positive real number α, an absolute moment of order α relative to the distribution exists (as a finite number) if and only if this is so for a certain truncated version of the corresponding Lévy measure. A generalized version of this result in the case of multivariate infinitely divisible distributions, involving the concept of g-moments, was given by Sato (1999) [6, Theorem 25.3]. We extend Ramachandran’s theorem to the multivariate case, keeping in mind the immediate requirements under appropriate assumptions of cumulant studies of the distributions referred to; the format of Sato’s theorem just referred to obviously varies from ours and seems to have a different agenda. Also, appealing to a further criterion based on the Lévy measure, we identify in a certain class of multivariate infinitely divisible distributions the distributions that are self-decomposable; this throws new light on structural aspects of certain multivariate distributions such as the multivariate generalized hyperbolic distributions studied by Barndorff-Nielsen (1977) [12] and others. Various points relevant to the study are also addressed through specific examples.     

A contribution to multivariate L-moments: L-comoment matrices

Multivariate statistical analysis relies heavily on moment assumptions of second order and higher. With increasing interest in heavy-tailed distributions, however, it is desirable to describe dispersion, skewness, and kurtosis under merely first order moment assumptions. Here, the univariate L-moments of Hosking [L-moments: analysis and estimation of distributions using linear combinations of order statistics, J. Roy. Statist. Soc. Ser. B 52 (1990) 105-124] are extended to “L-comoments” analogous to covariance. For certain models, the second order case yields correlational analysis coherent with classical correlation but also meaningful under just first moment assumptions. We develop properties and estimators for L-comoments, illustrate for several multivariate models, examine behavior of sample multivariate L-moments with heavy-tailed data, and discuss applications to financial risk analysis and regional frequency analysis.     

Third-order power comparisons for a class of tests for multivariate linear hypothesis under general distributions

The purpose of this paper is, in multivariate linear regression model (Part I) and GMANOVA model (Part II), to investigate the effect of nonnormality upon the nonnull distributions of some multivariate test statistics under normality. It is shown that whatever the underlying distributions, the difference of local powers up to order N⁻¹ after either Bartlett’s type adjustment or Cornish-Fisher’s type size adjustment under nonnormality coincides with that in Anderson [An Introduction to Multivariate Statistical Analysis, 2nd ed. and 3rd ed., Wiley, New York, 1984, 2003] under normality. The derivation of asymptotic expansions is based on the differential operator associated with the multivariate linear regression model under general distributions. The performance of higher-order results in finite samples, including monotone Bartlett’s type adjustment and monotone Cornish-Fisher’s type size adjustment, is examined using simulation studies.     

Canonical representation of conditionally specified multivariate discrete distributions

Most work on conditionally specified distributions has focused on approaches that operate on the probability space, and the constraints on the probability space often make the study of their properties challenging. We propose decomposing both the joint and conditional discrete distributions into characterizing sets of canonical interactions, and we prove that certain interactions of a joint distribution are shared with its conditional distributions. This invariance opens the door for checking the compatibility between conditional distributions involving the same set of variables. We formulate necessary and sufficient conditions for the existence and uniqueness of discrete conditional models, and we show how a joint distribution can be easily computed from the pool of interactions collected from the conditional distributions. Hence, the methods can be used to calculate the exact distribution of a Gibbs sampler. Furthermore, issues such as how near compatibility can be reconciled are also discussed. Using mixed parametrization, we show that the proposed approach is based on the canonical parameters, while the conventional approaches are based on the mean parameters. Our advantage is partly due to the invariance that holds only for the canonical parameters.     

A new class of bivariate distributions and its mixture

A new class of bivariate distributions is presented in this paper. The procedure used in this paper is based on a latent random variable with exponential distribution. The model introduced here is of Marshall-Olkin type. A mixture of the proposed bivariate distributions is also discussed. The results obtained here generalize those of the bivariate exponential distribution present in the literature.     

A new test for multivariate normality

We propose a new class of rotation invariant and consistent goodness-of-fit tests for multivariate distributions based on Euclidean distance between sample elements. The proposed test applies to any multivariate distribution with finite second moments. In this article we apply the new method for testing multivariate normality when parameters are estimated. The resulting test is affine invariant and consistent against all fixed alternatives. A comparative Monte Carlo study suggests that our test is a powerful competitor to existing tests, and is very sensitive against heavy tailed alternatives.     

A weighted multivariate signed-rank test for cluster-correlated data

A weighted multivariate signed-rank test is introduced for an analysis of multivariate clustered data. Observations in different clusters may then get different weights. The test provides a robust and efficient alternative to normal theory based methods. Asymptotic theory is developed to find the approximate p-value as well as to calculate the limiting Pitman efficiency of the test. A conditionally distribution-free version of the test is also discussed. The finite-sample behavior of different versions of the test statistic is explored by simulations and the new test is compared to the unweighted and weighted versions of Hotelling’s T² test and the multivariate spatial sign test introduced in [D. Larocque, J. Nevalainen, H. Oja, A weighted multivariate sign test for cluster-correlated data, Biometrika 94 (2007) 267-283]. Finally, a real data example is used to illustrate the theory.     

Principal points of a multivariate mixture distribution

A set of n-principal points of a distribution is defined as a set of n points that optimally represent the distribution in terms of mean squared distance. It provides an optimal n-point-approximation of the distribution. However, it is in general difficult to find a set of principal points of a multivariate distribution. Tarpey et al. [T. Tarpey, L. Li, B. Flury, Principal points and self-consistent points of elliptical distributions, Ann. Statist. 23 (1995) 103-112] established a theorem which states that any set of n-principal points of an elliptically symmetric distribution is in the linear subspace spanned by some principal eigenvectors of the covariance matrix. This theorem, called a “principal subspace theorem”, is a strong tool for the calculation of principal points. In practice, we often come across distributions consisting of several subgroups. Hence it is of interest to know whether the principal subspace theorem remains valid even under such complex distributions. In this paper, we define a multivariate location mixture model. A theorem is established that clarifies a linear subspace in which n-principal points exist.     

A note on the maximum correlation for Baker’s bivariate distributions with fixed marginals

We investigate Baker’s bivariate distributions with fixed marginals which are based on order statistics, and find conditions under which the correlation converges to the maximum for Fréchet-Hoeffding upper bound as the sample size tends to infinity. The convergence rate of the correlation is also investigated for some specific cases.     

Multivariate dependence of spacings of generalized order statistics

Multivariate dependence of spacings of generalized order statistics is studied. It is shown that spacings of generalized order statistics from DFR (IFR) distributions have the CIS (CDS) property. By restricting the choice of the model parameters and strengthening the assumptions on the underlying distribution, stronger dependence relations are established. For instance, if the model parameters are decreasingly ordered and the underlying distribution has a log-convex decreasing (log-concave) hazard rate, then the spacings satisfy the MTP₂ (S- MRR₂) property. Some consequences of the results are given. In particular, conditions for non-negativity of the best linear unbiased estimator of the scale parameter in a location-scale family are obtained. By applying a result for dual generalized order statistics, we show that in the particular situation of usual order statistics the assumptions can be weakened.     

Wishart-Laplace distributions associated with matrix quadratic forms

For a normal random matrix Y with mean zero, necessary and sufficient conditions are obtained for Y^′W_kY to be Wishart-Laplace distributed and {Y^′W_kY} to be independent, where each W_k is assumed to be symmetric rather than nonnegative definite.     

Stochastic comparisons in multivariate mixed model of proportional reversed hazard rate with applications

This paper studies the multivariate mixed proportional reversed hazard rate model having dependent mixing variables. Stochastic comparison as well as aging properties in this model are investigated, and stochastic monotone properties of the population vector with respect to the mixing vector are also discussed. Moreover, MTP₂ dependence among the mixing vectors is proved to imply the increasingness of the reversed hazard rate with respect to the baseline one. Finally, some interesting applications are presented as well.     

A decomposition for a stochastic matrix with an application to MANOVA

The aim of this paper is to propose a simple method in order to evaluate the (approximate) distribution of matrix quadratic forms when Wishartness conditions do not hold. The method is based upon a factorization of a general Gaussian stochastic matrix as a special linear combination of nonstochastic matrices with the standard Gaussian matrix. An application of previous result is proposed for matrix quadratic forms arising in MANOVA for a multivariate split-plot design with circular dependence structure.     

Eigenanalysis on a bivariate covariance kernel

Certain constructions of copulas can be interpreted as an eigendecomposition of a kernel. We study some properties of the eigenfunctions and their integrals of a covariance kernel related to a bivariate distribution. The covariance between functions of random variables in terms of the cumulative distribution function is used. Some bounds for the trace of the kernel and some inequalities for a continuous random variable concerning a function and its derivative are obtained. We also obtain relations to diagonal expansions and canonical correlation analysis and, as a by-product, series of constants for some particular distributions.     

Robustness of A-optimal designs

Suppose that Y=(Y_i) is a normal random vector with mean Xb and covariance σ²I_n, where b is a p-dimensional vector (b_j),X=(X_ij) is an n×p matrix. A-optimal designs X are chosen from the traditional set D of A-optimal designs for ρ=0 such that X is still A-optimal in D when the components Y_i are dependent, i.e., for i≠i′, the covariance of Y_i,Y_i^′ is ρ with ρ≠0. Such designs depend on the sign of ρ. The general results are applied to X=(X_ij), where X_ij∈{-1,1}; this corresponds to a factorial design with -1,1 representing low level or high level respectively, or corresponds to a weighing design with -1,1 representing an object j with weight b_j being weighed on the left and right of a chemical balance respectively.     

On Stein’s lemma, dependent covariates and functional monotonicity in multi-dimensional modeling

Tracking the correct directions of monotonicity in multi-dimensional modeling plays an important role in interpreting functional associations. In the presence of multiple predictors, we provide empirical evidence that the observed monotone directions via parametric, nonparametric or semiparametric fit of commonly used multi-dimensional models may entirely violate the actual directions of monotonicity. This breakdown is caused primarily by the dependence structure of covariates, with negligible influence from the bias of function estimation. To examine the linkage between the dependent covariates and monotone directions, we first generalize Stein’s Lemma for random variables which are mutually independent Gaussian to two important cases: dependent Gaussian, and independent non-Gaussian. We show that in both two cases, there is an explicit one-to-one correspondence between the monotone directions of a multi-dimensional function and the signs of a deterministic surrogate vector. Moreover, we demonstrate that the second case can be extended to accommodate a class of dependent covariates. This generalization further enables us to develop a de-correlation transform for arbitrarily dependent covariates. The transformed covariates preserve modeling interpretability with little loss in modeling efficiency. The simplicity and effectiveness of the proposed method are illustrated via simulation studies and real data application.