Erwarteter Wert der Log-Determinante einer Wishart-Matrix

16

Sei , dh verteilt nach einer D × D- dimensionalen Wishart-Verteilung mit dem Mittelwert ν Ψ und Freiheitsgraden ν . Ich möchte einen Ausdruck für E ( log | Λ | ) wo | Λ |ΛWD(ν,Ψ)D×DνΨνE(Log|Λ|)|Λ|ist die Determinante.

Ich habe ein bisschen nach der Antwort auf diese Frage gesucht und einige widersprüchliche Informationen erhalten. In diesem Artikel wird ausdrücklich angegeben, dass wobeiψ()die Digamma-Funktiond bezeichnet

E(Log|Λ|)=DLog2+Log|Ψ|+ich=1Dψ(ν-ich+12)
ψ() ; Soweit ich das beurteilen kann, gibt die Zeitung keine Quelle für diese Tatsache an. Dies ist auch die Formel, die auf derWikipedia-Seite für den Wishart verwendet wirdddxLogΓ(x) wird.

Auf der anderen Seite hat Google diese Diskussion mit einem verlinkten Artikel aufgegriffen , der besagt, dass Sie schließen mit der Aussage, dass E ( log | Λ | ) = D log 2 - D log ν + log | Ψ | + D i = 1 ψ ( ν - i + 1

νD|Λ||Ψ|χν2χν-12χν-D+12.()
die sich aus der Tatsache ergibt, dassE(logχ 2 ν )=log(2)+ψ(ν/2). Ich habe diese Berechnung ab()und es scheintOrdnung, aber wir haben eine extra-Dlogν.
E(Log|Λ|)=DLog2-DLogν+Log|Ψ|+ich=1Dψ(ν-ich+12)
E(Logχν2)=Log(2)+ψ(ν/2)()-DLogν
Kerl
quelle

Antworten:

9

Als ich mich darauf vorbereitete, dies zu posten, konnte ich meine eigene Frage beantworten. In Übereinstimmung mit der allgemeinen StackExchange-Etikette habe ich beschlossen, sie trotzdem zu veröffentlichen, in der Hoffnung, dass jemand anderes, der auf dieses Problem stößt, dies in Zukunft möglicherweise findet, nachdem er möglicherweise auf dieselben Probleme mit Quellen gestoßen ist, die ich getan habe. Ich habe beschlossen, es sofort zu beantworten, damit niemand Zeit damit verschwendet, da die Lösung nicht interessant ist.

() ist falsch, weil das in der Diskussion verlinkte Papier eine andere Parametrisierung des Wisharts verwendete; Dies wurde von den Diskutanten nicht bemerkt. Was wir eigentlich haben sollten, ist

|Λ||Ψ|χν2χν-12χν-D+12.()
Nach dieser Korrektur führen die beiden Formeln zu derselben Antwort.

Ich denke jedenfalls () ist eine interessante Beziehung.

BEARBEITEN:

Nach probabilisticlogics Rat können wir schreiben Λ=dΨ1/2LLTΨ1/2 wo niedriger dreieckig L hat N(0,1) Elemente aus der Diagonale und χν-ich+12,(ich=1,...,D)Elemente auf der Diagonale. Die Determinante beider Seiten nehmen gibt() sofort.

Kerl
quelle
2
Mir gefällt die Cholesky-Version besser - Sie haben die Quadratwurzel von Chi-Quadrat auf der Diagonale und die Normale auf dem unteren Dreieck.
Wahrscheinlichkeitslogik
@probabilityislogic Danke für den Tipp! Sich daran zu erinnern, scheint einfacher und nützlicher.
Kerl
Hey, ich versuche, die Erwartung des Protokolls Wishart (in Bishops Buch angegeben) abzuleiten, das kompliziert aussieht. Haben Sie eine Quelle gefunden, um das Ergebnis abzuleiten?
Avocado