On the singular gamma, Wishart, and beta matrix‐variate density functions
Arak M. Mathai,Serge B. Provost
DOI: https://doi.org/10.1002/cjs.11710
2022-07-09
Canadian Journal of Statistics
Abstract:When a p×p real positive definite matrix S follows a Wishart or, more generally, a matrix‐variate gamma distribution with shape parameter α and positive definite scale parameter matrix B, one can represent S as XX′ for some matrix X of dimension p×q. When p>q, S has a singular distribution whose properties can be studied via the density function of X. It will be shown that when X follows a matrix‐variate extended Gaussian distribution, the density function of the resulting singular gamma distribution can be obtained by making use of successive transformations and their associated Jacobians. The singular Wishart distribution will then be obtained as a particular case. The marginal and conditional density functions resulting from an arbitrary partitioning of X will be considered as well. The same technique will also be applied to the derivation of the density functions of real and complex singular type‐1 and type‐2 beta‐distributed matrices. It so happens that the proposed approach, which is based on manipulations involving the wedge products of certain differential elements, generally proves more efficient than the intricate procedures that have hitherto been employed in the literature. Résumé Lorsqu'une matrice réelle définie positive S (p×p) suit une loi de Wishart ou, plus généralement, une distribution gamma de paramètre de forme α et de paramètre d'échelle une matrice définie positive B, il est possible de représenter S sous forme XX′ où X est une matrice de dimension p×q. Lorsque p>q, S a une distribution singulière dont les propriétés peuvent être étudiées via la fonction de densité de X. Les auteurs de cet article montrent que lorsque X suit une distribution gaussienne étendue à variable matricielle, la fonction de densité de la distribution gamma singulière résultante peut être obtenue en utilisant des transformations successives et leurs Jacobiens associés. La distribution de Wishart singulière en découle comme un cas particulier. Le cas de fonctions de densité marginale et conditionnelle résultant d'un partitionnement arbitraire de X a également été traité. La même technique a été utilisée pour déterminer les fonctions de densité de matrices singulières réelles et complexes de lois beta de type 1 et de type 2. Ainsi l'approche proposée, basée sur des manipulations mettant en jeu les produits extérieurs de certains éléments différentiels, s'avère généralement plus efficace que les procédures complexes utilisées jusqu'à présent dans la littérature.