PyMc3：错误的后验预测结果

Question

尝试从后验分布中抽取100个预测样本：

<figure class='bkg'>
  <b></b><b></b><b></b>
</figure>

我收到以下错误：

data = norm(50,5).rvs(30)
with pm.Model() as model_x:
   mu = pm.Uniform("mu",30,90)
   sigma = pm.HalfNormal("sigma", sd=10)
   dat = pm.Normal("dat", mu=mu, sd=sigma, observed = data)
   trace_x = pm.sample(5000)
   pred = pm.sample_posterior_predictive(trace_x, 100, model_x, size=30)

我在其他数据或使用Metropolis时遇到相同的错误。 此错误的原因是什么？如何正确采样？蒂亚。

Answer 1

将NumPy至少更新为1.13，这是引入numpy.divmod的时候。