Connection between the Hadamard and matrix products with an application to matrix-variate Birnbaum-Saunders distributions

In this paper, we establish a connection between the Hadamard product and the usual matrix multiplication. In addition, we study some new properties of the Hadamard product and explore the inverse problem associated with the established connection, which facilitates diverse applications. Furthermore, we propose a matrix-variate generalized Birnbaum-Saunders (GBS) distribution. Three representations of the matrix-variate GBS density are provided, one of them by using the mentioned connection. The main motivation of this article is based on the fact that the representation of the matrix-variate GBS density based on element-by-element specification does not allow matrix transformations. Consequently, some statistical procedures based on this representation, such as multivariate data analysis and statistical shape theory, cannot be performed. For this reason, the primary goal of this work is to obtain a matrix representation of the matrix-variate GBS density that is useful for some statistical applications. When the GBS density is expressed by means of a matrix representation based on the Hadamard product, such a density is defined in terms of the original matrices, as is common for many matrix-variate distributions, allowing matrix transformations to be handled in a natural way and then suitable statistical procedures to be developed.     

On fundamental skew distributions

A new class of multivariate skew-normal distributions, fundamental skew-normal distributions and their canonical version, is developed. It contains the product of independent univariate skew-normal distributions as a special case. Stochastic representations and other main properties of the associated distribution theory of linear and quadratic forms are considered. A unified procedure for extending this class to other families of skew distributions such as the fundamental skew-symmetric, fundamental skew-elliptical, and fundamental skew-spherical class of distributions is also discussed.     

Characterization of Wishart-Laplace distributions via Jordan algebra homomorphisms

For a real, Hermitian, or quaternion normal random matrix Y with mean zero, necessary and sufficient conditions for a quadratic form Q(Y) to have a Wishart-Laplace distribution (the distribution of the difference of two independent central Wishart W_p(m_i,Σ) random matrices) are given in terms of a certain Jordan algebra homomorphism ρ. Further, it is shown that {Q_k(Y)} is independent Laplace-Wishart if and only if in addition to the aforementioned conditions, the images ρ_k(Σ⁺) of the Moore-Penrose inverse Σ⁺ of Σ are mutually orthogonal: ρ_k(Σ⁺)ρ_?(Σ⁺)=0 for k≠?.     

Distributions of order statistics and linear combinations of order statistics from an elliptical distribution as mixtures of unified skew-elliptical distributions

We consider here the distributions of order statistics and linear combinations of order statistics from an elliptical distribution. We show that these distributions can be expressed as mixtures of unified skew-elliptical distributions, and then use these mixture representations to derive their moment generating functions and moments, when they exist.     

Stein's Lemma for elliptical random vectors

For the family of multivariate normal distribution functions, Stein's Lemma presents a useful tool for calculating covariances between functions of the component random variables. Motivated by applications to corporate finance, we prove a generalization of Stein's Lemma to the family of elliptical distributions.     

On Pickands coordinates in arbitrary dimensions

Pickands coordinates were introduced as a crucial tool for the investigation of bivariate extreme value models. We extend their definition to arbitrary dimensions and, thus, we can generalize many known results for bivariate extreme value and generalized Pareto models to higher dimensions and arbitrary extreme value margins.In particular we characterize multivariate generalized Pareto distributions (GPs) and spectral δ-neighborhoods of GPs in terms of best attainable rates of convergence of extremes, which are well-known results in the univariate case. A sufficient univariate condition for a multivariate distribution function (df) to belong to the domain of attraction of an extreme value df is derived. Bounds for the variational distance in peaks-over-threshold models are established, which are based on Pickands coordinates.     

Hierarchical orbital decompositions and extended decomposable distributions

Elliptically contoured distributions can be considered to be the distributions for which the contours of the density functions are proportional ellipsoids. Kamiya, Takemura and Kuriki [Star-shaped distributions and their generalizations, J. Statist. Plann. Inference, 2006, available at 〈http://arxiv.org/abs/math.ST/0605600〉, to appear] generalized the elliptically contoured distributions to star-shaped distributions, for which the contours are allowed to be arbitrary proportional star-shaped sets. This was achieved by considering the so-called orbital decomposition of the sample space in the general framework of group invariance. In the present paper, we extend their results by conducting the orbital decompositions in steps and obtaining a further, hierarchical decomposition of the sample space. This allows us to construct probability models and distributions with further independence structures. The general results are applied to the star-shaped distributions with a certain symmetric structure, the distributions related to the two-sample Wishart problem and the distributions of preference rankings.     

On the spectrum of resolvable orthogonal arrays invariant under the Klein group K4

It is shown that there exists a resolvablen ² by 4 orthogonal array which is invariant under the Klein 4-groupK ₄ for all positive integersn congruent to 0 modulo 4 except possibly forn {12, 24, 156, 348}.     

Proper bayes minimax estimators for a multivariate normal mean with unknown common variance under a convex loss function

Summary Let X ∼ N_p(μ,σ²I_p) and let s/σ² ∼ χ _n ² , independent ofX, where μ and σ² are unknown. This paper considers the estimation of μ (by δ) relative to a convex loss function given by (δ−μ)′[(1−α)I_p/σ²+αQ](δ−μ)/[(1−α)p/σ²+α tr (Q)], whereQ is a knownp×p diagonal matrix and 0≦α≦1. Two classes of minimax estimators are obtained for μ whenp≦3; the first is a new result and the second is a generalization of a result of Strawderman (1973,Ann. Statist.,1, 1189–1194). A proper Bayes estimator is also obtained which is shown to satisfy the conditions of the second class of minimax estimators. The paper concludes by discussing the estimation of μ relative to another convex loss function. This work was supported by the Army, Navy and Air Force under Office of Naval Research Contract No. N00014-80-C-0093. Reproduction in whole or in part is permitted for any purpose of the United States Government.     

On the chernoff-savage theorem for dependent sequences

Summary Given a sequence of ϕ-mixing random variables not necessarily stationary, a Chernoff-Savage theorem for two-sample linear rank statistics is proved using the Pyke-Shorack [5] approach based on weak convergence properties of empirical processes in an extended metric. This result is a generalization of Fears and Mehra [4] in that the stationarity is not required and that the condition imposed on the mixing numbers is substantially relaxed. A similar result is shown to hold for strong mixing sequences under slightly stronger conditions on the mixing numbers. Research partially supported by the National Research Council of Canada under Grant No. A-3954.     

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.     

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.     

Multivariate skewness and kurtosis measures with an application in ICA

In this paper skewness and kurtosis characteristics of a multivariate p-dimensional distribution are introduced. The skewness measure is defined as a p-vector while the kurtosis is characterized by a p×p-matrix. The introduced notions are extensions of the corresponding measures of Mardia [K.V. Mardia, Measures of multivariate skewness and kurtosis with applications, Biometrika 57 (1970) 519–530] and Móri, Rohatgi & Székely [T.F. Móri, V.K. Rohatgi, G.J. Székely, On multivariate skewness and kurtosis, Theory Probab. Appl. 38 (1993) 547–551]. Basic properties of the characteristics are examined and compared with both the above-mentioned results in the literature. Expressions for the measures of skewness and kurtosis are derived for the multivariate Laplace distribution. The kurtosis matrix is used in Independent Component Analysis (ICA) where the solution of an eigenvalue problem of the kurtosis matrix determines the transformation matrix of interest [A. Hyvärinen, J. Karhunen, E. Oja, Independent Component Analysis, Wiley, New York, 2001].     

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.     

On Riesz and Wishart distributions associated with decomposable undirected graphs

Classical Wishart distributions on the open convex cone of positive definite matrices and their fundamental features are extended to generalized Riesz and Wishart distributions associated with decomposable undirected graphs using the basic theory of exponential families. The families of these distributions are parameterized by their expectations/natural parameter and multivariate shape parameter and have a non-trivial overlap with the generalized Wishart distributions defined in Andersson and Wojnar (2004) [4] and [8]. This work also extends the Wishart distributions of type I in Letac and Massam (2007) [7] and, more importantly, presents an alternative point of view on the latter paper.     

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.