Although it is common practice to fit a complex Bayesian model using Markov chain Monte Carlo (MCMC) methods, we provide an alternative sampling-based method to fit a two-stage hierarchical model in ...
When Gaussian errors are inappropriate in a multivariate linear regression setting, it is often assumed that the errors are iid from a distribution that is a scale mixture of multivariate normals.