问题标签 [rjags]

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 投票
1 回答
237 浏览

r - 如何覆盖观察节点的值?

JAGS 的新手并使用本文中提到的模型。看起来作者使用了更接近 BUGS 的早期版本的 JAGS,因为在某些时候,这出现在模型块中(第 28-29 行):

JAGS 抱怨以下错误:

检查 JAGS 手册,错误很明显。在 A.4 数据转换部分中,它指出 BUGS 允许在关系的左侧放置一个随机节点两次。作为一种解决方法,它提供了将数据转换包含在单独的data块中。它仍然失败。

有没有人试图复制这项工作并取得成功?有什么提示吗?


下面是完整的 JAGS 模型。可疑的任务是z[i,K] <- 0beta[j,K] <- 0

0 投票
1 回答
436 浏览

r - rjags 模型负二项似然和伽马先验

我读入了我的数据。我制作模型字符串。我把它交给JAGS。我得到“节点 y[1] 中的错误 - 节点与父节点不一致”。

首先,我不知道我的模型是否正确。我什至找不到 JAGS 的基本文档。我实际上很惭愧地承认这一点,因为这应该像互联网搜索一样简单,但我找不到任何文件可以告诉我 1)如何设置 JAGS 模型或 2)什么样的功能/分布/参数是在 JAGS 中可用。我之所以能走到这一步,是因为我发现有人在做类似的模型。如果有人知道 JAGS wiki 或文档,那就太好了。

编辑:如果有人能告诉我 dnegbin 的参数是什么,那将是一个巨大的帮助。l当我为and插入随机数时r,它“工作”,因为它为anddnegbin(l,r)绘制数字,但我不知道它是否意味着什么。lr

0 投票
1 回答
265 浏览

r - for循环中的锯齿先验

我试图通过使用 for 循环来最小化我的 jags 代码。我的原始代码是

如您所见,我手动运行每个先前的程序,而不是我想在 for 循环中运行它们我的尝试;

有没有办法在 rjags 中做到这一点,因为我不确定 JAGS 软件中是否有 if 语句。

0 投票
1 回答
128 浏览

r - JAGS/BUGS 错误 - 可能涉及以下部分或全部节点的定向循环

我的数据集是关于点球的,包含 106 行,特征是:

  • 踢球方向
  • 球员的脚劲
  • HomeAway 告诉我们玩家是在他的城镇还是在客场比赛

我想对这些数据执行多项逻辑回归,以便根据其他两个数据建立一个踢球方向的模型。我以 Aligator BUGS 为例:http ://www.openbugs.net/Examples/Aligators.html

我的 BUGS 文件如下:

我使用 rjags 并出现以下错误:

我做错了什么 ?

提前致谢

0 投票
1 回答
100 浏览

jags - 使用 JAGS for R 的分层模型显示错误的均值数

在此处输入图像描述

我有 2 列数据,y并且grp我正在尝试创建一个JAGS如上所示的模型。grp是组,我有 5 个组。以下代码来自这里。我使用这段代码是因为标题下的描述Model and Data看起来几乎像这个分层模型。

但是mu当我查看摘要时,我只得到一个。应该有 5mu's个,每组一个。有人可以更正代码吗?您也可以指出其他地方提供的类似示例,我可以尝试对其进行修改。我在代码中遗漏了一些东西,我相信代码可能就像问题一样,但是当我像这样修改它时,即使有 5 种方法,我似乎也没有得到适当的方法。

不确定这个问题是否属于数学堆栈交换。

0 投票
1 回答
246 浏览

gaussian - update.jags(model, n.iter, ...) 中的错误:节点 sd[1] 中的错误父值无效

我在节点 sd[1] 中有错误,它说编译器中的父值无效。我正在使用 R 的“MASS”p 包中的“星系”数据的高斯模型。

0 投票
2 回答
570 浏览

r - 为什么与 rjags 和 R2Jags 拟合的模型输出不同?

我正在使用组级预测器拟合多级逻辑回归模型。我通过 R 使用 JAGS。当我将模型与包相匹配时,我会得到不同的runjags行为R2Jags

我试图编写一个可重现的示例来显示该问题。下面,我模拟来自二项式模型的数据,将数据索引到 8 个图和 2 个块,然后拟合多级逻辑回归以恢复下面代码中的成功概率 (b1b2)。滚动到底部以查看两个配合的摘要。

我的问题是:

  1. 为什么这两个拟合的后验不同?我使用相同的数据,单个模型规范,并在每个之前设置随机数生成器。为什么后验的平均值不同,为什么Rhat值如此不同?

包版本信息:

这是 rjags 拟合的输出:

这是 R2jags fit 的输出:

这是来自 2 次拟合的 mu.alpha 的跟踪图。首先,从 rjags 适合:

rjags fit 中 mu.alpha 的跟踪图

二、来自R2jags的合身:

来自 R2Jags 拟合的 mu.alpha 的跟踪图

0 投票
0 回答
119 浏览

r - JAGS 无法识别存储在 R 全局环境中的值

我正在使用 R2jags 包在 Windows 10 中使用 R v3.6.1 和 Rstudio v1.2 运行 JAGS。JAGS 似乎没有找到我在 R 中为 MCMC 设置(例如 n.iter、n.burn-in 等)创建的存储值,因为我的代码:

产生错误

如果我用数字(n.burnin = 10000)替换存储值的名称,它将返回相同的错误,但对于不同的 MCMC 设置。我假设 JAGS 也将无法找到我存储在全局环境中的列表和数据框。我认为我系统的防火墙阻止了 JAGS 访问全局环境,所以我禁用了我的防病毒措施,并以管理员身份运行了 R 和 Rstudio。这并没有解决问题。

0 投票
1 回答
165 浏览

bayesian - 如何使用 Rjags 随机初始化链?

我正在尝试随机初始化贝叶斯模型的链,今天我手动插入值,如下所示:

我尝试通过以下方式执行它,但导致错误:

其中p->3(模型的参数个数)

跟随错误:

错误:使用 rjags 编译和调整模型时出现以下错误:setParameters(init.values[[i]], i) 中的错误:运行时错误:为 alpha 提供的值中的尺寸不匹配

0 投票
1 回答
49 浏览

mcmc - 为先验定义约束时,JAGS 中出现错误消息尝试重新定义节点

我正在尝试从先前的分布中进行采样并将某些总和和差异限制为大于 0,但是,当我运行通常在 OpenBugs/WinBugs 中运行的此代码时,我收到以下错误“尝试重新定义节点”。我想知道是否有任何可能的方法来定义我的约束以防止发生此错误。任何帮助将不胜感激。

谢谢