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
:
![{\displaystyle \Gamma _{p}(a)=\int _{S>0}\left|S\right|^{a-(p+1)/2}\ \mathrm {e} ^{-{\rm {trace(S)}}}\mathrm {d} S,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3fef94413bc77f2f32ef6e683c873c1500050a8a)
où S > 0 signifie que S est une matrice définie positive.
En pratique, on utilise
![{\displaystyle \Gamma _{p}(a)=\pi ^{p(p-1)/4}\prod _{j=1}^{p}\Gamma \left[a+(1-j)/2\right].}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4dc9bcc58164ece0de7aef7b0ba2b1a0925a1a04)
Le calcul est facilité par les relations de récurrence :
![{\displaystyle {\begin{aligned}\Gamma _{p}(a)&=\pi ^{(p-1)/2}\Gamma (a)\Gamma _{p-1}(a-{\tfrac {1}{2}})\\&=\pi ^{(p-1)/2}\Gamma _{p-1}(a)\Gamma \left(a+{\tfrac {1-p}{2}}\right).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f317d0b65051c899237dbb6faf3710dedb31fc5e)
Ainsi,
![{\displaystyle \Gamma _{1}(a)=\Gamma (a)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e357ec080e62e5edd491f28ae54bd8d80cb0f82f)
![{\displaystyle \Gamma _{2}(a)=\pi ^{1/2}\Gamma (a)\Gamma (a-{\tfrac {1}{2}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1d679eeb6fb4029605c29f10f7d37bd268fe8a8a)
![{\displaystyle \Gamma _{3}(a)=\pi ^{3/2}\Gamma (a)\Gamma (a-{\tfrac {1}{2}})\Gamma (a-1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cd3e4895c1ff93c0f74dd98a1a3ffb0e7c365944)
etc.
Dérivation
On définit la fonction digamma multivariée :
![{\displaystyle \psi _{p}(a)={\frac {\partial \log \Gamma _{p}(a)}{\partial a}}=\sum _{i=1}^{p}\psi \left(a+{\frac {1-i}{2}}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/75473d58b841d97234a34eb6d1f3589fc51f4e57)
et la fonction polygamma généralisée :
![{\displaystyle \psi _{p}^{(n)}(a)={\frac {\partial ^{n}\log \Gamma _{p}(a)}{\partial a^{n}}}=\sum _{i=1}^{p}\psi ^{(n)}\left(a+{\frac {1-i}{2}}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ae0ec26c033c706e815b4356ded6873e1abe61fb)
Démonstration
Vu que
![{\displaystyle \Gamma _{p}(a)=\pi ^{p(p-1)/4}\prod _{j=1}^{p}\Gamma \left(a+{\frac {1-j}{2}}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a494f794a92067af72b2af98e467dbc8d9d5b96)
on tire
![{\displaystyle {\frac {\partial \Gamma _{p}(a)}{\partial a}}=\pi ^{p(p-1)/4}\sum _{i=1}^{p}{\frac {\partial \Gamma \left(a+{\frac {1-i}{2}}\right)}{\partial a}}\prod _{j=1,j\neq i}^{p}\Gamma \left(a+{\frac {1-j}{2}}\right).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aa54823ac1e5ecb10ffa7529e3b654a7bb6b891e)
Par définition de la fonction digamma ψ,
![{\displaystyle {\frac {\partial \Gamma (a+{\frac {1-i}{2}})}{\partial a}}=\psi \left(a+{\frac {1-i}{2}}\right)\Gamma \left(a+{\frac {1-i}{2}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a7d843df5028cf63c3b1d78f288df1884a939b9)
il s'ensuit que
![{\displaystyle {\begin{aligned}{\frac {\partial \Gamma _{p}(a)}{\partial a}}&=\pi ^{p(p-1)/4}\prod _{j=1}^{p}\Gamma \left(a+{\frac {1-j}{2}}\right)\sum _{i=1}^{p}\psi \left(a+{\frac {1-i}{2}}\right)\\&=\Gamma _{p}(a)\sum _{i=1}^{p}\psi \left(a+{\frac {1-i}{2}}\right).\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/93708ce635b31db65c604dcddca30fc486c77b33)
Références
Portail des probabilités et de la statistique