pymc

    5热度

    1回答

    我更新了一些计算,我曾经pymc2到pymc3,我有一些问题,采样的行为,当我在我的模型的一些离散的随机变量。作为一个例子,考虑下面的模型中使用pymc2: import pymc as pm N = 100 data = 10 p = pm.Beta('p', alpha=1.0, beta=1.0) q = pm.Beta('q', alpha=1.0, beta=1.0) A

    0热度

    2回答

    我需要建模并估计从资产类别方差 - 协方差矩阵返回,所以我一直在寻找在股票收益例子在第6章给出的https://github.com/CamDavidsonPilon/Probabilistic-Programming-and-Bayesian-Methods-for-Hackers 这是我简单的实现方式,其中我从使用已知均值和方差 - 协方差矩阵的多变量正态分布开始。然后我尝试使用非信息性的pr

    0热度

    1回答

    我阅读以下纸(http://www3.stat.sinica.edu.tw/statistica/oldpdf/A10n416.pdf),其中它们的方差 - 协方差矩阵Σ建模为: Σ= DIAG(S)* R * DIAG(S)...(式1中(S)是对角元素S的对角矩阵,R是k×k相关矩阵。 我该如何使用PyMC来实现这个功能? 下面是一些初步的代码,我写道: import numpy as np

    3热度

    1回答

    我无法使用pymc以及psycopg2。从教程下面这个简单的代码片段: import pymc as pm with pm.Model() as model: x = pm.Normal('x', mu=0., sd=1) 导致以下错误: Exception: The environment variable 'DYLD_FALLBACK_LIBRARY_PATH' does no

    4热度

    1回答

    想我会开始下面这个例子: http://www.databozo.com/2014/01/17/Exploring_PyMC3.html 但是,当我按照精确使用pymc 2.3的例子中,我得到一个出口,并告诉该API已经改变 UserWarning: MCMC()语法已弃用。请通过M = MCMC(输入)明确传入节点。 'MCMC()语法已弃用。请通过M = MCMC(输入)明确地传入节点。')但

    1热度

    1回答

    继续使用this example做相当简单的贝叶斯线性回归使用PYMC3(学习,我希望)我得到的初始示例运行,但然后尝试使用我自己的数据,并得到: ValueError: Optimization error: max, logp or dlogp at max have non-finite values. Some values may be outside of distribution

    0热度

    1回答

    我一直在使用PyMC分析一些高能物理数据。它一直在努力完善,分析已经完成,我们正在撰写论文。 但是我有一个小问题。我使用RAM数据库后端运行采样器。这些痕迹现在已经在IPython内核进程的内存中呆了几个月。问题是工作站支持人员想要执行内核升级并重启该工作站。这会让我失去痕迹。我想保留这些痕迹(而不是仅仅产生新的痕迹),因为它们就是我所做的所有情节。我还想包含部分痕迹(仅限感兴趣的参数)作为出版物

    1热度

    2回答

    在我写的模型中,我必须计算一个数量的误差函数。我正在试图做的是这样的: from math import erf import numpy as np import pymc as pm sig = pm.Exponential('sig', beta=0.1, size=10) x = erf(sig ** 2) 因为erf不会对数组工作,就会失败。我想: @pm.determin

    0热度

    1回答

    我在pymc2以下模型: import pymc from scipy.stats import gamma alpha = pymc.Uniform('alpha', 0.01, 2.0) scale = pymc.Uniform('scale', 1.0, 4.0) @pymc.deterministic(plot=False) def beta(scale=scale):

    5热度

    1回答

    我试图端口的简单生存模式from here (the first one in introduction)形式PyMC 2至3 PyMC但是,我没有找到“发现”任何同等装饰和我试图写一个新的分发失败。有人可以提供一个例子,这是如何在PyMC 3中完成的?