Scalable Rejection Sampling for Bayesian Hierarchical Models
Description
Functions for implementing the Braun and Damien (2015) rejection
sampling algorithm for Bayesian hierarchical models. The algorithm generates
posterior samples in parallel, and is scalable when the individual units are
conditionally independent.