python - 在 pymc 中两次执行 MAP 给出不同的值

Question

我想了解 MAP 优化 im pymc。lambda我使用 mcmc 采样 后得到以下后验分布

显然，后验在 lambda = 0.20 处最大，95% 的区间为 [0.17, 0.24]（如果我错了，请纠正我）

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

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

使用 MAP 0.200091865615 之前 使用 MAP 0.197584715205 之后

使用 MAP 1.28960939539 之前 使用 MAP 2.70871770586 之后

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

Answer 1

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

MAP(method='fmin_powell')