我想期望一个函数,它是几个随机变量的非线性函数。有关PyMC3的所有文档似乎都与发现后验分布有关,但在我看来,这应该是可能的。这是一个点头的例子:
import numpy as np
import pymc3 as pm
def non_linear_fun(a):
return np.abs(a, 1.)
with pm.Model() as model:
a = pm.Exponential('a', lam=7.0)
b = pm.Deterministic('b', non_linear_fun(a))
trace = pm.sample(10000)
这失败了,因为我认为我应该用输入和输出类型来修饰该方法。
1)它们应该是什么? 2)给定用户定义的功能,我是否需要使用特定的采样器?
为幼稚的问题致歉。