I assume the distributions in Figure 1 are calculated according to Eqn. 2.19 and Section 2.5, from the data they fitted to their multivariate linear model according to their procedure. After a few minutes looking through the paper (available here: http://www.isds.duke.edu/courses/Fall05/sta395/joelucas1.pdf) they also mention that ‘the posterior distribution of ß_j is a mixture of scaled and shifted t-distributions.’
So it would seem that you will have to find (or simulate) some multivariate data, scale it as they describe, and fit it according to their procedure. Then use their methods to calculate the distribution. It doesn't appear (in my brief reading of the paper, anyway) to have an analytic expression.
It's an approach to multivariate regression ('Bayesian' as compared to the more common 'Frequentist') that I wasn't aware of before (since I don't often do multivariate regression). I learned something.
