×

通过扩散极限最小化MCMC方差,并应用于模拟回火。 (英语) Zbl 1298.60078号

如果马尔可夫链(X_n)具有平稳分布(pi),那么对于适当大的(n),(int h(X)pi(dx))通常可以通过(n^{-1}sum_{i=1}^n h(X_i)来估计。估计的效率可以用渐近方差(mathrm{Var}(h,P)=lim{n\to-infty}n^{-1}mathrm}Var}(sum{i=1}^nh(X_i))来度量,其中(P)是马尔可夫链的转移核。给定两个马尔可夫链核(P_1,P_2),它们都具有相同的不变测度(pi),我们说,对于所有可容许的(h),如果(mathrm{Var}(h,P_1)leq\mathrm}Var}(h,P2)),则(P_1\)支配(P_2)。
本文发展了这种方法,用于比较朗之万扩散的渐近方差。给出了充分条件。然后将结果应用于高维空间中的模拟回火算法。证明了它们的极限描述是朗之万扩散,并描述了最有效的马尔可夫链蒙特卡罗算法。

MSC公司:

60J22型 马尔可夫链中的计算方法
65立方厘米 马尔可夫链的数值分析或方法
62米05 马尔可夫过程:估计;隐马尔可夫模型
10层62层 点估计
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Atchadé,Y.F.、Roberts,G.O.和Rosenthal,J.S.(2011年)。都市耦合马尔可夫链蒙特卡罗的最优尺度。统计计算。21 555-568. ·Zbl 1223.65001号 ·doi:10.1007/s11222-010-9192-1
[2] 贝达德,M.(2007)。非身份识别目标分布Metropolis算法的弱收敛性。附录申请。普罗巴伯。17 1222-1244. ·Zbl 1144.60016号 ·doi:10.1214/1050516070000096
[3] Bédard,M.(2008)。Metropolis算法的最佳接受率:超过0.234。随机过程。申请。118 2198-2222. ·Zbl 1155.60310号 ·doi:10.1016/j.spa.2007.12.005
[4] Bédard,M.和Rosenthal,J.S.(2008)。Metropolis算法的最佳缩放:朝向一般目标分布。加拿大。J.统计。36 483-503. ·Zbl 1167.60344号 ·doi:10.1002/cjs.5550360401
[5] Beskos,A.、Roberts,G.和Stuart,A.(2009年)。高维非产品目标上本地Metropolis-Hastings链的最佳缩放。附录申请。普罗巴伯。19 863-898. ·Zbl 1172.60328号 ·doi:10.1214/08-AAP563
[6] Bhattacharya,R.N.和Waymire,E.C.(1990年)。随机过程及其应用。威利,纽约·Zbl 0715.60046号
[7] Brooks,S.、Gelman,A.、Jones,G.L.和Meng,X.-L.编辑(2011年)。马尔可夫链蒙特卡罗手册。查普曼和霍尔/CRC,佛罗里达州博卡拉顿·Zbl 1218.65001号
[8] Ethier,S.N.和Kurtz,T.G.(1986年)。马尔可夫过程:特征和收敛。纽约威利·Zbl 0592.60049号
[9] Fort,G.、Moulines,E.和Priouret,P.(2011年)。自适应和交互式马尔可夫链蒙特卡罗算法的收敛性。安。统计师。39 3262-3289. ·Zbl 1246.65003号
[10] Geyer,C.(1992)。实用马尔可夫链蒙特卡罗法。统计师。科学。7 473-483.
[11] Kofke,D.A.(2002年)。关于重复交换蒙特卡罗试验的接受概率。化学杂志。物理。117 6911. 勘误表:J.Chem。物理。120 10852.
[12] Leisen,F.和Mira,A.(2008年)。将Peskun和Tierney序推广到连续时间马尔可夫链。统计师。Sinica 18 1641-1651年·Zbl 1252.60074号
[13] Liu,J.S.、Wong,W.H.和Kong,A.(1994年)。吉布斯采样器的协方差结构及其在估计量和增强方案比较中的应用。生物特征81 27-40·Zbl 0811.62080号 ·doi:10.1093/biomet/81.1.27
[14] Marinari,E.和Parisi,G.(1992年)。模拟回火:一种新的蒙特卡罗方案。欧罗普提斯。莱特。19 451-458。
[15] 米拉,A.(2001)。排序和改进蒙特卡罗-马尔可夫链的性能。统计师。科学。16 340-350. ·Zbl 1127.60312号 ·doi:10.1214/ss/1015346319
[16] Mira,A.和Geyer,C.J.(2000年)。关于不可逆马尔可夫链。蒙特卡洛方法(多伦多,安大略省,1998年)。菲尔德学院通信26 95-110。阿默尔。数学。佛罗里达州普罗维登斯Soc·Zbl 0969.60071号
[17] Mira,A.和Leisen,F.(2009年)。离散和连续时间马尔可夫链的协方差排序。统计师。中国19 651-666·Zbl 1168.62073号
[18] Peskun,P.H.(1973)。使用马尔可夫链的最佳蒙特卡罗抽样。生物特征60 607-612·Zbl 0271.62041号 ·doi:10.1093/biomet/60.3.607
[19] Predescu,C.、Predescu,M.和Ciobanu,C.V.(2004)。经典正则系统平行回火抽样的不完全β函数定律。化学杂志。物理。120 4119-4128.
[20] Roberts,G.O.、Gelman,A.和Gilks,W.R.(1997年)。随机行走Metropolis算法的弱收敛性和最优尺度。附录申请。普罗巴伯。7 110-120. ·Zbl 0876.60015号 ·doi:10.1214/aoap/1034625254
[21] Roberts,G.O.和Rosenthal,J.S.(1998年)。朗之万扩散离散近似的最佳缩放。J.R.统计社会服务。B统计方法。60 255-268. ·兹比尔0913.60060 ·doi:10.111/1467-9868.00123
[22] Roberts,G.O.和Rosenthal,J.S.(2001年)。各种Metropolis-Hastings算法的最佳缩放。统计师。科学。16 351-367. ·Zbl 1127.65305号 ·doi:10.1214/ss/1015346320
[23] Sinclair,A.(1992年)。改进了马尔可夫链和多商品流混合速率的界。组合概率。计算。1 351-370. ·Zbl 0801.90039号 ·网址:10.1017/S096354830000390
[24] Tierney,L.(1994)。用于探索后验分布的马尔可夫链。安。统计师。22 1701-1762. ·Zbl 0829.62080号 ·doi:10.1214/aos/1176325750
[25] Tierney,L.(1998)。关于一般状态空间的Metropolis-Hastings核的注记。附录申请。普罗巴伯。8 1-9. ·Zbl 0935.60053号 ·doi:10.1214/oap/1027961031
[26] Ward,A.R.和Glynn,P.W.(2003)。一类带有叛逆的马尔可夫队列的扩散近似。排队系统。43 103-128. ·Zbl 1054.60100号 ·doi:10.1023/A:1021804515162
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。