TruncatedNormal bug还是我的错误？

Question

这个简单的例子：

import numpy as np
import pymc as pm
init = 0.5; minv = 0.; maxv = 1.
A = pm.Uniform('A', value=init, lower=minv, upper=maxv)
B = pm.TruncatedNormal('B', value=init, mu=A, tau=1., a=minv, b=maxv)
C = pm.Normal('C', value=np.random.normal(init,1,10), mu=B, tau=1., observed=True)
model = pm.MCMC(pm.Model([A, B, C]))
model.sample(1000)

几乎总是抛出以下错误（pymc2.3.3）：

pymc.Node.ZeroProbability: Stochastic B's value is outside its support,
 or it forbids its parents' current values.

我相信这个例子很小，因为如果我删除A或C它似乎有用（至少，它的工作频率更高）。

我错过了什么吗？如果没有，我将发布到github问题跟踪器。

Answer 1

它连续5次为我工作，但我还在使用PyMC版本2.3.2。我认为你应该发布问题跟踪器。