-
13 votesanswersviews
使用sklearn.datasets进行PyMC3贝叶斯线性回归预测
我一直在尝试使用来自 sklearn.datasets 中数据集的真实数据 PyMC3 来实现 Bayesian Linear Regression 模型（即不是来自线性函数高斯噪声） . 我选择了具有最小数量属性的回归数据集（即 load_diabetes() ），其形状为 (442, 10) ;也就是 442 samples 和 10 attributes . 我相信我的模型正常工作，后面看... -
0 votesanswersviews
使用pymc3和Theano导入错误
我一直在尝试使用Spyder IDE的Pymc3 . 我在https://anaconda.org/anaconda/pymc3中按照命令安装了pymc3 . 安装看起来很好 . 包装theano也包括在内，所以我没有单独安装 . 但是当我试图导入pymc3和theano时，我遇到了以下错误： for pymc3 ImportError：无法导入名称时间序列 for theano Import... -
2 votesanswersviews
用PyMC3预测贝叶斯线性回归中的新数据后验
使用PyMC3执行贝叶斯线性回归 . 我构建了我的模型，我想使用相同的模型预测新X值的后验 . 我一直在尝试按照文档网站上的说明进行操作：https://pymc-devs.github.io/pymc3/notebooks/posterior_predictive.html（参见预测） . 这包括在分析之前使您的X值成为theano共享变量，然后在模型构建之后更改值，并运行 run_ppc() ... -
1 votesanswersviews
PyMC 3中的确定性变量和Fortran Scipy函数
我正在尝试构建一个简单的PyMC 3模型，其中我估计两个切割点和潜在双变量高斯密度中的相关参数，产生（多项式）计数向量的四个预测概率 . （我希望，这将成为一个更大模型的一部分，在这个模型中，对于许多潜在的多元高斯密度估计这些和其他参数 . ） 因此，我想将切割点cx和cy建模为正常的随机变量，并将相关参数rho建模为缩放的Beta随机变量（作为旁注，我希望听到更好的方法来处理rho - 做P... -
4 votesanswersviews
使用PYMC3回归
我在这里发布了一个IPython笔记本http://nbviewer.ipython.org/gist/dartdog/9008026 我通过标准Statsmodels OLS和PYMC3与Pandas提供的数据一起工作，顺便说一下，这部分工作得很好 . 我看不出如何从PYMC3中获得更多标准参数？这些示例似乎只是使用OLS来绘制基本回归线 . 看来PYMC3模型数据应该能够给出回归线的参数吗？除... -
0 votesanswersviews
PyMC3中的加权回归
我正在尝试找到适合二项式调查数据的函数的参数 . 每个数据点实际上是大小为n的样本 . 如果我使用具有贝叶斯预测的动态模型（如West＆Harrison的贝叶斯预测和动态模型第14.4章），我可以使用β二项式共轭先验方案来解释不同的样本大小 . 但是，我还希望以包含 data.n 列的方式在PyMC3中执行此操作 . 我已将模型设置如下： import pymc3 as pm import pan... -
4 votesanswersviews
使用PyMC3计算最大似然
有些情况下，我是Hessian . PyMC3具有执行此操作的功能，但 find_MAP 似乎以转换形式返回模型参数，具体取决于它们上的先前分布 . 有没有一种简单的方法可以从这些中获得未转换的值？ find_hessian 的输出对我来说更不清楚，但它最有可能也在转型空间中 . -
4 votesanswersviews
从pymc3中的推断参数预测
我试图从非贝叶斯背景中理解这一点 . 在线性回归或黑盒机器学习工具中，工作流程类似于以下内容 . 获取数据 准备数据 模型数据（从中学习[或部分，训练集]） 测试模型（通常在测试集上） 如果模型根据某个指标是好的，转到6，否则调查和修改工作 . 模型足够好;用它来预测/分类等 因此，假设我使用pymc3来了解广告支出与销售商品收入之间的关系 . 如果从1到5的所有阶段都... -
3 votesanswersviews
用PyMC3进行基本贝叶斯线性回归预测
当新数据可用时，我想使用我的PyMC3 LR模型获得预测变量 y 的80％HPD范围 . 因此，为我的原始数据集中的 x 的新值推断 y 的值的可靠分布 . 模型： with pm.Model() as model_tlr: alpha = pm.Normal('alpha', mu=0, sd=10) beta = pm.Normal('beta', mu=0, sd=10) ... -
0 votesanswersviews
在建模PyMC3中的序数预测变量时使用度量预测器
我试图在PyMC3中实现Doing Bayesian Data Analysis（Kruschke）第23.4章的有序概率回归模型 . 采样后，截距和斜率的后验分布与本书的结果不具有可比性 . 我认为模型定义存在一些基本问题，但我没有看到它 . 数据：X是度量预测器（标准化为zX），Y是序数结果（1-7） . nYlevels3 = df3.Y.nunique() # Setting the t... -
0 votesanswersviews
从PyMC3中的混合物中采样点时获取分布
我有一个带有 pm.NormalMixture() 的模型，当我从普通混合物中采样时，我也想知道哪些混合分布是从中采样的 . import numpy as np import pymc3 as pm obs = np.concatenate([np.random.normal(5,1,100), np.random.normal(10,2,200)... -
6 votesanswersviews
PyMC3中的链条是什么？
我正在学习PyMC3进行贝叶斯建模 . 您可以使用以下命令创建模型和示例： import pandas as pd import pymc3 as pm # obs is a DataFrame with a single column, containing # the observed values for variable height obs = pd.DataFrame(...) #... -
4 votesanswersviews
与pymc3和信念的线性回归
我试图用 pymc3 掌握Bayesain统计数据 我运行此代码进行简单的线性回归 #Generating data y=a+bx import pymc3 import numpy as np N=1000 alpha,beta, sigma = 2.0, 0.5, 1.0 np.random.seed(47) X = np.linspace(0, 1, N) Y = alpha + beta*...