Bernoulli

An adaptive Metropolis algorithm

Heikki Haario, Eero Saksman, and Johanna Tamminen
Source: Bernoulli Volume 7, Number 2 (2001), 223-242.

Abstract

A proper choice of a proposal distribution for Markov chain Monte Carlo methods, for example for the Metropolis-Hastings algorithm, is well known to be a crucial factor for the convergence of the algorithm. In this paper we introduce an adaptive Metropolis (AM) algorithm, where the Gaussian proposal distribution is updated along the process using the full information cumulated so far. Due to the adaptive nature of the process, the AM algorithm is non-Markovian, but we establish here that it has the correct ergodic properties. We also include the results of our numerical tests, which indicate that the AM algorithm competes well with traditional Metropolis-Hastings algorithms, and demonstrate that the AM algorithm is easy to use in practical computation.

First Page: Show Hide
Full-text: Open access
Links and Identifiers

Permanent link to this document: http://projecteuclid.org/euclid.bj/1080222083
Mathematical Reviews number (MathSciNet): MR1828504
Zentralblatt MATH identifier: 0989.65004

References

[1] Davidson, J. and de Jong, R. (1997) Strong laws of large numbers for dependent heterogeneous processes: a synthesis of recent and new results. Econometric Rev., 16, 251-279. Abstract can also be found in the ISI/STMA publication
Mathematical Reviews (MathSciNet): MR1463083
Zentralblatt MATH: 0903.60025
Digital Object Identifier: doi:10.1080/07474939708800387
[2] Dobrushin, R. (1956) Central limit theorems for non-stationary Markov chains II. Theory Probab. Appl., 1, 329-383.
Mathematical Reviews (MathSciNet): MR86436
[3] Evans, M. (1991) Chaining via annealing. Ann. Statist., 19, 382-393. Abstract can also be found in the ISI/STMA publication
Mathematical Reviews (MathSciNet): MR1091858
Zentralblatt MATH: 0739.65119
Digital Object Identifier: doi:10.1214/aos/1176347989
Project Euclid: euclid.aos/1176347989
[4] Fishman, G.S. (1996) Monte Carlo: Concepts, Algorithms and Applications. New York: Springer- Verlag.
Mathematical Reviews (MathSciNet): MR1392474
[5] Gelfand, A.E. and Sahu, S.K. (1994) On Markov chain Monte Carlo acceleration. J. Comput. Graph. Statist., 3, 261-276. Abstract can also be found in the ISI/STMA publication
Mathematical Reviews (MathSciNet): MR1292118
Digital Object Identifier: doi:10.2307/1390911
[6] Gelman, A.G., Roberts, G.O. and Gilks, W.R. (1996) Efficient Metropolis jumping rules. In J.M. Bernardo, J.O. Berger, A.F. David and A.F.M. Smith (eds), Bayesian Statistics V, pp. 599-608. Oxford: Oxford University Press.
Mathematical Reviews (MathSciNet): MR1425429
[7] Gilks, W.R. and Roberts, G.O. (1995) Strategies for improving MCMC. In W.R. Gilks, S. Richardson and D.J. Spiegelhalter (eds), Markov Chain Monte Carlo in Practice, pp. 75-88. London: Chapman & Hall.
[8] Gilks, W.R., Roberts, G.O. and George, E.I. (1994) Adaptive direction sampling. The Statistician, 43, 179-189. Abstract can also be found in the ISI/STMA publication
Digital Object Identifier: doi:10.2307/2348942
[9] Gilks, W.R., Richardson, S. and Spiegelhalter, D.J. (1995) Introducing Markov chain Monte Carlo. In W.R. Gilks, S. Richardson and D.J. Spiegelhalter (eds), Markov Chain Monte Carlo in Practice, pp. 1-19. London: Chapman & Hall.
[10] Gilks, W.R., Roberts, G.O. and Sahu, S.K. (1998) Adaptive Markov chain Monte Carlo. J. Amer. Statist. Assoc., 93, 1045-1054. Abstract can also be found in the ISI/STMA publication
Mathematical Reviews (MathSciNet): MR1649199
Zentralblatt MATH: 1064.65503
Digital Object Identifier: doi:10.2307/2669848
[11] Haario, H. and Saksman, E. (1991) Simulated annealing process in general state space. Adv. Appl. Probab., 23, 866-893. Abstract can also be found in the ISI/STMA publication
Mathematical Reviews (MathSciNet): MR1133733
Zentralblatt MATH: 0744.60104
Digital Object Identifier: doi:10.2307/1427681
[12] Haario, H., Saksman, E. and Tamminen, J. (1999) Adaptive proposal distribution for random walk Metropolis algorithm. Comput. Statist., 14, 375-395. Abstract can also be found in the ISI/STMA publication
Digital Object Identifier: doi:10.1007/s001800050022
[13] Hall, P. and Heyde, C.C. (1980) Martingale Limit Theory and Its Application. New York, Academic Press.
Mathematical Reviews (MathSciNet): MR624435
Zentralblatt MATH: 0462.60045
[14] Hastings, W.K. (1970) Monte Carlo sampling methods using Markov chains and their applications. Biometrika, 57, 97-109.
Zentralblatt MATH: 0219.65008
Digital Object Identifier: doi:10.1093/biomet/57.1.97
[15] McLeish, D.L. (1975) A maximal inequality and dependent strong laws. Ann. Probab., 3, 829-839.
Mathematical Reviews (MathSciNet): MR400382
Digital Object Identifier: doi:10.1214/aop/1176996269
[16] Metropolis, N., Rosenbluth, A.W., Rosenbluth, M.N., Teller, A.H. and Teller, E. (1953) Equations of state calculations by fast computing machines. J. Chem. Phys., 21, 1087-1091.
Digital Object Identifier: doi:10.1063/1.1699114
[17] Neveu, J. (1965) Mathematical Foundations of the Calculus of Probability. San Francisco: Holden-Day.
Mathematical Reviews (MathSciNet): MR198505
Zentralblatt MATH: 0137.11301
[18] Nummelin, E. (1984) General Irreducible Markov Chains and Non-negative Operators. Cambridge: Cambridge University Press.
Mathematical Reviews (MathSciNet): MR776608
[19] Roberts, G.O., Gelman, A. and Gilks, W.R. (1997) Weak convergence and optimal scaling of random walk Metropolis algorithms. Ann. Appl. Probab., 7, 110-120. Abstract can also be found in the ISI/STMA publication
Mathematical Reviews (MathSciNet): MR1428751
Zentralblatt MATH: 0876.60015
Digital Object Identifier: doi:10.1214/aoap/1034625254
Project Euclid: euclid.aoap/1034625254
[20] Sahu, S.K. and Zhigljavsky, A.A. (1999) Self regenerative Markov chain Monte Carlo with adaptation. Preprint. http://www.statslab.cam.ac.uk/ mcmc.
Mathematical Reviews (MathSciNet): MR1997490
Digital Object Identifier: doi:10.3150/bj/1065444811
Project Euclid: euclid.bj/1065444811
[21] Tierney, L. (1994) Markov chains for exploring posterior distributions (with discussion). Ann. Statist., 22, 1701-1762. Abstract can also be found in the ISI/STMA publication
Mathematical Reviews (MathSciNet): MR1329166
Zentralblatt MATH: 0829.62080
Digital Object Identifier: doi:10.1214/aos/1176325750
Project Euclid: euclid.aos/1176325750

2012 © Bernoulli Society for Mathematical Statistics and Probability

Bernoulli

Bernoulli

Turn MathJax Off
What is MathJax?