问题标签 [r2winbugs]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
811 浏览

r - R2WinBUGS - 如何再次将 WinBUGS 计算结果重新加载到错误对象

我使用 R 中的 R2WinBUGS 包在 WinBUGS(百万次迭代)中运行了一个 looong 计算:

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

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

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

0 投票
3 回答
6026 浏览

r - R2WinBUGS - 警告消息

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

代码

(请只考虑部分:### 5.4.使用WinBUGS分析)

我收到此错误消息:

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

谢谢。

基督教

PS:

这是R代码:

0 投票
1 回答
2707 浏览

r - R2WinBUGS - 使用模拟数据进行逻辑回归

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

谢谢。

基督教

PS:

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

0 投票
1 回答
1244 浏览

r - OpenBUGS:伯努利分布中的缺失值

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

我的模型是这样的:

R 代码如下所示:

0 投票
1 回答
1317 浏览

r - 从 R 调用 WinBUGS14

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

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

然后,我调用 WinBUGS:

我在 WinBUGS 中收到以下错误:

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

0 投票
0 回答
95 浏览

r - R2OpenBUGS:从中间结果建模最终治疗结果

有人在使用 R2OpenBUGS 吗?我应该使用 r2winbugs 吗?...

我正在尝试为具有(单个)中间(3 个月)结果的患者样本模拟最终(2 年)治疗结果(例如成功、死亡、违约或失败)。

R2OpenBUGS 在多项式节点上给了我一个奇怪的后验,其中两个结果是恒定的,另外两个结果是相等的,并且结果的总数大于队列大小。

我究竟做错了什么?提前谢谢了!代码输出如下。

以下是一些输出:

0 投票
3 回答
1605 浏览

r - 如何将 mcmc.list 转换为 bugs 对象?

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

0 投票
1 回答
168 浏览

r - 避免在 R 中舍入

我正在 R 中进行彻底的“R2Winbugs”包中的贝叶斯分析。在 R 输出中,估计的参数被四舍五入。我的问题是如何控制这个?

这是我的命令:

这是我输出的一部分

我想要的是将 beta1 视为 0.30000。我试过digits=5,但没有用。有谁知道该怎么做?

任何帮助将不胜感激

干杯

居纳尔

0 投票
3 回答
1153 浏览

r - 使用 WinBUGS/R2WinBUGS 的 Logit 回归未收敛到似然估计

我是 R 2.15.2 的 Win-7 用户

有人可以帮我为什么下面的模型不能很好地接近简单的 logit 模型估计吗?

已编辑

似然估计的模型代码:

使用R2WinBUGS的模型代码:(已编辑)

然后使用传统方法/不使用贝叶斯方法

请建议我是否需要使用不同的模型来更改先验或语法...

0 投票
1 回答
3992 浏览

r - R2WinBUGS 错误 - 陷阱 - 不兼容的副本

我试图从 R 调用 winBUGS 来估计逻辑回归。我正在使用以下代码:

WinBUGS 打开,但出现错误。“不兼容的副本”。我一直在寻找这个错误,发现很多人都遇到过这种情况。我已经尝试了建议的解决方案:model.bug 的文件路径更短,确保工作目录中没有更多的“数据”文件,使用数据变量名称的名称列表等。

但我仍然无法让它工作。我究竟做错了什么?