问题标签 [winbugs]

我希望能够指定 t-1 中的“alpha”有多少进入 t 的先验。我试过 alpha[t,1:k] <- theta[t-1,1:k]*alpha0*n[t-1]; alpha0<-.5，但它没有用。

我尝试在 sum(delta[j,1:5])=0 的条件下为每个“j”估计 delta[j,k]。但是当我编译代码时，软件输出是“节点 W[1] 的多个定义”。有人可以帮助我吗？

谢谢

我无法使用 PyMC 从提供的 BUGS 代码中复制结果。BUGS 模型是 Andersen-Gill 乘法强度 Cox PH 模型。

我使用以下初始值

我使用单链（目前）和 5000 次迭代进行老化。我运行了 10000 次额外迭代的估计，并收到了与论文中报告的相同的点估计。这些也接近于早期的常客估计。

我尝试使用以下代码在 PyMC 2.3.2 中复制它。完整的复制代码可在此处获得

但是，我并没有接近相同的点估计。我得到类似的东西

最令人担忧的是，迹象不同。我想这可能只是一个收敛问题，所以我一夜之间运行了 50,000 次迭代，没有太大变化。也许我的 PyMC 模型存在一些错误或差异，尤其是 dL0 规范？

我尝试过不同的起始值。我试过让模型运行多次迭代。我将先验集中在 BUGS 的点估计上。

我正在尝试使用 BUGS 编写狄利克雷多项式模型。基本上我有 18 个地区和每个地区 3 个类别。例如，区域 1：0.50 属于 Low，0.30 属于 Middle，0.20 属于 High。该列表当然以不同的比例继续到第 18 区。我得到的唯一代码是这个

例如，我首先将其缩短为 3 个区域。单击“gen inits”后，它会显示“Dirichlet36”。请帮我编写代码。

我在运行我的 OpenBUGS 程序时遇到问题。该模型在语法上是正确的，但是一旦加载了数据，我就无法编译，因为出现错误“数组索引大于 beta的数组上限”。我一直在检查所有内容，但似乎找不到错误。

我正在尝试学习 WINBUGS，并尝试构建一个小型模型，根据教科书中的示例进行调整（代码如下），假设隐藏的受感染携带者群体具有增长率（“R0”）和去除率（筛选和处理）随时间变化。但是，我倾向于收到一系列错误消息（“扫描的令牌无效或意外”、“无法执行初始化”等）。因此，有更多WINBUGS经验的人能否如此善意地关注我是否在我对WINBUGS的理解上犯了愚蠢的错误？我特别不确定是否可以在 WINBUGS 中完成人口的顺序更新（N.est[t+1] <- N.est[t] + newcases - obs）？提前谢谢了

我正在通过 R 运行一个 BUGS 模型，但我遇到了一个问题，BUGS 说我的一个数据集是一个“未定义的变量”。它遇到问题的数据集有很多 NA，但为什么会导致问题？当 NA 被 0 替换时它工作得很好，所以我知道使用 NA 的 BUGS 存在问题。使用 0 而不是 NA 不是可接受的解决方案，因为我需要记录转换结果数量。有没有人遇到类似的事情并找到了解决方案？

这是我正在努力解决的模型部分：

该模型还有很多内容，但其余部分运行良好。我已经剪掉了模型的其余部分，只包括了给我带来问题的部分。我已经运行了模型的这一部分，我得到了相同的“c.eff.unr 是一个未定义的变量”错误

我是贝叶斯统计和 MCMC 的绝对新手，所以我正在研究 John Kruschke 的“做贝叶斯数据分析：R 和 BUGS 教程”。为了测试我的理解，我正在尝试将他的示例从 BUGS 转换为 PyMC。

英寸。在图 8 中，他观察了两个（可能有偏差的）硬币中的每一个的一系列翻转，并试图估计它们的偏差的差异。低于每个硬币的偏差是theta，观察到的翻转是y：

当然，这两个硬币是独立的。如果我分别模拟它们，然后查看 thetas 的差异，我会得到与本书相同的答案。（它也与网格上的分析解决方案和集成一致。）

另一方面，如果我尝试同时模拟两个硬币，作为同一模型的一部分，我会得到一个完全不同的答案，与任何其他方法都不一致。

所以我的第一个问题是，为什么？根据我对该理论的了解很少，我的猜测是 MCMC 很难充分探索二参数模型。（但是，BUGS 似乎处理得很好。）

真正奇怪的是，我一直在 iPython 笔记本中做这一切，而且 PyMC 中似乎存在一个错误。如果我运行独立硬币模型，重新启动我的内核（Kernel | Restart or File | Close and Halt），然后运行联合硬币模型，联合硬币将产生与独立硬币相似（但不相同）的答案（平均dtheta ~ 0.23)。如果我以相反的顺序运行模型（在两者之间重新启动内核），它们都会从联合硬币模型中产生（不正确的）答案，平均 dtheta ~ 0.33。如果我在两者之间完全关闭我的 iPython 笔记本服务器，我只能让这两个模型产生不同的答案。由于这也会从内存中卸载所有共享库，我希望这意味着 PyMC 的 Fortran/C 部分正在将这些模型中的某些内容缓存在一个内存位置 s 在 Python 解释器实例之间共享。版本是 Numpy 1.8.2、PyMC 2.3.3、Python 2.7.8、iPython 2.1.0、Anaconda 2.0.0。

对于理解这里发生的事情，我真的很感激任何帮助。我意识到这些都是愚蠢的、微不足道的模型，但 PyMC 的奇怪行为目前并不能激发信心！

我到处寻找，但仍然找不到 R 中可以帮助我从 mcmc 对象编写 coda 文件的函数。我正在使用 R2OpenBUGS 包在 R 中运行 OpenBUGS：

output <- bugs(...)

然后经过一些操作，我留下了一个名为chain_matclass的对象mcmc。这个对象相当长（超过 200 万行），所以我想创建一个 coda 文件来保存绘图。我知道您可以直接从bugs函数创建 coda 文件，但我需要先进行一些操作，然后才能创建文件。

任何想法将不胜感激。

我对 Winbugs 比较陌生。当我尝试使用以下模型加载数据时出现“预期变量名称”错误。我无法弄清楚数据中的哪个变量名称不是模型中的变量名称。请帮忙。