2014-06-11 27 views
0

我想了解MAP优化im pymc。我得到的lambda以下后验分布采用抽样MCMC enter image description herepymc中MAP的两次执行给出了不同的值

显然后,后是最大的波长= 0.20和95%的区间为[0.17,0.24](请纠正我,如果我错了)

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

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

使用MAP 0.200091865615 使用MAP 0.197584715205

使用MAP 1.28960939539 使用MAP 2.70871770586后

有人能解释一下是怎么回事,如何摆脱这一问题的前后前?

回答

0

可能是优化程序在第二个起始值时的收敛问题。你可以尝试一个不同的优化器吗?例如:

MAP(method='fmin_powell') 
+0

但是,我如何事先知道我将不得不使用其他优化器?我已经有了后门,因此我可以在这里核实,但这可能不是每次都是这样。 – turing

+0

你不知道。优化背后的艺术正在能够选择合理的初始值。你可以通过设置'value ='参数在PyMC中做到这一点(否则它使用随机抽取)。 –

相关问题