我想在 PyMC3 中实现基本的 model 捕获和重新捕获（捕获 100 只动物并标记它们，然后释放它们并在它们混合后重新捕获其中的 100 只并注释标记了多少）。 这是我的代码： 但是我在plot_ppc中得到错误'var names: "[\'likelihood\'] are not pre ...
运行代码 az.plot_trace(result_final); 面临以下错误TypeError: expected dtype object, got 'numpy.dtype[float64]' 上述异常是以下异常的直接原因： SystemError: CPUDispatcher(<0 ...
我有一个来自 PyMC3 的样本，我试图从中获得累积概率，例如 P(X < 0)。 我目前使用这个： 有没有更好的方法来做到这一点，或者使用来自 Arviz 或 XArray 的一些帮助程序 function？ 我还没有找到任何.cdf()方法或类似的东西。 奇怪的是缺少这些基本功能，但有更高 ...
以明显的方式（ conda install -c conda-forge arviz ）安装ArviZ后，出现以下故障。 如果有人能阐明如何解决这个问题，那就太好了。 按照@merv 的建议，我得到以下内容，它解释了问题（它是一个不同的 package，称为 snappy，用于计算拓扑）。 但是我该 ...
我正在尝试将非对称高斯拟合到我的数据中。 我的数据只是一个名为 wave (x) 的 numpy 数组和一个名为 spec (y) 的 numpy 数组，看起来像一个非对称高斯。 这是带有非对称高斯拟合曲线拟合数据的图像（这也有连续统，但这现在并不重要。 这是函数： 我正在使用此代码来拟合参 ...
我正在使用 python 版本 3.8.5 的 Spyder IDE，我使用的是 pymc3（不记得确切版本），它工作正常，升级 scipy、spyder 和 python 后，一切都变得混乱。 我花了过去 2 天的时间试图解决这个问题，每次我安装 pymc3 并运行程序时，我都会收到关于 thea ...
我正在使用 python 中的 pyMC3 库估计 model。 在我的“真实”model 中，有四个参数 arrays，其中两个参数超过 170,000 个。 总结这组参数在我的计算机上的计算量太大。 我一直在试图弄清楚 arviz 中的汇总 function 是否允许我仅汇总数组中的一个（或少量 ...
我是 python 和 pymc3 的新手。 如果这个问题相当愚蠢，请原谅我。 我有一个名为 toy_data.csv 的数据集： 批号 batch_id 积分 pred1 pred2 12-150 1 1 70.26 10.766 12-150 1 2 65.72 10.51 ...
我正在尝试生成我的 model 的跟踪 plot 但它显示module 'pymc3' has no attribute 'traceplot'错误。 我的代码是： 然后给出以下错误： 我在windows10上，我已经下载了带有pip的pymc3，因为它不包含在我下载的anaconda中。 ...
通常，我发现我们使用pm.binomial()等等。 并使用 pm.plot_posterior 我们得到离散图，我的问题是为什么我必须得到的图是曲线？ 有人可以启发我吗？ 我想得到： 但我得到了图表： ...
我正在使用 PyMC3 学习统计反思课程。 在第 4 章的末尾，他们要求提供不在原始 (.Kung) 数据集中的数据点的单个值的 HDI？ 在 PyMC3 中可以这样做吗？ 在 scikit-learn 中，您有fit()和predict() ，并且您能够预测 output 以获得全新的输入。 使 ...
我有一篇旧博文，我正在训练 PyMC3 模型。 您可以在此处找到博客文章，但模型的要点如下所示。 with pm.Model() as model: mu_intercept = pm.Normal('mu_intercept', mu=40, sd=5) mu_slope = pm ...
我有一个数据集，我们称它为 A，当我打印时它看起来像这样： 但现在我想提取 158.0 和 725.2 并将它们保存在一个单独的数组中，但我不知道如何获取这些值。 我试过 Ax、A.get() 等。谢谢！ ...
假设我使用 pymc3 为伽马分布生成了一些样本数据： 现在，我将创建一个模型来拟合该数据的伽马分布： 这将正确获得创建原始假数据的值和分布。 现在，我想绘制 pdf（但我不知道该怎么做！）。 我看到了一个这样做的例子： 它创建了一个矩阵，其中包含来自估计 pdf 的样本。 我用以下 ...
我正在玩 PyMC3，试图在 PyMC3 文档中适应采矿灾难开关点 model的修改版本。 假设您有两个煤矿（mine1 和 mine2），每个煤矿在相同的年份范围内具有相似的灾难计数。 然而，mine1 在实施降低灾难次数的安全程序变更方面晚了 5 年： 我有兴趣测试更好的 model 是否适合 ...
我正在尝试更改我正在使用 Arviz 创建的 plot 的标题。 通常我已经使用backend_kwargs完成了此操作，但它似乎不适用于某些地块：在本例中为plot_density 。 这是我正在使用的代码： 它识别出kwargs的高度变化，所以它不像没有注意到它们在那里。 我还尝试了其他变体，例 ...
我正在尝试使用 pymc3 和 arviz 获得 2D mcmc 链的有效样本大小 上面的代码适用于 1D，但不适用于 2D，我看到有一个 az.convert_to_dataset 可能会有所帮助，但我不知道如何使用它？ 样本将是一个 N x 2 数组，它应该只给出一个数字，如 output 谢谢 ...
我正在使用 python 研究贝叶斯分析一书，并尝试在 python 中实现所有代码，为此我使用 google colab，但是当我运行书中的以下代码时 它给出了这样的错误： 在我尝试实现这些库之前 但仍然没有成功。 如何修复此错误？ ...
我正在尝试为多个组可视化具有最高后验密度（hpd）的简单线性回归。 但是，我在为每个条件应用 hpd 时遇到问题。 每当我运行这段代码时，我都会为每个条件提取相同的后验密度。 我想可视化与其条件相对应的后密度。 如何为每个组提供 plot hpd？ 编辑：问题已在PyMC3 话语中得到解决 ...
我需要做些什么才能在 PyCharm 中显示我的 PyMC3 plots.traceplot结果吗？ 我刚刚发现我需要一个pm.summary()周围的print()语句。 所以我想知道这些情节是否还有另一个 function 。 谢谢你！ ...