从MCMC采样器输出中查找最可能的值

时间:2013-09-10 13:38:36

标签: python statistics mcmc pymc

我对PyMC比较陌生,我对MCMC采样器的输出有一个快速的问题。我希望找到MCMC采样器找到的最可能的变量值(后验的最大值)。有快速的方法吗?据推测,MCMC采样器发现的后验最大值的变量值可能与PyMC的MAP方法发现的变量值大不相同。

感谢开发人员提供PyMC。这对我的工作非常有用。此问题也已发布到PyMC google小组。

1 个答案:

答案 0 :(得分:1)

编辑:正如Ahmed Fasih指出的那样,我误读了这个问题。他在问题评论中的建议是如何做到这一点=)


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

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

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