Valeur attendue du log-déterminant d'une matrice de Wishart

16

Soit ΛWD(ν,Ψ) , c'est-à-dire distribué selon une distribution de Wishart dimensionnelle avec la moyenne et les degrés de liberté . Je voudrais une expression pour oùest le déterminant.ν Ψ ν E ( log | Λ | ) | Λ |D×DνΨνE(Journal|Λ|)|Λ|

J'ai recherché un peu sur Google la réponse à cette question et j'ai obtenu des informations contradictoires. Ce document indique explicitement que où désigne la fonction digamma ; le papier ne donne pas une source pour ce fait pour autant que je puisse dire. C'est également la formule utilisée sur la page wikipedia pour le Wishart , qui contient le texte de reconnaissance des motifs de Bishop.ψ()d

E(Journal|Λ|)=Journal2+Journal|Ψ|+je=1ψ(ν-je+12)
ψ()XJournalΓ(X)

D'un autre côté, google a lancé cette discussion avec un document lié qui indique que Ils concluent en déclarant que qui est dérivé en utilisant le fait que . J'ai vérifié ce calcul à partir de (\ dagger) et cela semble correct, mais nous avons un -D \ log \ nu supplémentaire .E ( log | Λ | ) = D log 2 - D log ν + log | Ψ | + D i = 1 ψ ( ν - i + 1

ν|Λ||Ψ|χν2χν-12χν-+12.()
E(Journal|Λ|)=Journal2-Journalν+Journal|Ψ|+je=1ψ(ν-je+12)
E(logχν2)=log(2)+ψ(ν/2)()Dlogν
gars
la source

Réponses:

9

Alors que je m'apprêtais à poster ceci, j'ai pu répondre à ma propre question. Conformément à l'étiquette générale de StackExchange, j'ai décidé de le publier de toute façon dans l'espoir que quelqu'un d'autre qui rencontre ce problème puisse le trouver à l'avenir, peut-être après avoir rencontré les mêmes problèmes avec des sources que moi. J'ai décidé d'y répondre immédiatement afin que personne n'y perde de temps car la solution n'est pas intéressante.

est faux, car le document lié à la discussion utilisait une paramétrisation différente du Wishart; cela n'a pas été remarqué par les présentateurs. Ce que nous devrions réellement avoir, c'est |() Après cette correction, les deux formules conduisent à la même réponse.

|Λ||Ψ|χν2χν12χνD+12.()

En tout cas, je pense que est une relation intéressante.()

ÉDITER:

À la suite de l'avis de probabilisticlogic on peut écrire où triangulaire inférieure L a N ( 0 , 1 ) des éléments hors de la diagonale et Λ=dΨ1/2LLTΨ1/2LN(0,1)éléments sur la diagonale. Prendre le déterminant des deux côtés donne()immédiatement.χνi+12,(i=1,...,D)()

gars
la source
2
J'aime mieux la version Cholesky - vous avez la racine carrée du chi carré sur la diagonale et la normale standard sur le triangle inférieur.
probabilityislogic
@probabilityislogic Merci pour le conseil! Se souvenir de cela comme ça semble plus facile et plus utile.
gars
Hé, j'essaie de dériver l'attente du journal Wishart (indiqué dans le livre de Bishop), qui semble compliqué, avez-vous trouvé une source pour dériver le résultat?
avocat