近似 Metropolis 算法 - 有意义吗?

机器算法验证 模拟 马尔可夫链蒙特卡罗 累积分布函数 大都会黑斯廷斯
2022-04-08 20:03:11

前段时间西安问MCMC for pdfs的cdfs相当于什么?天真的答案是在形式上使用“近似” Metropolis 算法

给定 1. 生成 2. 取X(t)=x(t)
Yq(y|x(t))

X(t+1)={Y with probability min(F(Y+ε)F(Yε)F(x(t)+ε)F(x(t)ε),1)x(t) otherwise.

其中F是一个目标 CDF,而ε是一些小的常数。这使我们能够将 Metropolis 算法与 CDF 一起使用。

问题是:有什么理由说明这实际上可能是个坏主意吗?

1个回答

不,我不明白为什么这是一个坏主意。在我看来,从 CDF 中抽取样本是一种自然(且有趣)的扩展。

但是,我认为接受应该是

min(F(Y+ε)F(Y)F(x(t)+ε)F(x(t)),1)

因为根据定义

limε0F(x+ε)F(x)ε=PDF(x)

不过,这是我第一次看到这种情况。有了从非平凡 CDF 中采样的有力案例,这可能会成为一个有趣的出版物。