使用PyMC3进行增量式模型更新
在pyMC3中是否可以逐步更新模型？我目前找不到相关信息。所有文档都是使用先验已知数据进行操作的。 但是据我理解，贝叶斯模型也意味着能够更新信念。这在pyMC3中可能吗？在哪里可以找到相关信息？ 谢谢 :)
使用PyMC3计算最大似然
有时我对Bayesian推断的完整后验分布并不感兴趣，而只关心最大似然(或最大后验概率，如果先验设定合适)，以及可能的Hessian矩阵。PyMC3有相关函数可以实现，但find_MAP函数返回的模型参数可能会因为先验分布而被转换。是否有一种简单的方法可以得到这些未经变换的值呢？find_he...
用Pymc3编写自定义似然函数
我正在努力使用pymc3实现具有自定义似然函数的线性回归。 我之前在CrossValidated上发布了这个问题，并被建议在这里发布，因为这个问题更加代码化（关闭帖子here） 假设你有两个独立变量x1、x2和一个目标变量y，以及一个名为delta的指示变量。 当delta为0时，似然函数是...
如何制作截断正态先验：将Pymc2转换为Pymc3
在pymc3中如何配置截断正态先验分布？在pymc2中，这很简单（如下所示），但在pymc3中似乎不再提供截断正态分布。 Pymc2: TruncatedNormal('gamma_own_%i_' % i, mu=go, tau=v_gamma_inv, value=0, a=-np....
PYMC3 季节性变量
我相对较新于PYMC3，正在尝试实现没有回归变量的贝叶斯结构时间序列(BSTS)模型，例如R中的这里所拟合的模型。该模型如下： 我可以使用GaussianRandomWalk实现局部线性趋势模型，如下所示： delta = pymc3.GaussianRandomWalk('delta...
使用PyMC3创建的模型如何绘制成图形
我正在使用以下代码来创建一个简单的PyMC3模型： import pymc3 as pm import theano.tensor as tt with pm.Model() as model: p = pm.Uniform("freq_cheating", 0, 1) p...
在pymc3中自定义似然函数
我如何在PyMC3中定义自定义似然函数？在PyMC2中，我可以使用@pymc.potential。我尝试在PyMC3中使用pymc.Potential，但是似乎无法对参数应用布尔运算（当我这样做时会收到类似于此类的错误）。例如，以下代码无法运行： from pymc import * wi...
如何使用Pymc3适配属于实例的方法？
我无法使用PyMc3将属于类实例的方法作为确定性函数进行拟合。您能向我展示如何做到这一点吗？ 为了简单起见，以下是一个简单示例总结了我的情况。实际上，我的限制是一切都通过GUI完成，像“find_MAP”这样的操作应该在与pyqt按钮链接的方法内部。 我想要将“FunctionIWantT...
PYMC3贝叶斯预测锥
我仍在学习PYMC3，但在文档中找不到有关以下问题的任何信息。考虑没有季节性的贝叶斯结构时间序列（BSTS）模型this question。这可以通过以下方式在PYMC3中建模： import pymc3, numpy, matplotlib.pyplot # generate some t...
pymc3: 多个观测变量的分层模型
我有一个简单的分层模型，其中有很多个体，我从正态分布中获取了小样本。这些分布的均值也遵循正态分布。 import numpy as np n_individuals = 200 points_per_individual = 10 means = np.random.normal(30, ...