Articulo de referencia

Función gamma multivariada

En matemáticas , la función gamma multivariada Γ p es una generalización de la función gamma . Es útil en estadística multivariada , apareciendo en la función de densidad de pro...

En matemáticas , la función gamma multivariada Γ p es una generalización de la función gamma . Es útil en estadística multivariada , apareciendo en la función de densidad de probabilidad de las distribuciones de Wishart e inversa de Wishart , y en la distribución beta matricial . [ 1 ]

Tiene dos definiciones equivalentes. Una se da como la siguiente integral sobre lapag×pag{\displaystyle p\times p}matrices reales definidas positivas :

Γpag(a)=S>0exp(tr(S))|S|apag+12dS,{\displaystyle \Gamma _{p}(a)=\int _{S>0}\exp \left(-{\rm {tr}}(S)\right)\,\left|S\right|^{a-{\frac {p+1}{2}}}dS,}

dónde|S|{\displaystyle |S|}denota el determinante deS{\displaystyle S}La otra, más útil para obtener un resultado numérico, es:

Γpag(a)=πpag(pag1)/4j=1pagΓ(a+(1j)/2).{\displaystyle \Gamma _ {p}(a)=\pi ^{p(p-1)/4}\prod _ {j=1}^{p}\Gamma (a+(1-j)/2).}

En ambas definiciones,a{\displaystyle a}es un número complejo cuya parte real satisface(a)>(pag1)/2{\displaystyle \Re (a)>(p-1)/2}. Tenga en cuenta queΓ1(a){\displaystyle \Gamma _{1}(a)}se reduce a la función gamma ordinaria. La segunda de las definiciones anteriores permite obtener directamente las relaciones recursivas parapag2{\displaystyle p\geq 2}:

Γpag(a)=π(pag1)/2Γ(a)Γpag1(a12)=π(pag1)/2Γpag1(a)Γ(a+(1pag)/2).{\displaystyle \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 (a+(1-p)/2).}

De este modo

  • Γ2(a)=π1/2Γ(a)Γ(a1/2){\displaystyle \Gamma _{2}(a)=\pi ^{1/2}\Gamma (a)\Gamma (a-1/2)}
  • Γ3(a)=π3/2Γ(a)Γ(a1/2)Γ(a1){\displaystyle \Gamma _{3}(a)=\pi ^{3/2}\Gamma (a)\Gamma (a-1/2)\Gamma (a-1)}

etcétera.

Esto también se puede extender a valores no enteros depag{\displaystyle p}con la expresión:

Γpag(a)=πpag(pag1)/4GRAMO(a+12)GRAMO(a+1)GRAMO(a+1pag2)GRAMO(a+1pag2){\displaystyle \Gamma _{p}(a)=\pi ^{p(p-1)/4}{\frac {G(a+{\frac {1}{2}})G(a+1)}{G(a+{\frac {1-p}{2}})G(a+1-{\frac {p}{2}})}}}

Donde G es la función G de Barnes , el producto indefinido de la función Gamma .

La función es derivada por Anderson [ 2 ] a partir de primeros principios, quien también cita trabajos anteriores de Wishart , Mahalanobis y otros.

También existe una versión de la función gamma multivariada que, en lugar de un único número complejo, toma unpag{\displaystyle p}vector -dimensional de números complejos como su argumento. Generaliza la función gamma multivariada definida anteriormente en la medida en que esta última se obtiene mediante una elección particular del argumento multivariado de la primera. [ 3 ]

Derivados

Podemos definir la función digamma multivariada como

ψpag(a)=registroΓpag(a)a=i=1pagψ(a+(1i)/2),{\displaystyle \psi _{p}(a)={\frac {\partial \log \Gamma _{p}(a)}{\partial a}}=\sum _{i=1}^{p}\psi (a+(1-i)/2),}

y la función poligamma general como

ψpag(norte)(a)=norteregistroΓpag(a)anorte=i=1pagψ(norte)(a+(1i)/2).{\displaystyle \psi _{p}^{(n)}(a)={\frac {\partial ^{n}\log \Gamma _{p}(a)}{\partial a^{n}}}=\sum _{i=1}^{p}\psi ^{(n)}(a+(1-i)/2).}

Pasos del cálculo

  • Desde
Γpag(a)=πpag(pag1)/4j=1pagΓ(a+1j2),{\displaystyle \Gamma _{p}(a)=\pi ^{p(p-1)/4}\prod _{j=1}^{p}\Gamma \left(a+{\frac {1-j}{2}}\right),}
resulta que
Γpag(a)a=πpag(pag1)/4i=1pagΓ(a+1i2)aj=1,jipagΓ(a+1j2).{\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).}
Γ(a+(1i)/2)a=ψ(a+(i1)/2)Γ(a+(i1)/2){\displaystyle {\frac {\partial \Gamma (a+(1-i)/2)}{\partial a}}=\psi (a+(i-1)/2)\Gamma (a+(i-1)/2)}
resulta que
Γpag(a)a=πpag(pag1)/4j=1pagΓ(a+(1j)/2)i=1pagψ(a+(1i)/2)=Γpag(a)i=1pagψ(a+(1i)/2).{\displaystyle {\begin{aligned}{\frac {\partial \Gamma _{p}(a)}{\partial a}}&=\pi ^{p(p-1)/4}\prod _{j=1}^{p}\Gamma (a+(1-j)/2)\sum _{i=1}^{p}\psi (a+(1-i)/2)\\[4pt]&=\Gamma _{p}(a)\sum _{i=1}^{p}\psi (a+(1-i)/2).\end{aligned}}}

Referencias

  1. James, Alan T. (junio de 1964). "Distribuciones de variables matriciales y raíces latentes derivadas de muestras normales" . The Annals of Mathematical Statistics . 35 (2): 475– 501. doi : 10.1214/aoms/1177703550 . ISSN 0003-4851 . 
  2. Anderson, TW (1984). Introducción al análisis estadístico multivariante . Nueva York: John Wiley and Sons. págs. Cap. 7. ISBN  0-471-88987-3.
  3. D. St. P. Richards (s.f.). "Capítulo 35 Funciones de argumento matricial" . Biblioteca digital de funciones matemáticas . Consultado el 23 de mayo de 2022 .
  • 1. James, A. (1964). "Distribuciones de variables matriciales y raíces latentes derivadas de muestras normales" . Annals of Mathematical Statistics . 35 (2): 475– 501. doi : 10.1214 /aoms/1177703550 . MR 0181057. Zbl 0121.36605 .  
  • 2. AK Gupta y DK Nagar 1999. "Distribuciones matriciales variables". Chapman and Hall.
Obtenido de " https://en.wikipedia.org/w/index.php?title=Multivariate_gamma_function&oldid=1353126675 "