我对PyMC比较陌生,我有一个关于MCMC采样器输出的快速问题。我希望找到MCMC采样器发现的变量的最大可能值(后验最大值)。有什么快速的方法可以做到吗?假设MCMC采样器所发现的最大后变量值与PyMC的MAP方法发现的变量有很大的不同。
感谢开发人员提供PyMC。这对我的工作非常有用。这个问题也被发布到了PyMC谷歌集团。
最佳答案
编辑:正如艾哈迈德法西所指出的,我误解了这个问题。他在评论这个问题时的建议是怎样做
该功能内置于PyMC中。主pymc命名空间中的类MAP
接受pymc变量数组(或Model
类),并公开fit
方法。
map = mc.MAP(model) #or [var1, var2, .. ]
map.fit()
然后,所有的PYMC变量将被设置为最大顶点。