pymc 相关问题

PyMc是一个Python模块,用于提供贝叶斯统计模型,算法和估计。目前广泛使用两种版本:2和3,它们是显着不同的。版本2不再受支持,但版本3与以前的代码不完全兼容,将V2代码转换为V3并不总是直截了当。如果你有一个特定于V3的PyMc问题,除了[pymc]标签外,还要考虑使用[pymc3]标签。

如何在使用自定义对数似然的 PyMC V5 模型上应用 Blackjax 采样器

我正在使用 PyMC v5 在模型中执行哈密顿蒙特卡罗。我可以运行下面的代码,但即使有多个核心,它也非常慢。我有一个 applyMCMC 函数用于此目的: # 定义一个

回答 1 投票 0

贝叶斯 MMM 用 Pytensor 代替 Theano

我有一段代码,可以使用 Theano 执行几何广告库存衰减。 这是一段旧代码,我需要使用最新版本的 PyTensor 来更新它。 有人可以帮我转换吗

回答 1 投票 0

MAP 估计与后验平均值不一致

MAP 估计值不应该接近迹线样本直方图的中心吗? 在对我的模型进行建模时,我发现 find_MAP 估计值与跟踪变量的直方图之间非常一致

回答 1 投票 0

如何从 pymc.LKJCorr 构造相关矩阵?

我想使用 pymc.LKJCorr 分布类显式构造一个相关矩阵,但我不相信我对 pymc.expand_packed_triangle 的理解。这是一个最小的工作

回答 2 投票 0

如何在 PyMC (5.10.0) 中参数化模型中每个自变量的滞后参数?

我是 PyMC 新手,我使用的是 5.10.0 版本。我正在运行一个简单的媒体混合模型,我需要参数化独立(媒体活动)和依赖(sa...

回答 1 投票 0

无法使用 Pymc 3.11.4 加载之前保存的跟踪以供将来推理

我正在尝试实现一个简单的分层贝叶斯推理。我使用 Pymc 推断简单线性模型 y = mx 的参数“m”。我想保存跟踪以便稍后加载。

回答 1 投票 0

添加新发行版时 FreeRV 类型出现问题

我正在尝试通过包装 Agner Fogs c++ 版本(https://www.agner.org/random/)向 PyMC3(Wallenius 非中心超几何)添加新的离散分布。 我已经成功把

回答 1 投票 0

具有 MVN 先验的分层变化效应模型

我正在努力做什么 我已经在 pymc 中处理过多元先验(我使用的是 4.0.1),但我无法在分层模型中使用它们。在我的示例中,我正在建模回归专业...

回答 1 投票 0

应用贝叶斯框架估计线性回归参数

我是贝叶斯框架的新手,我从Kruschke的书和BDA课程开始了解一些理论来解决我的问题。 我的问题很常见:我收集了一些数据 x,y 并且我......

回答 0 投票 0

pymc.Model 的示例子类中的 super().__init__() 中的模型参数是什么?

pymc.Model 文档显示了这个例子: 自定义模型类(模型): # 1) 覆盖初始化 def __init__(self, mean=0, sigma=1, name=''): # 2) 首先调用 super 的 init,传递模型和名称...

回答 1 投票 0

使用 PyMC 进行贝叶斯逻辑回归并收到错误:“您的数据似乎只有一个值或没有有限值”?

我正在尝试按照本教程在 Python 中进行贝叶斯逻辑回归。 在创建我的轨迹图时,我总是会遇到同样的错误:“您的数据似乎只有一个

回答 0 投票 0

批量矩阵 RV 的维度规范(例如 Wishart dist)

我目前正在构建一个具有多元正态似然和基于截断定义的狄利克雷过程混合模型。具体来说,这是 $K$ 组件的混合,每个组件都与...

回答 0 投票 0

什么是 pymc.MutableData?

文档并没有真正解释什么是 pymc.MutableData。从这个例子来看,这似乎是一种将数据放入模型的方法,但我不清楚这与使用 observed=data i... 有何不同...

回答 0 投票 0

在 PyMC 中,“瘦”论点似乎不再使用。可以用什么代替?

试图通过示例理解 PyMC,我制作了一个小模型,参考一组鸟类体重(观察到的数据:y20),据说来自正常种群 N(mu, sigma)。 将 pymc 导入为 pm ...

回答 1 投票 0

如何决定在PyMC3中对参数使用什么priors分布?

我正在寻找pymc3包,我有兴趣在一个场景中实现包,我有几个不同的信号,每个信号有不同的振幅。然而,我是...

回答 1 投票 0

贝叶斯IRT Pymc3 - 参数推理

我想用PyMC3来估计IRT模型,我生成的数据有以下分布:alpha_fix = 4 beta_fix = 100 theta= np.random.normal(100,15,1000) prob = np.exp(alpha_fix*(theta-. 我生成的数据具有以下分布:alpha_fix = 4 beta_fix = 100 theta= np.random.normal(100,15,1000) prob = np.exp(alpha_fix*(theta-......)

回答 1 投票 0

在Python中实现贝叶斯分层建模。

我正试图在不同的数据集上实现贝叶斯层次建模,我通过互联网搜索,发现了Pymc3的这个文档。https:/docs.pymc.ionotebooksGLM-hierarchical...。

回答 1 投票 -1

使用pyMC3的贝叶斯网络中给定数据的参数的计算概率

我是pyMC3的新手,我想知道是否可以使用它来解决以下问题:我有一个贝叶斯网络(我的BN的图像:我的问题的贝叶斯网络),但我不知道...

回答 1 投票 0

从Dirichlet分布中进行采样的PyMC中的FloatingPointError

[使用修饰符定义“指数随机变量的对数”的随机对象失败后,我决定使用...]手动编写用于此新发行版的代码...

回答 2 投票 21


© www.soinside.com 2019 - 2024. All rights reserved.