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. denotes a matrix with
rows and columns; is the transpose matrix, and if
has an inverse, it will be denoted by . is a symmetric matrix if . If all their eigenvalues are positive then is a positive definite
matrix, a fact denoted as . An identity matrix will be denoted by
. To specify the size of the identity, we will use . denotes the trace of matrix . If then by we mean the vector formed by stacking the columns of
under each other; that is, if , where for ,
The Frobenius norm of a matrix will be denoted as . Typically the
Frobenius norm is denoted as , to differentiate it from other matrix norms. Since
we will use only the Frobenius norm, it just be denoted as . It is
defined by
Finally, denotes the Stiefel manifold, the space of all matrices
() with orthogonal columns, that is, . In addition, if
defines an invariant measure on the Stiefel manifold
, from Theorem 2.1.15, p. 70 in Muirhead (2005),
(1)
where denotes the multivariate Gamma function, see Muirhead (2005, Definiton 2.1.10, p. 61)
Now, let random matrix with a matrix variate elliptical
distribution with respect to the Lebesgue measure , see Gupta and Varga (1993). Therefore its density function is
given by
(2)
The location parameter is ; and the scale parameters
and , are positive definite
matrices. The distribution shall be denoted by , and indexed by the kernel function , where .
When , and a special case of a matrix variate elliptical distribution appears, in this case it is
said that has a matrix variate spherical distribution.
Note that for constant , the substitution and in Fang et al. (1990, Equation 2.21, p. 26) provides that
(3)
For some applications we will require James zonal polynomials, notated by , see Muirhead (2005). They are an orthogonal base expansion in . Here is positive definite and the summation is over all ordered partitions , into not more than parts, such that .
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):
(4)
where the function is independent of and , with , and the and are complex numbers. If , then the well known hypergeometric series is obtained, see for example Muirhead (2005).
We also introduce the notation
(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 has a matrix variate
spherical distribution, with , , . Define and , .
The termed multimatrix variate Pearson type VII is the marginal density
of and is given by
(6)
where , and .
Now define , .
The multimatrix variate beta type II distribution
is given by
(7)
Lemma 2
Suppose that has a matrix variate
spherical distribution, with , , . Define and , . The multimatrix variate Pearson type II distribution is
(8)
Moreover, Assuming that and with , the multimatrix variate beta type I distribution is written as
(9)
Finally
Lemma 3
Assume that has a matrix variate spherical distribution, with , , . Define and
, .
Then, the joint density
is given by
(10)
where . This distribution shall be termed multimatrix variate generalised Gamma - generalised Wishart
distribution.
The joint density function
is not invariant under spherical elliptical functions but it can be written in following closed form:
(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 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 by , by and by . Where the 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 , where 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 , because the partitions are trivially truncated by parts, but exact expressions for with arbitrary 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 of arbitrary 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 in the expansion of is just the number of partitions, just imagine the problem of providing an exact formula for . 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 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 ; 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 () and beta () distribution. The Wishart lower bound probability and the integral for matrix beta function back to Constantine (1963) (See Arias et al. (2021) for the ). However, as Constantine (1963) states: “the complementary probabilities (, and ) 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 (). 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 is a positive definite matrix for all constant positive definite matrices , .
Theorem 1
Assume the hypothesis of Lemma 1 and suitable parameters for convergence of series (4). For positive definite matrices , , the probability is given by:
where the nested summation is given
by
(12)
Here (12) involves sums of type (4) and depends on the indexes and matrices ,
where ,, , .
Proof 1
Let ,
, and , . Consider backwards integration starting in . To avoid the appearance of generalized binomial expansions of zonal polynomials of sums of matrices, the integrand is written as
Performing a convergent power expansion of the form and a variable substitution , with Jacobian , lead to the application of Muirhead (2005, Th. 7.2.10) and (4).
Thus,
(13)
Repeat a similar procedure to (13) for the next integral on by the corresponding splitting of the term indexed by ; thus the nested series is obtained. Explicitly, partial integration on takes the form
Induction on , by a similar integration procedure leading to (13), provides the nested series for the probability on . Finally, the last integral in is given by:
As usual in Wishart type probabilities, if is the largest latent root of , then Theorem 1 provides the distribution function of , , by taking 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 . 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 , then
Excepting the extra summation, the procedure for Theorem 1 can be reproduced next for probabilities on positive definite matrices , .
Theorem 2
Consider the hypothesis of Lemma 3 and positive definite matrices , . Then we have the following probabilities:
1) Gaussian (independent) case:
(14)
2) General elliptical model, (Caro-Lopera et al. (2016)):
(15)
3) General elliptical model, for all :
(16)
Proof 2
(14) is derived by using the generator function
in (11), and hence the integrals follows by 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 . It requires a recurrent use of the binomial theorem, the representation (5) and the known integral
Lemma 4
In the setting of Lemma 3 and the cone probabilities here derived, the only feasible lower probability is reached for and is given by:
(17)
For the probability turns in terms of invariant polynomials and it cannot be computable.
Proof 3
Let and , such that and . The integrand in (9) simplifies to . Then, the result (17) follows after binomial theorem and application of
For 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 , where , stands for Real; , for Complex; , for Quaternionic; and , for Octonion. Other notations for the algebras are given by , 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 be the linear space of all matrices of rank over a real finite-dimensional normed division algebra with distinct positive singular values. Let be the set of all matrices over . The dimension of
over is . Let , then denotes the usual conjugate
transpose.
The set of matrices such that
is a manifold denoted , is termed the Stiefel manifold ( is also known as
semi-orthogonal (), semi-unitary (), semi-symplectic
() and semi-exceptional type () matrices, see Dray and Manogue (1999)). The
dimension of over is . In
particular, with dimension over , , is
the maximal compact subgroup of and
consists of all matrices such that
. Therefore, is the
real orthogonal group (), the unitary group
(), compact symplectic group ()
or exceptional type matrices (), for ,
, or , respectively.
We denote by the real vector space of all such that . Let
be the cone of positive definite matrices ; then is an open subset of . Over , consist of symmetric matrices;
over , Hermitian matrices; over , quaternionic
Hermitian matrices (also termed self-dual matrices) and over ,
octonionic Hermitian matrices. Generically, the elements of
are termed
Hermitian matrices, irrespective of the nature of . The
dimension of over is .
Let be the diagonal subgroup of
consisting of all , .
For any matrix , denotes the
matrix of differentials . Finally, we define the measure or volume element
when ,
or .
If then (the Lebesgue measure in
) denotes the exterior product of the functionally
independent variables
If (or )
then (the Lebesgue measure in or in
) denotes the exterior product of the functionally
independent variables (or denotes the exterior product of the functionally
independent variables, if for all )
The context generally establishes the conditions on the elements of , that is, if
, , or . It is
considered that
Observe, too, that for the Lebesgue measure defined thus, it is required that
, that is, must be a non singular
Hermitian matrix (Hermitian positive definite matrix).
If then (the Legesgue
measure in ) denotes the exterior product of the
functionally independent variables
If then
where . It can be
proved that this differential form does not depend on the choice of the
matrix. When ; defines the unit sphere in
. This is, of course, an - dimensional surface in
. When and denoting by ,
is termed the Haar measure on .
The surface area or volume of the Stiefel manifold is
(18)
where denotes the multivariate Gamma function for the space
, and is defined by
where , denotes the determinant and , see Gross and Richards (1987). If then by we mean the vector formed by stacking the columns of
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 have the same affine shape if , for some translation and . Then the configuration matrix compressing the meaningful geometrical information of the original matrix , is obtained in the sequence of filtering geometrical stages . Here and 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 has the following multivariate elliptically contoured distribution for real normed division algebras respect to the Lebesgue measure:
(19)
where , , , and generator function , satisfies . This fact will be denoted by Finally, for a convergent Taylor series of , if , the affine shape density of is given by
(20)
with
(21)
and
(22)
, , .
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 general derivative of the kernel function . Inspired by Theorem 2, we can avoid integration and non null evaluation of the derivative in (22) by a simple computation of . 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 for Gaussian, Pearson type VII and II models is given by
(23)
where:
1.
Gaussian or Hermitian:
with .
2.
T type I or Pearson type VII:
with
3.
Gegenbauer type I or Pearson type II:
with
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 and . Finally, the dependent sample joint distribution was estimated as:
where , , , for matrices 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 the corresponding symmetrized matrix of the nearest atoms in the protein.
The latent roots of are , a fact that reflects an embracing neighborhood cavity. Meanwhile, the more stable ligand (highly flat) is represented by the symmetric matrix with latent roots .
Thus, by Theorem 1, the probability that a symmetrized ligand is near to the symmetrized neighborhood in terms of definite positivity is given by:
In this case , thus the can be written as
and 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 of a variable , as a decreasing distance from , respectively. Explicitly, we take triples extracted from the sequences from to by , . Then the upper bound definite matrix , representing the bent or deformed ligand, have the latent roots , , respectively. Finally, the 1000 results are depicted in Figure 2, they are increasing probabilities from , which is the rigid optimized ligand in the last time (56) of the dependent sample, to with the specified latent roots corresponding to the target neighborhood in the protein.
The first probability is which is less than the associated probability .
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 as weights of the corresponding coefficient of the 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 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