问题标签 [hierarchical-bayesian]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
42 浏览

machine-learning - 使用数据驱动方法进行信息检索的贝叶斯网络建模?

我正在探索使用数据驱动方法构建用于信息检索的贝叶斯网络的方法和代码。我确实发现数据集或代码不可用的非常旧的论文。我是新手,正在探索这个领域。请,如果有人可以提供代码链接或最新论文的建议,可以帮助为贝叶斯网络建设提供实施。

0 投票
1 回答
76 浏览

r - 分层狄利克雷回归(锯齿)中的随机截距

我有以下数据结构:

  • y:3 列,这些列是多年来观察到的死亡比例。
  • x1:GDP - 与每年相关的连续变量
  • x2:与死亡有关的年龄

这里的型号规格:

模型

这里模拟数据:

锯齿模型

我收到了这个错误:

对模型规范有什么建议吗?

0 投票
0 回答
45 浏览

bayesian - 为 mcmc 样本绘制散点图矩阵的原因是什么?

我发现一些贝叶斯论文尝试绘制参数的散点图矩阵。

我只是想知道绘制这个散点图矩阵的目标是什么?

如果我看到样本之间存在一些线性或非线性关系,这意味着什么。

谢谢 在此处输入图像描述

0 投票
0 回答
132 浏览

projection - WinBugs 中的“此初始值不对应于随机节点”

我给了所有随机节点的初始值,但 WinBUGS 仍然给我这样的信息

此初始值不对应于随机节点

我在这里缺少什么节点?

0 投票
0 回答
8 浏览

jupyter - Windows 10 更新:是什么导致我的分层建模代码运行缓慢而痛苦?

我最近在夜间强制更新了 Windows 10,当我打开 Jupyter Notebook 运行一些贝叶斯建模时,它变得非常缓慢。我正在使用马尔可夫链蒙特卡洛,每个样本需要 3 秒以上。我以前没有这个问题。

没有成功关闭并重新启动内核。我认为这可能是 Jupyter Notebook 的问题,所以我尝试在 Jupyter Lab 和 Spyder 中运行代码,他们遇到了同样的问题。

假设它与更新有关,我试图恢复到旧的还原点,但不知何故,这些都没有被记录,所以我没有想法了。

我将 py36 用于特定包(HDDM)。还有其他人有这个问题吗?

0 投票
1 回答
72 浏览

python - 从具有 2 个预测变量的分层模型中绘制多个散点图

0 投票
0 回答
33 浏览

bayesian - 如何避免过度分散的泊松回归过拟合?

我有一个数据集,包括三个变量,包括公司 id(有 96 家公司)、专家 id(有 38 位专家)和专家给公司的分数。点是从 0 到 100 的离散值。我尝试将过度分散的泊松拟合到专家给出的模型点。但是我不知道为什么模型会过度拟合,尽管我使用的是线性似然。这是我的 JAGS 代码:

任何人都知道为什么这个模型过拟合以及如何修复它?

0 投票
1 回答
54 浏览

r - 有没有办法在 JAGS 中获取和存储矩阵元素的位置?

我正在用 JAGS 中的 BUGS 代码在 R 中开发贝叶斯层次模型。

在我的模型中,我有两个矩阵,它们在相同的矩阵位置包含彼此相关的信息。我的信息按行排列。我按行对第一个矩阵 Distmat 应用数学运算:

我有兴趣在新向量中记录每行 diffmat 中每个最小值的列位置,然后将此向量应用于第二个矩阵。这在使用函数“which”或“which.min”的常规 R 代码中相对容易:

然后将向量“a”应用于第二个矩阵(terrmat)以获得与 Distmat 位置相关的值:

但是,显然 BUGS 代码无法识别 which 或 which.min(),我正在努力寻找一种将这些矩阵行位置存储在向量中的方法。也许有一个非常简单的解决方案,但我真的被困在那里。希望我的信息足够清楚。

任何建议将不胜感激。谢谢你的时间!

0 投票
1 回答
88 浏览

r - 分层狄利克雷回归(锯齿)...过拟合

早上好,我需要社区帮助以了解编写此模型时出现的一些问题。我的目标是使用“log_GDP”(对数尺度的国内生产总值)和“log_h”(对数尺度每 1000 人的医院床位)作为预测因子来建模死亡比例的原因

  • y:3 列,这些列是多年来观察到的死亡比例。
  • x1:“log_GDP”(对数刻度的国内生产总值)
  • x2:“log_h”(每 1,000 人的对数病床数)

从上图中的估计结果可以看出,我得到了很高的噪声水平。在我只使用一个协变量即 log_GDP 工作的地方,我得到了平滑的结果

这里的型号规格:

在此处输入图像描述

这里模拟数据:

锯齿模型

收集系数并估计比例

预测和观察。价值数据库

阴谋

在此处输入图像描述

0 投票
1 回答
179 浏览

r - “变量名中不允许有空槽”(OpenBUGS,R2OpenBUGS)

我是 OpenBUGS 的初学者,我通过 R2OpenBUGS R 包使用它。我尝试设置状态空间模型来识别非常嘈杂的数据中的对数正态信号。经过多次试验和错误,我设法获得了此代码,但我仍然收到以下错误消息:“变量名错误 pos 664 中不允许空槽”,我不明白。谁能知道代码有什么问题?

免责声明:

  • alt = 测量高度
  • true_alt = 我尝试评估的内容
  • nbird = 个体数量
  • nobs = 观察次数(这个数字对于每只鸟都不一样)
  • nstate = '飞行状态',这是鸟类的行为方式(nstate = 3,因为有 3 种不同的行为)

我尝试确定每个状态的 true_alt 的对数正态分布。

谢谢!!!