ค่าที่คาดหวังของตัวกำหนดล็อกของเมทริกซ์ Wishart


16

ให้คือกระจายตามD × Dกระจายมิติริชาร์ตที่มีค่าเฉลี่ยν Ψและองศาอิสระν ฉันต้องการนิพจน์สำหรับE ( log | Λ | )โดยที่| Λ | เป็นตัวกำหนดΛWD(ν,Ψ)D×DνΨνE(log|Λ|)|Λ|

ฉันขอคำตอบเล็กน้อยจาก google และได้รับข้อมูลที่ขัดแย้งกันบ้าง บทความนี้ระบุว่า โดยที่ψ()หมายถึงฟังก์ชั่น digammad

E(log|Λ|)=Dlog2+log|Ψ|+i=1Dψ(νi+12)
ψ()ddxlogΓ(x) ; กระดาษไม่ได้ให้แหล่งข้อมูลสำหรับข้อเท็จจริงนี้เท่าที่ฉันสามารถบอกได้ นี่เป็นสูตรที่ใช้ในหน้าวิกิพีเดียสำหรับ Wishartซึ่งเป็นเว็บไซต์ที่ข้อความการจดจำรูปแบบของ Bishop

ในทางกลับกัน Google ก็เปิดการสนทนานี้ด้วยเอกสารที่เชื่อมโยงซึ่งระบุว่า

νD|Λ||Ψ|χν2χν12χνD+12.()
E(log|Λ|)=Dlog2Dlogν+log|Ψ|+i=1Dψ(νi+12)
E(logχν2)=log(2)+ψ(ν/2)()Dlogν

คำตอบ:


9

As I was getting ready to post this, I was able to answer my own question. In accordance with general StackExchange etiquette I've decided to post it anyways in hopes that someone else who runs into this problem might find this in the future, possibly after running into the same issues with sources that I did. I've decided to answer it immediately so that no one wastes time on it since the solution isn't interesting.

() is wrong, because the paper linked to in the discussion was using a different parametrization of the Wishart; this wasn't noticed by the discussants. What we should actually have is

|Λ||Ψ|χν2χν12χνD+12.()
After this correction, the two formulas lead to the same answer.

At any rate, I think think () is an interesting relationship.

EDIT:

Following probabilisticlogic's advice we can write Λ=dΨ1/2LLTΨ1/2 where lower triangular L has N(0,1) elements off the diagonal and χνi+12,(i=1,...,D) elements on the diagonal. Taking the determinent of both sides gives () immediately.


2
I like the Cholesky version better - you have square root of chi-square on the diagonal and standard normal on the lower triangle.
probabilityislogic

@probabilityislogic Thanks for the tip! Remembering it like that seems easier and more useful.
guy

Hey I'm trying to derive the expectation of the log Wishart (stated in Bishop's book), that looks complicated, did you find any source to derive the result?
avocado
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.