问题标签 [winbugs]

我有一个 OpenBUGS 模型，它使用随时间 (x.values) 观察到的数据 (y.values) 来模拟多次运行 (~100000)，并为每次运行提供新的 y 值 (y.est) 估计值。观察到的数据显示出从最大值明显下降。

我想跟踪每次运行从最大丰度 (T.max) 下降到最大丰度的 10% (T.10%) 所需的时间长度。因为最大丰度值随运行而变化，该最大值的 10% 也将随运行而变化，因此 T.10% 将随运行而变化。

设置一个参数来存储 T.max 很容易，它不会因运行而异，因为最大值远远大于任何其他值。

我想不通的是如何存储 y-est 值和 T.10% 的交集。

我的第一次尝试是使用以下函数确定每个 y-est 值是高于还是低于 T.10% step()：

这会为每个 y.est 值生成一串一和零（例如，0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1 , 1, 1, 1, 1, 1 等）如果每次运行只是从最大值连续下降到最小值，我可以使用该rank()函数来确定有多少above.below[i]值出现在 T.10% 以上：

在此示例中，decline.length将等于上述字符串中“0”的数量，即 9。不幸的是，y-est 值偶尔会在下降到 T.10% 以下后显示增长期。因此，above.below值向量可以如下所示： 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 , 1, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1 等。因此，decline.length给定向量中的后续 0，将等于 14 而不是 9。

我想要做的是弄清楚如何在above.below第一个“1”之前只存储“0”的数量；above.below[1:10]而不是above.below[1:N]. above.below不幸的是，第一个“1”并不总是出现在第 10 个时间步长，所以我需要在模拟过程中使每次运行的变化范围达到最大。

我正在努力在 OpenBUGS 中实现这一点，因为它是一种非过程语言，但我认为它可以做到，我只是不知道这样做的诀窍。我希望更熟悉step()和rank()功能的人可以提供一些专家建议。

非常感谢任何指导！

我正在尝试使用在 Windows 7 中运行的 R 和 R2WinBUGS 在 WinBUGS 中运行一些模拟。

在我准备好所有文件 a 并启动函数 bugs() 后，WinBUGS 窗口会打开，但不会启动模拟。旧版本的 R 和 OpenBUGS 也会发生同样的情况。是Windows的问题吗？

这是代码：

任何帮助，将不胜感激！干杯马里奥

如果我在 OpenBUGS 中运行以下代码：

我希望item为 0，因为没有任何组件variable可以低于 1。这是基于rankOpenBUGS 用户手册中的函数定义：

" rank(v, s):v小于或等于s"的组件数

然而，for的后验item实际上是：

这意味着大约 25 个分量variable小于或等于 1，显然没有一个分量variable可以小于 1——由variable值的后验分布证实：

这里发生了什么？

我使用了 OpenBUGS，它产生了 MCMC 输出的 coda 文件。要计算和绘制 Gelman Rubin 和 Geweke 诊断，我需要将此 coda.odc 文件转换为 R 中的 mcmc 对象吗？有没有办法做到这一点？或者你是否推荐我一些其他的方法来做这个分析？

谢谢

我给了随机节点的所有值，但 WinBUGS 仍然给我一条链包含未初始化变量的消息。当我尝试生成它时，BUGS 给了我一个错误undefined real result。我在这里缺少什么节点？

我有附加的模型（由 Rasmus Bååth 提供），我正在努力改进它以增加我对 MCMC 和 R 的理解。模型中输入的是足球比赛，每行包含（HomeTeam、AwayTeam、Season、HomeGoals、AwayGoals、MatchResult ，家庭表格，离开表格）。该模型用于预测两支球队之间的比赛中的足球结果，其中 HomeGoals 和 AwayGoals 被建模为泊松分布。过去几天我一直在努力完成的是使用这些关于模型中团队当前形式的数据，但我想不出一种方法来实现它。我希望就如何做到这一点提出一些建议。

一些澄清：

• 赛季是比赛所在的赛季
• 如果客队获胜，MatchResult 为 -1，如果平局，则为 0，如果主队获胜，则为 1。
• HomeForm 和 AwayForm 分别是量化主队和客队当前状态的值（-1 到 1，球队最后一场比赛的滚动平均比赛结果，窗口为 5，如果球队踢客场则为负比赛结果） .
• 文本中的列在传递给模型（因子）之前被转换为整数。

谢谢！

模型：

数据（部分）：

我想构建一个模型

我收到一个错误：Attempt to redefine node f[1]

有什么方法可以在 JAGS/WinBUGS 中编写矢量文字？我在 JAGS 中尝试了这些，都导致错误：

我不想一一分配。

在这里或手册中实际上没有找到任何有用的东西。

我有一个简单的贝叶斯分层模型（带有随机截距的线性混合模型），应该很容易运行。问题是在成功加载模型和数据后，我尝试编译模型时出现以下错误，“array index is greater than array upper bound for a.” 看起来这应该很容易解决，但我反复检查了索引和数据，但找不到问题。我已经在不同的示例数据集上进行了实验，并且没有任何问题，这让我认为除了索引之外的其他问题是问题所在。

任何建议将不胜感激！

假设我们有以下语句：

让我们这么说吧Y[1] = 5。jags 会评估所有pi[1,1:100]节点，还是唯一需要的节点，即pi[1,5]？

根据我的经验，似乎 JAGS 评估所有父节点的效率很低，因为在我摆脱 dcat 后，我​​的模型加速了 3 倍。我必须使用多个 for 循环来处理不同的结果Y[i]。

现在我意识到，dcat在 JAGS 中实际上不需要sum(pi[]) = 1，并且dcat它将归一化pi[]以使其总和为 1。这意味着它必须评估所有节点。

这是非常可悲的。是否有任何智能等价物dcat只会评估唯一需要的一个父节点？WinBUGS 和斯坦呢？