Valor esperado del log-determinante de una matriz Wishart


16

Sea , es decir, distribuido según una distribución dimensional de Wishart con media y grados de libertad . Me gustaría una expresión para dondeEs el determinante.E ( log | Λ | ) | Λ |ΛWre(ν,Ψ)re×reνΨνmi(Iniciar sesiónEl |ΛEl |)El |ΛEl |

Busqué en Google un poco la respuesta a esto y obtuve información contradictoria. Este documento declara explícitamente que donde denota la función digamma ; el documento no proporciona una fuente para este hecho por lo que puedo decir. Esta es también la fórmula utilizada en la página de Wikipedia para Wishart , que ubica el texto de Bishop's Pattern Recognition.

mi(Iniciar sesiónEl |ΛEl |)=reIniciar sesión2+Iniciar sesiónEl |ΨEl |+yo=1reψ(ν-yo+12)
ψ()rereXIniciar sesiónΓ(X)

Por otro lado, Google presentó esta discusión con un documento vinculado que dice que Concluyen afirmando que que se deriva del hecho de que . Verifiqué este cálculo a partir de y parece estar bien, pero tenemos un adicional .E ( log | Λ | ) = D log 2 - D log ν + log | Ψ | + D i = 1 ψ ( ν - i + 1

νreEl |ΛEl |El |ΨEl |χν2χν-12χν-re+12.()
mi(Iniciar sesiónEl |ΛEl |)=reIniciar sesión2-reIniciar sesiónν+Iniciar sesiónEl |ΨEl |+yo=1reψ(ν-yo+12)
mi(Iniciar sesiónχν2)=Iniciar sesión(2)+ψ(ν/ /2)()-reIniciar sesiónν

Respuestas:


9

Cuando me estaba preparando para publicar esto, pude responder mi propia pregunta. De acuerdo con la etiqueta general de StackExchange, he decidido publicarlo de todos modos con la esperanza de que otra persona que se encuentre con este problema pueda encontrarlo en el futuro, posiblemente después de encontrarse con los mismos problemas con las fuentes que hice. He decidido responderlo de inmediato para que nadie pierda el tiempo, ya que la solución no es interesante.

() está mal, porque el artículo vinculado en la discusión estaba usando una parametrización diferente de Wishart; los comentaristas no lo notaron. Lo que deberíamos tener es Después de esta corrección, las dos fórmulas conducen a la misma respuesta.

El |ΛEl |El |ΨEl |χν2χν-12χν-re+12.()

En cualquier caso, creo que es una relación interesante.()

EDITAR:

Siguiendo el consejo de probabilisticlogic podemos escribir donde triangular inferior tiene elementos fuera de la diagonal y elementos en la diagonal. Tomar el determinante de ambos lados da inmediatamente.Λ=reΨ1/ /2LLTΨ1/ /2Lnorte(0 0,1)χν-yo+12,(yo=1,...,re)()


2
Me gusta más la versión Cholesky: tienes una raíz cuadrada de chi-cuadrado en la diagonal y normal estándar en el triángulo inferior.
probabilidadislogic

@probabilityislogic ¡Gracias por el consejo! Recordarlo así parece más fácil y más útil.
chico

Hola, estoy tratando de deducir la expectativa del registro Wishart (indicado en el libro de Bishop), que parece complicado, ¿encontraste alguna fuente para obtener el resultado?
aguacate
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.