我需要做一个模拟来评估一个 3 参数函数的积分,我们说,它有一个非常复杂的公式。要求使用 MCMC 方法对其进行计算并实现 Metropolis-Hastings 算法以生成分布为的值,并建议使用 3 变量正态作为提议分布。阅读一些关于它的例子,我发现其中一些使用具有固定参数的法线,而一些使用变量均值,其中是最后接受的值根据分布。我对这两种方法都有一些疑问:
1)选择最后一个接受的值作为我们提案分布的新均值是什么意思?我的直觉说它应该保证我们的值将更接近分布为的值,并且接受的机会会更大。但它不是过于集中我们的样本吗?可以保证,如果我获得更多样本,链条将变得静止?
2)选择固定参数(因为真的很难分析)会不会真的很困难并且取决于我们需要选择启动算法的第一个样本?在这种情况下,找出哪个更好的最佳方法是什么?
这些方法中的一种是否比另一种更好,或者这取决于具体情况?
我希望我的疑问是明确的,如果可以提供一些文献我会很高兴(我已经阅读了一些关于这个主题的论文,但越多越好!)
提前致谢!