问题标签 [r2winbugs]

我是 OpenBUGS 的新手，在使用 logit() 函数拟合模型时遇到了一些问题。

阅读周围我发现一个可能的解决方案是明确指定 logit 函数而不使用 WinBUGS 自己的 logit 函数：

在更复杂的模型中，我们在使用 WinBUGS 自己的 logit 函数时经常遇到问题，例如实现收敛（实际上，即使使用相当简单的模型也可能出现问题。）。因此，通常最好通过 logit.p[i] <- log(p[i] / (1 – p[i])), p[i] <- exp(logit.p[i ]) / (1 + exp(logit.p[i])) 或 p[i] <- 1 / (1 + exp(- logit.p[i]))。

（更多信息在这里http://www.mbr-pwrc.usgs.gov/software/kerybook/AppendixA_list_of_WinBUGS_tricks.pdf第 14 点。）。

问题是我不明白该怎么做，假设我使用 WinBUGS 集成 logit 函数的原始似然函数是：

我如何明确地写出来？

非常感谢。

文森佐

我正在编写使用 NETWORK META-ANALYSIS 计算不同治疗效果大小的指南，如示例 6.a 中所做的那样。 这里

它在 winBugs 中运行良好，但我想使用 R2winugs 在 R 中进行分析，以便我可以自动化数据输入。

它根据日志很好地读取了模型，但是在读取数据时它被挂断了。

EG display(log) check(C:/Users/Temp User/.../Plaque_Psoriasis_Project/RE_Psoriasis.bug.txt) 模型是语法正确的数据(C:/Users/Temp User/.../Plaque_Psoriasis_Project/data.txt )

这是程序挂起的地方。陷阱屏幕显示

不兼容的副本

BugsCmds.TextError [000003A1H]

…………

我尝试将数据读取为：

数据 <- 列表（t=t，C=C，r=r，n=n，na=na，nc=nc，ns=ns，nt=nt，Cmax=Cmax，meanA=meanA，precA=precA）

并作为

数据<-列表（“t”，“C”，“r”，“n”，“na”，“nc”，“ns”，“nt”，“Cmax”，“meanA”，“precA”）

两者都不起作用。

我让 R2winbugs 在文档中做玩具学校的例子，这样它就可以工作了。

有什么想法吗？

我正在尝试使用大型数据集（590 万行）运行模型。我得到错误made use of undefined node ...

我已经尝试了很多样本​​数据（来自我更大的数据集），这些样本数据很好地通过了模型，直到样本数据集大小达到 ~8900 行，此时错误发生。

如有必要，我可以尝试生成一些示例数据来演示模型，但是不知道它是否仅与我正在使用的数据量有关？我对错误消息感到有些困惑，因为有问题的变量（节点）是在其他地方定义的，并且可以与示例数据一起正常工作。

该模型相当长，所以我一直在期待内存问题等。

亲切的问候

我正在尝试在 JAGS 中实现为 WinBUGS 编写的以下模型：

我认为在 JAGS 中使用它会很简单：

只需将 替换I()为T(). 这会产生coda.samples()错误：

如果我忽略监控weibmed然后工作weibmean但coda.samples()参数估计：

与我在使用 WinBUGS 时得到的那些没有可比性：

有什么想法或建议吗？

我们希望使用 R 和 WinBugs 运行给定的贝叶斯模型进行决策。我们正在使用 RStudio 和 R2WinBugs。模型由下式给出：

底层数据由一个 matlab 文件给出，我们通过 readMat() 方法读入该文件。

所有数据 <-readMat("file.mat")。

这实际上是有效的，因为我们可以在 RStudio 中看到这个文件给出的数据。该文件包含数组 x、v、p、y、m。这些数组存储在列表中

数据<-列表（所有数据$x，所有数据$v，所有数据$p，所有数据$y，所有数据$m）。

目前我们使用 NULL 作为 inits 值，这意味着 WinBugs 生成数据。但是我们不确定是否需要 WinBugs 来生成数据，或者我们是否需要自己指定：

样本 <- 错误（数据，inits=NULL，参数，model.file = "Path/SearchStop.txt"，bugs.directory = bugsdir，debug = TRUE）

这会带来以下错误：

在 R2WinBugs 中没有关于哪个变量未定义的确切信息，但是当我们在 Matlab 中尝试它时，错误消息说 z1 未定义。

全部！

我正在使用 winbugs 进行简单的线性回归。但是，系统总是给出错误消息，应为逗号。

这是我的模型声明：

它有什么问题？谢谢你。

我的数据文件如下所示：

假设这个文件被保存为“data.txt”。

如果我在“R2OpenBUGS”中工作，它可以让我毫无问题地将数据作为文件传递：

但是在 JAGS 中，如果我将数据作为“data.txt”传递，它会说：“数据必须是列表或环境”。这里有什么问题？另外，如果没有办法解决它，有没有办法可以将数据作为 R 中的列表读取？

我的模型是：

我在 R 的“R2OpenBUGS”中工作。我有一些 mcmc 链：

我的问题是当codaPkg = TRUE时如何打印 DIC 。如果 codaPkg = FaLSE，那么如果我只执行 print(mcmc)，它会在最后打印 DIC 值。

但是当“codaPkf=TRUE”时，它只会打印如下内容：

那么，我现在如何打印 DIC 值呢？

第一次在这里发帖，所以我希望我能够勾选所有要求的框。我正在使用 NICE 技术支持文档通过 R 和 R2WinBUGS 提供的代码运行网络元分析。使用均值差的模型效果很好，但是当我转换为对冲时，GI 突然得到一个未定义的实际误差，这似乎是非正定协方差矩阵的结果。我确定我遗漏了一些明显的东西，但如果有任何帮助，我将不胜感激。请看下面的代码。

我在弄清楚如何使用 R2WinBUGS 运行一些模拟研究时遇到问题。目的是模拟 n 个数据集（以 1000 个为目标，但从 10 个开始），并将它们全部作为矩阵放入 R2WinBUGS 代码中，以便当它移植到 WinBUGS 时，它将运行生成 n 个数据集的估计值。这是我目前拥有的：

该模型：

我使用的错误代码是：

其中 Ndat 是数据集的数量，py 是一个 13 xn 矩阵，初始值是：

有什么帮助吗？