问题标签 [r2winbugs]

我目前正在对几项临床试验的生存数据进行荟萃分析。

为此，我使用相同方法从已发布的分析中获取代码。但是，当使用已发布分析中的数据运行此代码时，我无法复制他们的结果。事实上，结果无法收敛到任何合理的估计。

代码本身（不包括数据）应该是正确的，因为它直接来自作者。我认为问题必须与如何运行采样的初始值或参数有关，但是在播放了许多初始值之后，老化的长度，变薄等......我没有得到有意义的结果。

对于如何运行它（初始值等）以使其正常运行的任何建议，我将不胜感激。或者，如果代码中存在问题，或者数据的设置方式与代码不匹配，那么了解这将很有用。

作为旁注，我正在使用 R2WinBUGs 进行分析，尽管我单独使用 WinBUGs 时遇到了同样的问题。

该方法的一些额外背景：

其工作方式是通过使用随机效应估计多个研究的处理之间重新参数化的 Weibull 分布的形状和尺度参数的差异。

Weibull 分布被重新参数化，使得危险率的对数为 a+b*log(t)，其中 a 是尺度参数，b 是形状参数。由此，您可以计算给定数量的患者在一个时间间隔内给定数量的失败的似然函数。

不幸的是，这篇文章是公开的，但如果你可以访问这里是链接： http ://onlinelibrary.wiley.com/doi/10.1002/jrsm.25/abstract;jsessionid=2BA8F0D9BEF9A33F84975618D33F8DD9.f03t03?userIsAuthenticated=false&deniedAccessCustomisedMessage=

输入模型的变量的快速摘要：

NT：包括的单独治疗的数量。

N：主数据集中的行数。NS：研究数量

s：数据行对应的研究（编号为1:6）

r：此治疗/研究间隔失败的患者人数

n：此治疗/研究间隔开始时有风险的患者人数

t：这行数据对应的处理（编号为 1:3）

b：指示哪个处理是与其他处理进行比较的基线（每行设置为 1）。

bs：作为本研究控制组的治疗

bt：治疗是本研究的研究部门

WinBUGS 代码（包括数据）：

在尝试使用 调整工作WinBUGS模型并将其缓解为 RR2WinBUGS时，我收到了几条错误消息。我相信与 的规范有关inits，但一直无法解决问题。

第一个是错误消息正在使用inits1：

在阅读了 Uwe Liggers “列表包含列表解决方案”的建议修复后，我将 inits1 修改为inits2：

并收到错误：

我还尝试了 AndyC 在这篇文章“getting-winbugs-leuk-example-to-work-from-r-using-r2winbugs”中建议的修复。通过更改inits为：

并收到错误：

这是我对R2WinBUGS. 请注意，我inits在代码中包含了几个不起作用的代码：

这是工作WinBUGS代码，改编自WinBUGS布里斯托尔和莱斯特大学举办的课程：

问：谁能告诉我一个 WinBUGS 数据文件中的最大数据行数？

我如何发现有一个最大值。 在 R2WinBUGS 中构建和测试 WinBUGS 模型时，我构建了一个虚拟数据集，其中包含 6 个变量（18750 行）的所有可能数据组合。通过在 R 中运行代码，WinBUGS 构建模型并加载数据然后崩溃。直接在 WinBUGS 中重复此操作后，我在加载数据时收到消息“预期关键字 END”，尽管在我的数据文件末尾有“END”（包括之后所需的空行）。我花了一些时间才发现削减大量数据行可以解决问题。我在网上找到的 WinBUGS 手册或教程中没有任何信息（据我所知）。

更新：发现372行太长，216行还可以。我猜254是最大长度。

回答

偶然发现了这个。问题不在于数据长度，而是存在 X.Xe-4 而不是 X.XE-4（手册中明确说明，只是我看到错误时没有想到 'expected key word END'。限制数据长度只是限制了在我的数据集中获得小数字的可能性。

我想用命令 plot() 来解释 bug 模型中的收敛。输出示例如下图所示

我不确定我能不能很好地阅读这个输出，谢谢大家:)

当我使用库（R2WINBUGS）运行上述代码时，winbugs 面板显示以下消息：

预期的集合运算符 c compile(2) inits(1,C:/Users/molasey/AppData/Local/Temp/RtmpmAA9bi/inits1.txt) 命令#Bugs:inits 无法执行（灰显）

请帮帮我。

我真的希望这不是一个愚蠢的问题，但我正在尝试模拟一个非常简单的空间贝叶斯模型，并且在 WinBUGS 窗口上不断收到某种“无法复制”的消息。

我正在模拟数据：

这是模型：

任何线索为什么这件事不起作用？如果这真的很简单，我很抱歉。非常感谢！

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

模型如下。

在日志中我得到：

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

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

这是我的 R 代码数据：

谁能帮我解决错误？

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

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

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

因此，我一直在研究这段代码一段时间，为人群运行危险模型，而对于我的一个初始参数，它无法弄清楚如何修复它。这是给我带来问题的代码区域。我在 Notebook++ 中编码，在 R 中运行，它在 WinBUGS 中打开以运行它。Hz.Scale 是我正在努力解决的参数，并且出现了错误。