An adaptive Metropolis algorithm
Haario, Heikki ; Saksman, Eero ; Tamminen, Johanna
Bernoulli, Tome 7 (2001) no. 6, p. 223-242 / Harvested from Project Euclid
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.
Publié le : 2001-04-14
Classification:  adaptive Markov chain Monte Carlo,  comparison,  convergence,  ergodicity,  Markov chain Monte Carlo,  Metropolis-Hastings algorithm
@article{1080222083,
     author = {Haario, Heikki and Saksman, Eero and Tamminen, Johanna},
     title = {An adaptive Metropolis algorithm},
     journal = {Bernoulli},
     volume = {7},
     number = {6},
     year = {2001},
     pages = { 223-242},
     language = {en},
     url = {http://dml.mathdoc.fr/item/1080222083}
}
Haario, Heikki; Saksman, Eero; Tamminen, Johanna. An adaptive Metropolis algorithm. Bernoulli, Tome 7 (2001) no. 6, pp.  223-242. http://gdmltest.u-ga.fr/item/1080222083/