问题标签 [winbugs]

我正在尝试将 ANOVA 模型拟合到 rjags 中。模型是这样的

我需要适应

这是此模型的零约束之和。但是，rjags 给了我信息

如果我删除约束部分，模型编译得很好但不会收敛。在 BUGS 中，模型被接受。

如何在 rjags 中实现这些约束？

我一直在尝试使用贝叶斯方法解决问题，但我不知道如何编写模型。

我试图找到值 mu 使得 x = mu + Error, Error ~ Gamma(a, b) （或其他分布）。

基本上我希望能够说类似的话：

但是，这不起作用，因为 x[i] 需要有一个分布，我想不出如何做到这一点（我试过让 x[i] ~ dgamma(0.001, 0.001)，但这没有't a. 将 x 限制为 > mu，并且 b. 不能帮助我估计 mu。

如果您能提供帮助，将不胜感激。

谢谢！

由于我是 winbugs 新手，我需要在 R 中进行贝叶斯分析。为此，我需要将下面的 winbugs 代码转换为 R 代码：

这是我使用的 R 代码

这是上面使用的“reg.txt”文件

最后这是winbugs发起的错误：

显然，我没有定义变量 n。有谁知道出了什么问题以及如何解决。任何帮助将不胜感激。

提前谢谢了

居纳尔

我正在 winbugs 中进行贝叶斯分析。这是我的模型：

其中 b[i] 是第 i 个随机效应。我想知道如何指定 tau、alpha 和 betas 的先前分布。考虑哪些点？任何帮助将不胜感激。

干杯

我正在使用 R 中的 Winbugs 进行贝叶斯分析。我需要将两个 Winbugs 脚本合并为一个：但是，我收到一条错误消息 ( Variable x2 is not defined in model or in data set)。这是winbugs代码：

其中x2是 1002*2 矩阵并且y是向量

这是定义数据的 R 代码：

有谁知道出了什么问题？

我一直在尝试运行 R 代码来使用 OpenBugs。直到昨天它运行良好，我不知道我做了什么，它开始显示一些错误。然后我决定看看一些简单的例子是否有效。所以我尝试了R2WinBUGSvignettes中的学校示例。

文件schools.txt中的模型是

运行程序其余部分的代码如下：

不幸的是，我遇到了和以前一样的错误。

我也尝试过指定program=c('OpenBugs') 但是鉴于 OpenBugs 在路径中，这给了我另一个奇怪的错误

编辑

traceback()给...

编辑-2

我应该如何删除sd我定义的那个坐在里面.GlobalEnv？

EDIT-3

问题解决了。sd我从工作区中删除了对象，一切正常。

我通过从 R 运行 Winbugs 进行了贝叶斯分析，并得出了拟合值及其贝叶斯区间。这是相关的 Winbugs 输出，其中 mu[i] 是第 i 个拟合值。

我想要做的是找到这 5 个拟合值的平均值的贝叶斯区间。知道怎么做吗？

我正在 R 的 winbugs 中进行 logit 回归分析。我必须强制该模型的所有系数为正。因此，我对所有系数都使用了统一的先验，但 winbugs 对此并不满意：它会生成一个愚蠢的错误窗口。当我将dnorm(0.0,1.0E-4))所有系数用作先验时，问题就解决了。在下面给出的这个模型中，可以做些什么来获得正贝塔？

我进行了贝叶斯分析并在winbugs中获得了一些图。有谁如何将历史图放在Latex的WinBUGS中？我无法将绘图保存为 pdf。

我需要计算每个状态的平均值。在 R 中，我通过使用

tapply(y,state,mean).

这是一个例子：

其中 y 表示拟合值的向量，state 表示作为因子的因变量。

我想知道如何在 Winbugs 中做到这一点。有人对此有任何想法吗？