A Level-Set Hit-and-Run Sampler for Quasi-Concave Distributions
Abstract
We develop a new sampling strategy that uses the hit-and-run algorithm within level sets of the target density. Our method can be applied to any quasi-concave density, which covers a broad class of models. Our sampler performs well in high-dimensional settings, which we illustrate with a comparison to Gibbs sampling on a spike-and-slab mixture model. We also extend our method to exponentially-tilted quasi-concave densities, which arise often in Bayesian models consisting of a log-concave likelihood and quasi-concave prior density. Within this class of models, our method is effective at sampling from posterior distributions with high dependence between parameters, which we illustrate with a simple multivariate normal example. We also implement our level-set sampler on a Cauchy-normal model where we demonstrate the ability of our level set sampler to handle multi-modal posterior distributions.
- Publication:
-
arXiv e-prints
- Pub Date:
- February 2012
- arXiv:
- arXiv:1202.4094
- Bibcode:
- 2012arXiv1202.4094F
- Keywords:
-
- Statistics - Computation