[论文解读] Optimal Proposal Design for Random Walk Type Metropolis Algorithms with Gaussian Random Field Priors
本文提出了一种带有 Ornstein-Uhlenbeck 提议的可逆 Metropolis-Hastings 算法,用于从希尔伯特空间上相对于高斯先验测度绝对连续的测度中进行采样。通过分析由此产生的马尔可夫链的扩散极限,证明其收敛于由具有与高斯先验匹配的空间相关性的维纳过程驱动的噪声梯度流 SDE,从而建立了该方法在目标测度下的渐近可逆性与稳定性。
Consider a probability measure on a Hilbert space defined via its density with respect to a Gaussian. The purpose of this paper is to demonstrate that an appropriately defined Markov chain, which is reversible with respect to the measure in question, exhibits a diffusion limit to a noisy gradient flow, also reversible with respect to the same measure. The Markov chain is defined by applying a Metropolis-Hastings accept-reject mechanism to an Ornstein-Uhlenbeck proposal which is itself reversible with respect to the underlying Gaussian measure. The resulting noisy gradient flow is a stochastic partial differential equation driven by a Wiener process with spatial correlation given by the underlying Gaussian structure.
研究动机与目标
- 开发一种马尔可夫链蒙特卡洛方法,用于从相对于高斯先验测度绝对连续的希尔伯特空间测度中进行采样。
- 确保马尔可夫链相对于目标测度是可逆的,以保持细致平衡。
- 建立该链的扩散极限,使其收敛于由具有空间相关性的维纳过程驱动的随机偏微分方程(噪声梯度流)。
- 分析该算法在目标测度下的渐近行为,确保其长期稳定性和正确性。
提出的方法
- 定义一种使用与底层高斯测度可逆的 Ornstein-Uhlenbeck 提议的 Metropolis-Hastings 算法。
- 对 OU 提议应用接受-拒绝机制,以在希尔伯特空间上构建一个可逆的马尔可夫链。
- 在扩散缩放下分析马尔可夫链的极限,证明其收敛于随机偏微分方程。
- 将极限 SDE 推导为一个噪声梯度流,其中噪声的空间相关性由高斯先验的协方差结构决定。
- 证明极限 SDE 也相对于目标测度是可逆的,从而在极限中保持细致平衡。
- 使用泛函中心极限定理技术,以证明马尔可夫链收敛于随机 PDE。
实验结果
研究问题
- RQ1所提出的带有 OU 提议的 Metropolis-Hastings 算法在扩散极限下是否收敛于一个定义良好的随机 PDE?
- RQ2极限 SDE 是否相对于相对于高斯先验密度定义的目标测度是可逆的?
- RQ3极限 SDE 中噪声的空间相关性如何与底层高斯测度的协方差结构相关联?
- RQ4在扩散缩放下,马尔可夫链的渐近行为是什么?其在极限中是否保持细致平衡?
- RQ5该算法能否被合理地视为无限维空间中基于梯度采样方法的有效近似?
主要发现
- 在扩散缩放下,由 Metropolis-Hastings 算法与 OU 提议定义的马尔可夫链以分布收敛于一个随机 PDE。
- 极限 SDE 是一个相对于目标测度可逆的噪声梯度流。
- 极限 SDE 中的噪声根据底层高斯先验的协方差算子具有空间相关性。
- 收敛结果确保了该算法在无限维极限下保持细致平衡。
- 该方法为在具有高斯先验的希尔伯特空间中使用随机游走型 Metropolis 算法提供了严格的理论依据。
- 该分析证实了该算法在渐近区域的稳定性和正确性,支持其在高维与无限维贝叶斯推断中的应用。
更好的研究,从现在开始
从论文设计到论文写作,大幅缩短您的研究时间。
无需绑定信用卡
本解读由 AI 生成,并经人工编辑审核。