103
v1v2v3 (latest)

Notes on the runtime of A* sampling

Abstract

The challenge of simulating random variables is a central problem in Statistics and Machine Learning. Given a tractable proposal distribution PP, from which we can draw exact samples, and a target distribution QQ which is absolutely continuous with respect to PP, the A* sampling algorithm allows simulating exact samples from QQ, provided we can evaluate the Radon-Nikodym derivative of QQ with respect to PP. Maddison et al. originally showed that for a target distribution QQ and proposal distribution PP, the runtime of A* sampling is upper bounded by O(exp(D[QP]))\mathcal{O}(\exp(D_{\infty}[Q||P])) where D[QP]D_{\infty}[Q||P] is the Renyi divergence from QQ to PP. This runtime can be prohibitively large for many cases of practical interest. Here, we show that with additional restrictive assumptions on QQ and PP, we can achieve much faster runtimes. Specifically, we show that if QQ and PP are distributions on R\mathbb{R} and their Radon-Nikodym derivative is unimodal, the runtime of A* sampling is O(D[QP])\mathcal{O}(D_{\infty}[Q||P]), which is exponentially faster than A* sampling without assumptions.

View on arXiv
Comments on this paper