Metropolis-Hastings algorithms with adaptive proposals

被引:30
|
作者
Cai, Bo [1 ]
Meyer, Renate [2 ]
Perron, Francois [3 ]
机构
[1] Univ S Carolina, Dept Epidemiol & Biostat, Columbia, SC 29208 USA
[2] Univ Auckland, Dept Stat, Auckland 1, New Zealand
[3] Univ Montreal, Dept Math & Stat, Montreal, PQ H3C 3J7, Canada
基金
加拿大自然科学与工程研究理事会;
关键词
Adaptive rejection Metropolis sampling; Bayesian inference; Markov chain Monte Carlo; Non-conjugate distribution; State-space model;
D O I
10.1007/s11222-008-9051-5
中图分类号
TP301 [理论、方法];
学科分类号
081202 ;
摘要
Different strategies have been proposed to improve mixing and convergence properties of Markov Chain Monte Carlo algorithms. These are mainly concerned with customizing the proposal density in the Metropolis-Hastings algorithm to the specific target density and require a detailed exploratory analysis of the stationary distribution and/or some preliminary experiments to determine an efficient proposal. Various Metropolis-Hastings algorithms have been suggested that make use of previously sampled states in defining an adaptive proposal density. Here we propose a general class of adaptive Metropolis-Hastings algorithms based on Metropolis-Hastings-within-Gibbs sampling. For the case of a one-dimensional target distribution, we present two novel algorithms using mixtures of triangular and trapezoidal densities. These can also be seen as improved versions of the all-purpose adaptive rejection Metropolis sampling (ARMS) algorithm to sample from non-logconcave univariate densities. Using various different examples, we demonstrate their properties and efficiencies and point out their advantages over ARMS and other adaptive alternatives such as the Normal Kernel Coupler.
引用
收藏
页码:421 / 433
页数:13
相关论文
共 50 条