我正在复制艾伦·唐尼( Allen )在“想想贝斯”(Think)中向pymc3展示的一些例子。
他的伟大著作为我们提供了一些关于贝叶斯方法的介绍性例子,并且是使用Allen自己的库完成的。我如何定义自定义的可能性是pymc3?我用DensityDist来创建我自己的似然函数,但是我正在复制的这个函数依赖于从1到N的假设(假设N= 100),而在pymc3中,我找不到从张量得到X的方法。
我正在尝试使用pymc3定义一个复杂的自定义似然函数。似然函数涉及大量迭代,因此我尝试使用theano的扫描方法直接在theano中定义迭代。这里有一个非常简单的例子,它说明了我所面临的挑战。我试图定义的(伪)似然函数是两个pymc3随机变量p和θ的和。当然,我可以简单地返回p+theta,但我尝试编写的实际似然函数更复杂,而且我认为我需要使用theano.scan,因为它涉及很多迭代。import pymc3 as pmimpo