问题标签 [winbugs]

我正在尝试编写一个 Winbugs/Jags 模型来建模多粒度主题模型（正是这篇论文-> http://www.ryanmcd.com/papers/mg_lda.pdf）

在这里，我想根据特定值选择不同的分布。例如：我想做类似的事情

这可以在 Winbugs/JAGS 中完成吗？

我有一个输入需求模型，数据不平衡。因变量是yijt，其中i表示生产函数的投入(i=1,2,3)，j表示公司(j=1,..,21)，t表示时间(ti=1,..,Ti)。所以，因为面板是不平衡的，有三个方程，在 winbugs 中使用它的正确格式是什么？

我尝试以 STATA 长格式 3*107 的形式使用它，其中索引有点不同。行 ( i) 表示输入，列 ( j) 表示面板i*t。所以，没有明确的时间索引。我对列使用嵌套索引，将数据形成为面板。一小段代码：

其中a[i]是一个特定于输入的自由参数，u是一个期望的特定于公司的误差分量，应该对公司进行估计（21）。为此，我使用嵌套索引，其中模型中 u 的表达式超过列u[firm[j]]，而先验 foru超过公司u[k]。

这是关于嵌套索引的问题。如果模型包含许多协变量，这些协变量也是时间相关的，我是否应该对所有这些协变量使用嵌套索引，否则 Winbugs 不会将其作为面板？

我尝试对所有多维变量使用多维数组。为了平衡很多 NA-s 被包括在内。例如，因变量 yijt 是一个 3 维数组y[ , , ]，（使用 list 命令），但由于某种原因，Winbugs 无法识别 NA-s？要点是后验分布不收敛，数据格式化可能是原因之一

我尝试使用WinBUGSfrom Rvia BRugsand R2WinBUGS，代码如下：

结果没有出来，我发现了部分log.txt的WinBUGS：

很明显erroris multiple definitions of node y[1]，但它是什么意思？我不认为y[1]有多个定义，因为我使用y[i]但y不在loop.

代码如下：

我不知道为什么，当我替换为时程序会正确theta[j] <- b*x[j-1]，theta[j] <- 1*x[j-1]但我已经定义了b ~ dunif(-1, 1). 确实，我需要theta[j] <- a - b*x[j-1]在最终模型中进行设置，当我尝试添加a并b进入它时结果是错误的。有谁发现问题出在哪里？

我用Winbugs了from RviaR2WinBUGS和BRugspackage，我写了模型并尝试运行它，但是结果并没有出来，我看了WinBUGS日志，似乎没有任何问题。这是模型：

模型检查通过，日志在WinBUGS这里：

WinBUGS没有继续，并且没有显示错误消息，所以我找不到任何错误。

我需要从 Excel 将一些数据加载到 WinBugs 中，由于某种原因，Winbugs 无法识别它。我尝试了很多事情（所有这些都在下面列出），但这些都不起作用。

1. 我试图从 Excel 复制并将“选择性粘贴”作为“纯文本”复制到 WinBugs 中。我[]在每个变量之后添加，甚至'END'在最后添加。错误消息如下：

   LoadNumeric Data抱歉，模块中的程序出了点问题BugRectData

2. 我复制并粘贴到记事本中，[]在每个变量之后添加，甚至'END'在最后添加，然后删除行之间的所有空格，每个数据之间只有一个空格。当我将它拖到 Winbugs 中时，它分散在行之间的双倍空间中，数据之间的间隙很大，并且出现错误消息（尝试加载数据时）：

   LoadNumeric Data抱歉，模块中的程序出了点问题BugRectData

3. 我做的最后一件事是将 .txt 文件打开到 Winbugs 中，它变得一团糟，当我尝试加载数据时，它反而突出显示了一些NA2.0应该是的数据，NA 2.0并说它“期望一个数字或 NA 或结尾”。在纠正了这些错误之后，数据似乎很好，直到它说的最后，它'END'说的是同样的事情，只是突出显示E了END.

如果有人能提供任何见解，我将不胜感激。

我需要帮助来查找我的 WINBUGS 代码 (v. 1.4.3) 中的错误。

在“模型规范”步骤中，模型看起来语法正确。但是，在我尝试加载数据时，我收到了这个错误：

数组索引大于 phi3 的数组上限

有人可以帮我吗？我的模型提供如下：

我pnorm在 R 中使用标准正态分布。我想知道我可以在 Winbugs 中使用什么代码来执行此操作，因为我将使用贝叶斯概率。

例如，我将 p = Phi(Xbeta) 计算为 pnorm(Xbeta)。

我的问题是如何在 Winbugs 中做到这一点。

我找不到我的 BUGS 模型的问题。模型语法正确，数据加载正常，但在加载我的初始化时出现以下错误消息：

初始值的数量不等于向量的长度

我的代码在下面提供以供您欣赏。

一切顺利！

• 卡米拉

  /li>

如果等号括号中的元素之一是 NA，请有人告诉我 WinBUGS 中的等号函数返回什么？

具体来说，我想做转换：

并且想要保留 NA，因此 oldvar[i] 中带有 NA 的记录在 newvar[i] 中也有 NA，但我担心这种转换会将 NA 变成 0。

非常感谢任何和所有建议。

董里