Expected value of the log-determinant of a Wishart matrix

Let ΛWD(ν,Ψ), i.e. distributed according to a D×D dimensional Wishart distribution with mean νΨ and degrees of freedom ν. I would like an expression for E(log|Λ|) where |Λ| is the determinant.

I’ve google’d a bit for the answer to this and have gotten some conflicting information. This paper explicitly states that
E(log|Λ|)=Dlog2+log|Ψ|+Di=1ψ(νi+12)
where ψ() denotes the digamma function ddxlogΓ(x); the paper does not give a source for this fact as far as I can tell. This is also the formula used on the wikipedia page for the Wishart, which sites Bishop’s Pattern Recognition text.

On the other hand, google turned up this discussion with a linked paper that states that
νD|Λ||Ψ|χ2νχ2ν1χ2νD+1.()
They conclude by stating that
E(log|Λ|)=Dlog2Dlogν+log|Ψ|+Di=1ψ(νi+12)
which is derived using the fact that E(logχ2ν)=log(2)+ψ(ν/2). I checked this calculation starting from () and it seems okay, but we have an extra Dlogν.

Answer

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νχ2ν1χ2νD+1.()
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 χ2νi+1,(i=1,...,D) elements on the diagonal. Taking the determinent of both sides gives () immediately.

Attribution
Source : Link , Question Author : guy , Answer Author : guy

Leave a Comment