问题标签 [jags]

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

r - Attach() -> 找不到对象错误 (R2Jags)

我正在使用 R2Jags 包。

我首先获得模型如下:

但是当我使用以下方法获得 X1 的后valeus 时:

我得到:

错误:找不到对象“X1”

我的问题是:我该怎么做才能解决这个问题?

提前致谢。

0 投票
0 回答
416 浏览

r - 为什么我的 beta-binomial 模型使用 jags 的参数估计不同于最大似然估计

我有一个像这样的 beta-binomial 模型 在此处输入图像描述

其中 $B$ 是 beta 函数。

我想估计参数 $\theta_1,\theta_2,\ldots,\theta_5$。

我使用了最大似然法:

然后,我使用 MCMC 如下:

MCMC 估计与最大似然估计一致。

最后,我使用 jags 估计参数如下:

这种方式的估计很奇怪,不像以前的值。我想知道 jags 功能有什么问题!?感谢您提前提出任何意见或建议。

https://ehc.ac/p/mcmc-jags/discussion/610037/thread/dc35eac1/#4823

0 投票
1 回答
1412 浏览

r - 响应为比例时的逻辑回归(使用 JAGS)

我正在尝试在 JAGS 中拟合逻辑回归模型,但我的数据形式为 (#success y, #尝试 n),而不是二进制变量。在 R 中,可以通过使用带有“权重”参数的 glm(y/n ~) 将模型拟合到诸如此类的数据,但我不确定如何在 JAGS 中拟合它。

这是一个简单的例子,我希望能解决我想要问的问题。请注意,我使用的是 rjags 包。谢谢你的帮助!

0 投票
1 回答
1271 浏览

survival-analysis - JAGS 中的逐步右删失生存分析

这是对 SE 上一篇较早帖子的一种跟进:https ://stats.stackexchange.com/questions/70858/right-censored-survival-fit-with-jags

但是在这里,我希望看到一个完整的 R 脚本(从头到尾)对 JAGS 中的右删失数据进行生存分析。我发现的所有网站都需要非常熟练的 JAGS,所以我很难理解如何从一行代码转换到另一行代码。我知道这有很多问题要问...

无论如何,这里有一些示例生存数据。组是 t1、t2、t3。NA 指的是右删失数据(审查截断 = 3)。

我知道有很多问题要问,但是我花了几天时间试图拼凑一些东西,但我一直迷路/困惑。我知道现在有一些包可以运行这种分析,但我真的很想学习如何从头开始自己构建它!谢谢各位读者!

0 投票
1 回答
1748 浏览

r - 在使用 jags 创建的模型上使用 update() 的原因

我正在使用 rjags 在 R 中进行分析(基于这篇博文:http: //www.sumsar.net/blog/2013/08/bayesian-estimation-of-correlation/)我有一个问题. 使用 update() 更新具有 500 个样本的模型与仅将原始模型创建中的 n.adapt 参数设置为 1000 而不是 500 之间有什么区别吗?

换句话说:当我执行以下操作时,run1 和 run2 之间是否存在重要区别:

0 投票
1 回答
735 浏览

r - JAGS 运行时错误:无法将节点插入 X[]。尺寸不匹配

我正在尝试向数据增强捕获重新捕获模型添加一些代码,并且遇到了一些我以前从未遇到过的错误。简而言之,我想估计一系列生存阶段,每个阶段都持续一个以上的时间间隔。我希望该模型估计每个幸存阶段的长度,并使用它来改进捕获重新捕获模型。我尝试了几种不同的方法,但都失败了,现在我尝试使用生存阶段的切换状态数组来实现这一点:

因此,这将创建一个 [5,t,5] 数组,其中幸存状态只能切换到后续状态而不能向后切换(例如,1 到 2、4 到 5,但不能从 4 到 3)。现在我创建一个定义生存状态的向量:

我们总是从状态 1 开始,然后在每个时间步 't' 进行分类抽签,以保持当前状态或在给定数组内概率的情况下进入下一个状态。我想要最多 5 个状态(假设模型通过估计从状态 3 移动到 4 并接近 0 的概率,或者使后续状态的生存值相同或相似,如果它们属于到现实中相同的生存值)。所以我创建了 5 个分层生存概率:

现在下一步是错误开始的地方。现在我已经为我的 PhiState 向量分配了值 1-5,它应该看起来像这样:

或者可能

我现在想将 mean.phi[] 分配给我的实际 phi[] 项,它会输入模型:

但是,当我尝试运行它时,我收到以下错误:

值得注意的是,当我使用以下 phi[] 确定时,模型工作得很好:

或者

或者

我确实读过这篇文章:https ://sourceforge.net/p/mcmc-jags/discussion/610037/thread/36c48f25/

但我看不到在这种情况下我在哪里重新定义变量......任何解决这个问题的帮助或关于更好方法的建议都将受到欢迎!

非常感谢,乔希

0 投票
0 回答
181 浏览

r - 使用 jags.parallel() 时出现 Jags 错误:“所有连接都在使用中”

我正在使用该R2Jags软件包,并且我想使用该jags.parallel函数来加快计算速度。我使用以下代码:

我得到的错误是Error in file(con, "w") : all connections are in use.

有谁知道这可能是什么原因?提前致谢。

0 投票
1 回答
496 浏览

bayesian - JAGS 贝叶斯状态空间建模

我正在尝试使用状态空间模型来估计人口统计数据(繁殖力、存活率、人口增长、人口规模)。我们有 4 个不同的年龄状态。

这是错误代码:

不幸的是,当我运行代码时出现以下错误。

有人对我们为什么会收到此错误有任何建议吗?不知道为什么 J0 的分布是错误的。

0 投票
0 回答
773 浏览

r - JAGS - 无法找到合适的采样器

我正在尝试在 JAGS 中开发分层 Dirichlet-多项式过程隐藏马尔可夫模型,以根据民意调查结果估计多方、主要投票意图。我还使用主要投票估计来计算澳大利亚优先投票制度下的两党优先投票份额。

dmulti() 多项分布失败并显示运行时错误消息:无法找到合适的采样器。我有一个使用一系列二项式分布和总和到 N 约束的解决方法。从理论上讲,这应该会产生相同的结果,但会在模型中造成空间和时间效率低下。

我的问题是我是否可以在下面模型的隐藏的时间部分中做一些事情来使多项分布工作。

模型(和周围的 R 代码)如下:

六个月内模型的输入数据如下。

下面是输出的比较(与我拥有的其他模型相比)。下一张图表中的红线是从上面生成的。

在此处输入图像描述

0 投票
1 回答
503 浏览

r - 尝试将 dmulti 与 rjags 一起使用时尺寸不匹配

我正在尝试使用 rjags 拟合以下模型,但出现尺寸不匹配错误。该模型在 WinBUGS 中运行良好,我不知道如何更改代码。提前感谢您的帮助。

数据:

链的初始值

模型

错误

jags.model 中的错误(文件 =“model.txt”,数据 = 数据集,inits = inits,:运行时错误:无法将节点插入 y1[1:2,1:2]。尺寸不匹配