解决“统计反思”中的问题时-我正在使用pymc3解决“萨拉曼德斯”问题10H4。(下面的代码)
采样器总是失败,并显示以下错误:
RV
alpha
.ravel()[0]的导数为零。
但是,当我将init ='map'选项传递给采样器时,与使用R代码的其他人得到的结果相似。强烈建议不要在采样器的输出中使用“ map” init选项,但许多其他init选项也存在问题。
数据在这里
我曾尝试根据与错误消息相关的网络搜索中的建议更改先验条件。
with pm.Model() as hw_10_4:
alpha = pm.Normal('alpha',mu=0,sd=5)
beta = pm.Normal('beta',mu=0,sd=5)
l = pm.Deterministic('lambda',pm.math.exp(alpha + beta*pctcover))
s = pm.Poisson('salam',l,observed=salaman)
trace_10_4 = pm.sample(1000, tune=1000, init='map')
我得到了想要的结果,但是我有兴趣进一步了解为什么我首先遇到了问题。我是PYMC3的新手,我想完全了解发生了什么。
谢谢!
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句