So, I have:
f(mu,R|B,m,W,v) = Normal(mu|m,inv(B*R)).Wishart(R|W,v)
Using the property of the Wishart distribution that the inverse of R (the covariance matrix, C) has distribution:
C ~ Inverse-Wishart(V, v)
where V = W^-1
and the property of the inverse-Wishart distribution that if C and V are partitioned conformably with each other (in my case, partitioned into components for each of the 2 dimensions), then
C_11 ~ Inverse-Wishart(V_11, v-1)
Which for scalar C_11, V_11 reduces to
C_11 ~ Inverse-Gamma((v-1)/2, (V_11)/2)
Together with the property of the normal distribution that under such a partition (again, by component), if mu ~ Normal(m, C) then
mu_1 ~ Normal(m_1,C_11)
We obtain a distribution for mu_1, C_11:
f(mu_1,C_11|v,W,B,m) ~ Normal(mu_1|m_1, (1/B)*C_11).Inverse-Gamma(C_11|(v-1)/2, (V_11)/2)
I think...




Reply With Quote