问题标签 [winbugs14]

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 回答
404 浏览

bayesian - “多变量分布必须有多个组件” - WinBUGS 中的错误

在编译过程中,我收到一条错误消息——“多元分布必须有多个组件”,我无法获得有关此错误的太多信息。任何人有任何想法如何解决这个问题,请分享。谢谢。

0 投票
1 回答
96 浏览

winbugs - winbugs 的脚本命令

如果我通过脚本在 Winbugs14 中运行我的模型,那么细化或绘制箱线图的命令是什么?我知道一些基本命令,例如:

但我在 Winbugs14 文档中找不到用于细化、绘制模型拟合或绘制箱线图的命令。有任何想法吗?

0 投票
1 回答
151 浏览

bayesian - WinBUGS 示例第 1 卷,染料示例返回错误

目前正在浏览示例第 1 卷,并遇到了染料示例的错误。

当我尝试从示例中加载 inits 时,它返回“此链包含未初始化的变量。我不确定它的哪一部分是不正确的被排除在外。

我在帮助选项卡下直接使用示例第 1 卷中的代码。相同的错误发生在所有三种选择之间的先验上。

我真的很感激任何关于这个问题的建议。提前致谢。

下面是我直接从染料示例中复制的代码。

数据

初始值1

初始化2

初始化3

0 投票
1 回答
1231 浏览

statistics - 如何在 WinBUGS 中使用多元回归获得多项概率

在 WinBUGS 中,我指定了一个具有多项似然函数的模型,并且我需要确保多项概率都在 0 和 1 之间并且总和为 1。

这是指定可能性的代码部分:

这里,数组 P[] 指定多项分布的概率。

这些概率将根据我的数据(矩阵 e[])使用对一系列固定和随机效应的多重线性回归进行估计。例如,这里是用于预测 P[] 的元素之一的多元线性回归:

在编译时,模型会产生错误:

如果我理解正确,这意味着向量 P[k,i,1:9] 的元素(上面的多项似然函数中的概率向量)可能是非常大(或小)的数字。实际上,它们都需要在 0 和 1 之间,并且总和为 1。

我是 WinBUGS 的新手,但从周围的阅读来看,似乎以某种方式使用 beta 回归而不是多重线性回归可能是前进的方向。然而,虽然这将允许每个元素介于 0 和 1 之间,但它似乎并没有触及问题的核心,即 P[k,i,1:9] 的所有元素都必须为正且总和为 1。

可能响应变量可以非常简单地转换为比例。我已经尝试过将每个元素除以 P[k,i,1:9] 的总和,但到目前为止还没有成功。

任何提示将不胜感激!

(我已经提供了模型有问题的部分;整个内容相当长。)

0 投票
0 回答
792 浏览

r - 解码陷阱 - WinBUGS 中的“索引超出范围”消息

我的任务是使用更新的数据重新运行其他人的模型(我进行了数据更新和格式化)。我了解模型的工作原理但没有编写它,而且它很长而且很详细。最重要的是,我专门使用 JAGS,这是我第一次涉足 WinBUGS,所以界面和错误消息对我来说是新的。该模型需要两天时间才能编译(它确实编译),但随后我收到一个陷阱错误,显示索引超出范围,然后是几十行我无法理解的行。我已经检查了整个模型、for 循环和数据结构,但我看不出哪里有任何索引问题。而且,令人沮丧的是,与 JAGS 不同,WinBUGS 似乎不会告诉您哪个索引超出范围。

我非常感谢解码 Trap 错误消息的任何帮助,以防其中有比我意识到的更多信息来指定索引不匹配的位置。我将在下面发布模型以及陷阱消息,但我看不出发布实际数据的意义,因为没有人会花 2 天时间等待模型编译以重现错误. 相反,我将只发布数据结构,以便每个数据对象的维度都很明显。

一些注意事项:我使用的是以前使用的确切模型(字面意思是复制和粘贴),所以我知道它应该可以工作。之前的数据文件是矩形格式,而我创建和导出的数据文件是 S 格式。我的理解是 WinBUGS 应该能够同时处理这两者,但我想知道这是否可能是问题的一部分。我直接从 WinBUGS 运行模型,手动加载模型和数据文件。我之前尝试使用 r2winbugs 从 R 中执行此操作,但遇到了内存错误和崩溃。

模型:

输入数据(我在 R 中进行了所有数据格式化,因此我将其编译为一个列表,然后使用 r2winbugs 以 S/R 格式导出它,所以它应该可以在 WinBUGS 中工作):

数据结构(我在实际的文本文件中手动确认了这些):

最后是陷阱错误:

非常感谢您的帮助!-乔什

0 投票
1 回答
103 浏览

r - 从 R 上传 winbugs 中的数据集

我想知道如何将数据集从 R 包上传到 winbugs。

特别是 R 中的“LearnBayes”包有太多的数据集。我想在 Winbugs 中使用其中之一。

谁能帮我这个?

0 投票
0 回答
31 浏览

loops - 无法使用 R2WinBUGS 访问 R 中的 WinBUGS 代码

我在 WinBUGS 中有模型规范,并且在 R 控制台中指定了所有参数和首字母缩写。在使用 bugs() 函数时说,bugs(data, parameters, initial,model.file="c:/ProgramFiles/WinBUGS14/sample1",n.chain,n.iters,bugs.directory="c:/ProgramFiles/ WinBUGS14/") 我收到一条警告消息,说“c:/ProgramFiles/WinBUGS14/sample1”不存在。有人可以帮忙破解一下吗?

0 投票
0 回答
33 浏览

winbugs - 使用 R2Winbugs 时在哪里保存模型文件?

我正在学习 R2WinBUGS。我几乎没有疑问。谁能帮我解决它?当我运行此示例时,我收到一条错误消息,提示“ratsmodel1.txt”不存在。我想在哪里保存我的模型文件代码?我应该如何使用 R2WinBUGS 从 R 调用它?

WinBUGS 中的模型文件代码

R编码

0 投票
1 回答
79 浏览

bayesian - “多元正态 Y[1,1] 的平均值必须与 Y[1,1] 具有相同数量的分量” - WinBUGS 中的错误

如何处理编译错误:“多元正态 Y[1,1] 的平均值必须与 Y[1,1] 具有相同数量的分量”?

这是我在 WinBugs 中的代码和数据:

数据:

0 投票
0 回答
61 浏览

initialization - WinBUGS中变量的初始化

我目前正在使用 WinBUGS 研究具有有序分类变量的 SEM(结构方程建模)。我对 WinBUGS 比较陌生,所以我按照Lee 和 Song 的书5.2.4 节中的代码进行操作。我的目标是将数据拟合到指定的模型(参见下面的 WinBUGS 代码)并获得参数的估计值。

下面提供了我的 WinBUGS 代码。数据集太大(Sample size 1986 with 19 variables)无法直接包含在下面的代码中,所以我在这里附上了两个数据集z.dataxi.data

所有 19 个变量都是在各种尺度上测量的有序分类变量,范围从四点量表到十点量表。我遵循书中所述的常见做法,使用 R 找到阈值并将它们作为“thd”包含在数据集中。

我希望上述代码应该在 WinBUGS 中成功运行。模型检查、数据加载和模型编译步骤没有问题。但是,当我继续在 WinBUGS 中加载两组初始值时,两条链都显示“此链包含未初始化的变量”消息,但我不知道哪些变量未初始化。当我尝试使用“gen inits”选项生成初始值时,在初始化模型时,会显示一个陷阱日志消息“未定义的实际结果”和一长串错误。我还用Lee 和 Song 的书5.2.4 节中的代码进行了试运行,但也出现了类似的错误(代码和数据集可以在同一个站点下载)。

我想到的问题是:

  1. 模型中的哪些变量未初始化?我以为我已经指定了所有必要的参数。
  2. 当我尝试使用“gen inits”选项生成初始值时,是什么导致了“未定义的实际结果”问题以及如何解决这个问题?

您能否提一些建议?非常感谢!