yongsheng2021 2021-04-15 12:15 采纳率: 25%
浏览 209
已采纳

关于pymc3报错问题

import pymc3 as pm
import numpy as np
import scipy.stats as stats
import matplotlib.pyplot as plt
import seaborn as sns
palette = 'muted'
sns.set_palette(palette); sns.set_color_codes(palette)

np.random.seed(123)
n_experiments = 4
theta_real = 0.35  # unkwon value in a real experiment
data = stats.bernoulli.rvs(p=theta_real, size=n_experiments)

with pm.Model() as our_first_model:
    # a priori
    theta = pm.Beta('theta', alpha=1, beta=1)
    # likelihood
    y = pm.Bernoulli('y', p=theta, observed=data)
    #y = pm.Binomial('theta',n=n_experimentos, p=theta, observed=sum(datos))
    start = pm.find_MAP()
    step = pm.Metropolis()
    trace = pm.sample(1000, step=step, start=start,cores=1)

with our_first_model:
    step = pm.Metropolis()
    multi_trace = pm.sample(1000, step=step,cores=1)

if __name__ == '__main__':
    burnin = 0  # no burnin
    multi_chain = multi_trace[burnin:]
    pm.traceplot(multi_chain, lines={'theta':theta_real})
    plt.show()
    # print(pm.gelman_rubin(multi_chain))
    pm.forestplot(multi_chain, varnames=['theta'])
    plt.show()

总这样报错呢!?????????????????????????

TypeError                                 Traceback (most recent call last)
<ipython-input-1-88ae9a9ad243> in <module>
     34     plt.show()
     35     # print(pm.gelman_rubin(multi_chain))
---> 36     pm.forestplot(multi_chain, varnames=['theta'])
     37     plt.show()
     38 

~\anaconda3\lib\site-packages\pymc3\plots\__init__.py in wrapped(*args, **kwargs)
     37         if "varnames" in kwargs:
     38             raise DeprecationWarning(
---> 39                 f"The `varnames` kwarg was renamed to `var_names`.", stacklevel=2
     40             )
     41         original = func.__name__

TypeError: DeprecationWarning() takes no keyword arguments
  • 写回答

4条回答 默认 最新

  • CSDN专家-HGJ 2021-04-15 14:49
    关注

    源码中有一个弃用警告,就是说pymc3中的有关plot函数的关键字参数varnames已经弃用,改成了var_names,用var_names=['theta']试试。

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(3条)

报告相同问题?

悬赏问题

  • ¥15 进行基于用户的协同过滤推荐,使用准确率和召回率评价推荐模型的效果
  • ¥15 关于#cesium#的问题:cesium 根据经纬度、高度、方位角、俯仰角和距离求另一点高度
  • ¥15 word脚注设置不起空格符号
  • ¥15 为什么apriori的运行时间会比fp growth的运行时间短呢
  • ¥15 求解密json或者美化
  • ¥15 帮我解决一下膳食平衡的线性规划模型的数据实例
  • ¥40 万年历缺少农历,需要和阳历同时显示
  • ¥250 雷电模拟器内存穿透、寻基址和特征码的教学
  • ¥200 比特币ord程序wallet_constructor.rs文件支持一次性铸造1000个代币,并将它们分配到40个UTXO上(每个UTXO上分配25个代币),并设置找零地址
  • ¥15 关于Java的学习问题