我对PyMC比较陌生,我对MCMC采样器的输出有一个快速的问题。我希望找到MCMC采样器找到的最可能的变量值(后验的最大值)。有快速的方法吗?据推测,MCMC采样器发现的后验最大值的变量值可能与PyMC的MAP方法发现的变量值大不相同。
感谢开发人员提供PyMC。这对我的工作非常有用。此问题也已发布到PyMC google小组。
答案 0 :(得分:1)
编辑:正如Ahmed Fasih指出的那样,我误读了这个问题。他在问题评论中的建议是如何做到这一点=)
该功能内置于PyMC中。主pymc命名空间中的类MAP
接受一个pymc变量数组(或Model
类),并公开fit
方法。
map = mc.MAP(model) #or [var1, var2, .. ]
map.fit()
然后将所有pymc变量设置为最大值。