中文

Metropolis Hastings 算法瞬态阶段的最优尺度:长期行为

统计方法学 2014-10-22 v3

摘要

我们考虑在 Rn\mathbb{R}^n 上具有高斯提议且目标概率测度为一维分布的 nn 重积的 Random Walk Metropolis 算法。众所周知(见 Roberts 等 (Ann. Appl. Probab. 7 (1997) 110-120)),在 nn\to\infty 的极限下,从平衡态出发,并对方差和时间尺度作为维数 nn 的函数进行适当缩放,可获得 Markov 链各分量的扩散极限。在 Jourdain 等(《Random Walk Metropolis 算法瞬态阶段的最优缩放:平均场极限》(2012) 预印本)中,我们在初始分布非目标概率测量的情况下推广了这一结果。所得的扩散极限是在 McKean 意义下非线性的随机微分方程的解。在本文中,我们证明了该方程向平衡态的收敛。我们讨论了为了优化提议分布的方差以加速向平衡态收敛的实际对策。我们的分析证实了最初在 Roberts 等 (Ann. Appl. Probab. 7 (1997) 110-120) 中提出的常数接受率策略(接受率在 1/41/41/31/3 之间)的优势。我们还探讨了 Metropolis-Adjusted Langevin Algorithm 的缩放问题。当从平衡态出发时,Roberts 和 Rosenthal (J. R. Stat. Soc. Ser. B. Stat. Methodol. 60 (1998) 255-268) 获得了方差最优缩放的扩散极限。在瞬态情形下,我们在形式上发现,最优方差在 nn 中的缩放方式因某分布矩的符号不同而大相径庭,该分布在平衡态时消失。这表明,很难从此类渐近结果中为 MALA 推导出实用的建议。

关键词

引用

@article{arxiv.1212.5517,
  title  = {Optimal scaling for the transient phase of Metropolis Hastings algorithms: The longtime behavior},
  author = {Benjamin Jourdain and Tony Lelièvre and Błażej Miasojedow},
  journal= {arXiv preprint arXiv:1212.5517},
  year   = {2014}
}

备注

Published in at http://dx.doi.org/10.3150/13-BEJ546 the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm)