问题标签 [winbugs]

我使用 R 中的 R2WinBUGS 包在 WinBUGS（百万次迭代）中运行了一个 looong 计算：

但是R崩溃了。如何在不再次运行 winbugs 的情况下再次将 bugs.object 重新加载到 R 中？我试过这个（我有 3 条链）：

但是out数据结构与 bugs 对象完全不同（因为它是不可用的）。我尝试将其转换为as.bugs.array：

但它不起作用。请帮忙。谢谢。

我正在尝试进行分层贝叶斯分析，但在使用 R 和 WinBUGS 代码时遇到了一些问题。我没有平衡的数据，并且正在努力编码。我每天使用 iButtons（温度记录设备）在横断面上收集温度数据，并试图生成一个将其与遥感数据相关联的模型。不幸的是，每个样带都有不同数量的 iButton，因此在样带（j）中创建按钮（i）的 3D 矩阵，在第（t）天重复“采样”对我来说是一个问题。

最终，我的模型将类似于：

级别 1 Temp[ijk] ~ N(theta[ijk], tau) theta[ijk] = b0 + b1*x1 + . . . + bn*xn

级别 2 b0 = a00 + a01*y1 + 。. . an*yn b1 = a10 + a11*y1 ...

3 级（也许？） - 2 级随机拦截

通常我会做这样的事情： Wide <- reshape(Data1, idvar = c("iButton","block"), timevar = "julian", direction = "wide")

然后我可以有一个 3D 矩阵，我可以像这样在 WinBUGS 或 OpenBUGS 中循环：

无论如何，不​​要担心上面代码的细节，它只是作为其他分析的示例。我的主要问题是，当我没有平衡设计且每个样带具有相同数量的 iButton 时，如何进行此类分析？任何帮助将不胜感激。我显然是 R 和 WinBUGS 的新手，以前没有太多的计算机编码经验。

谢谢！

哦，这是长（堆叠）格式的数据：

我正在尝试从 R 调用 WinBUGS，但我并不是很成功。

我在 R 中输入了这段代码：

然后，我调用 WinBUGS：

我在 WinBUGS 中收到以下错误：

如果有人可以帮助我，我将不胜感激。我的文件在同一个

我可以在 WinBUGS14 中毫无问题地运行模型并获得结果，但是当我从 MatLab 运行相同的 WinBUGS 模型（如下所示）时出现问题。它看起来程序永远不会停止运行并且没有结果返回。

谁能帮我。任何建议将不胜感激。谢谢。

1) 我的 WinBUGS 代码 --- CHK_model.txt

2) 我的 MatLab 代码

正如这篇文章的标题所说，当我尝试使用BRugsFit(with coda=T) 从 R 运行在 WinBUGS 中正常工作的代码和数据时，我收到以下错误：

当我tail()在 coda 对象上做时，我一遍又一遍地得到相同的数字。另一方面，当我运行 WinBUGS、保存结尾并将其加载到 R 中时，我得到了一些我预期的随机变化，并且没有关于收敛的警告。

这是我的模型（它使用“一个技巧”来找到 Logistic-Makeham 分布参数的后验）。

这是我的数据的典型示例：

...和典型的初始化（我使用 4 个链，细化 20，burnin 2000、20000 次迭代）

OpenBUGS 是否只是简单地舍入到比 WinBUGS 更少的有效数字，如果是这样，也许我可以更改一个设置以使其停止这样做？

首先，有人知道 WinBUGS 的开发人员指南吗？该网站充满了 Doodles 的详细示例和模型语言的文档，但我还没有找到任何关于如何解释陷阱窗口的信息。

其次，有没有人找到任何方法来简化检查/加载/编译/初始化/监控/更新周期？我的意思是，似乎没有任何方式可以说“不要费心重新检查模型或将任何设置恢复为默认值（！！！），只需继续从这些文件加载​​数据，从这些文件，并为每个文件生成一个新的尾声”。甚至标准的 Windows 快捷方式在这里也被取消了，迫使用户不断地点击并用相同的值一遍又一遍地填充相同的字段。这似乎是一个小问题，但是当您一个接一个地进行许多类似的分析时，它很快就会变老。

我正准备使用 TRON.EXE 向程序发送虚假的鼠标点击，但在走到那个极端之前，我希望有一些本机和更优雅的方式来自动执行重复的 WinBUGS 任务。

我正在阅读这篇论文（旨在使用贝叶斯分层方法对纳米线 (NW) 生长进行建模）。在第 7 页，作者提出了一个模型来描述纳米线的生长。我正在尝试编写 WinBUGS 代码（请参见下文），但我不确定它是否正确。

我从论文中引用了这一点作为摘要：

案例研究将展示不确定性下 NW 增长过程的分层建模和估计过程。...其中的数据是在六种生长条件（T = 365°C、380°C、400°C、420°C、430° C，440°C）。前四个条件用于模型构建。我们不考虑这两种高温条件，因为在 180 秒和 900 秒时没有观察到。“</p>

所以我们要处理三个变量，包括时间（以秒为单位）、温度（以摄氏度为单位）和纳米线的长度（以厘米为单位）。有四个温度级别（365°C、380°C、400°C 和 420°C）和四个时间点（15 秒、30 秒、180 秒、900 秒）。我没有每个条件下的确切长度值，因此我通过查看论文中的图 6 生成了一些值。

这是 WinBUGS 代码。它运行完美，但我不确定它是否正确。我很感激任何意见。谢谢！

如何在 WinBUGS 14 的语法中为差异项功能(DIF) 添加数据？

我是 WinBUGS 和 OpenBUGS 的新手。我刚刚运行了一个示例模型，想知道是否可以得到 WinBUGS/OpenBUGS 生成的预测。如果没有，是否有任何方便的方法来实现这一点（例如，借助 R 等其他应用程序）？

我正在使用rjagsR 库。该函数coda.samples产生一个mcmc.list，例如（从example(coda.samples)）：

但是，我想使用plot.bugs需要一个bugs对象作为输入的函数。

是否可以将一个对象从一个对象转换mcmc.list为一个bugs对象，这样plot.bugs(LINE.out)？

请注意，在 stats.SE上有一个类似的问题，一个多月以来一直没有答案。这个问题有一个赏金，于 2012 年 8 月 29 日结束。

更多提示：

我发现 R2WinBUGS 包有一个函数“as.bugs.array”函数 - 但不清楚该函数如何应用于 mcmc.list。