问题标签 [winbugs14]

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

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

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

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

我正在尝试使用此示例使用 R2WinBUGS：

代码

（请只考虑部分：### 5.4.使用WinBUGS分析）

我收到此错误消息：

我不确定这对于正确的功能是否至关重要（其他一切看起来都不错）。有没有办法摆脱这个？

谢谢。

基督教

PS：

这是R代码：

我只是想知道是否有人有一些使用包 R2WinBUGS 运行逻辑回归的 R 代码 - 理想情况下使用模拟数据来生成“真相”和两个连续的协变量。

谢谢。

基督教

PS：

生成人工数据（一维案例）并通过 r2winbugs 运行 winbugs 的潜在代码（它还不能工作）。

我正在尝试通过 R（R2OpenBUGS）将观察“时间”建模为 OpenBUGS 的随机变量。如果所有观察时间都可用（没有 NA）一切正常，但如果我将其中一个时间设置为 NA，则什么也不会发生。我用 WinBUGS 测试了相同的代码，我得到陷阱错误'NIL dereference (read)'。所以我的问题是我的代码中是否真的有问题，或者我的模型对于错误来说太奇怪了？

我的模型是这样的：

R 代码如下所示：

正如这篇文章的标题所说，当我尝试使用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) 添加数据？

我正在使用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。

从stan 网页启动标准示例时，如下所示：

（这是从这里获得的）

然而，这只为我提供了参数后验的分位数。所以我的问题是：如何获得其他百分位数？我想它应该类似于错误（？）

备注：我试图介绍标签stan，但是我的名声太小了；）抱歉