Studying Convergence of Markov Chain Monte Carlo Algorithms Using Coupled Sample Paths
提出一种通过比较同一采样器在不同随机数种子下的耦合链路径,来诊断马尔可夫链蒙特卡洛算法收敛性的简单方法,并用于检测多峰后验分布的模式。
Abstract I describe a simple procedure for investigating the convergence properties of Markov chain Monte Carlo sampling schemes. The procedure uses coupled chains from the same sampler, obtained by using the same sequence of random deviates for each run. By examining the distribution of the iteration at which all sample paths couple, convergence properties for the system can be established. The procedure also provides a simple diagnostic for detecting modes in multimodal posteriors. Several examples of the procedure are provided. In Ising models, the relation between the correlation parameter and the convergence rate of rudimentary Gibbs samplers is investigated. In another example, the effects of multiple modes on the convergence of coupled paths are explored using mixtures of bivariate normal distributions. The technique is also used to evaluate the convergence properties of a Gibbs sampling scheme applied to a model for rat growth rates.