问题标签 [openbugs]

我最近开始尝试使用 OpenBUGS 和 R 进行一些贝叶斯建模R2OpenBUGS。

我关注了 David Eagles 的安装 wine 和 OpenBUGS 的一个很好的链接，那里发布的学校示例运行良好。

然后，我尝试运行我拥有的一些代码，并在the model is syntactically correct告诉我数据未正确加载后不断收到错误消息。

经过几天的故障排除后，似乎当数据文件超过一定长度时，最高位（您通常list在 OpenBUGS 中突出显示以加载数据）丢失并且 OpenBUGS 无法加载数据。

我的 OpenBUGS 程序的屏幕截图，其中包含已加载的列表部分已被切断的数据

这张图片显示，在我的 data.txt 文件的顶部，不再有一个list{基本上停止了我运行模型的所有尝试。

为了确认这一点，我运行了一个简单线性回归的示例，其中我将点数设置为 1000，然后设置为 100：

该模型运行良好，linedata <- list(Y = c(1:100), x = c(1:100), N = 100, xbar = 3)但在加载数据阶段时崩溃了linedata <- list(Y = c(1:1000), x = c(1:1000), N = 1000, xbar = 3)。

我已经用谷歌搜索了 wine 和 OpenBUGS 的数据文件大小限制，但没有发现任何有用的信息。

任何人有任何关于尝试什么/从哪里开始/以前体验过的建议？

我在 OpenBUGS 中初始化以下模型时遇到问题

这是我的数据

这是我姓名缩写的列表

模型检查，数据加载和编译步骤都可以。加载首字母时，OpenBUGS 表示已加载初始值，但链包含未初始化的变量。然后我也尝试初始化 theta0，但问题仍然存在。有人可以帮我解决这个问题吗？谢谢哈立德

我是 R2OpenBUGS 的新手，非常神秘的错误非常令人沮丧。

我尝试运行一个非常简单的模型。我之前曾成功运行过类似的模型。我的问题是因为我有一个二维数组（矩阵）吗？

我尝试简化模型但没有成功。

以下是错误：

这是我使用的代码

数据文件：

我对 JAGS/OpenBUGS 完全陌生，所以在指定我的模型时，我真的很感激朝着正确的方向前进。我正在使用由 103 个国家在 15 年内编制的不平衡纵向数据，在这种情况下选择 12 年。DV 是基尼系数，它不应该被建模为 log-Normal，而可能是 Beta，尽管现在的重点是了解如何在 JAGS 中编译模型。我暂时使用的是固定效应模型。

我正在运行的数据和代码：

而且我不断收到此运行时错误：

任何关于我应该做什么不同的建议将不胜感激！我知道有 3 个“技巧”可以应用于不平衡数据，但我仍然对所有这些工作原理感到有些困惑，即 JAGS 如何读取数据输入。

干杯

Ĵ

我正在尝试使用 R2OpenBUGS 库从 R 评估分层模型。我已经看过类似问题的答案，但似乎没有任何效果。如果有人可以提供帮助，我将非常高兴。

模型如下。

在日志中我得到：

当我尝试将先验均值向量priora[]和priord[]赋给a[i]和d[i]时，问题就开始了。如果我将这些手段设置为“0”，则代码有效。

我想计算 x1b1+x2b2+......+xpbp。但是，p不知道，我不能提前写下来。所以我想做以下事情。但在 BUGS 中，它表示多个定义错误。

有没有办法在错误模型中完成它？

我正在尝试在 OpenBUGS 中运行马尔可夫模型来估计过渡概率。但我收到一个错误“数组索引不是整数”这是我的模型代码：

这是我的 R 代码数据：

谁能帮我解决错误？

我正在尝试在 OpenBUGS 中运行模型。但我收到错误“未定义变量 TT”。这是我的错误代码：

这是我运行 BUGS 代码的数据和 R 代码：

任何人都可以帮助解决这个问题吗？提前致谢。

我正在使用 R 中的 BRugs 包运行模型以进行贝叶斯推理。该模型在语法上是正确的并且数据加载，但是 modelCompile 函数永远运行，我在 15 小时后停止了它，但我不知道问题可能是什么。

数据文件很大，但我可以根据要求提供...

我是 OpenBUGS 的新手，在使用 logit() 函数拟合模型时遇到了一些问题。

阅读周围我发现一个可能的解决方案是明确指定 logit 函数而不使用 WinBUGS 自己的 logit 函数：

在更复杂的模型中，我们在使用 WinBUGS 自己的 logit 函数时经常遇到问题，例如实现收敛（实际上，即使使用相当简单的模型也可能出现问题。）。因此，通常最好通过 logit.p[i] <- log(p[i] / (1 – p[i])), p[i] <- exp(logit.p[i ]) / (1 + exp(logit.p[i])) 或 p[i] <- 1 / (1 + exp(- logit.p[i]))。

（更多信息在这里http://www.mbr-pwrc.usgs.gov/software/kerybook/AppendixA_list_of_WinBUGS_tricks.pdf第 14 点。）。

问题是我不明白该怎么做，假设我使用 WinBUGS 集成 logit 函数的原始似然函数是：

我如何明确地写出来？

非常感谢。

文森佐