如何用PyMC3期待不可微函数

时间:2019-04-05 14:31:43

标签: python pymc3 expectations

我想期望一个函数,它是几个随机变量的非线性函数。有关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)给定用户定义的功能,我是否需要使用特定的采样器?

为幼稚的问题致歉。

0 个答案:

没有答案