从pymc3中推断的参数生成预测
我遇到了一个常见问题，想知道是否有人可以帮忙。我经常想以两种模式使用pymc3：训练（即实际运行参数推断）和评估（即使用推断的参数生成预测）。 通常情况下，我希望得到预测的后验概率分布，而不仅仅是点估计（这是贝叶斯框架的优势之一，不是吗？）。当您的训练数据固定时，通常通过添加与观察变量类似的...
PyMC3贝叶斯线性回归预测与sklearn.datasets
我一直在尝试使用PyMC3和sklearn.datasets数据集中的真实数据（即不是线性函数+高斯噪声）来实现Bayesian Linear Regression模型。我选择了属性最少的回归数据集（即load_diabetes()），其形状为(442,10)，即442个样本和10个属性。我相...
PyMC3中的“chain”是什么？
我正在学习用PyMC3进行贝叶斯建模。您可以使用以下代码创建一个模型并进行采样：import pandas as pd import pymc3 as pm # obs is a DataFrame with a single column, containing # the observe...
如何在PyMC3中从狄利克雷过程中提取无监督聚类？
我刚刚完成了Python中的贝叶斯分析一书，作者是Osvaldo Martin（这是一本非常好的书，可以帮助您理解贝叶斯概念和一些花式numpy索引）。 我真的想扩展我的理解，以便于使用贝叶斯混合模型对样本进行无监督聚类。我所有的谷歌搜索都导向了Austin Rochford的教程，它真的很...
如何加速PyMC马尔可夫模型？
有没有办法加速这个简单的PyMC模型？在20-40个数据点上，拟合需要大约5-11秒。 import pymc import time import numpy as np from collections import OrderedDict # prior probability of...
能否在不加载原始模型到内存中的情况下，加载PyMC3跟踪并访问其值？
我仍在学习如何使用PyMC3，希望文档中已经涵盖了这些基础知识。基本思路是组合模型、采样生成后验分布并保存链。如果我按照Backends页面的建议来读取链（如trace = pm.backends.text.load('test_txt')），则会出现TypeError: No context...
我们如何将一个Python上下文管理器与出现在其块中的变量"关联"起来？
据我所知，上下文管理器用于在Python中为对象定义初始化和终止代码片段（__enter__和__exit__）。然而，在PyMC3的教程中，他们展示了以下上下文管理器示例： basic_model = pm.Model() with basic_model: # Priors ...
pymc3: 多个观测值
我有一些观测数据需要估计参数，想着试试PYMC3。我的数据结构是一系列记录。每个记录包含一对与固定一小时期间相关的观测值。一个观测值是在给定小时内发生的事件总数，另一个观测值是该时间段内成功的数量。例如，一个数据点可能指定在给定的一个1小时周期内，总共有1000个事件，其中100个是成功的。在...
PyMC的并行化
有人能够给出如何并行化 PyMC MCMC 代码的一般指导吗？我正在尝试运行这里提供的示例中的 LASSO 回归。我在某处读到并行采样是默认完成的，但我是否仍然需要使用类似于 Parallel Python 的东西使其工作？ 下面是一些参考代码，我希望能够在我的机器上并行化运行。x1 = n...
PyMC3和Theano - 在导入pymc3后停止工作的Theano代码
当我引入pymc3时，一些完美运作的简单theano代码停止工作。 以下是重现错误所需的一些片段：#Initial Theano Code (this works) import theano.tensor as tsr x = tsr.dscalar('x') y = tsr.d...