在pymc中执行两次MAP会产生不同的值

Question

我想了解MAP优化im pymc。使用mcmc取样后，我得到了lambda的后验分布

显然，后验在λ= 0.20处是最大值，95％间隔是[0.17,0.24]（如果我错了请纠正我）

据我所知，MAP给出了一个点估计值（lambda的值具有最大后验概率），但是当我运行MAP两次时，我得到的值不同，不应该是这种情况。

我在使用map执行d相同程序后的nd之前打印lambda的值。

使用MAP之前 0.200091865615 使用MAP后 0.197584715205

使用MAP之前 1.28960939539 使用MAP后 2.70871770586

有人可以解释发生了什么以及如何摆脱这个问题吗？

Answer 1

可能是优化器在第二个起始值处出现收敛问题。你可以尝试不同的优化器吗？例如：

MAP(method='fmin_powell')