Scalable Rejection Sampling for Bayesian Hierarchical Models
提出一种无需MCMC的贝叶斯模型后验抽样方法，样本独立可并行计算，适用于大规模数据，还能计算边际似然。
Bayesian hierarchical modeling is a popular approach to capturing unobserved heterogeneity across individual units. However, standard estimation methods such as Markov chain Monte Carlo (MCMC) can be impracticable for modeling outcomes from a large number of units. We develop a new method to sample from posterior distributions of Bayesian models, without using MCMC. Samples are independent, so they can be collected in parallel, and we do not need to be concerned with issues like chain convergence and autocorrelation. The algorithm is scalable under the weak assumption that individual units are conditionally independent, making it applicable for large data sets. It can also be used to compute marginal likelihoods. Data, as supplemental material, are available at http://dx.doi.org/10.1287/mksc.2014.0901 .