Nilai yang diharapkan dari penentu log dari matriks Wishart

16

Mari , yaitu didistribusikan menurut D × D distribusi Wishart dimensi dengan mean ν Ψ dan derajat kebebasan ν . Saya ingin ungkapan untuk E ( log | Λ | ) di mana | Λ | adalah penentu.ΛWD(ν,Ψ)D×DνΨνE(catatan|Λ|)|Λ|

Saya telah sedikit Google 'untuk jawaban untuk ini dan mendapatkan beberapa informasi yang saling bertentangan. Makalah ini secara eksplisit menyatakan bahwa manaψ()menunjukkan fungsi digammad

E(catatan|Λ|)=Dcatatan2+catatan|Ψ|+saya=1Dψ(ν-saya+12)
ψ(); makalah tidak memberikan sumber untuk fakta ini sejauh yang saya tahu. Ini juga merupakan formula yang digunakan padahalaman wikipedia untuk Wishart, yang memuat teks Pengenalan Pola Uskup.ddxcatatanΓ(x)

Di sisi lain, google membuka diskusi ini dengan makalah tertaut yang menyatakan bahwa Mereka menyimpulkan dengan menyatakan bahwa E ( log | Λ | ) = D log 2 - D log ν + log | Ψ | + D i = 1 ψ ( ν - i + 1

νD|Λ||Ψ|χν2χν-12χν-D+12.()
yang diturunkan menggunakan fakta bahwaE(logχ 2 ν )=log(2)+ψ(ν/2). Saya memeriksa perhitungan ini mulai dari()dan sepertinya oke, tapi kami punya ekstra-Dlogν.
E(log|Λ|)=Dlog2Dcatatanν+catatan|Ψ|+saya=1Dψ(ν-saya+12)
E(catatanχν2)=catatan(2)+ψ(ν/2)()-Dcatatanν
orang
sumber

Jawaban:

9

Ketika saya bersiap-siap untuk memposting ini, saya dapat menjawab pertanyaan saya sendiri. Sesuai dengan etiket StackExchange umum, saya telah memutuskan untuk mempostingnya dengan harapan bahwa orang lain yang mengalami masalah ini dapat menemukan ini di masa depan, mungkin setelah mengalami masalah yang sama dengan sumber yang saya lakukan. Saya telah memutuskan untuk menjawabnya segera sehingga tidak ada yang membuang waktu karena solusinya tidak menarik.

()

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

()

EDIT:

Λ=dΨ1/2L.L.TΨ1/2L.N(0,1)χν-saya+12,(saya=1,...,D)elemen diagonal. Mengambil determinasi dari kedua belah pihak memberi() segera.

orang
sumber
2
Saya suka versi Cholesky lebih baik - Anda memiliki akar kuadrat dari chi-square pada diagonal dan standar normal pada segitiga bawah.
probabilityislogic
@probabilityislogic Terima kasih atas tipnya! Mengingatnya seperti itu tampaknya lebih mudah dan lebih bermanfaat.
pria
Hei, saya sedang mencoba untuk mendapatkan harapan dari log Wishart (dinyatakan dalam buku Bishop), yang terlihat rumit, apakah Anda menemukan sumber untuk mendapatkan hasilnya?
alpukat