Skip to main content
Log in

MALA with annealed proposals: a generalization of locally and globally balanced proposal distributions

  • Published:
Statistics and Computing Aims and scope Submit manuscript

Abstract

We introduce a generalized version of the Metropolis-adjusted Langevin algorithm (MALA). The informed proposal distribution of this new sampler features two tuning parameters: the usual step size parameter \(\sigma ^2\) and an interpolation parameter \(\gamma \) that may be adjusted to accommodate the dimension of the target distribution. We theoretically study the efficiency of the sampler by making use of the local- and global-balance concepts introduced in Zanella (JASA 115:852–865, 2020) and provide efficient tuning guidelines that work well with a variety of target distributions. Although the usual MALA (\(\gamma =1\)) is shown to be optimal for infinite-dimensional targets, in practice, the generalized MALA (\(1<\gamma \le 2\)) remains the most appealing option, even in high-dimensional contexts. Simulation studies and numerical experiments are presented to illustrate our findings. We apply the new sampler to a Bayesian logistic regression context and show that its efficiency compares favourably to competing algorithms.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6

Similar content being viewed by others

Data Availability Statement

Data sets used in §4 are publicly available at http://archive.ics.uci.edu (German Credit, Australian Credit, Heart) and in the package MASS on R (Pima Indian).

References

  • Boisvert-Beaudry, G.: Efficacité des distributions instrumentales en équilibre dans un algorithme de type Metropolis-Hastings. Université de Montréal, Thesis (2019)

  • Geyer, C.J.: Practical Markov chain Monte Carlo. Stat. Sci. 7(4), 473–483 (1992)

    Google Scholar 

  • Girolami, M., Calderhead, B.: Riemann manifold Langevin and Hamiltonian Monte Ccarlo methods. J. R. Stat. Soc. Ser. B 73(2), 123–214 (2011)

    Article  MathSciNet  Google Scholar 

  • Hastings, W.K.: Monte Carlo sampling methods using Markov chains and their applications. Biometrika 57, 97–109 (1970)

    Article  MathSciNet  Google Scholar 

  • Metropolis, N., Rosenbluth, A.W., Rosenbluth, M.N., Teller, A.H., Teller, E.: Equation of state calculations by fast computing machines. J. Chem. Phys. 21(6), 1087–1092 (1953)

    Article  Google Scholar 

  • Michie, D., Spiegelhalter, D.J., Taylor, C.C., Campbell, J. (eds.): Machine learning, neural and statistical classification. Ellis Horwood, Upper Saddle River, NJ, USA (1994)

    MATH  Google Scholar 

  • Peskun, P.H.: Optimum Monte-Carlo sampling using Markov chains. Biometrika 60(3), 607–612 (1973)

    Article  MathSciNet  Google Scholar 

  • Ripley, B.D., Hjort, N.: Pattern recognition and neural networks. Cambridge University Press, Cambridge (1996)

    Book  Google Scholar 

  • Roberts, G., Rosenthal, J.: Optimal scaling for various Metropolis-Hastings algorithms. Stat. Sci. 16(4), 351–367 (2001)

    Article  MathSciNet  Google Scholar 

  • Roberts, G.O., Gelman, A., Gilks, W.R.: Weak convergence and optimal scaling of random walk Metropolis algorithms. Ann. Appl. Probabil. 7, 110–120 (1997)

    MathSciNet  MATH  Google Scholar 

  • Roberts, G.O., Rosenthal, J.S.: Optimal scaling of discrete approximations to Langevin diffusions. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 60(1), 255–268 (1998)

    Article  MathSciNet  Google Scholar 

  • Tierney, L.: A note on Metropolis-Hastings kernels for general state spaces. Annals Appl. Probabil. , 1–9 (1998)

  • Zanella, G.: Informed proposals for local MCMC in discrete spaces. J. Am. Stat. Assoc. 115(530), 852–865 (2020)

    Article  MathSciNet  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mylène Bédard.

Additional information

The authors are grateful to the two anonymous reviewers whose suggestions and comments led to an improved manuscript. They also acknowledge the support of the Natural Sciences and Engineering Research Council of Canada

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Boisvert-Beaudry, G., Bédard, M. MALA with annealed proposals: a generalization of locally and globally balanced proposal distributions. Stat Comput 32, 5 (2022). https://doi.org/10.1007/s11222-021-10063-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1007/s11222-021-10063-1

Keywords

Mathematics Subject Classification

Navigation