Probabilities in multimatrix variate distributions: an application in SARS-CoV-2

Francisco J. Caro-Lopera [email protected] José A. Díaz-García [email protected] University of Medellin, Faculty of Basic Sciences, Carrera 87 No.30-65, Medellin, Colombia Universidad Autónoma de Chihuahua, Facultad de Zootecnia y Ecología, Periférico Francisco R. Almada Km 1, Zootecnia, 33820 Chihuahua, Chihuahua, México.
Abstract

Recently the termed multimatrix variate distributions were proposed in Díaz-García and Caro-Lopera (2024a) as an alternative for univariate and vector variate copulas. The distributions are based on sample probabilistic dependent elliptically countered models and most of them are also invariant under this family of laws. Despite a large of results on matrix variate distributions since the last 70 years, the spherical multimatrix distributions and the associated probabilities on hyper cones can be computable. The multiple probabilities are set in terms of recurrent integrations allowing several matrix computation a feasible task. An application of the emerging probabilities is placed into a dynamic molecular docking in the SARS-CoV-2 main protease. Finally, integration over multimatrix Wishart distribution provides a simplification of a complex kernel integral in elliptical models under real normed division algebras and the solution was applied in elliptical affine shape theory.

keywords:
Multimatrix variate distributions, matrix variate elliptical distributions, probabilities on cones, zonal polynomials, real normed division algebras, elliptical affine shape theory.
MSC:
60E05, 62E15, 15A23, 15B52

1 Introduction

The theory of matrix variate distributions has been developed profusely in the last century. It has emerged as a natural extension of the univariate case, covering all statistical and probabilistic disciplines.

However, the new extensions are becoming more complicated, giving rise to complex calculations that remained impossible for decades. Other generalizations took longer, specifically those related to joint distributions of dependent samples in non-Gaussian models and had usually been addressed by theories such as copulas among others.

Each new extension of the univariate or vector case brings with it not only computability, but increasingly demanding applications in support distributions, estimation and dependence. In response to such goal, alternatives to copula theory have appeared, which not only extend the results of the univariate and vector cases, but also consider all real normed division algebras (real, complex, quaternion, octonion).

Depending on the type of function that describes the relationships of the multiple matrices, two types of distributions have recently been proposed applied to probabilistically dependent samples and under robust elliptically countered models. These are called multimatricvariate and multimatrix variate distributions and the seminal works appear in Díaz-García and Caro-Lopera (2022), Díaz-García and Caro-Lopera (2024a) and Díaz-García and Caro-Lopera (2024b). The multivector variate case can be seen in Díaz-García et al. (2022). The importance of dependent samples, robust distributions beyond Gaussian, general calculations for all real normed division algebras and estimation based on non-independent likelihoods have been studied in detail in the aforementioned articles.

In this work we focus on the computational problem of distributions and on the calculation of multimatrix probabilities.

A multimatrix theory that solves the addressed previous problems but that cannot be computed would add to the hundreds of theoretical results in matrix analysis developed in the last seven decades. Since that time, the advent of the theory of integration on orthogonal groups due to A.T. James, allowed the calculation of the joint distribution of the latent roots of any function of a positive definite matrix, among many others applications, see for example James (1960). The solution created a set of orthogonal polynomials for the expansion of the power of the trace of a matrix and with it the central matrix variate statistics on Gaussian or elliptic models could be built in all real normed division algebras. Only a few years ago the results involving hypergeometric type series of zonal polynomials could be approximated numerically, so it is quite a challenge for new theories that can be expressed in computable distributions. The central and isotropic cases have such an opportunity, but the non-central case, useful for solving, for example, the distribution of the latent roots of a Wishart matrix, came about with the creation in 1979 of an extension of the James zonal polynomials and which A.W.Davis called invariant polynomials of several matrix arguments (Davis (1979)). The success of the calculation of zonal polynomials lay in their recurring construction using the Laplace-Beltrami operator (James (1968)), so Davis conjectured in 1979 that his polynomials could have such a construction. Unfortunately, in 2016 the impossibility of the method was demonstrated, leaving open the calculation problem for results involving Davis polynomials (Caro-Lopera (2016)). In this context, and for now, only results involving at most zonal polynomials can be computed and be useful for applications. A natural challenge for a multimatrix theory is that the joint distribution can be integrated to compute matrix event probability problems. This aspect has had little treatment in the literature and very few explicit calculations of probabilities in cones are known. The problem of calculating multiple probabilities in a joint distribution is precisely what we also hope to solve, in such a way that at most they are expressed in computable series of zonal polynomials.

The above discussion is placed in the present article by computing some probabilities in multimatrix variate distributions. Some preliminary integrals of interest as well as the definition of the matrix variate elliptical contoured distributions and generalized series of zonal polynomials are collected in Section 2. Then, Section 3 revises the problem of computation in matrix variate distribution theory and provides some results on multiple probabilities on cones written in terms of computable series of zonal polynomials. Some of the probabilities are invariant for the complete family of elliptical distributions, meanwhile the non invariant distributions promoted the simplification of kernel elliptical integrals indexed by general derivatives; that result is applied in affine shape theory under real normed division algebras. Section 4 is based on the dependent joint distribution of a sample derived from a molecular docking in a new cavity of SARS-CoV-2 main protease. Then a simulation of a dynamic molecular docking is set in terms of the probabilities arising from a deformation of the ligand into the corresponding binding site.

2 Preliminary results

Matrix notations, matrix variate elliptical contoured distributions, zonal polynomials and some multimatrix variate distributions are presented in this section. First, we start with notations and terminologies. 𝐀n×m𝐀superscript𝑛𝑚\mathbf{A}\in\Re^{n\times m}bold_A ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT denotes a matrix with n𝑛nitalic_n rows and m𝑚mitalic_m columns; 𝐀m×nsuperscript𝐀superscript𝑚𝑛\mathbf{A}^{\prime}\in\Re^{m\times n}bold_A start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_m × italic_n end_POSTSUPERSCRIPT is the transpose matrix, and if 𝐀n×n𝐀superscript𝑛𝑛\mathbf{A}\in\Re^{n\times n}bold_A ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × italic_n end_POSTSUPERSCRIPT has an inverse, it will be denoted by 𝐀1n×nsuperscript𝐀1superscript𝑛𝑛\mathbf{A}^{-1}\in\Re^{n\times n}bold_A start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × italic_n end_POSTSUPERSCRIPT. 𝐀n×n𝐀superscript𝑛𝑛\mathbf{A}\in\Re^{n\times n}bold_A ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × italic_n end_POSTSUPERSCRIPT is a symmetric matrix if 𝐀=𝐀𝐀superscript𝐀\mathbf{A}=\mathbf{A}^{\prime}bold_A = bold_A start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT. If all their eigenvalues are positive then 𝐀𝐀\mathbf{A}bold_A is a positive definite matrix, a fact denoted as 𝐀>𝟎𝐀0\mathbf{A}>\mathbf{0}bold_A > bold_0. An identity matrix will be denoted by 𝐈n×n𝐈superscript𝑛𝑛\mathbf{I}\in\Re^{n\times n}bold_I ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × italic_n end_POSTSUPERSCRIPT. To specify the size of the identity, we will use 𝐈nsubscript𝐈𝑛\mathbf{I}_{n}bold_I start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT. tr(𝐀)tr𝐀\mathop{\rm tr}\nolimits(\mathbf{A})roman_tr ( bold_A ) denotes the trace of matrix 𝐀m×m𝐀superscript𝑚𝑚\mathbf{A}\in\Re^{m\times m}bold_A ∈ roman_ℜ start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT. If 𝐀n×m𝐀superscript𝑛𝑚\mathbf{A}\in\Re^{n\times m}bold_A ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT then by vec(𝐀)vec𝐀\mathop{\rm vec}\nolimits(\mathbf{A})roman_vec ( bold_A ) we mean the mn×1𝑚𝑛1mn\times 1italic_m italic_n × 1 vector formed by stacking the columns of 𝐀𝐀\mathbf{A}bold_A under each other; that is, if 𝐀=[𝐚1𝐚2𝐚m]𝐀delimited-[]subscript𝐚1subscript𝐚2subscript𝐚𝑚\mathbf{A}=[\mathbf{a}_{1}\mathbf{a}_{2}\dots\mathbf{a}_{m}]bold_A = [ bold_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT bold_a start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT … bold_a start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ], where 𝐚jn×1subscript𝐚𝑗superscript𝑛1\mathbf{a}_{j}\in\Re^{n\times 1}bold_a start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × 1 end_POSTSUPERSCRIPT for j=1,2,,m𝑗12𝑚j=1,2,\dots,mitalic_j = 1 , 2 , … , italic_m, vec(𝐀)=[𝐚1𝐚2𝐚m].vec𝐀superscriptdelimited-[]superscriptsubscript𝐚1superscriptsubscript𝐚2superscriptsubscript𝐚𝑚\mathop{\rm vec}\nolimits(\mathbf{A})=\left[\mathbf{a}_{1}^{\prime}\\ \mathbf{a}_{2}^{\prime}\\ \cdots\\ \mathbf{a}_{m}^{\prime}\right]^{\prime}.roman_vec ( bold_A ) = [ bold_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT bold_a start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ⋯ bold_a start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT . The Frobenius norm of a matrix 𝐀𝐀\mathbf{A}bold_A will be denoted as 𝐀norm𝐀||\mathbf{A}||| | bold_A | |. Typically the Frobenius norm is denoted as 𝐀Fsubscriptnorm𝐀𝐹||\mathbf{A}||_{F}| | bold_A | | start_POSTSUBSCRIPT italic_F end_POSTSUBSCRIPT, to differentiate it from other matrix norms. Since we will use only the Frobenius norm, it just be denoted as 𝐀norm𝐀||\mathbf{A}||| | bold_A | |. It is defined by

𝐀=tr(𝐀𝐀)=vec(𝐀)vec(𝐀).norm𝐀trsuperscript𝐀𝐀superscriptvec𝐀vec𝐀||\mathbf{A}||=\sqrt{\mathop{\rm tr}\nolimits(\mathbf{A}^{\prime}\mathbf{A})}=% \sqrt{\mathop{\rm vec}\nolimits^{\prime}(\mathbf{A})\mathop{\rm vec}\nolimits(% \mathbf{A})}.| | bold_A | | = square-root start_ARG roman_tr ( bold_A start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT bold_A ) end_ARG = square-root start_ARG roman_vec start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( bold_A ) roman_vec ( bold_A ) end_ARG .

Finally, 𝒱m,nsubscript𝒱𝑚𝑛\mathcal{V}_{m,n}caligraphic_V start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT denotes the Stiefel manifold, the space of all matrices 𝐇1n×msubscript𝐇1superscript𝑛𝑚\mathbf{H}_{1}\in\Re^{n\times m}bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT (nm𝑛𝑚n\geq mitalic_n ≥ italic_m) with orthogonal columns, that is, 𝒱m,n={𝐇1n×m;𝐇1𝐇1=𝐈m}subscript𝒱𝑚𝑛formulae-sequencesubscript𝐇1superscript𝑛𝑚subscriptsuperscript𝐇1subscript𝐇1subscript𝐈𝑚\mathcal{V}_{m,n}=\{\mathbf{H}_{1}\in\Re^{n\times m};\mathbf{H}^{\prime}_{1}% \mathbf{H}_{1}=\mathbf{I}_{m}\}caligraphic_V start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT = { bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT ; bold_H start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = bold_I start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT }. In addition, if (𝐇1d𝐇1)subscriptsuperscript𝐇1𝑑subscript𝐇1(\mathbf{H}^{\prime}_{1}d\mathbf{H}_{1})( bold_H start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_d bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) defines an invariant measure on the Stiefel manifold 𝒱m,nsubscript𝒱𝑚𝑛\mathcal{V}_{m,n}caligraphic_V start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT, from Theorem 2.1.15, p. 70 in Muirhead (2005),

𝒱m,n(𝐇1d𝐇1)=2mπmn/2Γm[n/2].subscriptsubscript𝒱𝑚𝑛subscriptsuperscript𝐇1𝑑subscript𝐇1superscript2𝑚superscript𝜋𝑚𝑛2subscriptΓ𝑚delimited-[]𝑛2\int_{\mathcal{V}_{m,n}}(\mathbf{H}^{\prime}_{1}d\mathbf{H}_{1})=\frac{2^{m}% \pi^{mn/2}}{\Gamma_{m}[n/2]}.∫ start_POSTSUBSCRIPT caligraphic_V start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_H start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_d bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = divide start_ARG 2 start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_π start_POSTSUPERSCRIPT italic_m italic_n / 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_n / 2 ] end_ARG . (1)

where Γm[a]subscriptΓ𝑚delimited-[]𝑎\Gamma_{m}[a]roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a ] denotes the multivariate Gamma function, see Muirhead (2005, Definiton 2.1.10, p. 61)

Now, let 𝐕N×m𝐕superscript𝑁𝑚\mathbf{V}\in\Re^{N\times m}bold_V ∈ roman_ℜ start_POSTSUPERSCRIPT italic_N × italic_m end_POSTSUPERSCRIPT random matrix with a matrix variate elliptical distribution with respect to the Lebesgue measure (d𝐕)𝑑𝐕(d\mathbf{V})( italic_d bold_V ), see Gupta and Varga (1993). Therefore its density function is given by

dF𝐕(𝐕)=|𝚺|N/2|𝚯|m/2h{tr[(𝐕𝝁)T𝚺1(𝐕𝝁)𝚯1]}(d𝐕).𝑑subscript𝐹𝐕𝐕superscript𝚺𝑁2superscript𝚯𝑚2trdelimited-[]superscript𝐕𝝁𝑇superscript𝚺1𝐕𝝁superscript𝚯1𝑑𝐕dF_{{}_{\mathbf{V}}}(\mathbf{V})=|\mathbf{\Sigma}|^{-N/2}|\mathbf{\Theta}|^{-m% /2}h\left\{\mathop{\rm tr}\nolimits\left[(\mathbf{V}-\boldsymbol{\mu})^{T}% \mathbf{\Sigma}^{-1}(\mathbf{V}-\boldsymbol{\mu})\mathbf{\Theta}^{-1}\right]% \right\}(d\mathbf{V}).italic_d italic_F start_POSTSUBSCRIPT start_FLOATSUBSCRIPT bold_V end_FLOATSUBSCRIPT end_POSTSUBSCRIPT ( bold_V ) = | bold_Σ | start_POSTSUPERSCRIPT - italic_N / 2 end_POSTSUPERSCRIPT | bold_Θ | start_POSTSUPERSCRIPT - italic_m / 2 end_POSTSUPERSCRIPT italic_h { roman_tr [ ( bold_V - bold_italic_μ ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT bold_Σ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( bold_V - bold_italic_μ ) bold_Θ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ] } ( italic_d bold_V ) . (2)

The location parameter is 𝝁N×m𝝁superscript𝑁𝑚\boldsymbol{\mu}\in\Re^{N\times m}bold_italic_μ ∈ roman_ℜ start_POSTSUPERSCRIPT italic_N × italic_m end_POSTSUPERSCRIPT; and the scale parameters 𝚺N×N𝚺superscript𝑁𝑁\mathbf{\Sigma}\in\Re^{N\times N}bold_Σ ∈ roman_ℜ start_POSTSUPERSCRIPT italic_N × italic_N end_POSTSUPERSCRIPT and 𝚯m×m𝚯superscript𝑚𝑚\mathbf{\Theta}\in\Re^{m\times m}bold_Θ ∈ roman_ℜ start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT, are positive definite matrices. The distribution shall be denoted by 𝐕N×m(𝝁,𝚺,𝚯;h)similar-to𝐕subscript𝑁𝑚𝝁𝚺𝚯\mathbf{V}\sim\mathcal{E}_{N\times m}(\boldsymbol{\mu},\mathbf{\Sigma},\mathbf% {\Theta};h)bold_V ∼ caligraphic_E start_POSTSUBSCRIPT italic_N × italic_m end_POSTSUBSCRIPT ( bold_italic_μ , bold_Σ , bold_Θ ; italic_h ), and indexed by the kernel function h[0,)0h\mbox{: }\Re\to[0,\infty)italic_h : roman_ℜ → [ 0 , ∞ ), where 0uNm/21h(u)𝑑u<superscriptsubscript0superscript𝑢𝑁𝑚21𝑢differential-d𝑢\int_{0}^{\infty}u^{Nm/2-1}h(u)du<\infty∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_u start_POSTSUPERSCRIPT italic_N italic_m / 2 - 1 end_POSTSUPERSCRIPT italic_h ( italic_u ) italic_d italic_u < ∞.

When 𝝁=𝟎𝝁0\boldsymbol{\mu}=\mathbf{0}bold_italic_μ = bold_0, 𝚺=𝐈N𝚺subscript𝐈𝑁\mathbf{\Sigma}=\mathbf{I}_{N}bold_Σ = bold_I start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT and 𝚯=𝐈m𝚯subscript𝐈𝑚\mathbf{\Theta}=\mathbf{I}_{m}bold_Θ = bold_I start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT a special case of a matrix variate elliptical distribution appears, in this case it is said that 𝐕𝐕\mathbf{V}bold_V has a matrix variate spherical distribution.

Note that for constant a𝑎a\in\Reitalic_a ∈ roman_ℜ, the substitution v=u/a𝑣𝑢𝑎v=u/aitalic_v = italic_u / italic_a and (du)=a(dv)𝑑𝑢𝑎𝑑𝑣(du)=a(dv)( italic_d italic_u ) = italic_a ( italic_d italic_v ) in Fang et al. (1990, Equation 2.21, p. 26) provides that

v>0vNm/21h(av)(dv)=aNm/2Γ1[Nm/2]πNm/2subscript𝑣0superscript𝑣𝑁𝑚21𝑎𝑣𝑑𝑣superscript𝑎𝑁𝑚2subscriptΓ1delimited-[]𝑁𝑚2superscript𝜋𝑁𝑚2\int_{v>0}v^{Nm/2-1}h(av)(dv)=\frac{a^{-Nm/2}\Gamma_{1}[Nm/2]}{\pi^{Nm/2}}∫ start_POSTSUBSCRIPT italic_v > 0 end_POSTSUBSCRIPT italic_v start_POSTSUPERSCRIPT italic_N italic_m / 2 - 1 end_POSTSUPERSCRIPT italic_h ( italic_a italic_v ) ( italic_d italic_v ) = divide start_ARG italic_a start_POSTSUPERSCRIPT - italic_N italic_m / 2 end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_N italic_m / 2 ] end_ARG start_ARG italic_π start_POSTSUPERSCRIPT italic_N italic_m / 2 end_POSTSUPERSCRIPT end_ARG (3)

For some applications we will require James zonal polynomials, notated by Cρ(𝐀)subscript𝐶𝜌𝐀C_{\rho}(\mathbf{A})italic_C start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( bold_A ), see Muirhead (2005). They are an orthogonal base expansion in (tr𝐀)r=ρrCρ(𝐀)superscripttr𝐀𝑟subscript𝜌𝑟subscript𝐶𝜌𝐀(\mathop{\rm tr}\nolimits\mathbf{A})^{r}=\sum_{\rho\in r}C_{\rho}(\mathbf{A})( roman_tr bold_A ) start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT = ∑ start_POSTSUBSCRIPT italic_ρ ∈ italic_r end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( bold_A ). Here 𝐀m×m𝐀superscript𝑚𝑚\mathbf{A}\in\Re^{m\times m}bold_A ∈ roman_ℜ start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT is positive definite and the summation is over all ordered partitions ρ=(r1,,rm)𝜌subscript𝑟1subscript𝑟𝑚\rho=(r_{1},...,r_{m})italic_ρ = ( italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_r start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ), into not more than m𝑚mitalic_m parts, such that r1r2rm0subscript𝑟1subscript𝑟2subscript𝑟𝑚0r_{1}\geq r_{2}\geq\cdots\geq r_{m}\geq 0italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≥ italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ≥ ⋯ ≥ italic_r start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ≥ 0. Series of zonal polynomials in this work will required the following expression given by Caro-Lopera et al. (2010) and Caro-Lopera and Díaz-García (2012):

Pqpr[f(r,𝐗):a1,,ap;b1,,bq;𝐗]=r=0f(r,𝐗)r!ρr(a1)ρ(ap)ρ(b1)ρ(bq)ρCρ(𝐗){}_{p}^{r}P_{q}[f(r,\mathbf{X}):a_{1},\ldots,a_{p};b_{1},\ldots,b_{q};\mathbf{% X}]=\sum_{r=0}^{\infty}\frac{f(r,\mathbf{X})}{r!}\sum_{\rho\in r}\frac{(a_{1})% _{\rho}\cdots(a_{p})_{\rho}}{(b_{1})_{\rho}\cdots(b_{q})_{\rho}}C_{\rho}(% \mathbf{X})start_FLOATSUBSCRIPT italic_p end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ italic_f ( italic_r , bold_X ) : italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_a start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ; bold_X ] = ∑ start_POSTSUBSCRIPT italic_r = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG italic_f ( italic_r , bold_X ) end_ARG start_ARG italic_r ! end_ARG ∑ start_POSTSUBSCRIPT italic_ρ ∈ italic_r end_POSTSUBSCRIPT divide start_ARG ( italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ⋯ ( italic_a start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT end_ARG start_ARG ( italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ⋯ ( italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT end_ARG italic_C start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( bold_X ) (4)

where the function f(r,𝐗)𝑓𝑟𝐗f(r,\mathbf{X})italic_f ( italic_r , bold_X ) is independent of ρ𝜌\rhoitalic_ρ and (w)ρ=i=1m(w12(i1))risubscript𝑤𝜌superscriptsubscriptproduct𝑖1𝑚subscript𝑤12𝑖1subscript𝑟𝑖(w)_{\rho}=\prod_{i=1}^{m}(w-\frac{1}{2}(i-1))_{r_{i}}( italic_w ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT = ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( italic_w - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( italic_i - 1 ) ) start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT, with (w)0=1subscript𝑤01(w)_{0}=1( italic_w ) start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1, (w)t=w(w+1)(w+t1)subscript𝑤𝑡𝑤𝑤1𝑤𝑡1(w)_{t}=w(w+1)\cdots(w+t-1)( italic_w ) start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = italic_w ( italic_w + 1 ) ⋯ ( italic_w + italic_t - 1 ) and the aissuperscriptsubscript𝑎𝑖𝑠a_{i}^{\prime}sitalic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_s and bjssuperscriptsubscript𝑏𝑗𝑠b_{j}^{\prime}sitalic_b start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_s are complex numbers. If f(r,𝐗)=1𝑓𝑟𝐗1f(r,\mathbf{X})=1italic_f ( italic_r , bold_X ) = 1, then the well known hypergeometric series Fqp(a1,,ap;b1,,bq;𝐗)subscriptsubscript𝐹𝑞𝑝subscript𝑎1subscript𝑎𝑝subscript𝑏1subscript𝑏𝑞𝐗{}_{p}F_{q}(a_{1},\ldots,a_{p};b_{1},\ldots,b_{q};\mathbf{X})start_FLOATSUBSCRIPT italic_p end_FLOATSUBSCRIPT italic_F start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ( italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_a start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ; bold_X ) is obtained, see for example Muirhead (2005). We also introduce the notation

Qqpr[a1,,ap;b1,,bq;𝐗]=ρr(a1)ρ(ap)ρ(b1)ρ(bq)ρCρ(𝐗)subscriptsuperscriptsubscript𝑄𝑞𝑟𝑝subscript𝑎1subscript𝑎𝑝subscript𝑏1subscript𝑏𝑞𝐗subscript𝜌𝑟subscriptsubscript𝑎1𝜌subscriptsubscript𝑎𝑝𝜌subscriptsubscript𝑏1𝜌subscriptsubscript𝑏𝑞𝜌subscript𝐶𝜌𝐗{}_{p}^{r}Q_{q}[a_{1},\ldots,a_{p};b_{1},\ldots,b_{q};\mathbf{X}]=\sum_{\rho% \in r}\frac{(a_{1})_{\rho}\cdots(a_{p})_{\rho}}{(b_{1})_{\rho}\cdots(b_{q})_{% \rho}}C_{\rho}(\mathbf{X})start_FLOATSUBSCRIPT italic_p end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT [ italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_a start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ; bold_X ] = ∑ start_POSTSUBSCRIPT italic_ρ ∈ italic_r end_POSTSUBSCRIPT divide start_ARG ( italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ⋯ ( italic_a start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT end_ARG start_ARG ( italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ⋯ ( italic_b start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT end_ARG italic_C start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( bold_X ) (5)

Finally, we focus on some multimatrix variate distributions derived by Díaz-García and Caro-Lopera (2024a) which are invariant under family of elliptically countered distributions. i.e. they are not dependent of the kernel function; a preferable property for applications that avoids any prior knowledge of the underlying distribution.

Lemma 1

Assume that 𝐗=(𝐗0,,𝐗k)𝐗superscriptsubscriptsuperscript𝐗0subscriptsuperscript𝐗𝑘\mathbf{X}=\left(\mathbf{X}^{\prime}_{0},\dots,\mathbf{X}^{\prime}_{k}\right)^% {\prime}bold_X = ( bold_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , … , bold_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT has a matrix variate spherical distribution, with 𝐗ini×msubscript𝐗𝑖superscriptsubscript𝑛𝑖𝑚\mathbf{X}_{i}\in\Re^{n_{i}\times m}bold_X start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT × italic_m end_POSTSUPERSCRIPT, nimsubscript𝑛𝑖𝑚n_{i}\geq mitalic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≥ italic_m, i=0,1,,k𝑖01𝑘i=0,1,\dots,kitalic_i = 0 , 1 , … , italic_k. Define V=𝐗02𝑉superscriptnormsubscript𝐗02V=||\mathbf{X}_{0}||^{2}italic_V = | | bold_X start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and 𝐓i=V1/2𝐗isubscript𝐓𝑖superscript𝑉12subscript𝐗𝑖\mathbf{T}_{i}=V^{-1/2}\mathbf{X}_{i}bold_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_V start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT bold_X start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, i=1,,k𝑖1𝑘i=1,\dots,kitalic_i = 1 , … , italic_k. The termed multimatrix variate Pearson type VII is the marginal density dF𝐓1,,𝐓k(𝐓1,,𝐓k)𝑑subscript𝐹subscript𝐓1subscript𝐓𝑘subscript𝐓1subscript𝐓𝑘dF_{\mathbf{T}_{1},\dots,\mathbf{T}_{k}}(\mathbf{T}_{1},\dots,\mathbf{T}_{k})italic_d italic_F start_POSTSUBSCRIPT bold_T start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_T start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_T start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_T start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) of 𝐓1,,𝐓ksubscript𝐓1subscript𝐓𝑘\mathbf{T}_{1},\dots,\mathbf{T}_{k}bold_T start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_T start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT and is given by

Γ1[Nm/2]π(Nn0)m/2Γ1[n0m/2](1+i=1k𝐓i2)Nm/2i=1k(d𝐓i),subscriptΓ1delimited-[]𝑁𝑚2superscript𝜋𝑁subscript𝑛0𝑚2subscriptΓ1delimited-[]subscript𝑛0𝑚2superscript1superscriptsubscript𝑖1𝑘superscriptnormsubscript𝐓𝑖2𝑁𝑚2superscriptsubscript𝑖1𝑘𝑑subscript𝐓𝑖\frac{\Gamma_{1}[Nm/2]}{\pi^{(N-n_{0})m/2}\Gamma_{1}[n_{0}m/2]}\left(1+% \displaystyle\sum_{i=1}^{k}||\mathbf{T}_{i}||^{2}\right)^{-Nm/2}\bigwedge_{i=1% }^{k}\left(d\mathbf{T}_{i}\right),divide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_N italic_m / 2 ] end_ARG start_ARG italic_π start_POSTSUPERSCRIPT ( italic_N - italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_m / 2 end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 ] end_ARG ( 1 + ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT | | bold_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - italic_N italic_m / 2 end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_d bold_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , (6)

where 𝐓ini×msubscript𝐓𝑖superscriptsubscript𝑛𝑖𝑚\mathbf{T}_{i}\in\Re^{n_{i}\times m}bold_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT × italic_m end_POSTSUPERSCRIPT, nimsubscript𝑛𝑖𝑚n_{i}\geq mitalic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≥ italic_m and N=n0+n1++nk𝑁subscript𝑛0subscript𝑛1subscript𝑛𝑘N=n_{0}+n_{1}+\cdots+n_{k}italic_N = italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + ⋯ + italic_n start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT. Now define 𝐅i=𝐓i𝐓i>𝟎subscript𝐅𝑖subscriptsuperscript𝐓𝑖subscript𝐓𝑖0\mathbf{F}_{i}=\mathbf{T}^{\prime}_{i}\mathbf{T}_{i}>\mathbf{0}bold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = bold_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT bold_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT > bold_0, i=1,,k𝑖1𝑘i=1,\dots,kitalic_i = 1 , … , italic_k. The multimatrix variate beta type II distribution dF𝐅1,,𝐅k(𝐅1,,𝐅k)𝑑subscript𝐹subscript𝐅1subscript𝐅𝑘subscript𝐅1subscript𝐅𝑘dF_{\mathbf{F}_{1},\dots,\mathbf{F}_{k}}(\mathbf{F}_{1},\dots,\mathbf{F}_{k})italic_d italic_F start_POSTSUBSCRIPT bold_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is given by

Γ1[Nm/2]Γ1[n0m/2]i=1kΓm[ni/2]i=1k|𝐅i|(nim1)/2(1+i=1ktr𝐅i)Nm/2i=1k(d𝐅i).subscriptΓ1delimited-[]𝑁𝑚2subscriptΓ1delimited-[]subscript𝑛0𝑚2superscriptsubscriptproduct𝑖1𝑘subscriptΓ𝑚delimited-[]subscript𝑛𝑖2superscriptsubscriptproduct𝑖1𝑘superscriptsubscript𝐅𝑖subscript𝑛𝑖𝑚12superscript1superscriptsubscript𝑖1𝑘trsubscript𝐅𝑖𝑁𝑚2superscriptsubscript𝑖1𝑘𝑑subscript𝐅𝑖\frac{\Gamma_{1}[Nm/2]}{\Gamma_{1}[n_{0}m/2]\displaystyle\prod_{i=1}^{k}\Gamma% _{m}[n_{i}/2]}\prod_{i=1}^{k}|\mathbf{F}_{i}|^{(n_{i}-m-1)/2}\left(1+% \displaystyle\sum_{i=1}^{k}\mathop{\rm tr}\nolimits\mathbf{F}_{i}\right)^{-Nm/% 2}\bigwedge_{i=1}^{k}\left(d\mathbf{F}_{i}\right).divide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_N italic_m / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 ] ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 ] end_ARG ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT | bold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_m - 1 ) / 2 end_POSTSUPERSCRIPT ( 1 + ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_tr bold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - italic_N italic_m / 2 end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_d bold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) . (7)
Lemma 2

Suppose that 𝐗=(𝐗0,,𝐗k)𝐗superscriptsubscriptsuperscript𝐗0subscriptsuperscript𝐗𝑘\mathbf{X}=\left(\mathbf{X}^{\prime}_{0},\dots,\mathbf{X}^{\prime}_{k}\right)^% {\prime}bold_X = ( bold_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , … , bold_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT has a matrix variate spherical distribution, with 𝐗ini×msubscript𝐗𝑖superscriptsubscript𝑛𝑖𝑚\mathbf{X}_{i}\in\Re^{n_{i}\times m}bold_X start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT × italic_m end_POSTSUPERSCRIPT, nimsubscript𝑛𝑖𝑚n_{i}\geq mitalic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≥ italic_m, i=0,1,,k𝑖01𝑘i=0,1,\dots,kitalic_i = 0 , 1 , … , italic_k. Define V=𝐗02𝑉superscriptnormsubscript𝐗02V=||\mathbf{X}_{0}||^{2}italic_V = | | bold_X start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and 𝐑i=(V+𝐗i2)1/2𝐗isubscript𝐑𝑖superscript𝑉superscriptnormsubscript𝐗𝑖212subscript𝐗𝑖\mathbf{R}_{i}=\left(V+||\mathbf{X}_{i}||^{2}\right)^{-1/2}\mathbf{X}_{i}bold_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( italic_V + | | bold_X start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT bold_X start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, i=1,,k𝑖1𝑘i=1,\dots,kitalic_i = 1 , … , italic_k. The multimatrix variate Pearson type II distribution dF𝐑1,,𝐑k(𝐑1,,𝐑k)𝑑subscript𝐹subscript𝐑1subscript𝐑𝑘subscript𝐑1subscript𝐑𝑘dF_{\mathbf{R}_{1},\dots,\mathbf{R}_{k}}(\mathbf{R}_{1},\dots,\mathbf{R}_{k})italic_d italic_F start_POSTSUBSCRIPT bold_R start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_R start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_R start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_R start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is

Γ1[Nm/2]π(Nn0)m/2Γ1[n0m/2][1+i=1k𝐑i2(1𝐑i2)]Nm/2subscriptΓ1delimited-[]𝑁𝑚2superscript𝜋𝑁subscript𝑛0𝑚2subscriptΓ1delimited-[]subscript𝑛0𝑚2superscriptdelimited-[]1superscriptsubscript𝑖1𝑘superscriptnormsubscript𝐑𝑖21superscriptnormsubscript𝐑𝑖2𝑁𝑚2\frac{\Gamma_{1}[Nm/2]}{\pi^{(N-n_{0})m/2}\Gamma_{1}[n_{0}m/2]}\left[1+% \displaystyle\sum_{i=1}^{k}\frac{||\mathbf{R}_{i}||^{2}}{\left(1-||\mathbf{R}_% {i}||^{2}\right)}\right]^{-Nm/2}\hskip 56.9055ptdivide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_N italic_m / 2 ] end_ARG start_ARG italic_π start_POSTSUPERSCRIPT ( italic_N - italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_m / 2 end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 ] end_ARG [ 1 + ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT divide start_ARG | | bold_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ( 1 - | | bold_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG ] start_POSTSUPERSCRIPT - italic_N italic_m / 2 end_POSTSUPERSCRIPT
×i=ik(1||𝐑i||2)nim/21i=1k(d𝐑i),\hskip 85.35826pt\times\prod_{i=i}^{k}\left(1-||\mathbf{R}_{i}||^{2}\right)^{-% n_{i}m/2-1}\bigwedge_{i=1}^{k}\left(d\mathbf{R}_{i}\right),× ∏ start_POSTSUBSCRIPT italic_i = italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( 1 - | | bold_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_m / 2 - 1 end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_d bold_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , (8)

Moreover, Assuming that 𝐁i=𝐑i𝐑i>𝟎subscript𝐁𝑖subscriptsuperscript𝐑𝑖subscript𝐑𝑖0\mathbf{B}_{i}=\mathbf{R}^{\prime}_{i}\mathbf{R}_{i}>\mathbf{0}bold_B start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = bold_R start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT bold_R start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT > bold_0 and tr𝐁i1trsubscript𝐁𝑖1\mathop{\rm tr}\nolimits\mathbf{B}_{i}\leq 1roman_tr bold_B start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≤ 1 with i=1,,k𝑖1𝑘i=1,\dots,kitalic_i = 1 , … , italic_k, the multimatrix variate beta type I distribution dF𝐁1,,𝐁k(𝐁1,,𝐁k)𝑑subscript𝐹subscript𝐁1subscript𝐁𝑘subscript𝐁1subscript𝐁𝑘dF_{\mathbf{B}_{1},\dots,\mathbf{B}_{k}}(\mathbf{B}_{1},\dots,\mathbf{B}_{k})italic_d italic_F start_POSTSUBSCRIPT bold_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_B start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_B start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is written as

Γ1[Nm/2]Γ1[n0m/2]i=1k(|𝐁i|(nim1)/2Γm[ni/2])[1+i=1ktr𝐁i(1tr𝐁i)]Nm/2subscriptΓ1delimited-[]𝑁𝑚2subscriptΓ1delimited-[]subscript𝑛0𝑚2superscriptsubscriptproduct𝑖1𝑘superscriptsubscript𝐁𝑖subscript𝑛𝑖𝑚12subscriptΓ𝑚delimited-[]subscript𝑛𝑖2superscriptdelimited-[]1superscriptsubscript𝑖1𝑘trsubscript𝐁𝑖1trsubscript𝐁𝑖𝑁𝑚2\frac{\Gamma_{1}[Nm/2]}{\Gamma_{1}[n_{0}m/2]}\prod_{i=1}^{k}\left(\frac{|% \mathbf{B}_{i}|^{(n_{i}-m-1)/2}}{\Gamma_{m}[n_{i}/2]}\right)\left[1+% \displaystyle\sum_{i=1}^{k}\frac{\mathop{\rm tr}\nolimits\mathbf{B}_{i}}{\left% (1-\mathop{\rm tr}\nolimits\mathbf{B}_{i}\right)}\right]^{-Nm/2}\hskip 85.3582% 6ptdivide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_N italic_m / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 ] end_ARG ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( divide start_ARG | bold_B start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_m - 1 ) / 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 ] end_ARG ) [ 1 + ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT divide start_ARG roman_tr bold_B start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG ( 1 - roman_tr bold_B start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) end_ARG ] start_POSTSUPERSCRIPT - italic_N italic_m / 2 end_POSTSUPERSCRIPT
×i=ik(1tr𝐁i)nim/21i=1k(d𝐁i).\hskip 113.81102pt\times\prod_{i=i}^{k}\left(1-\mathop{\rm tr}\nolimits\mathbf% {B}_{i}\right)^{-n_{i}m/2-1}\bigwedge_{i=1}^{k}\left(d\mathbf{B}_{i}\right).× ∏ start_POSTSUBSCRIPT italic_i = italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( 1 - roman_tr bold_B start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_m / 2 - 1 end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_d bold_B start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) . (9)

Finally

Lemma 3

Assume that 𝐗=(𝐗0,,𝐗k)𝐗superscriptsubscriptsuperscript𝐗0subscriptsuperscript𝐗𝑘\mathbf{X}=\left(\mathbf{X}^{\prime}_{0},\dots,\mathbf{X}^{\prime}_{k}\right)^% {\prime}bold_X = ( bold_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , … , bold_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT has a matrix variate spherical distribution, with 𝐗ini×msubscript𝐗𝑖superscriptsubscript𝑛𝑖𝑚\mathbf{X}_{i}\in\Re^{n_{i}\times m}bold_X start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT × italic_m end_POSTSUPERSCRIPT, nimsubscript𝑛𝑖𝑚n_{i}\geq mitalic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ≥ italic_m, i=0,1,,k𝑖01𝑘i=0,1,\dots,kitalic_i = 0 , 1 , … , italic_k. Define V=𝐗02𝑉superscriptnormsubscript𝐗02V=||\mathbf{X}_{0}||^{2}italic_V = | | bold_X start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and 𝐖i=𝐗i𝐗i>𝟎subscript𝐖𝑖subscriptsuperscript𝐗𝑖subscript𝐗𝑖0\mathbf{W}_{i}=\mathbf{X}^{\prime}_{i}\mathbf{X}_{i}>\mathbf{0}bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = bold_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT bold_X start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT > bold_0, i=1,,k𝑖1𝑘i=1,\dots,kitalic_i = 1 , … , italic_k.
Then, the joint density dFV,𝐖1,,𝐖k(v,𝐗1,,𝐖k)𝑑subscript𝐹𝑉subscript𝐖1subscript𝐖𝑘𝑣subscript𝐗1subscript𝐖𝑘dF_{V,\mathbf{W}_{1},\dots,\mathbf{W}_{k}}(v,\mathbf{X}_{1},\dots,\mathbf{W}_{% k})italic_d italic_F start_POSTSUBSCRIPT italic_V , bold_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_W start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_v , bold_X start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_W start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is given by

πNm/2vNm/21Γ1[n0m/2]i=1k(|𝐖i|(nim1)/2Γm[ni/2])h[v+i=1ktr𝐖i](dv)i=1k(d𝐖i),superscript𝜋𝑁𝑚2superscript𝑣𝑁𝑚21subscriptΓ1delimited-[]subscript𝑛0𝑚2superscriptsubscriptproduct𝑖1𝑘superscriptsubscript𝐖𝑖subscript𝑛𝑖𝑚12subscriptΓ𝑚delimited-[]subscript𝑛𝑖2delimited-[]𝑣superscriptsubscript𝑖1𝑘trsubscript𝐖𝑖𝑑𝑣superscriptsubscript𝑖1𝑘𝑑subscript𝐖𝑖\frac{\pi^{Nm/2}v^{Nm/2-1}}{\Gamma_{1}[n_{0}m/2]}\prod_{i=1}^{k}\left(\frac{|% \mathbf{W}_{i}|^{(n_{i}-m-1)/2}}{\Gamma_{m}[n_{i}/2]}\right)h\left[v+% \displaystyle\sum_{i=1}^{k}\mathop{\rm tr}\nolimits\mathbf{W}_{i}\right](dv)% \bigwedge_{i=1}^{k}\left(d\mathbf{W}_{i}\right),divide start_ARG italic_π start_POSTSUPERSCRIPT italic_N italic_m / 2 end_POSTSUPERSCRIPT italic_v start_POSTSUPERSCRIPT italic_N italic_m / 2 - 1 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 ] end_ARG ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( divide start_ARG | bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_m - 1 ) / 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 ] end_ARG ) italic_h [ italic_v + ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_tr bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ] ( italic_d italic_v ) ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_d bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , (10)

where V>0𝑉0V>0italic_V > 0. This distribution shall be termed multimatrix variate generalised Gamma - generalised Wishart distribution. The joint density function dF𝐖1,,𝐖k(𝐗1,,𝐖k)𝑑subscript𝐹subscript𝐖1subscript𝐖𝑘subscript𝐗1subscript𝐖𝑘dF_{\mathbf{W}_{1},\dots,\mathbf{W}_{k}}(\mathbf{X}_{1},\dots,\mathbf{W}_{k})italic_d italic_F start_POSTSUBSCRIPT bold_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_W start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_X start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_W start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is not invariant under spherical elliptical functions but it can be written in following closed form:

π(Nn0)m/2i=1k(|𝐖i|(nim1)/2Γm[ni/2])h[i=1ktr𝐖i]i=1k(d𝐖i),superscript𝜋𝑁subscript𝑛0𝑚2superscriptsubscriptproduct𝑖1𝑘superscriptsubscript𝐖𝑖subscript𝑛𝑖𝑚12subscriptΓ𝑚delimited-[]subscript𝑛𝑖2delimited-[]superscriptsubscript𝑖1𝑘trsubscript𝐖𝑖superscriptsubscript𝑖1𝑘𝑑subscript𝐖𝑖\pi^{(N-n_{0})m/2}\prod_{i=1}^{k}\left(\frac{|\mathbf{W}_{i}|^{(n_{i}-m-1)/2}}% {\Gamma_{m}[n_{i}/2]}\right)h\left[\displaystyle\sum_{i=1}^{k}\mathop{\rm tr}% \nolimits\mathbf{W}_{i}\right]\bigwedge_{i=1}^{k}\left(d\mathbf{W}_{i}\right),italic_π start_POSTSUPERSCRIPT ( italic_N - italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_m / 2 end_POSTSUPERSCRIPT ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( divide start_ARG | bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_m - 1 ) / 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 ] end_ARG ) italic_h [ ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_tr bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ] ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT ( italic_d bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , (11)

This marginal distribution shall be termed multimatrix variate generalised Wishart distribution.

As in the preceding lemmas, result (10) was derived in Díaz-García and Caro-Lopera (2024a). However, proposing the joint marginal (11) by expansion of h[v+i=1ktr𝐖i]delimited-[]𝑣superscriptsubscript𝑖1𝑘trsubscript𝐖𝑖h\left[v+\displaystyle\sum_{i=1}^{k}\mathop{\rm tr}\nolimits\mathbf{W}_{i}\right]italic_h [ italic_v + ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_tr bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ] will lead a simplification of complex kernel integrals involving general derivatives of elliptical models.

Observe that the parameter domain of the multimatrix variate distributions can be extended to the complex or real fields. However their geometrical and/or statistical explication perhaps can be lost. These distributions are valid if we replace ni/2subscript𝑛𝑖2n_{i}/2italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 by aisubscript𝑎𝑖a_{i}italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, n0m/2subscript𝑛0𝑚2n_{0}m/2italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 by a0subscript𝑎0a_{0}italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and Nm/2𝑁𝑚2Nm/2italic_N italic_m / 2 by a𝑎aitalic_a. Where the assuperscript𝑎𝑠a^{\prime s}italic_a start_POSTSUPERSCRIPT ′ italic_s end_POSTSUPERSCRIPT are complex numbers with positive real part. From a practical point of view for parametric estimation, this domain extension allows the use of nonlinear optimisation rather integer nonlinear optimisation, among other possibilities.

3 Computation and Probabilities on cones

We devote a few lines to the problem of computation of certain distributions of type (4). Certainly, there are hundreds of papers about matrix variate distribution theory, however the explicit computability of similar results to (4) is not usually addressed or applied. For example, the excellent work of Muirhead (2005) has attempted certain exact distributions by some approximations, but the majority of the results in non-central models forces the appearance and computation of A.T. James zonal polynomials of one matrix argument (James (1964)) or A.W. Davis invariant polynomials of several matrix arguments (Davis (1980)). The last problem was closed recently in Caro-Lopera (2016) by refuting Davis (1979) conjecture about a recursion computation of Davis functions in the same way as James polynomials, leaving intractable a number of theoretical results in non-central distribution theory. But, the remarkable property of James polynomials (James (1968)) (also known Jack polynomials in real normed division algebras) opened the possibility of computing a number of old series of zonal polynomials. The classical central cases on positive definite matrices were expressed in terms of hypergeometric series; for example probabilities of Wisharts matrices bounded by positive definite matrices, and so on. The key fact consists of provide a numerical computation of (tr𝐀)r=ρrCρ(𝐀)superscripttr𝐀𝑟subscript𝜌𝑟subscript𝐶𝜌𝐀(\mathop{\rm tr}\nolimits\mathbf{A})^{r}=\sum_{\rho\in r}C_{\rho}(\mathbf{A})( roman_tr bold_A ) start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT = ∑ start_POSTSUBSCRIPT italic_ρ ∈ italic_r end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( bold_A ), where 𝐀m×m𝐀superscript𝑚𝑚\mathbf{A}\in\Re^{m\times m}bold_A ∈ roman_ℜ start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT is positive definite, see details in (4). The numerical solution appeared very late in Koev and Edelman (2006) after more than a half century of theoretical results. Given the addressed James recurrence construction proposed in 1968, those algorithms are sufficient for low values of m𝑚mitalic_m, because the partitions ρ𝜌\rhoitalic_ρ are trivially truncated by m𝑚mitalic_m parts, but exact expressions for Cρ(𝐀)subscript𝐶𝜌𝐀C_{\rho}(\mathbf{A})italic_C start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( bold_A ) with arbitrary m=r𝑚𝑟m=ritalic_m = italic_r are out of any knowledge. In that context, and referring to S. Ramanujan and G.H. Hardy, knowing an exact formula for the number of partitions ρ𝜌\rhoitalic_ρ of arbitrary r𝑟ritalic_r is one of the biggest problem of mathematics in all history. The problem is such challenging that is so far to be included into the reasonable famous list of the Millennium Prize Problems (Clay Mathematics Institute of Cambridge (2000)). If the number of Cρ(𝐀)subscript𝐶𝜌𝐀C_{\rho}(\mathbf{A})italic_C start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( bold_A ) in the expansion of (tr𝐀)rsuperscripttr𝐀𝑟(\mathop{\rm tr}\nolimits\mathbf{A})^{r}( roman_tr bold_A ) start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT is just the number of partitions, just imagine the problem of providing an exact formula for ρrCρ(𝐀)subscript𝜌𝑟subscript𝐶𝜌𝐀\sum_{\rho\in r}C_{\rho}(\mathbf{A})∑ start_POSTSUBSCRIPT italic_ρ ∈ italic_r end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_ρ end_POSTSUBSCRIPT ( bold_A ). Other apparitions of the number of partitions arrive in permanents (Caro-Lopera et al. (2013)), and the general derivatives of a composite function and a Kotz model generator (Caro-Lopera et al. (2010)), among many others. Moreover, computing generalized hypergeometric series of zonal polynomials (4) only can be achieved under truncation, low values of m𝑚mitalic_m and suitable functions depending on the series index. A number of such series were computed recently by modification of the algorithms of Koev and Edelman (2006) in the context of statistical shape theory, see Caro-Lopera and Díaz-García (2012) and the references therein. In other context, several families of computable polynomial distributions based on zonal functions can be also seen in Caro-Lopera (2018).

Under this point of revision, we can mention that all the multimatricvariate and multimatrix variate distributions presented in Díaz-García and Caro-Lopera (2022), Díaz-García and Caro-Lopera (2024a) and Díaz-García and Caro-Lopera (2024b) are completelly computable. Most of them are free of series representation, then the computation is straighforwardly. Such is the case of multimatrix variate Pearson and Beta type distributions given in (6), (8), (7) and (9). The highlighted distributions involve the important property of invariance under the spherical family with generator h()h(\cdot)italic_h ( ⋅ ); this is crucial for a researcher, because no previous knowledge of the underlying distribution is required. No fitting distribution test must be done, except that the general assumption on ellipticity should be held.

Now, if the computation of matrix variate distribution is problematic, we can imagine the issues involved in finding a matrix probability. It should be noted that even simple univariate pdfs and cdfs for small, large or any beta, F and Wishart latent roots themselves have involved important historical problems in the last century. Discarding the trivial probabilities emerging from integrals involving zonal polynomials of an specific partition, it seems that the only existing computable matrix bounded probability are those for central Wishart (𝐖𝐖\mathbf{W}bold_W) and beta (𝐁𝐁\mathbf{B}bold_B) distribution. The Wishart lower bound probability P(𝐖<𝛀)𝑃𝐖𝛀P(\mathbf{W}<\mathbf{\Omega})italic_P ( bold_W < bold_Ω ) and the integral for matrix beta function back to Constantine (1963) (See Arias et al. (2021) for the P(𝐁<𝛀)𝑃𝐁𝛀P(\mathbf{B}<\mathbf{\Omega})italic_P ( bold_B < bold_Ω )). However, as Constantine (1963) states: “the complementary probabilities (P(𝐖>𝛀)𝑃𝐖𝛀P(\mathbf{W}>\mathbf{\Omega})italic_P ( bold_W > bold_Ω ), and P(𝐁<𝛀)𝑃𝐁𝛀P(\mathbf{B}<\mathbf{\Omega})italic_P ( bold_B < bold_Ω )) seem difficult to evaluate”. The Beta probability appeared recently in Arias et al. (2021). For the Wishart probability, the explicit solution in the Gaussian case arrived in 1982 (Muirhead (2005)), and a revision of the underlying proof was given in Caro-Lopera et al. (2016) joint with a generalization to elliptically contoured distributions for both lower and upper probabilities indexed by kernel (h()h(\cdot)italic_h ( ⋅ )). However, probabilities for rectangular matrices and non symmetric square matrices are still open problems. In fact, excepting for the interest of small and large latent roots, applications of the existing probabilities on Gaussian or elliptical Wishart distributions and classical Beta matrices seem out of consideration in literature. We try to set here the interest for those probabilities, in particular emerging from a multimatrix context, which naturally arises in experiment of probabilistic dependent samples. In this aspect, once the multimatrix joint distributions involve simplicity, we expect that some related measures considers tractable series representation. A suitable next stage of application in multimatrix variate Pearson and Beta type distributions should consider the computation of probabilities. In particular, multimatrix variate beta type II distribution in (7) has a potential application that we will explore in Section 4. For example, the following result considers the probability that 𝐀i𝐅isubscript𝐀𝑖subscript𝐅𝑖\mathbf{A}_{i}-\mathbf{F}_{i}bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - bold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is a positive definite matrix for all constant positive definite matrices 𝐀isubscript𝐀𝑖\mathbf{A}_{i}bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, i=1,,k𝑖1𝑘i=1,\ldots,kitalic_i = 1 , … , italic_k.

Theorem 1

Assume the hypothesis of Lemma 1 and suitable parameters for convergence of series (4). For positive definite matrices 𝐀im×msubscript𝐀𝑖superscript𝑚𝑚\mathbf{A}_{i}\in\Re^{m\times m}bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT, i=1,,k𝑖1𝑘i=1,\ldots,kitalic_i = 1 , … , italic_k, the probability P(𝟎<𝐅1<𝐀1,,𝟎<𝐅k<𝐀k,)P(\mathbf{0}<\mathbf{F}_{1}<\mathbf{A}_{1},\ldots,\mathbf{0}<\mathbf{F}_{k}<% \mathbf{A}_{k},)italic_P ( bold_0 < bold_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_0 < bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , ) is given by:

Γ1[Nm/2]Γmk[(m+1)/2]Γ1[n0m/2]i=1kΓm[(ni+m+1)/2]i=1k|𝐀i|ni/2×S(r2,,rk,𝐀1,,𝐀k),subscriptΓ1delimited-[]𝑁𝑚2superscriptsubscriptΓ𝑚𝑘delimited-[]𝑚12subscriptΓ1delimited-[]subscript𝑛0𝑚2superscriptsubscriptproduct𝑖1𝑘subscriptΓ𝑚delimited-[]subscript𝑛𝑖𝑚12superscriptsubscriptproduct𝑖1𝑘superscriptsubscript𝐀𝑖subscript𝑛𝑖2𝑆subscript𝑟2subscript𝑟𝑘subscript𝐀1subscript𝐀𝑘\frac{\Gamma_{1}[Nm/2]\Gamma_{m}^{k}[(m+1)/2]}{\Gamma_{1}[n_{0}m/2]% \displaystyle\prod_{i=1}^{k}\Gamma_{m}[(n_{i}+m+1)/2]}\prod_{i=1}^{k}|\mathbf{% A}_{i}|^{n_{i}/2}\times S(r_{2},\ldots,r_{k},\mathbf{A}_{1},\ldots,\mathbf{A}_% {k})\hskip 199.16928pt,divide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_N italic_m / 2 ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT [ ( italic_m + 1 ) / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 ] ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + italic_m + 1 ) / 2 ] end_ARG ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT | bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 end_POSTSUPERSCRIPT × italic_S ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , … , italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) ,

where the nested summation S(r2,,rk,𝐀1,,𝐀k)𝑆subscript𝑟2subscript𝑟𝑘subscript𝐀1subscript𝐀𝑘S(r_{2},\ldots,r_{k},\mathbf{A}_{1},\ldots,\mathbf{A}_{k})italic_S ( italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , … , italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) is given by

P11rk[hkP11rk1[hk1P11rk2P11r2[h2P11r1[h1:a1;b1;𝐀1]:a2;b2;𝐀2]]:{}_{1}^{r_{k}}P_{1}\left[h_{k}{}_{1}^{r_{k-1}}P_{1}\left[h_{k-1}{}_{1}^{r_{k-2% }}P_{1}{}_{1}^{r_{2}}P_{1}\left[h_{2}{}_{1}^{r_{1}}P_{1}\left[h_{1}:a_{1};b_{1% };-\mathbf{A}_{1}\right]:a_{2};b_{2};-\mathbf{A}_{2}\right]\cdots\right]\right.:start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_h start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_h start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k - 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_h start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_h start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT : italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; - bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ] : italic_a start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ; - bold_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ] ⋯ ] :
:ak1;bk1;𝐀k1]:ak;bk;𝐀k].\tiny:\left.\left.a_{k-1};b_{k-1};-\mathbf{A}_{k-1}\right]:a_{k};b_{k};-% \mathbf{A}_{k}\right].: italic_a start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT ; - bold_A start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT ] : italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ; - bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] . (12)

Here (12) involves k𝑘kitalic_k sums of type (4) and depends on the indexes r1,,rksubscript𝑟1subscript𝑟𝑘r_{1},\ldots,r_{k}italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT and matrices 𝐀1,,𝐀ksubscript𝐀1subscript𝐀𝑘\mathbf{A}_{1},\ldots,\mathbf{A}_{k}bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT, where hk=(c)rksubscript𝑘subscript𝑐subscript𝑟𝑘h_{k}=(c)_{r_{k}}italic_h start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = ( italic_c ) start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT,hj=(c+i=j+1kri)rj,j=1,,k1formulae-sequencesubscript𝑗subscript𝑐superscriptsubscript𝑖𝑗1𝑘subscript𝑟𝑖subscript𝑟𝑗𝑗1𝑘1h_{j}=\left(c+\sum_{i=j+1}^{k}r_{i}\right)_{r_{j}},j=1,\ldots,k-1italic_h start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = ( italic_c + ∑ start_POSTSUBSCRIPT italic_i = italic_j + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT , italic_j = 1 , … , italic_k - 1, c=Nm/2𝑐𝑁𝑚2c=Nm/2italic_c = italic_N italic_m / 2, ai=ni/2,bi=(ni+m+1)/2,i=1,,kformulae-sequencesubscript𝑎𝑖subscript𝑛𝑖2formulae-sequencesubscript𝑏𝑖subscript𝑛𝑖𝑚12𝑖1𝑘a_{i}=n_{i}/2,b_{i}=(n_{i}+m+1)/2,i=1,\ldots,kitalic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 , italic_b start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + italic_m + 1 ) / 2 , italic_i = 1 , … , italic_k.

Proof 1

Let c=Nm/2𝑐𝑁𝑚2c=Nm/2italic_c = italic_N italic_m / 2,

M=Γ1[Nm/2]Γ1[n0m/2]i=1kΓm[ni/2],𝑀subscriptΓ1delimited-[]𝑁𝑚2subscriptΓ1delimited-[]subscript𝑛0𝑚2superscriptsubscriptproduct𝑖1𝑘subscriptΓ𝑚delimited-[]subscript𝑛𝑖2M=\frac{\Gamma_{1}[Nm/2]}{\Gamma_{1}[n_{0}m/2]\prod_{i=1}^{k}\Gamma_{m}[n_{i}/% 2]},italic_M = divide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_N italic_m / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 ] ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 ] end_ARG ,

ai=ni/2,gi=(nim1)/2formulae-sequencesubscript𝑎𝑖subscript𝑛𝑖2subscript𝑔𝑖subscript𝑛𝑖𝑚12a_{i}=n_{i}/2,g_{i}=(n_{i}-m-1)/2italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 , italic_g start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_m - 1 ) / 2, and bi=(ni+m+1)/2subscript𝑏𝑖subscript𝑛𝑖𝑚12b_{i}=(n_{i}+m+1)/2italic_b start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + italic_m + 1 ) / 2, i=1,,k𝑖1𝑘i=1,\dots,kitalic_i = 1 , … , italic_k. Consider backwards integration starting in 𝟎<𝐅k<𝐀k0subscript𝐅𝑘subscript𝐀𝑘\mathbf{0}<\mathbf{F}_{k}<\mathbf{A}_{k}bold_0 < bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT. To avoid the appearance of generalized binomial expansions of zonal polynomials of sums of matrices, the integrand is written as

j=1k1|𝐅j|gj|𝐅k|gk(1+j=1k1tr𝐅j+tr𝐅k)c.superscriptsubscriptproduct𝑗1𝑘1superscriptsubscript𝐅𝑗subscript𝑔𝑗superscriptsubscript𝐅𝑘subscript𝑔𝑘superscript1superscriptsubscript𝑗1𝑘1trsubscript𝐅𝑗trsubscript𝐅𝑘𝑐\prod_{j=1}^{k-1}|\mathbf{F}_{j}|^{g_{j}}|\mathbf{F}_{k}|^{g_{k}}\left(1+% \displaystyle\sum_{j=1}^{k-1}\mathop{\rm tr}\nolimits\mathbf{F}_{j}+\mathop{% \rm tr}\nolimits\mathbf{F}_{k}\right)^{-c}.∏ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k - 1 end_POSTSUPERSCRIPT | bold_F start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT | bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( 1 + ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k - 1 end_POSTSUPERSCRIPT roman_tr bold_F start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + roman_tr bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - italic_c end_POSTSUPERSCRIPT .

Performing a convergent power expansion of the form (w+x)p=i=0(p)iwpii!superscript𝑤𝑥𝑝superscriptsubscript𝑖0subscript𝑝𝑖superscript𝑤𝑝𝑖𝑖(w+x)^{-p}=\sum_{i=0}^{\infty}\frac{(p)_{i}w^{-p-i}}{i!}( italic_w + italic_x ) start_POSTSUPERSCRIPT - italic_p end_POSTSUPERSCRIPT = ∑ start_POSTSUBSCRIPT italic_i = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG ( italic_p ) start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_w start_POSTSUPERSCRIPT - italic_p - italic_i end_POSTSUPERSCRIPT end_ARG start_ARG italic_i ! end_ARG and a variable substitution 𝐅k=𝐀k1/2𝐙k𝐀k1/2subscript𝐅𝑘superscriptsubscript𝐀𝑘12subscript𝐙𝑘superscriptsubscript𝐀𝑘12\mathbf{F}_{k}=\mathbf{A}_{k}^{1/2}\mathbf{Z}_{k}\mathbf{A}_{k}^{1/2}bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT bold_Z start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT, with Jacobian (d𝐅k)=|𝐀k|(m+1)/2(d𝐙k)𝑑subscript𝐅𝑘superscriptsubscript𝐀𝑘𝑚12𝑑subscript𝐙𝑘(d\mathbf{F}_{k})=|\mathbf{A}_{k}|^{(m+1)/2}(d\mathbf{Z}_{k})( italic_d bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) = | bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT ( italic_m + 1 ) / 2 end_POSTSUPERSCRIPT ( italic_d bold_Z start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ), lead to the application of Muirhead (2005, Th. 7.2.10) and (4). Thus,

j=1k1|𝐅j|gj𝟎<𝐅k<𝐀k|𝐅k|gk(1+j=1ktr𝐅j)c(d𝐅k)=Γm[ak]Γm[(m+1)/2]Γm[bk]|𝐀k|aksuperscriptsubscriptproduct𝑗1𝑘1superscriptsubscript𝐅𝑗subscript𝑔𝑗subscript0subscript𝐅𝑘subscript𝐀𝑘superscriptsubscript𝐅𝑘subscript𝑔𝑘superscript1superscriptsubscript𝑗1𝑘trsubscript𝐅𝑗𝑐𝑑subscript𝐅𝑘subscriptΓ𝑚delimited-[]subscript𝑎𝑘subscriptΓ𝑚delimited-[]𝑚12subscriptΓ𝑚delimited-[]subscript𝑏𝑘superscriptsubscript𝐀𝑘subscript𝑎𝑘\prod_{j=1}^{k-1}|\mathbf{F}_{j}|^{g_{j}}\int_{\mathbf{0}<\mathbf{F}_{k}<% \mathbf{A}_{k}}|\mathbf{F}_{k}|^{g_{k}}\left(1+\displaystyle\sum_{j=1}^{k}% \mathop{\rm tr}\nolimits\mathbf{F}_{j}\right)^{-c}\left(d\mathbf{F}_{k}\right)% =\frac{\Gamma_{m}[a_{k}]\Gamma_{m}[(m+1)/2]}{\Gamma_{m}[b_{k}]|\mathbf{A}_{k}|% ^{-a_{k}}}∏ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k - 1 end_POSTSUPERSCRIPT | bold_F start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT bold_0 < bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT | bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( 1 + ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_tr bold_F start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - italic_c end_POSTSUPERSCRIPT ( italic_d bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) = divide start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_m + 1 ) / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_b start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] | bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT - italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG
×P11rk(j=1k1|𝐅j|gj(1+j=1k1tr𝐅j)crk(c)rk1:ak;bk;𝐀k).\hskip 28.45274pt\times\,{}_{1}^{r_{k}}P_{1}\left(\frac{\prod_{j=1}^{k-1}|% \mathbf{F}_{j}|^{g_{j}}\left(1+\displaystyle\sum_{j=1}^{k-1}\mathop{\rm tr}% \nolimits\mathbf{F}_{j}\right)^{-c-r_{k}}}{(c)_{r_{k}}^{-1}}:a_{k};b_{k};-% \mathbf{A}_{k}\right).× start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( divide start_ARG ∏ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k - 1 end_POSTSUPERSCRIPT | bold_F start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( 1 + ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k - 1 end_POSTSUPERSCRIPT roman_tr bold_F start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - italic_c - italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_c ) start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT end_ARG : italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ; - bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) . (13)

Repeat a similar procedure to (13) for the next integral on 𝟎<𝐅k1<𝐀k10subscript𝐅𝑘1subscript𝐀𝑘1\mathbf{0}<\mathbf{F}_{k-1}<\mathbf{A}_{k-1}bold_0 < bold_F start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT by the corresponding splitting of the term indexed by rksubscript𝑟𝑘r_{k}italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT; thus the nested series P11rk1[]subscriptsuperscriptsubscript𝑃1subscript𝑟𝑘11delimited-[]{}_{1}^{r_{k-1}}P_{1}\left[\cdot\right]start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ⋅ ] is obtained. Explicitly, partial integration on 𝟎<𝐅k1<𝐀k1,𝟎<𝐅k<𝐀kformulae-sequence0subscript𝐅𝑘1subscript𝐀𝑘10subscript𝐅𝑘subscript𝐀𝑘\mathbf{0}<\mathbf{F}_{k-1}<\mathbf{A}_{k-1},\mathbf{0}<\mathbf{F}_{k}<\mathbf% {A}_{k}bold_0 < bold_F start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT , bold_0 < bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT takes the form

i=k1kΓm[ai]Γm2[(m+1)/2]i=k1kΓm[bi]i=k1k|𝐀i|aisuperscriptsubscriptproduct𝑖𝑘1𝑘subscriptΓ𝑚delimited-[]subscript𝑎𝑖superscriptsubscriptΓ𝑚2delimited-[]𝑚12superscriptsubscriptproduct𝑖𝑘1𝑘subscriptΓ𝑚delimited-[]subscript𝑏𝑖superscriptsubscriptproduct𝑖𝑘1𝑘superscriptsubscript𝐀𝑖subscript𝑎𝑖\frac{\prod_{i=k-1}^{k}\Gamma_{m}[a_{i}]\Gamma_{m}^{2}[(m+1)/2]}{\prod_{i=k-1}% ^{k}\Gamma_{m}[b_{i}]}\prod_{i=k-1}^{k}|\mathbf{A}_{i}|^{a_{i}}\hskip 199.1692% 8ptdivide start_ARG ∏ start_POSTSUBSCRIPT italic_i = italic_k - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ ( italic_m + 1 ) / 2 ] end_ARG start_ARG ∏ start_POSTSUBSCRIPT italic_i = italic_k - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_b start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ] end_ARG ∏ start_POSTSUBSCRIPT italic_i = italic_k - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT | bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUPERSCRIPT
×P11rk(P11rk1(j=1k2|𝐅j|gj(1+j=1k2tr𝐅j)ci=k1kri(c+i=k1kri)rk11:ak1;bk1;𝐀k1)(c)rk1:\tiny\times\,{}_{1}^{r_{k}}P_{1}\left(\frac{{}_{1}^{r_{k-1}}P_{1}\left(\frac{% \prod_{j=1}^{k-2}|\mathbf{F}_{j}|^{g_{j}}\left(1+\displaystyle\sum_{j=1}^{k-2}% \mathop{\rm tr}\nolimits\mathbf{F}_{j}\right)^{-c-\sum_{i=k-1}^{k}r_{i}}}{% \left(c+\sum_{i=k-1}^{k}r_{i}\right)_{r_{k-1}}^{-1}}:a_{k-1};b_{k-1};-\mathbf{% A}_{k-1}\right)}{(c)_{r_{k}}^{-1}}:\right.× start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( divide start_ARG start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( divide start_ARG ∏ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k - 2 end_POSTSUPERSCRIPT | bold_F start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_g start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( 1 + ∑ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k - 2 end_POSTSUPERSCRIPT roman_tr bold_F start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - italic_c - ∑ start_POSTSUBSCRIPT italic_i = italic_k - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_c + ∑ start_POSTSUBSCRIPT italic_i = italic_k - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT end_ARG : italic_a start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT ; - bold_A start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT ) end_ARG start_ARG ( italic_c ) start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT end_ARG :
:ak;bk;𝐀k).\hskip 284.52756pt\left.:a_{k};b_{k};-\mathbf{A}_{k}\right).: italic_a start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ; - bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ) .

Induction on risubscript𝑟𝑖r_{i}italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, by a similar integration procedure leading to (13), provides the nested series P11rk2[],,P11r2[]subscriptsuperscriptsubscript𝑃1subscript𝑟𝑘21delimited-[]subscriptsuperscriptsubscript𝑃1subscript𝑟21delimited-[]{}_{1}^{r_{k-2}}P_{1}\left[\cdot\right],\ldots,{}_{1}^{r_{2}}P_{1}\left[\cdot\right]start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k - 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ⋅ ] , … , start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ⋅ ] for the probability on 𝟎<𝐅2<𝐀2,,𝟎<𝐅k1<𝐀k1,𝟎<𝐅k<𝐀kformulae-sequence0subscript𝐅2subscript𝐀20subscript𝐅𝑘1subscript𝐀𝑘10subscript𝐅𝑘subscript𝐀𝑘\mathbf{0}<\mathbf{F}_{2}<\mathbf{A}_{2},\ldots,\mathbf{0}<\mathbf{F}_{k-1}<% \mathbf{A}_{k-1},\mathbf{0}<\mathbf{F}_{k}<\mathbf{A}_{k}bold_0 < bold_F start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , … , bold_0 < bold_F start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT , bold_0 < bold_F start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT. Finally, the last integral in 𝟎<𝐅1<𝐀10subscript𝐅1subscript𝐀1\mathbf{0}<\mathbf{F}_{1}<\mathbf{A}_{1}bold_0 < bold_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is given by:

𝟎<𝐅1<𝐀1|𝐅1|a1(1+tr𝐅1)ci=2kri(d𝐅1)=subscript0subscript𝐅1subscript𝐀1superscriptsubscript𝐅1subscript𝑎1superscript1trsubscript𝐅1𝑐superscriptsubscript𝑖2𝑘subscript𝑟𝑖𝑑subscript𝐅1absent\int_{\mathbf{0}<\mathbf{F}_{1}<\mathbf{A}_{1}}|\mathbf{F}_{1}|^{a_{1}}\left(1% +\mathop{\rm tr}\nolimits\mathbf{F}_{1}\right)^{-c-\sum_{i=2}^{k}r_{i}}(d% \mathbf{F}_{1})=∫ start_POSTSUBSCRIPT bold_0 < bold_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT | bold_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( 1 + roman_tr bold_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - italic_c - ∑ start_POSTSUBSCRIPT italic_i = 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( italic_d bold_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) =
Γm[a1]Γm[(m+1)/2]Γm[b1]|𝐀1|a1P11r1((c+i=2kri)r1:a1;b1;𝐀1),\frac{\Gamma_{m}[a_{1}]\Gamma_{m}[(m+1)/2]}{\Gamma_{m}[b_{1}]|\mathbf{A}_{1}|^% {-a_{1}}}\quad{}_{1}^{r_{1}}P_{1}\left(\left(c+\sum_{i=2}^{k}r_{i}\right)_{r_{% 1}}:a_{1};b_{1};\mathbf{A}_{1}\right),divide start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_m + 1 ) / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ] | bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT - italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT end_ARG start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( ( italic_c + ∑ start_POSTSUBSCRIPT italic_i = 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT : italic_a start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ; bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) ,

then the required probability (12) is obtained.

As usual in Wishart type probabilities, if λisubscript𝜆𝑖\lambda_{i}italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is the largest latent root of 𝐅isubscript𝐅𝑖\mathbf{F}_{i}bold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, then Theorem 1 provides the distribution function of λisubscript𝜆𝑖\lambda_{i}italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, P(λi<x)𝑃subscript𝜆𝑖𝑥P(\lambda_{i}<x)italic_P ( italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT < italic_x ), by taking 𝐀i=x𝐈subscript𝐀𝑖𝑥𝐈\mathbf{A}_{i}=x\mathbf{I}bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_x bold_I in (12).

We now focus on the multimatrix variate generalised Wishart distribution (11). This is the only distribution presented here that is not invariant under the spherical family. First note that (11) is indexed by v>0vNm/21h(r)(v)𝑑vsubscript𝑣0superscript𝑣𝑁𝑚21superscript𝑟𝑣differential-d𝑣\int_{v>0}v^{Nm/2-1}h^{(r)}(v)dv∫ start_POSTSUBSCRIPT italic_v > 0 end_POSTSUBSCRIPT italic_v start_POSTSUPERSCRIPT italic_N italic_m / 2 - 1 end_POSTSUPERSCRIPT italic_h start_POSTSUPERSCRIPT ( italic_r ) end_POSTSUPERSCRIPT ( italic_v ) italic_d italic_v. It is a simple integral once the general derivatives are obtained. They can be found in Caro-Lopera et al. (2010) for the classical elliptical generators of Pearson, Kotz, Bessel, Jensen-Logistic and so on. For example, the Gaussian case of the multimatrix variate generalised Wishart distribution is obtained by taking the generator function h(v)=(2π)Nm/2ev/2𝑣superscript2𝜋𝑁𝑚2superscript𝑒𝑣2h(v)=(2\pi)^{-Nm/2}e^{-v/2}italic_h ( italic_v ) = ( 2 italic_π ) start_POSTSUPERSCRIPT - italic_N italic_m / 2 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_v / 2 end_POSTSUPERSCRIPT, then

v>0vNm/21h(t)(v)𝑑v=πNm/2Γ[Nm/2](12)t.subscript𝑣0superscript𝑣𝑁𝑚21superscript𝑡𝑣differential-d𝑣superscript𝜋𝑁𝑚2Γdelimited-[]𝑁𝑚2superscript12𝑡\int_{v>0}v^{Nm/2-1}h^{(t)}(v)dv=\pi^{-Nm/2}\Gamma[Nm/2]\left(-\frac{1}{2}% \right)^{t}.∫ start_POSTSUBSCRIPT italic_v > 0 end_POSTSUBSCRIPT italic_v start_POSTSUPERSCRIPT italic_N italic_m / 2 - 1 end_POSTSUPERSCRIPT italic_h start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT ( italic_v ) italic_d italic_v = italic_π start_POSTSUPERSCRIPT - italic_N italic_m / 2 end_POSTSUPERSCRIPT roman_Γ [ italic_N italic_m / 2 ] ( - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ) start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT .

Excepting the extra summation, the procedure for Theorem 1 can be reproduced next for probabilities on m×m𝑚𝑚m\times mitalic_m × italic_m positive definite matrices 𝐖isubscript𝐖𝑖\mathbf{W}_{i}bold_W start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, i=1,,k𝑖1𝑘i=1,\ldots,kitalic_i = 1 , … , italic_k.

Theorem 2

Consider the hypothesis of Lemma 3 and positive definite matrices 𝐀im×msubscript𝐀𝑖superscript𝑚𝑚\mathbf{A}_{i}\in\Re^{m\times m}bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ roman_ℜ start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT, i=1,,k𝑖1𝑘i=1,\ldots,kitalic_i = 1 , … , italic_k. Then we have the following probabilities:

1) Gaussian (independent) case:

P(𝟎<𝐖1<𝐀1,,𝟎<𝐖k<𝐀k,)=P(\mathbf{0}<\mathbf{W}_{1}<\mathbf{A}_{1},\ldots,\mathbf{0}<\mathbf{W}_{k}<% \mathbf{A}_{k},)=\hskip 199.16928ptitalic_P ( bold_0 < bold_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_0 < bold_W start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , ) =
Γmk[(m+1)/2]i=1k|12𝐀i|ni/2i=1kΓm[(ni+m+1)/2]i=1kF11[ni2;ni+m+12;12𝐀i]superscriptsubscriptΓ𝑚𝑘delimited-[]𝑚12superscriptsubscriptproduct𝑖1𝑘superscript12subscript𝐀𝑖subscript𝑛𝑖2superscriptsubscriptproduct𝑖1𝑘subscriptΓ𝑚delimited-[]subscript𝑛𝑖𝑚12superscriptsubscriptproduct𝑖1𝑘subscriptsubscript𝐹11subscript𝑛𝑖2subscript𝑛𝑖𝑚1212subscript𝐀𝑖\frac{\Gamma_{m}^{k}[(m+1)/2]\displaystyle\prod_{i=1}^{k}\left|\frac{1}{2}% \mathbf{A}_{i}\right|^{n_{i}/2}}{\displaystyle\prod_{i=1}^{k}\Gamma_{m}[(n_{i}% +m+1)/2]}\prod_{i=1}^{k}{}_{1}F_{1}\left[\frac{n_{i}}{2};\frac{n_{i}+m+1}{2};-% \frac{1}{2}\mathbf{A}_{i}\right]divide start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT [ ( italic_m + 1 ) / 2 ] ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT | divide start_ARG 1 end_ARG start_ARG 2 end_ARG bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 end_POSTSUPERSCRIPT end_ARG start_ARG ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + italic_m + 1 ) / 2 ] end_ARG ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT italic_F start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ divide start_ARG italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ; divide start_ARG italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + italic_m + 1 end_ARG start_ARG 2 end_ARG ; - divide start_ARG 1 end_ARG start_ARG 2 end_ARG bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ] (14)

2) General elliptical model, k=1𝑘1k=1italic_k = 1 (Caro-Lopera et al. (2016)):

P(𝟎<𝐖1<𝐀1)=Γm[m+12]|𝐀1|n12Γm[n1+m+12]πn1m2t=0h(t)(0)t!Q11t[n12;n1+m+12;𝐀1]𝑃0subscript𝐖1subscript𝐀1subscriptΓ𝑚delimited-[]𝑚12superscriptsubscript𝐀1subscript𝑛12subscriptΓ𝑚delimited-[]subscript𝑛1𝑚12superscript𝜋subscript𝑛1𝑚2superscriptsubscript𝑡0superscript𝑡0𝑡subscriptsuperscriptsubscript𝑄1𝑡1subscript𝑛12subscript𝑛1𝑚12subscript𝐀1P(\mathbf{0}<\mathbf{W}_{1}<\mathbf{A}_{1})=\frac{\Gamma_{m}[\frac{m+1}{2}]|% \mathbf{A}_{1}|^{\frac{n_{1}}{2}}}{\Gamma_{m}[\frac{n_{1}+m+1}{2}]\pi^{-\frac{% n_{1}m}{2}}}\sum_{t=0}^{\infty}\frac{h^{(t)}(0)}{t!}\ {}_{1}^{t}Q_{1}\left[% \frac{n_{1}}{2};\frac{n_{1}+m+1}{2};\mathbf{A}_{1}\right]italic_P ( bold_0 < bold_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = divide start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ divide start_ARG italic_m + 1 end_ARG start_ARG 2 end_ARG ] | bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT divide start_ARG italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ divide start_ARG italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_m + 1 end_ARG start_ARG 2 end_ARG ] italic_π start_POSTSUPERSCRIPT - divide start_ARG italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_m end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_t = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG italic_h start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT ( 0 ) end_ARG start_ARG italic_t ! end_ARG start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ divide start_ARG italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ; divide start_ARG italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_m + 1 end_ARG start_ARG 2 end_ARG ; bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ] (15)

3) General elliptical model, for all k>1𝑘1k>1italic_k > 1:

P(𝟎<𝐖1<𝐀1,,𝟎<𝐖k<𝐀k,)=P(\mathbf{0}<\mathbf{W}_{1}<\mathbf{A}_{1},\ldots,\mathbf{0}<\mathbf{W}_{k}<% \mathbf{A}_{k},)=\hskip 199.16928ptitalic_P ( bold_0 < bold_W start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_0 < bold_W start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT , ) =
π(Nn0)m/2Γmk[(m+1)/2]i=1k|𝐀i|ni/2i=1kΓm[(ni+m+1)/2]t=0h(t)(0)t!superscript𝜋𝑁subscript𝑛0𝑚2superscriptsubscriptΓ𝑚𝑘delimited-[]𝑚12superscriptsubscriptproduct𝑖1𝑘superscriptsubscript𝐀𝑖subscript𝑛𝑖2superscriptsubscriptproduct𝑖1𝑘subscriptΓ𝑚delimited-[]subscript𝑛𝑖𝑚12superscriptsubscript𝑡0superscript𝑡0𝑡\frac{\pi^{(N-n_{0})m/2}\Gamma_{m}^{k}[(m+1)/2]\displaystyle\prod_{i=1}^{k}|% \mathbf{A}_{i}|^{n_{i}/2}}{\displaystyle\prod_{i=1}^{k}\Gamma_{m}[(n_{i}+m+1)/% 2]}\sum_{t=0}^{\infty}\frac{h^{(t)}(0)}{t!}\hskip 170.71652ptdivide start_ARG italic_π start_POSTSUPERSCRIPT ( italic_N - italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_m / 2 end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT [ ( italic_m + 1 ) / 2 ] ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT | bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT / 2 end_POSTSUPERSCRIPT end_ARG start_ARG ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT + italic_m + 1 ) / 2 ] end_ARG ∑ start_POSTSUBSCRIPT italic_t = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG italic_h start_POSTSUPERSCRIPT ( italic_t ) end_POSTSUPERSCRIPT ( 0 ) end_ARG start_ARG italic_t ! end_ARG
×r1=0t(tr1)Q11r1[n12;n1+m+12;𝐀1]\times\,\sum_{r_{1}=0}^{t}\binom{t}{r_{1}}\,{}_{1}^{r_{1}}Q_{1}\left[\frac{n_{% 1}}{2};\frac{n_{1}+m+1}{2};\mathbf{A}_{1}\right]× ∑ start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT ( FRACOP start_ARG italic_t end_ARG start_ARG italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG ) start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ divide start_ARG italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ; divide start_ARG italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_m + 1 end_ARG start_ARG 2 end_ARG ; bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ]
×j=1k2{rj+1=0ti=1jri(ti=1jrirj+1)Q11rj+1[nj+12;nj+1+m+12;𝐀j+1]}\times\prod_{j=1}^{k-2}\left\{\sum_{r_{j+1}=0}^{t-\sum_{i=1}^{j}r_{i}}\binom{t% -\sum_{i=1}^{j}r_{i}}{r_{j+1}}\,{}_{1}^{r_{j+1}}Q_{1}\left[\frac{n_{j+1}}{2};% \frac{n_{j+1}+m+1}{2};\mathbf{A}_{j+1}\right]\right\}× ∏ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k - 2 end_POSTSUPERSCRIPT { ∑ start_POSTSUBSCRIPT italic_r start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_t - ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( FRACOP start_ARG italic_t - ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_j end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_r start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT end_ARG ) start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ divide start_ARG italic_n start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ; divide start_ARG italic_n start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT + italic_m + 1 end_ARG start_ARG 2 end_ARG ; bold_A start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT ] }
×Q11rk[nk2;nk+m+12;𝐀k]absentsubscriptsuperscriptsubscript𝑄1subscript𝑟𝑘1subscript𝑛𝑘2subscript𝑛𝑘𝑚12subscript𝐀𝑘\times\,{}_{1}^{r_{k}}Q_{1}\left[\frac{n_{k}}{2};\frac{n_{k}+m+1}{2};\mathbf{A% }_{k}\right]× start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ divide start_ARG italic_n start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ; divide start_ARG italic_n start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT + italic_m + 1 end_ARG start_ARG 2 end_ARG ; bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ] (16)
Proof 2

(14) is derived by using the generator function h(v)=(2π)Nm/2ev/2𝑣superscript2𝜋𝑁𝑚2superscript𝑒𝑣2h(v)=(2\pi)^{-Nm/2}e^{-v/2}italic_h ( italic_v ) = ( 2 italic_π ) start_POSTSUPERSCRIPT - italic_N italic_m / 2 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_v / 2 end_POSTSUPERSCRIPT in (11), and hence the integrals follows by k𝑘kitalic_k independent application of the isotropic version of Caro-Lopera et al. (2016, Cor. 7) or Muirhead (2005, Th. 9.7.1 ). (15) was derived in Caro-Lopera et al. (2016, th.6) as a generalization of the Gaussian case given in Muirhead (2005, Th. 9.7.1 ).

Finally, (16) follows by simple induction on k𝑘kitalic_k. It requires a recurrent use of the binomial theorem, the representation (5) and the known integral

𝟎<𝐖<𝐀|𝐖|am+12trr(𝐖)(d𝐖)subscript0𝐖𝐀superscript𝐖𝑎𝑚12superscripttr𝑟𝐖𝑑𝐖\int_{\mathbf{0}<\mathbf{W}<\mathbf{A}}|\mathbf{W}|^{a-\frac{m+1}{2}}\mathop{% \rm tr}\nolimits^{r}(\mathbf{W})(d\mathbf{W})\hskip 170.71652pt∫ start_POSTSUBSCRIPT bold_0 < bold_W < bold_A end_POSTSUBSCRIPT | bold_W | start_POSTSUPERSCRIPT italic_a - divide start_ARG italic_m + 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT roman_tr start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT ( bold_W ) ( italic_d bold_W )
=Γm[a]Γm[(m+1)/2]Γm[a+(m+1)/2]|𝐀|aQ11r(a;a+m+12;𝐀).absentsubscriptΓ𝑚delimited-[]𝑎subscriptΓ𝑚delimited-[]𝑚12subscriptΓ𝑚delimited-[]𝑎𝑚12superscript𝐀𝑎subscriptsuperscriptsubscript𝑄1𝑟1𝑎𝑎𝑚12𝐀\hskip 113.81102pt=\frac{\Gamma_{m}[a]\Gamma_{m}[(m+1)/2]}{\Gamma_{m}[a+(m+1)/% 2]|\mathbf{A}|^{-a}}\,{}_{1}^{r}Q_{1}\left(a;a+\frac{m+1}{2};\mathbf{A}\right).= divide start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_m + 1 ) / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a + ( italic_m + 1 ) / 2 ] | bold_A | start_POSTSUPERSCRIPT - italic_a end_POSTSUPERSCRIPT end_ARG start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_a ; italic_a + divide start_ARG italic_m + 1 end_ARG start_ARG 2 end_ARG ; bold_A ) .
Lemma 4

In the setting of Lemma 3 and the cone probabilities here derived, the only feasible lower probability is reached for k=1𝑘1k=1italic_k = 1 and is given by:

P(𝟎<𝐁1<𝐀1<𝐈)=Γm[n1/2]Γm[(m+1))/2]Γm[(n1+m+1))/2]P(\mathbf{0}<\mathbf{B}_{1}<\mathbf{A}_{1}<\mathbf{I})=\frac{\Gamma_{m}[n_{1}/% 2]\Gamma_{m}[(m+1))/2]}{\Gamma_{m}[(n_{1}+m+1))/2]}\hskip 170.71652ptitalic_P ( bold_0 < bold_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_I ) = divide start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT / 2 ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_m + 1 ) ) / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_m + 1 ) ) / 2 ] end_ARG
×|𝐀1|n1/2t=0(n0m/21t)Q11t[n1/2;(n1+m+1)/2;𝐀1].absentsuperscriptsubscript𝐀1subscript𝑛12superscriptsubscript𝑡0binomialsubscript𝑛0𝑚21𝑡subscriptsuperscriptsubscript𝑄1𝑡1subscript𝑛12subscript𝑛1𝑚12subscript𝐀1\times|\mathbf{A}_{1}|^{n_{1}/2}\sum_{t=0}^{\infty}\binom{n_{0}m/2-1}{t}{}_{1}% ^{t}Q_{1}\left[n_{1}/2;(n_{1}+m+1)/2;-\mathbf{A}_{1}\right].× | bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT / 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_t = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ( FRACOP start_ARG italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 - 1 end_ARG start_ARG italic_t end_ARG ) start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_t end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT / 2 ; ( italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_m + 1 ) / 2 ; - bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ] . (17)

For k>1𝑘1k>1italic_k > 1 the probability P(𝟎<𝐁i<𝐀i<𝐈,,𝟎<𝐁k<𝐀k<𝐈)P(\mathbf{0}<\mathbf{B}_{i}<\mathbf{A}_{i}<\mathbf{I},\ldots,\mathbf{0}<% \mathbf{B}_{k}<\mathbf{A}_{k}<\mathbf{I})italic_P ( bold_0 < bold_B start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT < bold_I , … , bold_0 < bold_B start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT < bold_I ) turns in terms of invariant polynomials and it cannot be computable.

Proof 3

Let k=1𝑘1k=1italic_k = 1 and 𝐀1subscript𝐀1\mathbf{A}_{1}bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, such that 𝟎<𝐁1<𝐀1<𝐈0subscript𝐁1subscript𝐀1𝐈\mathbf{0}<\mathbf{B}_{1}<\mathbf{A}_{1}<\mathbf{I}bold_0 < bold_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < bold_I and tr𝐀11trsubscript𝐀11\mathop{\rm tr}\nolimits\mathbf{A}_{1}\leq 1roman_tr bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≤ 1. The integrand in (9) simplifies to |𝐁1|(n1m1)/2(1tr𝐁1)n0m/21superscriptsubscript𝐁1subscript𝑛1𝑚12superscript1trsubscript𝐁1subscript𝑛0𝑚21|\mathbf{B}_{1}|^{(n_{1}-m-1)/2}\left(1-\mathop{\rm tr}\nolimits\mathbf{B}_{1}% \right)^{n_{0}m/2-1}| bold_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT ( italic_n start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_m - 1 ) / 2 end_POSTSUPERSCRIPT ( 1 - roman_tr bold_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_n start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m / 2 - 1 end_POSTSUPERSCRIPT. Then, the result (17) follows after binomial theorem and application of

𝟎<𝐘<𝐗|𝐘|am+12(tr𝐘)r(d𝐘)subscript0𝐘𝐗superscript𝐘𝑎𝑚12superscripttr𝐘𝑟𝑑𝐘\int_{\mathbf{0}<\mathbf{Y}<\mathbf{X}}|\mathbf{Y}|^{a-\frac{m+1}{2}}(\mathop{% \rm tr}\nolimits\mathbf{Y})^{r}(d\mathbf{Y})\hskip 170.71652pt∫ start_POSTSUBSCRIPT bold_0 < bold_Y < bold_X end_POSTSUBSCRIPT | bold_Y | start_POSTSUPERSCRIPT italic_a - divide start_ARG italic_m + 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( roman_tr bold_Y ) start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT ( italic_d bold_Y )
=Γm[a]Γm[(m+1)/2]Γm[a+(m+1)/2]|𝐗|aQ11r(a;a+m+12;𝐗).absentsubscriptΓ𝑚delimited-[]𝑎subscriptΓ𝑚delimited-[]𝑚12subscriptΓ𝑚delimited-[]𝑎𝑚12superscript𝐗𝑎subscriptsuperscriptsubscript𝑄1𝑟1𝑎𝑎𝑚12𝐗\hskip 113.81102pt=\frac{\Gamma_{m}[a]\Gamma_{m}[(m+1)/2]}{\Gamma_{m}[a+(m+1)/% 2]|\mathbf{X}|^{-a}}\,{}_{1}^{r}Q_{1}\left(a;a+\frac{m+1}{2};\mathbf{X}\right).= divide start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_m + 1 ) / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a + ( italic_m + 1 ) / 2 ] | bold_X | start_POSTSUPERSCRIPT - italic_a end_POSTSUPERSCRIPT end_ARG start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_a ; italic_a + divide start_ARG italic_m + 1 end_ARG start_ARG 2 end_ARG ; bold_X ) .

For k>1𝑘1k>1\,italic_k > 1 appear product of powers of traces arising invariant polynomials in the multiple integration, then the computation turns impossible by a similar Laplace Beltrami operator computation of zonal polynomials (Caro-Lopera, 2016).

We end this section in the context of real normed division algebras. We refer to Baez (2002) and Díaz-García, and Gutiérrez-Jáimez (2013) and the references therein for a complete exposition of the topic. For our purposes we index the four normed division algebras by the real dimension β𝛽\betaitalic_β, where β=1𝛽1\beta=1italic_β = 1, stands for Real; β=2𝛽2\beta=2italic_β = 2, for Complex; β=4𝛽4\beta=4italic_β = 4, for Quaternionic; and β=8𝛽8\beta=8italic_β = 8, for Octonion. Other notations for the algebras are given by α=2/β𝛼2𝛽\alpha=2/\betaitalic_α = 2 / italic_β, see Edelman and Rao (2005).

For the sequel some notations and definitions are required. For an understandable comparison with the real case in Section 2, we provide a similar complete exposition for the real normed division algebras. Let m,nβsubscriptsuperscript𝛽𝑚𝑛{\mathcal{L}}^{\beta}_{m,n}caligraphic_L start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT be the linear space of all n×m𝑛𝑚n\times mitalic_n × italic_m matrices of rank mn𝑚𝑛m\leq nitalic_m ≤ italic_n over a real finite-dimensional normed division algebra 𝔉𝔉\mathfrak{F}fraktur_F with m𝑚mitalic_m distinct positive singular values. Let 𝔉n×msuperscript𝔉𝑛𝑚\mathfrak{F}^{n\times m}fraktur_F start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT be the set of all n×m𝑛𝑚n\times mitalic_n × italic_m matrices over 𝔉𝔉\mathfrak{F}fraktur_F. The dimension of 𝔉n×msuperscript𝔉𝑛𝑚\mathfrak{F}^{n\times m}fraktur_F start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT over \Reroman_ℜ is βmn𝛽𝑚𝑛\beta mnitalic_β italic_m italic_n. Let 𝐀𝔉n×m𝐀superscript𝔉𝑛𝑚\mathbf{A}\in\mathfrak{F}^{n\times m}bold_A ∈ fraktur_F start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT, then 𝐀H=𝐀¯Tsuperscript𝐀𝐻superscript¯𝐀𝑇\mathbf{A}^{H}=\overline{\mathbf{A}}^{T}bold_A start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT = over¯ start_ARG bold_A end_ARG start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT denotes the usual conjugate transpose.

The set of matrices 𝐇1𝔉n×msubscript𝐇1superscript𝔉𝑛𝑚\mathbf{H}_{1}\in\mathfrak{F}^{n\times m}bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∈ fraktur_F start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT such that 𝐇1H𝐇1=𝐈msuperscriptsubscript𝐇1𝐻subscript𝐇1subscript𝐈𝑚\mathbf{H}_{1}^{H}\mathbf{H}_{1}=\mathbf{I}_{m}bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = bold_I start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT is a manifold denoted 𝒱m,nβsuperscriptsubscript𝒱𝑚𝑛𝛽{\mathcal{V}}_{m,n}^{\beta}caligraphic_V start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT, is termed the Stiefel manifold (𝐇1subscript𝐇1\mathbf{H}_{1}bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is also known as semi-orthogonal (β=1𝛽1\beta=1italic_β = 1), semi-unitary (β=2𝛽2\beta=2italic_β = 2), semi-symplectic (β=4𝛽4\beta=4italic_β = 4) and semi-exceptional type (β=8𝛽8\beta=8italic_β = 8) matrices, see Dray and Manogue (1999)). The dimension of 𝒱m,nβsuperscriptsubscript𝒱𝑚𝑛𝛽\mathcal{V}_{m,n}^{\beta}caligraphic_V start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT over \Reroman_ℜ is [βmnm(m1)β/2m]delimited-[]𝛽𝑚𝑛𝑚𝑚1𝛽2𝑚[\beta mn-m(m-1)\beta/2-m][ italic_β italic_m italic_n - italic_m ( italic_m - 1 ) italic_β / 2 - italic_m ]. In particular, 𝒱m,mβsuperscriptsubscript𝒱𝑚𝑚𝛽{\mathcal{V}}_{m,m}^{\beta}caligraphic_V start_POSTSUBSCRIPT italic_m , italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT with dimension over \Reroman_ℜ, [m(m+1)β/2m]delimited-[]𝑚𝑚1𝛽2𝑚[m(m+1)\beta/2-m][ italic_m ( italic_m + 1 ) italic_β / 2 - italic_m ], is the maximal compact subgroup 𝔘β(m)superscript𝔘𝛽𝑚\mathfrak{U}^{\beta}(m)fraktur_U start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( italic_m ) of m,mβsubscriptsuperscript𝛽𝑚𝑚{\mathcal{L}}^{\beta}_{m,m}caligraphic_L start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m , italic_m end_POSTSUBSCRIPT and consists of all matrices 𝐇𝔉m×m𝐇superscript𝔉𝑚𝑚\mathbf{H}\in\mathfrak{F}^{m\times m}bold_H ∈ fraktur_F start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT such that 𝐇H𝐇=𝐈msuperscript𝐇𝐻𝐇subscript𝐈𝑚\mathbf{H}^{H}\mathbf{H}=\mathbf{I}_{m}bold_H start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_H = bold_I start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT. Therefore, 𝔘β(m)superscript𝔘𝛽𝑚\mathfrak{U}^{\beta}(m)fraktur_U start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( italic_m ) is the real orthogonal group 𝒪(m)𝒪𝑚\mathcal{O}(m)caligraphic_O ( italic_m ) (β=1𝛽1\beta=1italic_β = 1), the unitary group 𝒰(m)𝒰𝑚\mathcal{U}(m)caligraphic_U ( italic_m ) (β=2𝛽2\beta=2italic_β = 2), compact symplectic group 𝒮p(m)𝒮𝑝𝑚\mathcal{S}p(m)caligraphic_S italic_p ( italic_m ) (β=4𝛽4\beta=4italic_β = 4) or exceptional type matrices 𝒪o(m)𝒪𝑜𝑚\mathcal{O}o(m)caligraphic_O italic_o ( italic_m ) (β=8𝛽8\beta=8italic_β = 8), for 𝔉=𝔉\mathfrak{F}=\Refraktur_F = roman_ℜ, \mathfrak{C}fraktur_C, \mathfrak{H}fraktur_H or 𝔒𝔒\mathfrak{O}fraktur_O, respectively.

We denote by 𝔖mβsuperscriptsubscript𝔖𝑚𝛽{\mathfrak{S}}_{m}^{\beta}fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT the real vector space of all 𝐒𝔉m×m𝐒superscript𝔉𝑚𝑚\mathbf{S}\in\mathfrak{F}^{m\times m}bold_S ∈ fraktur_F start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT such that 𝐒=𝐒H𝐒superscript𝐒𝐻\mathbf{S}=\mathbf{S}^{H}bold_S = bold_S start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT. Let 𝔓mβsuperscriptsubscript𝔓𝑚𝛽\mathfrak{P}_{m}^{\beta}fraktur_P start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT be the cone of positive definite matrices 𝐒𝔉m×m𝐒superscript𝔉𝑚𝑚\mathbf{S}\in\mathfrak{F}^{m\times m}bold_S ∈ fraktur_F start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT; then 𝔓mβsuperscriptsubscript𝔓𝑚𝛽\mathfrak{P}_{m}^{\beta}fraktur_P start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT is an open subset of 𝔖mβsuperscriptsubscript𝔖𝑚𝛽{\mathfrak{S}}_{m}^{\beta}fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT. Over \Reroman_ℜ, 𝔖mβsuperscriptsubscript𝔖𝑚𝛽{\mathfrak{S}}_{m}^{\beta}fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT consist of symmetric matrices; over \mathfrak{C}fraktur_C, Hermitian matrices; over \mathfrak{H}fraktur_H, quaternionic Hermitian matrices (also termed self-dual matrices) and over 𝔒𝔒\mathfrak{O}fraktur_O, octonionic Hermitian matrices. Generically, the elements of 𝔖mβsuperscriptsubscript𝔖𝑚𝛽\mathfrak{S}_{m}^{\beta}fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT are termed Hermitian matrices, irrespective of the nature of 𝔉𝔉\mathfrak{F}fraktur_F. The dimension of 𝔖mβsuperscriptsubscript𝔖𝑚𝛽\mathfrak{S}_{m}^{\beta}fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT over \Reroman_ℜ is [m(m1)β+2m]/2delimited-[]𝑚𝑚1𝛽2𝑚2[m(m-1)\beta+2m]/2[ italic_m ( italic_m - 1 ) italic_β + 2 italic_m ] / 2.

Let 𝔇mβsuperscriptsubscript𝔇𝑚𝛽\mathfrak{D}_{m}^{\beta}fraktur_D start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT be the diagonal subgroup of m,mβsuperscriptsubscript𝑚𝑚𝛽\mathcal{L}_{m,m}^{\beta}caligraphic_L start_POSTSUBSCRIPT italic_m , italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT consisting of all 𝐃𝔉m×m𝐃superscript𝔉𝑚𝑚\mathbf{D}\in\mathfrak{F}^{m\times m}bold_D ∈ fraktur_F start_POSTSUPERSCRIPT italic_m × italic_m end_POSTSUPERSCRIPT, 𝐃=diag(d1,,dm)𝐃diagsubscript𝑑1subscript𝑑𝑚\mathbf{D}=\mathop{\rm diag}\nolimits(d_{1},\dots,d_{m})bold_D = roman_diag ( italic_d start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , italic_d start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ).

For any matrix 𝐗𝔉n×m𝐗superscript𝔉𝑛𝑚\mathbf{X}\in\mathfrak{F}^{n\times m}bold_X ∈ fraktur_F start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT, d𝐗𝑑𝐗d\mathbf{X}italic_d bold_X denotes the matrix of differentials (dxij)𝑑subscript𝑥𝑖𝑗(dx_{ij})( italic_d italic_x start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ). Finally, we define the measure or volume element (d𝐗)𝑑𝐗(d\mathbf{X})( italic_d bold_X ) when 𝐗𝔉m×n,𝔖mβ𝐗superscript𝔉𝑚𝑛superscriptsubscript𝔖𝑚𝛽\mathbf{X}\in\mathfrak{F}^{m\times n},\mathfrak{S}_{m}^{\beta}bold_X ∈ fraktur_F start_POSTSUPERSCRIPT italic_m × italic_n end_POSTSUPERSCRIPT , fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT, 𝔇mβsuperscriptsubscript𝔇𝑚𝛽\mathfrak{D}_{m}^{\beta}fraktur_D start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT or 𝒱m,nβsuperscriptsubscript𝒱𝑚𝑛𝛽\mathcal{V}_{m,n}^{\beta}caligraphic_V start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT.

If 𝐗𝔉n×m𝐗superscript𝔉𝑛𝑚\mathbf{X}\in\mathfrak{F}^{n\times m}bold_X ∈ fraktur_F start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT then (d𝐗)𝑑𝐗(d\mathbf{X})( italic_d bold_X ) (the Lebesgue measure in 𝔉n×msuperscript𝔉𝑛𝑚\mathfrak{F}^{n\times m}fraktur_F start_POSTSUPERSCRIPT italic_n × italic_m end_POSTSUPERSCRIPT) denotes the exterior product of the βmn𝛽𝑚𝑛\beta mnitalic_β italic_m italic_n functionally independent variables

(d𝐗)=i=1nj=1mdxij where dxij=r=1βdxij(r).formulae-sequence𝑑𝐗superscriptsubscript𝑖1𝑛superscriptsubscript𝑗1𝑚𝑑subscript𝑥𝑖𝑗 where 𝑑subscript𝑥𝑖𝑗superscriptsubscript𝑟1𝛽𝑑superscriptsubscript𝑥𝑖𝑗𝑟(d\mathbf{X})=\bigwedge_{i=1}^{n}\bigwedge_{j=1}^{m}dx_{ij}\quad\mbox{ where }% \quad dx_{ij}=\bigwedge_{r=1}^{\beta}dx_{ij}^{(r)}.( italic_d bold_X ) = ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_j = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_d italic_x start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT where italic_d italic_x start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = ⋀ start_POSTSUBSCRIPT italic_r = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT italic_d italic_x start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_r ) end_POSTSUPERSCRIPT .

If 𝐒𝔖mβ𝐒superscriptsubscript𝔖𝑚𝛽\mathbf{S}\in\mathfrak{S}_{m}^{\beta}bold_S ∈ fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT (or 𝐒𝔗Lβ(m)𝐒superscriptsubscript𝔗𝐿𝛽𝑚\mathbf{S}\in\mathfrak{T}_{L}^{\beta}(m)bold_S ∈ fraktur_T start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( italic_m )) then (d𝐒)𝑑𝐒(d\mathbf{S})( italic_d bold_S ) (the Lebesgue measure in 𝔖mβsuperscriptsubscript𝔖𝑚𝛽\mathfrak{S}_{m}^{\beta}fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT or in 𝔗Lβ(m)superscriptsubscript𝔗𝐿𝛽𝑚\mathfrak{T}_{L}^{\beta}(m)fraktur_T start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( italic_m )) denotes the exterior product of the m(m+1)β/2𝑚𝑚1𝛽2m(m+1)\beta/2italic_m ( italic_m + 1 ) italic_β / 2 functionally independent variables (or denotes the exterior product of the m(m1)β/2+m𝑚𝑚1𝛽2𝑚m(m-1)\beta/2+mitalic_m ( italic_m - 1 ) italic_β / 2 + italic_m functionally independent variables, if siisubscript𝑠𝑖𝑖s_{ii}\in\Reitalic_s start_POSTSUBSCRIPT italic_i italic_i end_POSTSUBSCRIPT ∈ roman_ℜ for all i=1,,m𝑖1𝑚i=1,\dots,mitalic_i = 1 , … , italic_m)

(d𝐒)={ijmr=1βdsij(r),i=1mdsiii<jmr=1βdsij(r),if sii.𝑑𝐒casessuperscriptsubscript𝑖𝑗𝑚superscriptsubscript𝑟1𝛽𝑑superscriptsubscript𝑠𝑖𝑗𝑟missing-subexpressionsuperscriptsubscript𝑖1𝑚𝑑subscript𝑠𝑖𝑖superscriptsubscript𝑖𝑗𝑚superscriptsubscript𝑟1𝛽𝑑superscriptsubscript𝑠𝑖𝑗𝑟if subscript𝑠𝑖𝑖(d\mathbf{S})=\left\{\begin{array}[]{ll}\displaystyle\bigwedge_{i\leq j}^{m}% \bigwedge_{r=1}^{\beta}ds_{ij}^{(r)},&\\ \displaystyle\bigwedge_{i=1}^{m}ds_{ii}\bigwedge_{i<j}^{m}\bigwedge_{r=1}^{% \beta}ds_{ij}^{(r)},&\hbox{if }s_{ii}\in\Re.\end{array}\right.( italic_d bold_S ) = { start_ARRAY start_ROW start_CELL ⋀ start_POSTSUBSCRIPT italic_i ≤ italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_r = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT italic_d italic_s start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_r ) end_POSTSUPERSCRIPT , end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_d italic_s start_POSTSUBSCRIPT italic_i italic_i end_POSTSUBSCRIPT ⋀ start_POSTSUBSCRIPT italic_i < italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_r = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT italic_d italic_s start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_r ) end_POSTSUPERSCRIPT , end_CELL start_CELL if italic_s start_POSTSUBSCRIPT italic_i italic_i end_POSTSUBSCRIPT ∈ roman_ℜ . end_CELL end_ROW end_ARRAY

The context generally establishes the conditions on the elements of 𝐒𝐒\mathbf{S}bold_S, that is, if sijsubscript𝑠𝑖𝑗s_{ij}\in\Reitalic_s start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ∈ roman_ℜ, absent\in\mathfrak{C}∈ fraktur_C, absent\in\mathfrak{H}∈ fraktur_H or 𝔒absent𝔒\in\mathfrak{O}∈ fraktur_O. It is considered that

(d𝐒)=ijmr=1βdsij(r)i=1mdsiii<jmr=1βdsij(r).𝑑𝐒superscriptsubscript𝑖𝑗𝑚superscriptsubscript𝑟1𝛽𝑑superscriptsubscript𝑠𝑖𝑗𝑟superscriptsubscript𝑖1𝑚𝑑subscript𝑠𝑖𝑖superscriptsubscript𝑖𝑗𝑚superscriptsubscript𝑟1𝛽𝑑superscriptsubscript𝑠𝑖𝑗𝑟(d\mathbf{S})=\bigwedge_{i\leq j}^{m}\bigwedge_{r=1}^{\beta}ds_{ij}^{(r)}% \equiv\bigwedge_{i=1}^{m}ds_{ii}\bigwedge_{i<j}^{m}\bigwedge_{r=1}^{\beta}ds_{% ij}^{(r)}.( italic_d bold_S ) = ⋀ start_POSTSUBSCRIPT italic_i ≤ italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_r = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT italic_d italic_s start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_r ) end_POSTSUPERSCRIPT ≡ ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_d italic_s start_POSTSUBSCRIPT italic_i italic_i end_POSTSUBSCRIPT ⋀ start_POSTSUBSCRIPT italic_i < italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_r = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT italic_d italic_s start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_r ) end_POSTSUPERSCRIPT .

Observe, too, that for the Lebesgue measure (d𝐒)𝑑𝐒(d\mathbf{S})( italic_d bold_S ) defined thus, it is required that 𝐒𝔓mβ𝐒superscriptsubscript𝔓𝑚𝛽\mathbf{S}\in\mathfrak{P}_{m}^{\beta}bold_S ∈ fraktur_P start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT, that is, 𝐒𝐒\mathbf{S}bold_S must be a non singular Hermitian matrix (Hermitian positive definite matrix).

If 𝚲𝔇mβ𝚲superscriptsubscript𝔇𝑚𝛽\mathbf{\Lambda}\in\mathfrak{D}_{m}^{\beta}bold_Λ ∈ fraktur_D start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT then (d𝚲)𝑑𝚲(d\mathbf{\Lambda})( italic_d bold_Λ ) (the Legesgue measure in 𝔇mβsuperscriptsubscript𝔇𝑚𝛽\mathfrak{D}_{m}^{\beta}fraktur_D start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT) denotes the exterior product of the βm𝛽𝑚\beta mitalic_β italic_m functionally independent variables

(d𝚲)=i=1nr=1βdλi(r).𝑑𝚲superscriptsubscript𝑖1𝑛superscriptsubscript𝑟1𝛽𝑑superscriptsubscript𝜆𝑖𝑟(d\mathbf{\Lambda})=\bigwedge_{i=1}^{n}\bigwedge_{r=1}^{\beta}d\lambda_{i}^{(r% )}.( italic_d bold_Λ ) = ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_r = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT italic_d italic_λ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_r ) end_POSTSUPERSCRIPT .

If 𝐇1𝒱m,nβsubscript𝐇1superscriptsubscript𝒱𝑚𝑛𝛽\mathbf{H}_{1}\in\mathcal{V}_{m,n}^{\beta}bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∈ caligraphic_V start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT then

(𝐇1Hd𝐇1)=i=1nj=i+1m𝐡jHd𝐡i.subscriptsuperscript𝐇𝐻1𝑑subscript𝐇1superscriptsubscript𝑖1𝑛superscriptsubscript𝑗𝑖1𝑚superscriptsubscript𝐡𝑗𝐻𝑑subscript𝐡𝑖(\mathbf{H}^{H}_{1}d\mathbf{H}_{1})=\bigwedge_{i=1}^{n}\bigwedge_{j=i+1}^{m}% \mathbf{h}_{j}^{H}d\mathbf{h}_{i}.( bold_H start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_d bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = ⋀ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ⋀ start_POSTSUBSCRIPT italic_j = italic_i + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT bold_h start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT italic_d bold_h start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT .

where 𝐇=(𝐇1|𝐇2)=(𝐡1,,𝐡m|𝐡m+1,,𝐡n)𝔘β(m)𝐇conditionalsubscript𝐇1subscript𝐇2subscript𝐡1conditionalsubscript𝐡𝑚subscript𝐡𝑚1subscript𝐡𝑛superscript𝔘𝛽𝑚\mathbf{H}=(\mathbf{H}_{1}|\mathbf{H}_{2})=(\mathbf{h}_{1},\dots,\mathbf{h}_{m% }|\mathbf{h}_{m+1},\dots,\mathbf{h}_{n})\in\mathfrak{U}^{\beta}(m)bold_H = ( bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | bold_H start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) = ( bold_h start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , … , bold_h start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT | bold_h start_POSTSUBSCRIPT italic_m + 1 end_POSTSUBSCRIPT , … , bold_h start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ) ∈ fraktur_U start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( italic_m ). It can be proved that this differential form does not depend on the choice of the 𝐇2subscript𝐇2\mathbf{H}_{2}bold_H start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT matrix. When m=1𝑚1m=1italic_m = 1; 𝒱1,nβsubscriptsuperscript𝒱𝛽1𝑛\mathcal{V}^{\beta}_{1,n}caligraphic_V start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 , italic_n end_POSTSUBSCRIPT defines the unit sphere in 𝔉nsuperscript𝔉𝑛\mathfrak{F}^{n}fraktur_F start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT. This is, of course, an (n1)β𝑛1𝛽(n-1)\beta( italic_n - 1 ) italic_β- dimensional surface in 𝔉nsuperscript𝔉𝑛\mathfrak{F}^{n}fraktur_F start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT. When m=n𝑚𝑛m=nitalic_m = italic_n and denoting 𝐇1subscript𝐇1\mathbf{H}_{1}bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT by 𝐇𝐇\mathbf{H}bold_H, (𝐇Hd𝐇)superscript𝐇𝐻𝑑𝐇(\mathbf{H}^{H}d\mathbf{H})( bold_H start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT italic_d bold_H ) is termed the Haar measure on 𝔘β(m)superscript𝔘𝛽𝑚\mathfrak{U}^{\beta}(m)fraktur_U start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( italic_m ).

The surface area or volume of the Stiefel manifold 𝒱m,nβsubscriptsuperscript𝒱𝛽𝑚𝑛\mathcal{V}^{\beta}_{m,n}caligraphic_V start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT is

Vol(𝒱m,nβ)=𝐇1𝒱m,nβ(𝐇1Hd𝐇1)=2mπmnβ/2Γmβ[nβ/2],Volsubscriptsuperscript𝒱𝛽𝑚𝑛subscriptsubscript𝐇1subscriptsuperscript𝒱𝛽𝑚𝑛subscriptsuperscript𝐇𝐻1𝑑subscript𝐇1superscript2𝑚superscript𝜋𝑚𝑛𝛽2subscriptsuperscriptΓ𝛽𝑚delimited-[]𝑛𝛽2\mathop{\rm Vol}\nolimits(\mathcal{V}^{\beta}_{m,n})=\int_{\mathbf{H}_{1}\in% \mathcal{V}^{\beta}_{m,n}}(\mathbf{H}^{H}_{1}d\mathbf{H}_{1})=\frac{2^{m}\pi^{% mn\beta/2}}{\Gamma^{\beta}_{m}[n\beta/2]},roman_Vol ( caligraphic_V start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT ) = ∫ start_POSTSUBSCRIPT bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ∈ caligraphic_V start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( bold_H start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_d bold_H start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = divide start_ARG 2 start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT italic_π start_POSTSUPERSCRIPT italic_m italic_n italic_β / 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Γ start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_n italic_β / 2 ] end_ARG , (18)

where Γmβ[a]subscriptsuperscriptΓ𝛽𝑚delimited-[]𝑎\Gamma^{\beta}_{m}[a]roman_Γ start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a ] denotes the multivariate Gamma function for the space 𝔖mβsuperscriptsubscript𝔖𝑚𝛽\mathfrak{S}_{m}^{\beta}fraktur_S start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT, and is defined by

Γmβ[a]superscriptsubscriptΓ𝑚𝛽delimited-[]𝑎\displaystyle\Gamma_{m}^{\beta}[a]roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT [ italic_a ] =\displaystyle== 𝐀𝔓mβetr{𝐀}|𝐀|a(m1)β/21(d𝐀)subscript𝐀superscriptsubscript𝔓𝑚𝛽etr𝐀superscript𝐀𝑎𝑚1𝛽21𝑑𝐀\displaystyle\displaystyle\int_{\mathbf{A}\in\mathfrak{P}_{m}^{\beta}}\mathop{% \rm etr}\nolimits\{-\mathbf{A}\}|\mathbf{A}|^{a-(m-1)\beta/2-1}(d\mathbf{A})∫ start_POSTSUBSCRIPT bold_A ∈ fraktur_P start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT end_POSTSUBSCRIPT roman_etr { - bold_A } | bold_A | start_POSTSUPERSCRIPT italic_a - ( italic_m - 1 ) italic_β / 2 - 1 end_POSTSUPERSCRIPT ( italic_d bold_A )
=\displaystyle== πm(m1)β/4i=1mΓ[a(i1)β/2],superscript𝜋𝑚𝑚1𝛽4superscriptsubscriptproduct𝑖1𝑚Γdelimited-[]𝑎𝑖1𝛽2\displaystyle\pi^{m(m-1)\beta/4}\displaystyle\prod_{i=1}^{m}\Gamma[a-(i-1)% \beta/2],italic_π start_POSTSUPERSCRIPT italic_m ( italic_m - 1 ) italic_β / 4 end_POSTSUPERSCRIPT ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT roman_Γ [ italic_a - ( italic_i - 1 ) italic_β / 2 ] ,

where etr()=exp(tr())etrtr\mathop{\rm etr}\nolimits(\cdot)=\exp(\mathop{\rm tr}\nolimits(\cdot))roman_etr ( ⋅ ) = roman_exp ( roman_tr ( ⋅ ) ), |||\cdot|| ⋅ | denotes the determinant and Re(a)>(m1)β/2Re𝑎𝑚1𝛽2\mathop{\rm Re}\nolimits(a)>(m-1)\beta/2roman_Re ( italic_a ) > ( italic_m - 1 ) italic_β / 2, see Gross and Richards (1987). If 𝐀m,nβ𝐀superscriptsubscript𝑚𝑛𝛽\mathbf{A}\in\mathcal{L}_{m,n}^{\beta}bold_A ∈ caligraphic_L start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT then by vec(𝐀)vec𝐀\mathop{\rm vec}\nolimits(\mathbf{A})roman_vec ( bold_A ) we mean the mn×1𝑚𝑛1mn\times 1italic_m italic_n × 1 vector formed by stacking the columns of 𝐀𝐀\mathbf{A}bold_A under each other.

Now, generalized statistical theory of shape has been developed by the authors in a number of settings: SVD, polar, PseudoWishart, QR, affine, forms, Eulerian, etc.. In particular, the real configuration or affine density was set in the addressed families of elliptically countored distributions in Caro-Lopera et al. (2010) as a revision and generalization of the Gaussian case given by Goodall and Mardia (1993). Then Díaz-García and Caro-Lopera (2016) proposed the affine shape theory under the general approach for real normed division algebras. The configuration or affine shape filters are interested in removing geometrical information about translation, scaling, rotation, reflection and/or uniform shearing of random objects summarized by matrices in the addressed four real normed divison algebras. Explicitly, two figures 𝐗K,Nβ𝐗subscriptsuperscript𝛽𝐾𝑁\mathbf{X}\in{\mathcal{L}}^{\beta}_{K,N}bold_X ∈ caligraphic_L start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_K , italic_N end_POSTSUBSCRIPT have the same affine shape if 𝐗1=𝐗𝐄+𝟏Nesubscript𝐗1𝐗𝐄subscript1𝑁superscript𝑒\mathbf{X}_{1}=\mathbf{X}\mathbf{E}+\mathbf{1}_{N}e^{*}bold_X start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = bold_XE + bold_1 start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT, for some translation e1,Kβ𝑒subscriptsuperscript𝛽1𝐾e\in{\mathcal{L}}^{\beta}_{1,K}italic_e ∈ caligraphic_L start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 , italic_K end_POSTSUBSCRIPT and 𝐄K,Kβ𝐄subscriptsuperscript𝛽𝐾𝐾\mathbf{E}\in{\mathcal{L}}^{\beta}_{K,K}bold_E ∈ caligraphic_L start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_K , italic_K end_POSTSUBSCRIPT. Then the (N1)×K𝑁1𝐾(N-1)\times K( italic_N - 1 ) × italic_K configuration matrix 𝐔=(𝐈|(𝐘2𝐘11))𝐔superscriptconditional𝐈superscriptsubscript𝐘2superscriptsubscript𝐘11\mathbf{U}=(\mathbf{I}|(\mathbf{Y}_{2}\mathbf{Y}_{1}^{-1})^{\prime})^{\prime}bold_U = ( bold_I | ( bold_Y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT bold_Y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT compressing the meaningful geometrical information of the original N×K𝑁𝐾N\times Kitalic_N × italic_K matrix 𝐙𝐙\mathbf{Z}bold_Z, is obtained in the sequence of filtering geometrical stages 𝐋𝐙=𝐘=𝐔𝐄𝐋𝐙𝐘𝐔𝐄\mathbf{L}\mathbf{Z}=\mathbf{Y}=\mathbf{U}\mathbf{E}bold_LZ = bold_Y = bold_UE. Here 𝐘=(𝐘1|𝐘2)𝐘superscriptconditionalsuperscriptsubscript𝐘1superscriptsubscript𝐘2\mathbf{Y}=(\mathbf{Y}_{1}^{\prime}|\mathbf{Y}_{2}^{\prime})^{\prime}bold_Y = ( bold_Y start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT | bold_Y start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT and 𝐋𝐋\mathbf{L}bold_L is a subHelmert matrix, see Díaz-García and Caro-Lopera (2016), for details. Now, a similar definition to (2) emerges for real normed division algebras. We say that 𝐗m,nβ𝐗subscriptsuperscript𝛽𝑚𝑛\mathbf{X}\in{\mathcal{L}}^{\beta}_{m,n}bold_X ∈ caligraphic_L start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT has the following multivariate elliptically contoured distribution for real normed division algebras respect to the Lebesgue measure:

F𝐗(𝐗)=|𝚺|βm/2|𝚯|βn/2h{tr[𝚯1(𝐗𝝁)𝚺1(𝐗𝝁)]},subscript𝐹𝐗𝐗superscript𝚺𝛽𝑚2superscript𝚯𝛽𝑛2trdelimited-[]superscript𝚯1superscript𝐗𝝁superscript𝚺1𝐗𝝁F_{{}_{\mathbf{X}}}(\mathbf{X})=|\mathbf{\Sigma}|^{-\beta m/2}|\mathbf{\Theta}% |^{-\beta n/2}h\left\{\mathop{\rm tr}\nolimits\left[\mathbf{\Theta}^{-1}(% \mathbf{X}-\boldsymbol{\mu})^{*}\mathbf{\Sigma}^{-1}(\mathbf{X}-\boldsymbol{% \mu})\right]\right\},italic_F start_POSTSUBSCRIPT start_FLOATSUBSCRIPT bold_X end_FLOATSUBSCRIPT end_POSTSUBSCRIPT ( bold_X ) = | bold_Σ | start_POSTSUPERSCRIPT - italic_β italic_m / 2 end_POSTSUPERSCRIPT | bold_Θ | start_POSTSUPERSCRIPT - italic_β italic_n / 2 end_POSTSUPERSCRIPT italic_h { roman_tr [ bold_Θ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( bold_X - bold_italic_μ ) start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT bold_Σ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( bold_X - bold_italic_μ ) ] } , (19)

where 𝝁m,nβ𝝁subscriptsuperscript𝛽𝑚𝑛\boldsymbol{\mu}\in{\mathcal{L}}^{\beta}_{m,n}bold_italic_μ ∈ caligraphic_L start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m , italic_n end_POSTSUBSCRIPT, 𝚺𝔓nβ𝚺superscriptsubscript𝔓𝑛𝛽\boldsymbol{\Sigma}\in\mathfrak{P}_{n}^{\beta}bold_Σ ∈ fraktur_P start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT, 𝚯𝔓nβ𝚯superscriptsubscript𝔓𝑛𝛽\boldsymbol{\Theta}\in\mathfrak{P}_{n}^{\beta}bold_Θ ∈ fraktur_P start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT, and generator function h𝔉[0,)𝔉0h\mbox{: }\mathfrak{F}\to[0,\infty)italic_h : fraktur_F → [ 0 , ∞ ), satisfies 𝔓1βuβnm1h(u2)𝑑u<subscriptsuperscriptsubscript𝔓1𝛽superscript𝑢𝛽𝑛𝑚1superscript𝑢2differential-d𝑢\int_{\mathfrak{P}_{1}^{\beta}}u^{\beta nm-1}h(u^{2})du<\infty∫ start_POSTSUBSCRIPT fraktur_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_β italic_n italic_m - 1 end_POSTSUPERSCRIPT italic_h ( italic_u start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_d italic_u < ∞. This fact will be denoted by 𝐗n×mβ(𝝁,𝚺𝚯;h)similar-to𝐗superscriptsubscript𝑛𝑚𝛽𝝁tensor-product𝚺𝚯\mathbf{X}\sim\mathcal{E}_{n\times m}^{\beta}(\boldsymbol{\mu},\boldsymbol{% \Sigma}\otimes\boldsymbol{\Theta};h)bold_X ∼ caligraphic_E start_POSTSUBSCRIPT italic_n × italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( bold_italic_μ , bold_Σ ⊗ bold_Θ ; italic_h ) Finally, for a convergent Taylor series of h()h(\cdot)italic_h ( ⋅ ), if 𝐘N1×Kβ(𝝁𝚯1/2,𝚺𝐈K,h)similar-to𝐘superscriptsubscript𝑁1𝐾𝛽𝝁superscript𝚯12tensor-product𝚺subscript𝐈𝐾\mathbf{Y}\sim\mathcal{E}_{N-1\times K}^{\beta}(\boldsymbol{\mu}\boldsymbol{% \Theta}^{-1/2},\boldsymbol{\Sigma}\otimes\mathbf{I}_{K},h)bold_Y ∼ caligraphic_E start_POSTSUBSCRIPT italic_N - 1 × italic_K end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( bold_italic_μ bold_Θ start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT , bold_Σ ⊗ bold_I start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT , italic_h ), the affine shape density of 𝐔𝐔\mathbf{U}bold_U is given by

g(𝐔,t,r)γt,r,𝑔𝐔𝑡𝑟subscript𝛾𝑡𝑟g(\mathbf{U},t,r)\gamma_{t,r},italic_g ( bold_U , italic_t , italic_r ) italic_γ start_POSTSUBSCRIPT italic_t , italic_r end_POSTSUBSCRIPT , (20)

with

g(𝐔,t,r)=πβK2/2ΓKβ[β(N1)/2]ΓKβ[βK/2]|𝚺|βK/2|𝐔𝚺1𝐔|β(N1)/2t=01t!Γ[K(N1)/2+t]𝑔𝐔𝑡𝑟superscript𝜋𝛽superscript𝐾22superscriptsubscriptΓ𝐾𝛽delimited-[]𝛽𝑁12superscriptsubscriptΓ𝐾𝛽delimited-[]𝛽𝐾2superscript𝚺𝛽𝐾2superscriptsuperscript𝐔superscript𝚺1𝐔𝛽𝑁12superscriptsubscript𝑡01𝑡Γdelimited-[]𝐾𝑁12𝑡g(\mathbf{U},t,r)=\frac{\pi^{\beta K^{2}/2}\Gamma_{K}^{\beta}[\beta(N-1)/2]}{% \Gamma_{K}^{\beta}[\beta K/2]|\boldsymbol{\Sigma}|^{\beta K/2}|\mathbf{U}^{*}% \boldsymbol{\Sigma}^{-1}\mathbf{U}|^{\beta(N-1)/2}}\sum_{t=0}^{\infty}\frac{1}% {t!\Gamma[K(N-1)/2+t]}\hskip 170.71652ptitalic_g ( bold_U , italic_t , italic_r ) = divide start_ARG italic_π start_POSTSUPERSCRIPT italic_β italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / 2 end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT [ italic_β ( italic_N - 1 ) / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT [ italic_β italic_K / 2 ] | bold_Σ | start_POSTSUPERSCRIPT italic_β italic_K / 2 end_POSTSUPERSCRIPT | bold_U start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT bold_Σ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_U | start_POSTSUPERSCRIPT italic_β ( italic_N - 1 ) / 2 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_t = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG italic_t ! roman_Γ [ italic_K ( italic_N - 1 ) / 2 + italic_t ] end_ARG
×r=0trr𝛀r!τ(β(N1)/2)τβ(βK/2)τβCτβ(𝐔Ω𝚺1𝐔(𝐔𝚺1𝐔)1),\times\sum_{r=0}^{\infty}\frac{\mathop{\rm tr}\nolimits^{r}\boldsymbol{\Omega}% }{r!}\sum_{\tau}\frac{\left(\beta(N-1)/2\right)_{\tau}^{\beta}}{\left(\beta K/% 2\right)_{\tau}^{\beta}}C_{\tau}^{\beta}\left(\mathbf{U}^{*}\Omega\boldsymbol{% \Sigma}^{-1}\mathbf{U}(\mathbf{U}^{*}\boldsymbol{\Sigma}^{-1}\mathbf{U})^{-1}% \right),× ∑ start_POSTSUBSCRIPT italic_r = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG roman_tr start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT bold_Ω end_ARG start_ARG italic_r ! end_ARG ∑ start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT divide start_ARG ( italic_β ( italic_N - 1 ) / 2 ) start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT end_ARG start_ARG ( italic_β italic_K / 2 ) start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT end_ARG italic_C start_POSTSUBSCRIPT italic_τ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT ( bold_U start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT roman_Ω bold_Σ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_U ( bold_U start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT bold_Σ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_U ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) , (21)

and

γt,r=𝔳1βvβK(N1)/2+t1h(2t+r)(v)𝑑v<,subscript𝛾𝑡𝑟subscriptsuperscriptsubscript𝔳1𝛽superscript𝑣𝛽𝐾𝑁12𝑡1superscript2𝑡𝑟𝑣differential-d𝑣\gamma_{t,r}=\int_{\mathfrak{v}_{1}^{\beta}}v^{\beta K(N-1)/2+t-1}h^{(2t+r)}(v% )dv<\infty,italic_γ start_POSTSUBSCRIPT italic_t , italic_r end_POSTSUBSCRIPT = ∫ start_POSTSUBSCRIPT fraktur_v start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_v start_POSTSUPERSCRIPT italic_β italic_K ( italic_N - 1 ) / 2 + italic_t - 1 end_POSTSUPERSCRIPT italic_h start_POSTSUPERSCRIPT ( 2 italic_t + italic_r ) end_POSTSUPERSCRIPT ( italic_v ) italic_d italic_v < ∞ , (22)

𝚺=𝐋𝚺𝐗𝐋𝚺𝐋subscript𝚺𝐗superscript𝐋\boldsymbol{\Sigma}=\mathbf{L}\boldsymbol{\Sigma}_{\mathbf{X}}\mathbf{L}^{*}bold_Σ = bold_L bold_Σ start_POSTSUBSCRIPT bold_X end_POSTSUBSCRIPT bold_L start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT, 𝝁=𝐋𝝁𝐗𝝁𝐋subscript𝝁𝐗\boldsymbol{\mu}=\mathbf{L}\boldsymbol{\mu}_{\mathbf{X}}bold_italic_μ = bold_L bold_italic_μ start_POSTSUBSCRIPT bold_X end_POSTSUBSCRIPT, 𝛀=𝚺1𝝁𝚯𝝁𝛀superscript𝚺1𝝁𝚯superscript𝝁\boldsymbol{\Omega}=\boldsymbol{\Sigma}^{-1}\boldsymbol{\mu}\boldsymbol{\Theta% }\boldsymbol{\mu}^{*}bold_Ω = bold_Σ start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT bold_italic_μ bold_Θ bold_italic_μ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT.

Now, in the computational context of this Section, the distribution (20) is completely feasible; Caro-Lopera et al. (2010) and Díaz-García and Caro-Lopera (2016) have provided a number of applications. However, the density requires the integral (22) which is set in terms of the 2t+r2𝑡𝑟2t+r2 italic_t + italic_r general derivative of the kernel function h()h(\cdot)italic_h ( ⋅ ). Inspired by Theorem 2, we can avoid integration and non null evaluation of the derivative in (22) by a simple computation of h(2t+r)(0)superscript2𝑡𝑟0h^{(2t+r)}(0)italic_h start_POSTSUPERSCRIPT ( 2 italic_t + italic_r ) end_POSTSUPERSCRIPT ( 0 ). Next result provides the solution for the three models given by Díaz-García, and Gutiérrez-Jáimez (2013).

Theorem 3

The affine shape density of 𝐔𝐔\mathbf{U}bold_U for Gaussian, Pearson type VII and II models is given by

g(𝐔,t,r)ζt,r,𝑔𝐔𝑡𝑟subscript𝜁𝑡𝑟g(\mathbf{U},t,r)\zeta_{t,r},italic_g ( bold_U , italic_t , italic_r ) italic_ζ start_POSTSUBSCRIPT italic_t , italic_r end_POSTSUBSCRIPT , (23)

where:

  • 1.

    Gaussian or Hermitian:

    ζt,r=(2β1)t+βK(N1)/2Γ[t+βK(N1)/2]h(2t+r)(0),subscript𝜁𝑡𝑟superscript2superscript𝛽1𝑡𝛽𝐾𝑁12Γdelimited-[]𝑡𝛽𝐾𝑁12superscript2𝑡𝑟0\zeta_{t,r}=(2\beta^{-1})^{t+\beta K(N-1)/2}\Gamma[t+\beta K(N-1)/2]h^{(2t+r)}% (0),italic_ζ start_POSTSUBSCRIPT italic_t , italic_r end_POSTSUBSCRIPT = ( 2 italic_β start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_t + italic_β italic_K ( italic_N - 1 ) / 2 end_POSTSUPERSCRIPT roman_Γ [ italic_t + italic_β italic_K ( italic_N - 1 ) / 2 ] italic_h start_POSTSUPERSCRIPT ( 2 italic_t + italic_r ) end_POSTSUPERSCRIPT ( 0 ) ,

    with h(y)=(2πβ1)βK(N1)/2exp(βy/2)𝑦superscript2𝜋superscript𝛽1𝛽𝐾𝑁12𝛽𝑦2h(y)=(2\pi\beta^{-1})^{-\beta K(N-1)/2}\exp(-\beta y/2)italic_h ( italic_y ) = ( 2 italic_π italic_β start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT - italic_β italic_K ( italic_N - 1 ) / 2 end_POSTSUPERSCRIPT roman_exp ( - italic_β italic_y / 2 ).

  • 2.

    T type I or Pearson type VII:

    ζt,r=βtβK(N1)/2Γ[t+βK(N1)/2]Γ[t+r+βν/2]Γ[2t+r+β(K(N1)+ν)/2]h(2t+r)(0),subscript𝜁𝑡𝑟superscript𝛽𝑡𝛽𝐾𝑁12Γdelimited-[]𝑡𝛽𝐾𝑁12Γdelimited-[]𝑡𝑟𝛽𝜈2Γdelimited-[]2𝑡𝑟𝛽𝐾𝑁1𝜈2superscript2𝑡𝑟0\zeta_{t,r}=\frac{\beta^{-t-\beta K(N-1)/2}\Gamma[t+\beta K(N-1)/2]\Gamma[t+r+% \beta\nu/2]}{\Gamma[2t+r+\beta(K(N-1)+\nu)/2]}h^{(2t+r)}(0),italic_ζ start_POSTSUBSCRIPT italic_t , italic_r end_POSTSUBSCRIPT = divide start_ARG italic_β start_POSTSUPERSCRIPT - italic_t - italic_β italic_K ( italic_N - 1 ) / 2 end_POSTSUPERSCRIPT roman_Γ [ italic_t + italic_β italic_K ( italic_N - 1 ) / 2 ] roman_Γ [ italic_t + italic_r + italic_β italic_ν / 2 ] end_ARG start_ARG roman_Γ [ 2 italic_t + italic_r + italic_β ( italic_K ( italic_N - 1 ) + italic_ν ) / 2 ] end_ARG italic_h start_POSTSUPERSCRIPT ( 2 italic_t + italic_r ) end_POSTSUPERSCRIPT ( 0 ) ,

    with

    h(y)=Γ1β[β(K(N1)+ν)/2](πβ1)βK(N1)/2Γ1β[βν/2](1+βy)β(K(N1)+ν)/2.𝑦superscriptsubscriptΓ1𝛽delimited-[]𝛽𝐾𝑁1𝜈2superscript𝜋superscript𝛽1𝛽𝐾𝑁12superscriptsubscriptΓ1𝛽delimited-[]𝛽𝜈2superscript1𝛽𝑦𝛽𝐾𝑁1𝜈2h(y)=\frac{\Gamma_{1}^{\beta}[\beta(K(N-1)+\nu)/2]}{(\pi\beta^{-1})^{\beta K(N% -1)/2}\Gamma_{1}^{\beta}[\beta\nu/2]}\left(1+\beta y\right)^{-\beta(K(N-1)+\nu% )/2}.italic_h ( italic_y ) = divide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT [ italic_β ( italic_K ( italic_N - 1 ) + italic_ν ) / 2 ] end_ARG start_ARG ( italic_π italic_β start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_β italic_K ( italic_N - 1 ) / 2 end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT [ italic_β italic_ν / 2 ] end_ARG ( 1 + italic_β italic_y ) start_POSTSUPERSCRIPT - italic_β ( italic_K ( italic_N - 1 ) + italic_ν ) / 2 end_POSTSUPERSCRIPT .
  • 3.

    Gegenbauer type I or Pearson type II:

    ζt,r=Γ[t+βK(N1)/2]Γ[r+tβK(N1)/2βq](β)t+βK(N1)/2Γ[2t+rβq]h(2t+r)(0),subscript𝜁𝑡𝑟Γdelimited-[]𝑡𝛽𝐾𝑁12Γdelimited-[]𝑟𝑡𝛽𝐾𝑁12𝛽𝑞superscript𝛽𝑡𝛽𝐾𝑁12Γdelimited-[]2𝑡𝑟𝛽𝑞superscript2𝑡𝑟0\zeta_{t,r}=\frac{\Gamma[t+\beta K(N-1)/2]\Gamma[r+t-\beta K(N-1)/2-\beta q]}{% (-\beta)^{t+\beta K(N-1)/2}\Gamma[2t+r-\beta q]}h^{(2t+r)}(0),italic_ζ start_POSTSUBSCRIPT italic_t , italic_r end_POSTSUBSCRIPT = divide start_ARG roman_Γ [ italic_t + italic_β italic_K ( italic_N - 1 ) / 2 ] roman_Γ [ italic_r + italic_t - italic_β italic_K ( italic_N - 1 ) / 2 - italic_β italic_q ] end_ARG start_ARG ( - italic_β ) start_POSTSUPERSCRIPT italic_t + italic_β italic_K ( italic_N - 1 ) / 2 end_POSTSUPERSCRIPT roman_Γ [ 2 italic_t + italic_r - italic_β italic_q ] end_ARG italic_h start_POSTSUPERSCRIPT ( 2 italic_t + italic_r ) end_POSTSUPERSCRIPT ( 0 ) ,

    with

    h(y)=Γ1β[βK(N1)/2+βq+1](πβ1)βK(N1)/2Γ1β[βq+1](1βy)βq.𝑦superscriptsubscriptΓ1𝛽delimited-[]𝛽𝐾𝑁12𝛽𝑞1superscript𝜋superscript𝛽1𝛽𝐾𝑁12superscriptsubscriptΓ1𝛽delimited-[]𝛽𝑞1superscript1𝛽𝑦𝛽𝑞h(y)=\frac{\Gamma_{1}^{\beta}[\beta K(N-1)/2+\beta q+1]}{(\pi\beta^{-1})^{% \beta K(N-1)/2}\Gamma_{1}^{\beta}[\beta q+1]}\left(1-\beta y\right)^{\beta q}.italic_h ( italic_y ) = divide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT [ italic_β italic_K ( italic_N - 1 ) / 2 + italic_β italic_q + 1 ] end_ARG start_ARG ( italic_π italic_β start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT italic_β italic_K ( italic_N - 1 ) / 2 end_POSTSUPERSCRIPT roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_β end_POSTSUPERSCRIPT [ italic_β italic_q + 1 ] end_ARG ( 1 - italic_β italic_y ) start_POSTSUPERSCRIPT italic_β italic_q end_POSTSUPERSCRIPT .

Finally, Díaz-García and Caro-Lopera (2024b) established a paralellism between the multimatrix variate distributions, studied here, and the multimatricvariate distributions. They share the simplicity that we have addressed for computation of multiple probabilities. Instead of traces, the corresponding determinants can be expanded in terms of zonal polynomials and then a similar computation can be performed for the multimatricvariate distributions termed as Pearson VII, Pearson type II and beta type II in Díaz-García and Caro-Lopera (2022).

4 Application in dynamic molecular docking in SARS-CoV-2

Recently, Díaz-García and Caro-Lopera (2024a) applied the multimatrix variate distributions in a problem of molecular docking, by finding a new cavity of 241 atoms inside the SARS-CoV-2 main protease for placing the inhibitor N3 of 21 atoms. The algorithm of searching was based on a theorem provided by Ramirez et al. (2022) and the main protease of 2387 atoms and the ligand emerged from Jin et al. (2020) and PDB (2020). A multimatrix setting models dependent sample experiments by a realistic estimation based on non independent likelihoods. In this situation, 56 movements of the rigid ligand N3 was recorded into the new pocket meanwhile the equilibrium is reached by a decreasing Lennard-Jones potential type 6126126-126 - 12 and 6106106-106 - 10. Finally, the dependent sample joint distribution was estimated as:

Γ1[(a0+ka)m]Γ1[a0m]Γmk[a]i=1k|𝐅i|am+12(1+i=1ktr𝐅i)(a0+ka)m,subscriptΓ1delimited-[]subscript𝑎0𝑘𝑎𝑚subscriptΓ1delimited-[]subscript𝑎0𝑚superscriptsubscriptΓ𝑚𝑘delimited-[]𝑎superscriptsubscriptproduct𝑖1𝑘superscriptsubscript𝐅𝑖𝑎𝑚12superscript1superscriptsubscript𝑖1𝑘trsubscript𝐅𝑖subscript𝑎0𝑘𝑎𝑚\frac{\Gamma_{1}[(a_{0}+ka)m]}{\Gamma_{1}[a_{0}m]\Gamma_{m}^{k}[a]}\prod_{i=1}% ^{k}|\mathbf{F}_{i}|^{a-\frac{m+1}{2}}\left(1+\displaystyle\sum_{i=1}^{k}% \mathop{\rm tr}\nolimits\mathbf{F}_{i}\right)^{-(a_{0}+ka)m},divide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ( italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_k italic_a ) italic_m ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT [ italic_a ] end_ARG ∏ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT | bold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_a - divide start_ARG italic_m + 1 end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ( 1 + ∑ start_POSTSUBSCRIPT italic_i = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_k end_POSTSUPERSCRIPT roman_tr bold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - ( italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_k italic_a ) italic_m end_POSTSUPERSCRIPT ,

where a0=0.34397,a=0.19735formulae-sequencesubscript𝑎00.34397𝑎0.19735a_{0}=0.34397,a=0.19735italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.34397 , italic_a = 0.19735, m=3,k=56formulae-sequence𝑚3𝑘56m=3,k=56italic_m = 3 , italic_k = 56, 𝐅i=𝐓i𝐓i,i=1,,kformulae-sequencesubscript𝐅𝑖superscriptsubscript𝐓𝑖subscript𝐓𝑖𝑖1𝑘\mathbf{F}_{i}=\mathbf{T}_{i}^{\prime}\mathbf{T}_{i},i=1,\ldots,kbold_F start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = bold_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT bold_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_i = 1 , … , italic_k, for 21×321321\times 321 × 3 matrices 𝐓isubscript𝐓𝑖\mathbf{T}_{i}bold_T start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT registering the spatial coordinates of the ligand inside the protease (see Díaz-García and Caro-Lopera (2024a) for details). We also highlight that the estimation was completed under the full invariant family of spherical distributions. A fact that eliminates the complex problem of previous knowledge and/or fitting of the hidden law. The dependent joint distribution is also crucial here because the molecular docking experiment demands a time dependent calibration which can not reached by the classical (probabilistic independent) likelihood estimation.

The addressed example of Díaz-García and Caro-Lopera (2024a) showed that the joint distribution functions based on spherical multimatrix distributions are easily computable and applied to real data. We now advance into a more complex example involving series of zonal polynomials.

The theorem in Ramirez et al. (2022, Th. Sec. 3.) holds only for rigid molecular docking, it means that the ligand is optimally placed in the pocket by rotations and without distortions. However, molecular biology states that the interaction of one inhibitor inside an active site is dynamic, forcing that the ligand could change its shape by suitable deformations. One way of emulating that molecular dynamic consists of computing the probability that the ligand changes its shape according a desirable coupling structure in the protein. The probability can be modeled by the latent roots of certain fixed positive definite matrix representing the active site near the ligand. For the sake of a simple mathematical illustration, and without any biological and expert study, consider the nearest 21 atoms to the optimized ligand (time 56). Figure 1 shows both groups of atoms. Let 𝐀𝐀\mathbf{A}bold_A the corresponding 3×3333\times 33 × 3 symmetrized matrix of the nearest atoms in the protein.

Refer to caption
Figure 1: Initial setting of the lowest Lennard-Jones potential rigid ligand (last locus 56, green) and the first 21 nearest atoms (magenta) in the selected cavity on SARS-CoV-2 main protease (grey). The symmetrized optimal ligand is apart 0.24884 in positive definite probability from the symmetrized surrounded neighborhood. It is expected to improve the affinity by deforming the ligand in such way that the P(𝟎<𝐅<𝐀)𝑃0𝐅𝐀P(\mathbf{0}<\mathbf{F}<\mathbf{A})italic_P ( bold_0 < bold_F < bold_A ) is greater, then the latent roots information about the ligand shape is near to the geometry of the symmetrized target.

The latent roots of 𝐀𝐀\mathbf{A}bold_A are A1=0.5864,A2=0.2351,A3=0.1785formulae-sequencesubscript𝐴10.5864formulae-sequencesubscript𝐴20.2351subscript𝐴30.1785A_{1}=0.5864,A_{2}=0.2351,A_{3}=0.1785italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0.5864 , italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0.2351 , italic_A start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = 0.1785, a fact that reflects an embracing neighborhood cavity. Meanwhile, the more stable ligand (highly flat) is represented by the symmetric matrix with latent roots L1=0.8663,L2=0.0991,L3=0.0346formulae-sequencesubscript𝐿10.8663formulae-sequencesubscript𝐿20.0991subscript𝐿30.0346L_{1}=0.8663,L_{2}=0.0991,L_{3}=0.0346italic_L start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0.8663 , italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0.0991 , italic_L start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT = 0.0346. Thus, by Theorem 1, the probability that a symmetrized ligand 𝐅𝐅\mathbf{F}bold_F is near to the symmetrized neighborhood 𝐀𝐀\mathbf{A}bold_A in terms of definite positivity is given by:

P(𝟎<𝐅<𝐀)=𝑃0𝐅𝐀absentP(\mathbf{0}<\mathbf{F}<\mathbf{A})=\hskip 284.52756ptitalic_P ( bold_0 < bold_F < bold_A ) =
Γ1[(a0+ka)m]Γm[(m+1)/2]Γ1[a0m]Γm[a+(m+1)/2]|𝐀1|aP11r[((a0+ka)m)r:a;a+(m+1);𝐀]\frac{\Gamma_{1}[(a_{0}+ka)m]\Gamma_{m}[(m+1)/2]}{\Gamma_{1}[a_{0}m]\Gamma_{m}% [a+(m+1)/2]}|\mathbf{A}_{1}|^{a}\,{}_{1}^{r}P_{1}\left[\left((a_{0}+ka)m\right% )_{r}:a;a+(m+1);-\mathbf{A}\right]divide start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ( italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_k italic_a ) italic_m ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ ( italic_m + 1 ) / 2 ] end_ARG start_ARG roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_m ] roman_Γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT [ italic_a + ( italic_m + 1 ) / 2 ] end_ARG | bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ( ( italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_k italic_a ) italic_m ) start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT : italic_a ; italic_a + ( italic_m + 1 ) ; - bold_A ]

In this case m=3𝑚3m=3italic_m = 3, thus the P11r[]subscriptsuperscriptsubscript𝑃1𝑟1delimited-[]{}_{1}^{r}P_{1}[\cdot]start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ⋅ ] can be written as

r=0((a0+ka)m)rQ11r[((a0+ka)m)r:a;a+(m+1);𝐀]r!,\sum_{r=0}^{\infty}\frac{\left((a_{0}+ka)m\right)_{r}\,{}_{1}^{r}Q_{1}\left[% \left((a_{0}+ka)m\right)_{r}:a;a+(m+1);-\mathbf{A}\right]}{r!},∑ start_POSTSUBSCRIPT italic_r = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT divide start_ARG ( ( italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_k italic_a ) italic_m ) start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ( ( italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_k italic_a ) italic_m ) start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT : italic_a ; italic_a + ( italic_m + 1 ) ; - bold_A ] end_ARG start_ARG italic_r ! end_ARG ,

and Q11r[]subscriptsuperscriptsubscript𝑄1𝑟1delimited-[]{}_{1}^{r}Q_{1}\left[\cdot\right]start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [ ⋅ ] is just a summation over partitions of 3 parts which can be easily write down, then the computation of the zonal polynomials of only 3 latent roots can be computed by using Gupta and Richards (1978) combined with the recurrence method of James (1968); or directly by the referred modification of the algorithms for hypergeometric series given by Koev and Edelman (2006). For tracing a path in the probability model, we just sweep the three latent roots δ1,δ2,δ3subscript𝛿1subscript𝛿2subscript𝛿3\delta_{1},\delta_{2},\delta_{3}italic_δ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_δ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_δ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT of a variable 𝐀𝐀\mathbf{A}bold_A, as a decreasing distance from L1,L2,L2subscript𝐿1subscript𝐿2subscript𝐿2L_{1},L_{2},L_{2}italic_L start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, respectively. Explicitly, we take n=1000𝑛1000n=1000italic_n = 1000 triples extracted from the sequences λi,jsubscript𝜆𝑖𝑗\lambda_{i,j}italic_λ start_POSTSUBSCRIPT italic_i , italic_j end_POSTSUBSCRIPT from 00 to LiAisubscript𝐿𝑖subscript𝐴𝑖L_{i}-A_{i}italic_L start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT by (LiAi)/nsubscript𝐿𝑖subscript𝐴𝑖𝑛(L_{i}-A_{i})/n( italic_L start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_A start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) / italic_n, i=1,2,3;j=1,,nformulae-sequence𝑖123𝑗1𝑛i=1,2,3;j=1,\ldots,nitalic_i = 1 , 2 , 3 ; italic_j = 1 , … , italic_n. Then the upper bound definite matrix 𝐀j,j=1,,nformulae-sequencesubscript𝐀𝑗𝑗1𝑛\mathbf{A}_{j},j=1,\ldots,nbold_A start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_j = 1 , … , italic_n, representing the bent or deformed ligand, have the latent roots A1,j=L1λ1,jsubscript𝐴1𝑗subscript𝐿1subscript𝜆1𝑗A_{1,j}=L_{1}-\lambda_{1,j}italic_A start_POSTSUBSCRIPT 1 , italic_j end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - italic_λ start_POSTSUBSCRIPT 1 , italic_j end_POSTSUBSCRIPT, A2,j=L2+λ2,jsubscript𝐴2𝑗subscript𝐿2subscript𝜆2𝑗A_{2,j}=L_{2}+\lambda_{2,j}italic_A start_POSTSUBSCRIPT 2 , italic_j end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT 2 , italic_j end_POSTSUBSCRIPT, A3,j=L3+λ3,j,subscript𝐴3𝑗subscript𝐿3subscript𝜆3𝑗A_{3,j}=L_{3}+\lambda_{3,j},italic_A start_POSTSUBSCRIPT 3 , italic_j end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + italic_λ start_POSTSUBSCRIPT 3 , italic_j end_POSTSUBSCRIPT , respectively. Finally, the 1000 results are depicted in Figure 2, they are increasing probabilities from 𝐀𝟏=𝐅𝟓𝟔subscript𝐀1subscript𝐅56\mathbf{A_{1}}=\mathbf{F_{56}}bold_A start_POSTSUBSCRIPT bold_1 end_POSTSUBSCRIPT = bold_F start_POSTSUBSCRIPT bold_56 end_POSTSUBSCRIPT, which is the rigid optimized ligand in the last time (56) of the dependent sample, to 𝐀𝟏𝟎𝟎𝟎subscript𝐀1000\mathbf{A_{1000}}bold_A start_POSTSUBSCRIPT bold_1000 end_POSTSUBSCRIPT with the specified latent roots A1,A2,A3subscript𝐴1subscript𝐴2subscript𝐴3A_{1},A_{2},A_{3}italic_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_A start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT , italic_A start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT corresponding to the target neighborhood in the protein. The first probability is P(𝟎<𝐅<𝐀1)=0.24884𝑃0𝐅subscript𝐀10.24884P(\mathbf{0}<\mathbf{F}<\mathbf{A}_{1})=0.24884italic_P ( bold_0 < bold_F < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = 0.24884 which is less than the associated probability P(𝟎<𝐅<𝐀1000)=0.37276𝑃0𝐅subscript𝐀10000.37276P(\mathbf{0}<\mathbf{F}<\mathbf{A}_{1000})=0.37276italic_P ( bold_0 < bold_F < bold_A start_POSTSUBSCRIPT 1000 end_POSTSUBSCRIPT ) = 0.37276.

Refer to caption
Figure 2: The probability path of 1000 transitions from the symmetrized optmized rigid and flat ligand to a positive definite close distance of symmetrized 21 surrounding atoms in a cavity found on SARS-CoV-2 main protease. The trayectory starts in P(𝟎<𝐅<𝐀1)=0.24884𝑃0𝐅subscript𝐀10.24884P(\mathbf{0}<\mathbf{F}<\mathbf{A}_{1})=0.24884italic_P ( bold_0 < bold_F < bold_A start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = 0.24884, where the ligand is flat and non-deformed, then by suitable increments of its second and third symmetrized latent roots, and decrement of the first, the probability is increased until P(𝟎<𝐅<𝐀1000)=0.37276𝑃0𝐅subscript𝐀10000.37276P(\mathbf{0}<\mathbf{F}<\mathbf{A}_{1000})=0.37276italic_P ( bold_0 < bold_F < bold_A start_POSTSUBSCRIPT 1000 end_POSTSUBSCRIPT ) = 0.37276. In the last stage the symetrized ligand reaches approximate near latent roots of the symmetrized neigborhood with greater probability, a fact that describes a plausible molecular docking by a non flat rigid coupling.

Then a distorted symmetrized ligand getting more similar in latent roots to the symmetrized target in the pocket appear with greater probability, a fact that can be studied in future as a punctuation of an effective molecular docking under non rigid movements and a plausible calibration with the Lennard-Jones potential.

Finally, the above technique can be extended to two o more probabilities on Theorem 1, it just needs to compute Qi1r[]subscriptsuperscriptsubscript𝑄𝑖𝑟1delimited-[]{}_{1}^{r}Q_{i}\left[\cdot\right]start_FLOATSUBSCRIPT 1 end_FLOATSUBSCRIPT start_POSTSUPERSCRIPT italic_r end_POSTSUPERSCRIPT italic_Q start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT [ ⋅ ] as weights of the corresponding coefficient of the ith𝑖𝑡i-thitalic_i - italic_t italic_h summation.

5 Conclusions

This work revises the problem of probability computations on cones in matrix variate distributions and applies the discussion in a several situations. In particular, the probabilities of multimatrix variate distributions are set in terms of computable series of zonal polynomials. Some of the distributions here applied are invariant under the complete family of elliptically countered distributions, they include the termed Pearson II, Pearson VII. The non invariant case is also considered in the generalised multimatrix variate Wishart distributions. This case promotes a simplification of a classical kernel integral in elliptically contoured distributions which was applied in statistical shape affine distributions under real normed division algebras. Finally, the study of probabilities on cones can be applied in some meaningful situations by understanding the positive definiteness probability in a context of dynamic molecular docking in SARS-CoV-2.

References

  • Arias et al. (2021) Arias-Serna, M.A., Caro-Lopera, F.J., Loubes, J.M. 2021. Risk measures a generalization from the univariate to the matrix variate. Journal of Risk, 23(4), 1–20.
  • Baez (2002) Baez, J. C. 2002. The octonions, Bull. Amer. Math. Soc. 39, 145–205.
  • Caro-Lopera (2016) Caro-Lopera, F. J., 2016. The impossibility of a recurrence construction of the invariant polynomials by using the Laplace-Beltrami operator. Far East J. of Mathematical Sciences 100 (8), 1265–1288.
  • Caro-Lopera (2018) Caro-Lopera, F. J., 2018. Families of computable matrix-variate polynomial distributions and applications. Far East J. of Mathematical Sciences 108 (2), 285–325.
  • Caro-Lopera et al. (2010) Caro-Lopera, F. J., Díaz-García, J.A., González-Farías, 2010. Noncentral elliptical configuration density. J. of Multivariate Analysis 101 (1), 32–43.
  • Caro-Lopera et al. (2013) Caro-Lopera, F. J., González-Farías, Balakrishnan, N. 2013. Determinants, permanents and some applications to statistical shape theory. J. of Multivariate Analysis 114, 29–39.
  • Caro-Lopera et al. (2016) Caro-Lopera, F. J., González-Farías, Balakrishnan, N. 2016. Matrix-variate distribution theory under elliptical models-4: Joint distribution of latent roots of covariance matrix and the largest and smallest latent roots. J. of Multivariate Analysis 145, 224–235.
  • Caro-Lopera and Díaz-García (2012) Caro-Lopera, F. J., Díaz-García, J.A., 2012. Matrix Kummer-Pearson VII Relation and Polynomial Pearson VII Configuration Density. J. of the Iranian Statist. Soc. 11, 217–230.
  • Clay Mathematics Institute of Cambridge (2000) Clay Mathematics Institute of Cambridge. The Millennium Prize Problems. 2000. https://www.claymath.org/millennium-problems/
  • Constantine (1963) Constantine, A.G., 1963. Some non-central distribution problems in multivariate analysis. The Annals of Mathematical Statistics 34, 1270–1285.
  • Davis (1979) Davis, A.W., 1979. Invariant polynomials with two matrix arguments extending the zonal polynomials: Applications to multivariate distribution theory. Ann. Inst. Stat. Math. 31, 465–485.
  • Davis (1980) Davis, A. W., 1980. Invariant polynomials with two matrix arguments, extending the zonal polynomials,Multivariate Analysis—V (ed. P. R. Krishnaiah), 287–299.
  • Díaz-García, and Gutiérrez-Jáimez (2013) Díaz-García, J. A., Gutiérrez-Jáimez, R. 2013. Spherical ensembles. Linear Algebra Appl. 438, 3174-3201.
  • Díaz-García and Caro-Lopera (2016) Díaz-García, J. A., Caro-Lopera, F. J. 2016. Elliptical affine shape distributions for real normed division algebras. J. Multiva. Anal. 144, 139-149.
  • Díaz-García and Caro-Lopera (2022) Díaz-García, J. A., Caro-Lopera, F. J. 2022. Multimatricvariate distribution under elliptical models. J. Stat. Plann. Infer. 2022, 109-117.
  • Díaz-García and Caro-Lopera (2024a) Díaz-García, J. A., Caro-Lopera, F. J. 2024a. Multimatrix variate distributions. arXiv:2405.02498. Also submitted.
  • Díaz-García and Caro-Lopera (2024b) Díaz-García, J. A., Caro-Lopera, F. J. 2024b. Multimatricvariate and multimatrix variate distributions based on elliptically contoured laws under real normed division algebras. arXiv:2405.06905. Also submitted.
  • Díaz-García et al. (2022) Díaz-García, J. A., Caro-Lopera, F. J., Pérez Ramírez, F. O. 2022. Multivector variate distributions: An application in Finance. Sankhyā 84-A,Part 2, 534-555.
  • Dray and Manogue (1999) Dray, T., Manogue, C. A. 1999. The exceptional Jordan eigenvalue problem, Inter. J. Theo. Phys. 38(11), 2901–2916.
  • Edelman and Rao (2005) Edelman, A. Rao, R. R: 2005. Random matrix theory, Acta Numer. 14, 233–297.
  • Fang et al. (1990) Fang, K. T., Zhang, Y. T., Ng, K. W. Symmetric Multivariate and realted distributions. Springer-Science+Business Media, B. V., New Delhi, 1990.
  • Goodall and Mardia (1993) Goodall, K. I. Mardia, K. V., 1993. Multivariate aspects of shape theory, Ann. Statist. 21, 848-866.
  • Gupta and Richards (1978) Gupta, R. D., Richards, D., 1978. Calculation of zonal polynomials of 3×3333\times 33 × 3 positive definite symmetric matrices. Ann. Inst. Statist. Math. 31, Part A, 207–213.
  • Gross and Richards (1987) Gross, K. I., Richards, D. ST, P. 1987, Special functions of matrix argument I: Algebraic induction zonal polynomials and hypergeometric functions, Trans. Amer. Math. Soc. 301(2), 475–501.
  • Gupta and Varga (1993) Gupta, A. K., Varga, T., 1993 Elliptically Contoured Models in Statistics, Kluwer Academic Publishers, Dordrecht.
  • James (1960) James, A. T., 1960. The Distribution of the latent roots of the covariance matrix, Ann. Math. Statist. 31(1): 151-158.
  • James (1964) James, A. T., 1964. Distributions of matrix variates and latent roots derived from normal samples, Ann. Math. Statist.,35, 475–501.
  • James (1968) James, A.T., 1968. Calculation of Zonal Polynomial Coefficients by Use of the Laplace-Beltrami Operator. The Annals of Mathematical Statistics, 39 (5), 1711–1718.
  • Jin et al. (2020) Jin, Z., Du, X., Xu, Y. et al. Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors. Nature 582, 289–293 (2020). https://doi.org/10.1038/s41586-020-2223-y
  • Koev and Edelman (2006) Koev, E., Edelman, A., 2006. The efficient evaluation of the hypergeometric function of a matrix argument. Mathematics of Computation, 75 (254), 833–846.
  • Muirhead (2005) Muirhead, R. J., 2005. Aspects of Multivariate Statistical Theory. John Wiley & Sons, New York.
  • PDB (2020) PDB - 6LU7. Liu, X., Zhang, B., Jin, Z., Yang, H., Rao, Z. The crystal structure of COVID-19 main protease in complex with an inhibitor N3. PDB DOI: https://doi.org/10.2210/pdb6lu7/pdb Primary publication DOI: https://doi.org/10.1038/s41586-020-2223-y
  • Ramirez et al. (2022) Ramirez-Velasquez, I., Bedoya-Calle, A., Velez, E., Caro-Lopera, F., 2022. Shape Theory Applied to Molecular Docking and Automatic Localization of Ligand Binding Pockets in Large Proteins. ACS Omega, 7(50), 45991-46002 doi: 10.1021/acsomega.2c02227