×

仿射不变集合采样器的高维“拉伸运动”特性。 (英语) Zbl 1521.62008年

总结:我们介绍了仿射不变集合采样器马尔可夫链蒙特卡罗方法的理论和实际性质。在高维情况下,取样器的“拉伸运动”具有不寻常和不理想的特性。我们用一个具有已知均值和协方差结构的多维相关高斯玩具问题以及Rosenbrock问题的多元版本来证明这一点。痕迹图的目视检查表明老化期很短。经过仔细检查,我们发现目标分布的平均值和方差与已知值不匹配,并且链需要很长时间才能收敛。当(n)增加到50以上时,这个问题变得更严重。我们还应用了适用于集成方法的不同诊断方法,以确定是否存在收敛不足。诊断包括Gelman-Rubin方法、海德堡-韦尔奇检验、综合自相关和接受率。单个步行者的轨迹图似乎也很有用。因此,我们得出结论,拉伸动作在中高维度中应谨慎使用。我们还提供了一些启发式结果来解释这种行为。

MSC公司:

62-08 统计问题的计算方法
65二氧化碳 蒙特卡罗方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Brooks,S.&Gelman,A.(1998年)。监测迭代模拟收敛性的一般方法。计算与图形统计杂志7,434-455。
[2] Christen,J.A.和Fox,C.(2010年)。一种用于连续分布(t‐walk)的通用采样算法。贝叶斯分析5263-281·Zbl 1330.62007
[3] Cotter,S.L.、Roberts,G.O.、Stuart,A.M.和White,D.(2013)。函数的MCMC方法:修改旧算法使其更快。arXiv.org.28,424-446·Zbl 1331.62132号
[4] Crossfield,I.J.M.,Petigura,E.,Schlieder,J.E.等人(2015)。附近的一颗M星,由K2发现,有三颗凌日超地球。天体物理学杂志804,10。
[5] Dixon,L.C.W.&Mills,D.J.(1994)。舍入误差对可变公制方法的影响。最优化理论与应用杂志80,175-179·Zbl 0797.90091号
[6] Ethier,S.N.和Kurtz,T.G.(1986年)。马尔可夫过程:特征和收敛。《概率与数理统计威利系列》,纽约:威利出版社·Zbl 0592.60049号
[7] Foreman‐Mackey,D.,Hogg,D.W.,Lang,D.&Goodman,J.(2013)。Emcee:MCMC锤子。太平洋天文学会出版物125,306-312。
[8] Foreman‐Mackey,D.C.(2012年)。主持MCMC锤子。可从以下URL获得:https://emcee.readthedocs.io。
[9] Gelman,A.和Rubin,D.(1992年)。使用多序列的迭代模拟推断。统计科学7,457-472·Zbl 1386.65060号
[10] Goodman,J.&Sokal,A.(1989)。多重网格蒙特卡罗方法。概念基础。物理审查。D、 《粒子与场》402035-2071。
[11] Goodman,J.和Weare,J.(2010年)。用仿射不变性集成采样器。应用数学与计算科学传播5,65-80·Zbl 1189.65014号
[12] 黑斯廷斯,W.K.(1970)。使用马尔可夫链的蒙特卡罗抽样方法及其应用。生物医学57,97-109·Zbl 0219.65008号
[13] Lampart,T.(2012)。仿射不变性集合采样器的实现与性能比较。技术报告,MOSAIC集团,苏黎世理工大学计算机科学系理论计算机科学研究所。
[14] Metropolis,N.、Rosenbluth,A.W.、Rosenbruth,M.N.、Teller,A.H.和Teller等人(1953年)。快速计算机器的状态方程计算。《化学物理杂志》21,1087-1092·兹比尔1431.65006
[15] Murray,I.、Adams,R.和MacKay,D.(2010年)。椭圆切片采样。《第十三届国际人工智能和统计会议记录》,第9卷。第541-548页。
[16] Neal,R.(2011)。MCMC使用哈密顿动力学。《马尔可夫链蒙特卡罗手册》,第5章。《现代统计方法手册》,纽约,第113-162页,美国纽约:查普曼和霍尔/CRC·Zbl 1229.65018号
[17] Neal,R.M.(2003年)。切片取样。统计年鉴31705-741·Zbl 1051.65007号
[18] Plummer,M.、Best,N.、Cowles,K.和Vines,K.(2006年)。Coda:MCMC的收敛诊断和输出分析。R新闻6,7-11。
[19] 罗森布罗克,H.H.(1960)。一种自动查找函数最大值或最小值的方法。《计算机杂志》3175-184。
[20] Sinharay,S.(2003)。评估马尔可夫链蒙特卡罗算法的收敛性:综述。ETS研究报告系列1,i-52。
[21] Ter Braak,C.&Vrugt,J.(2008)。具有斯诺克更新器和更少链的差分进化马尔可夫链。统计与计算18435-446。
[22] Ter Braak,C.J.F.(2006年)。遗传算法差分进化的马尔可夫链蒙特卡罗版本:实际参数空间的简单贝叶斯计算。统计与计算16,239-249。
[23] Vanderburg,A.、Monte,B.T.、Johnson,J.A.(2015)。K2行星发现的特征:一个超地球经过明亮的K矮星HIP 116454。天体物理学杂志800,59。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。