scispace - formally typeset
Open AccessJournal ArticleDOI

Robust adaptive Metropolis algorithm with coerced acceptance rate

Matti Vihola
- 01 Sep 2012 - 
- Vol. 22, Iss: 5, pp 997-1008
Reads0
Chats0
TLDR
A new robust adaptive Metropolis algorithm estimating the shape of the target distribution and simultaneously coercing the acceptance rate and showing promising behaviour in an example with Student target distribution having no finite second moment.
Abstract
The adaptive Metropolis (AM) algorithm of Haario, Saksman and Tamminen (Bernoulli 7(2):223---242, 2001) uses the estimated covariance of the target distribution in the proposal distribution. This paper introduces a new robust adaptive Metropolis algorithm estimating the shape of the target distribution and simultaneously coercing the acceptance rate. The adaptation rule is computationally simple adding no extra cost compared with the AM algorithm. The adaptation strategy can be seen as a multidimensional extension of the previously proposed method adapting the scale of the proposal distribution in order to attain a given acceptance rate. The empirical results show promising behaviour of the new algorithm in an example with Student target distribution having no finite second moment, where the AM covariance estimate is unstable. In the examples with finite second moments, the performance of the new approach seems to be competitive with the AM algorithm combined with scale adaptation.

read more

Content maybe subject to copyright    Report

Citations
More filters

Bayesian Filtering and Smoothing

Simo Särkkä
TL;DR: This compact, informal introduction for graduate students and advanced undergraduates presents the current state-of-the-art filtering and smoothing methods in a unified Bayesian framework and learns what non-linear Kalman filters and particle filters are, how they are related, and their relative advantages and disadvantages.

贝叶斯滤波与平滑 (Bayesian filtering and smoothing)

Simo Särkkä
TL;DR: This compact, informal introduction for graduate students and advanced undergraduates presents the current state-of-the-art filtering and smoothing methods in a unified Bayesian framework and learns what non-linear Kalman filters and particle filters are, how they are related, and their relative advantages and disadvantages.
Book

Bayesian Filtering and Smoothing

TL;DR: This compact, informal introduction for graduate students and advanced undergraduates presents the current state-of-the-art filtering and smoothing methods in a unified Bayesian framework, learning what non-linear Kalman filters and particle filters are, how they are related, and their relative advantages and disadvantages.
Journal ArticleDOI

Flexible and scalable methods for quantifying stochastic variability in the era of massive time-domain astronomical data sets

TL;DR: It is shown that the PSD of a CARMA model can be expressed as a sum of Lorentzian functions, which makes them extremely flexible and able to model a broad range of PSDs.
Journal ArticleDOI

Bayesian computation: a summary of the current state, and samples backwards and forwards

TL;DR: The difficulties of modelling and then handling ever more complex datasets most likely call for a new type of tool for computational inference that dramatically reduces the dimension and size of the raw data while capturing its essential aspects.
References
More filters
BookDOI

Markov Chain Monte Carlo in Practice

TL;DR: The Markov Chain Monte Carlo Implementation Results Summary and Discussion MEDICAL MONITORING Introduction Modelling Medical Monitoring Computing Posterior Distributions Forecasting Model Criticism Illustrative Application Discussion MCMC for NONLINEAR HIERARCHICAL MODELS.
Book

Monte Carlo Statistical Methods

TL;DR: This new edition contains five completely new chapters covering new developments and has sold 4300 copies worldwide of the first edition (1999).
Journal ArticleDOI

An adaptive Metropolis algorithm

TL;DR: An adaptive Metropolis (AM) algorithm, where the Gaussian proposal distribution is updated along the process using the full information cumulated so far, which establishes here that it has the correct ergodic properties.
Related Papers (5)