Fonction gamma multidimensionnelle

En analyse, la fonction gamma multivariée, Γ_p(·), est une généralisation de la fonction gamma. En statistique multivariée, elle apparait dans la fonction de densité de la loi de Wishart et de la loi de Wishart inverse.

Définitions

La définition formelle de la fonction gamma multivariée est, pour tout complexe a tel que $\mathrm{\Re }e(a)>\frac{p+1}{2}$:

${\mathrm{\Gamma }}_p(a)=\underset{S>0}{\int }{\left|S\right|}^{a-(p+1)/2}{\mathrm{e}}^{-\mathrm{t}\mathrm{r}\mathrm{a}\mathrm{c}\mathrm{e}\mathrm{(}\mathrm{S}\mathrm{)}}\mathrm{d}S,$

où Modèle:Nobr signifie que S est une matrice définie positive.

En pratique, on utilise

${\mathrm{\Gamma }}_p(a)={\pi }^{p(p-1)/4}\prod _{j=1}^p\mathrm{\Gamma }[a+(1-j)/2].$

Le calcul est facilité par les relations de récurrence :

$\begin{array}{cc}{\mathrm{\Gamma }}_p(a)& ={\pi }^{(p-1)/2}\mathrm{\Gamma }(a){\mathrm{\Gamma }}_{p-1}(a-\frac{1}{2})\\ & ={\pi }^{(p-1)/2}{\mathrm{\Gamma }}_{p-1}(a)\mathrm{\Gamma }(a+\frac{1-p}{2}).\end{array}$

Ainsi,

• ${\mathrm{\Gamma }}_1(a)=\mathrm{\Gamma }(a)$
• ${\mathrm{\Gamma }}_2(a)={\pi }^{1/2}\mathrm{\Gamma }(a)\mathrm{\Gamma }(a-\frac{1}{2})$
• ${\mathrm{\Gamma }}_3(a)={\pi }^{3/2}\mathrm{\Gamma }(a)\mathrm{\Gamma }(a-\frac{1}{2})\mathrm{\Gamma }(a-1)$

etc.

Dérivation

On définit la fonction digamma multivariée  :

${\psi }_p(a)=\frac{\partial \log {\mathrm{\Gamma }}_p(a)}{\partial a}=\sum _{i=1}^p\psi (a+\frac{1-i}{2}),$

et la fonction polygamma généralisée :

${\psi }_p^{(n)}(a)=\frac{{\partial }^n\log {\mathrm{\Gamma }}_p(a)}{\partial a^n}=\sum _{i=1}^p{\psi }^{(n)}(a+\frac{1-i}{2}).$

Références

• Modèle:Article
• Modèle:DLMF