Abstract
We derive explicit bounds for the computation of normalizing constants Z for log-concave densities π = e−U/Z w.r.t. the Lebesgue measure on Rd. Our approach relies on a Gaussian annealing combined with recent and precise bounds on the Unadjusted Langevin Algorithm [15]. Polynomial bounds in the dimension d are obtained with an exponent that depends on the assumptions made on U. The algorithm also provides a theoretically grounded choice of the annealing sequence of variances. A numerical experiment supports our findings. Results of independent interest on the mean squared error of the empirical average of locally Lipschitz functions are established.
| Original language | English |
|---|---|
| Pages (from-to) | 851-889 |
| Number of pages | 39 |
| Journal | Electronic Journal of Statistics |
| Volume | 12 |
| Issue number | 1 |
| DOIs | |
| Publication status | Published - 1 Jan 2018 |
Keywords
- Annealed importance sampling
- Bayes factor
- Normalizing constants
- Unadjusted langevin algorithm