Giá trị mong đợi của yếu tố xác định log của ma trận Wishart


16

Hãy , tức là được phân phối theo một D × D phân phối chiều Wishart với trung bình ν Ψ và bậc tự do ν . Tôi muốn một biểu thức cho E ( log | Λ | ) trong đó | Λ | là yếu tố quyết định.ΛWD(ν,Ψ)D×DνΨνE(log|Λ|)|Λ|

Tôi đã google một chút cho câu trả lời cho điều này và đã nhận được một số thông tin mâu thuẫn. Bài viết này nói rõ rằng trong đóbiểu thị hàm digamma; bài báo không cung cấp một nguồn cho sự thật này như tôi có thể nói. Đây cũng là công thức được sử dụng trêntrang wikipedia cho Wishart, trang web có văn bản Nhận dạng mẫu của Giám mục.

E(log|Λ|)=Dlog2+log|Ψ|+i=1Dψ(νi+12)
dψ()ddxlogΓ(x)

Mặt khác, google đã đưa ra cuộc thảo luận này với một bài viết được liên kết nói rằng Họ kết luận bằng cách nói rằng có nguồn gốc từ thực tế là . Tôi đã kiểm tra tính toán này bắt đầu từ và có vẻ ổn, nhưng chúng tôi có thêm .E ( log | Λ | ) = D log 2 - D log ν + log | Ψ | + D Σ i = 1 ψ ( ν - i + 1

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

Câu trả lời:


9

Khi tôi đã sẵn sàng để đăng bài này, tôi đã có thể trả lời câu hỏi của riêng tôi. Theo nguyên tắc chung StackExchange, tôi đã quyết định đăng nó bằng mọi cách với hy vọng rằng ai đó gặp phải vấn đề này có thể tìm thấy vấn đề này trong tương lai, có thể sau khi gặp vấn đề tương tự với các nguồn mà tôi đã làm. Tôi đã quyết định trả lời ngay lập tức để không ai lãng phí thời gian cho nó vì giải pháp không thú vị.

()

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

()

BIÊN TẬP:

Λ= =dΨ1/2LLTΨ1/2LN(0,1)χν-Tôi+12,(Tôi= =1,...,D)()


2
Tôi thích phiên bản Cholesky hơn - bạn có căn bậc hai của chi-vuông trên đường chéo và tiêu chuẩn bình thường ở tam giác dưới.
xác suất

@probabilityislogic Cảm ơn vì tiền boa! Ghi nhớ nó như thế có vẻ dễ dàng và hữu ích hơn.
anh chàng

Này, tôi đang cố gắng rút ra kỳ vọng về nhật ký Wishart (được nêu trong cuốn sách của Giám mục), điều đó có vẻ phức tạp, bạn có tìm thấy nguồn nào để rút ra kết quả không?
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.