Python PyMc3：后验预测检查结果出错

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

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)

我得到以下错误：

    AttributeError Traceback (most recent call last)
    <ipython-input-172-9812ddb325e6> in <module>()
    ----> 6     pred = pm.sample_posterior_predictive(trace_x, 100, model_x, size=30)

   .../python2.7/site-packages/pymc3/sampling.pyc in sample_posterior_predictive(trace, samples, model, vars, size, random_seed, progressbar)
       1128         for idx in indices:
       1129             if nchain > 1:
    -> 1130                 chain_idx, point_idx = np.divmod(idx, len_trace)
       1131                 param = trace._straces[chain_idx % nchain].point(point_idx)
       1132             else:

    AttributeError: 'module' object has no attribute 'divmod'

AttributeError回溯（最近一次调用）
在（）
---->6 pred=pm.样本后验预测（轨迹x，100，模型x，大小=30）
…/python2.7/site-packages/pymc3/sampling.pyc in sample\u posterior\u predictive（跟踪、样本、模型、变量、大小、随机种子、进度条）
1128对于指数中的idx：
1129如果nchain>1：
->1130 chain_idx，point_idx=np.divmod（idx，len_trace）
1131 param=trace.\u straces[chain\u idx%nchain]。点（point\u idx）
1132其他：
AttributeError:“模块”对象没有属性“divmod”

我在其他数据或使用Metropolis时也会遇到同样的错误。

---

这个错误的原因是什么？如何正确取样？Tia。

将NumPy更新到至少1.13，这是在推出时