mcmc - 如何制作截断的法线先验 : converting pymc2 to pymc3

Question

在 pymc3 中如何配置截断法线先验？在 pymc2 中，它非常简单(如下)，但在 pymc3 中，似乎不再有可用的截断正态分布。

Pymc2:

TruncatedNormal('gamma_own_%i_' % i, mu=go, tau=v_gamma_inv, value=0, a=-np.inf, b=0)

Pymc3:？

Answer 1

在 PyMC3 中，您可以使用 Bound 截断任何分布。首先，您必须构造有界分布(此处称为 BoundedNormal)，然后创建一个变量，在其中输入基础分布的常用参数:

with pm.Model() as model:
    BoundedNormal = pm.Bound(pm.Normal, lower=0, upper=1)
    n = BoundedNormal('n', mu=0, tau=10)
    tr = pm.sample(2000, pm.NUTS())

生成的分布如下所示: