an:06835781
Zbl 1417.65032
Cowles, Mary Kathryn; Bonett, Stephen; Seedorff, Michael
Independent sampling for Bayesian normal conditional autoregressive models with OpenCL acceleration
EN
Comput. Stat. 33, No. 1, 159-177 (2018).
0943-4062 1613-9658
2018
j
65C60 62F15 62P10
CAR; \texttt{CARrampsOcl}; Dirichlet density; Kronecker sum; rejection sampling; separability
Summary: A new computational strategy produces independent samples from the joint posterior distribution for a broad class of Bayesian spatial and spatiotemporal conditional autoregressive models. The method is based on reparameterization and marginalization of the posterior distribution and massive parallelization of rejection sampling using graphical processing units (GPUs) or other accelerators. It enables very fast sampling for small to moderate-sized datasets (up to approximately 10,000 observations) and feasible sampling for much larger datasets. Even using a mid-range GPU and a high-end CPU, the GPU-based implementation is up to 30 times faster than the same algorithm run serially on a single CPU, and the numbers of effective samples per second are orders of magnitude higher than those obtained with popular Markov chain Monte Carlo software. The method has been implemented in the \texttt{R} package \texttt{CARrampsOcl}. This work provides both a practical computing strategy for fitting a popular class of Bayesian models and a proof of concept that GPU acceleration can make independent sampling from Bayesian joint posterior densities feasible.