python - 从 MCMC 采样器输出中找到最可能的值

Question

我对 PyMC 比较陌生，我有一个关于 MCMC 采样器输出的快速问题。我希望找到 MCMC 采样器发现的变量的最可能值（后验的最大值）。有没有快速的方法来做到这一点？据推测，MCMC 采样器发现的后验最大值处的变量值可能与 PyMC 的 MAP 方法发现的变量值大不相同。

感谢开发者提供 PyMC。它对我的工作非常有用。这个问题也已经发布到 PyMC google group。

Answer 1

编辑：正如 Ahmed Fasih 指出的那样，我误读了这个问题。他在问题评论中的建议是这样做的方法=）

---

该功能内置于 PyMC 中。主 pymc 命名空间中的类MAP接受 pymc 变量数组（或Model类），并公开一个fit方法。

map = mc.MAP(model) #or [var1, var2, .. ]
map.fit()

然后所有 pymc 变量将设置为最大后验。