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 ...