- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我很难找出这个简单的 Rjags 模型中错误的位置:
dt=
[1] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 2.304157e-06 1.014568e-05
[8] 1.798720e-05 2.582872e-05 3.367025e-05 4.151177e-05 4.935329e-05 5.719481e-05 6.503634e-05
[15] 7.287786e-05 8.071938e-05 8.856090e-05 9.640242e-05 1.042439e-04 1.120855e-04 1.200225e-04
[22] 1.274917e-04 1.332419e-04 1.426773e-04 1.682728e-04 2.018543e-04 2.302676e-04 2.810248e-04
[29] 3.408895e-04 3.942806e-04 4.818033e-04 5.682188e-04 6.298504e-04 7.320266e-04 9.604787e-04
[36] 1.315336e-03 1.695445e-03 1.945434e-03 2.047065e-03 2.192386e-03 2.404184e-03 2.629989e-03
[43] 2.931693e-03 3.414607e-03 4.085636e-03 5.009995e-03 5.897174e-03 6.532098e-03 7.205807e-03
[50] 8.125645e-03 9.097041e-03 9.857559e-03 1.021277e-02 1.026859e-02 1.053671e-02 1.125948e-02
[57] 1.207536e-02 1.250217e-02 1.297049e-02 1.372405e-02 1.463063e-02 1.552144e-02 1.604347e-02
[64] 1.628257e-02 1.701820e-02 1.843134e-02 2.045894e-02 2.294199e-02 2.517307e-02 2.666305e-02
[71] 2.762321e-02 2.890740e-02 3.095278e-02 3.290344e-02 3.429784e-02 3.502680e-02 3.525143e-02
[78] 3.534536e-02 3.508330e-02 3.431497e-02 3.302332e-02 3.136818e-02 2.979437e-02 2.850674e-02
[85] 2.713331e-02 2.539821e-02 2.301579e-02 2.002319e-02 1.702893e-02 1.447840e-02 1.217091e-02
[92] 1.006765e-02 8.119336e-03 6.326128e-03 4.835775e-03 3.611105e-03 2.514688e-03 1.618690e-03
[99] 1.006864e-03 6.681004e-04 4.693963e-04 3.256632e-04
data = list('dt')
library(rjags);library(foreign);library(coda)
cat('model{
for(k in 1:102){ dt[k] ~ dgamma(a, 0.5)}
a=ax+ay
X ~ dgamma(ax, 0.5)
Y ~ dgamma(ay, 0.5)
ax ~ dgamma(3, 1/10)
ay ~ dgamma(3, 1/10)
}',
file='model1.bug')
inits1 = list('ax'=15, 'ay'=15)
jags_mod = jags.model('model1.bug', data=data, inits=inits1, n.chains=1, n.adapt=2000)
update(jags_mod, n.iter=2000)
posts=coda.samples(model=jags_mod, variable.names=c('ax','ay','a', 'X', 'Y'), n.iter=1000000, thin=2000)
主要思想应该是变量 dt
是一个 Gamma,即具有相同速率参数的两个 Gamma 随机变量的总和。我已经通过在两个比例参数上也分配 Gamma 先验来确保参数不能有负值,所以我不明白为什么 Rjags 在这里提示。我想推断两个参数 ax
、ay
(然后是 a
)。
感谢您提供的所有帮助!最好的,伊曼纽尔
最佳答案
gamma 分布严格为正,因此不能为零(或负)。因此,您数据的前 5 个观察结果与您的模型不一致,因此出现错误。
可能的解决方案是将其建模为删失,质量为零的混合模型,或者可能只是向数据添加一个常量,但在不知道数据代表什么的情况下无法说出哪个是最好的。
关于jags - Rjags 错误 : node inconsistent with parents,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/39901137/
我正在尝试将 ANOVA 模型放入 rjags。模型是这样的 for (r in 1:nE){ for ( j in 1:nP){ for ( i in 1:nA){ logi
与很多其他人一样,我在加载和安装 rjags 时遇到问题。 我得到错误: library(rjags) Error : .onLoad failed in loadNamespace() for 'r
我正在使用 RJAGS 修改现有模型。我想并行运行链,偶尔检查 Gelman-Rubin 收敛诊断,看看我是否需要继续运行。问题是,如果我需要根据诊断值恢复运行,重新编译的链会从第一个初始化的先前值而
我第一次尝试用 rjags 拟合分层模型;它在纸上看起来很简单,但我得到一个“ 解析模型文件时出错:第 5 行靠近“[”“的语法错误,我完全无法解释。 你能帮我告诉我我做错了什么吗? data = l
S.O.:Linux Ubuntu 14.04 LTS R: R version 3.2.1 (2015-06-18) -- 《世界著名宇航员》 版权所有 (C) 2015 统计计算 R 基金会 平台
我正在尝试使用 rjags 拟合多项逻辑回归模型因为结果是一个具有 3 个水平的分类(名义)变量( 结果 ),解释变量是 年龄 (连续)和 群 (分类为 3 个级别)。这样做时,我想获得 的后验均值和
在建立模型并使用 Gibbs Sampling 进行训练后,我得到了所有隐藏值预测的结果: jags <- jags.model('example.bug', data
R2WinBUGS 包有一个名为write.model() 的函数。 R 包 rjags 没有我所知道的这样的功能。 write.model 创建一个临时文本文件,WinBUGS 可以将其作为模型读取
我使用这两个包来进行贝叶斯分析,但有一些我不明白的差异: 首先包rjags允许适应阶段,使用 jags.model功能,而包r2jags没有这个阶段,并带有功能jags (或 jags.paralle
我正在尝试在 Ubuntu 上安装 Rjags,我首先安装了这个 sudo apt-get install JAGS 那我试过了 R install.packages("rjags") 我有以下错误
我很难找出这个简单的 Rjags 模型中错误的位置: dt= [1] 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+0
在 R 中安装 Jags 4.0 和 rjags 包时,Jags4.0 和“rjags”包之间的匹配库似乎有问题。 我首先将 Jags 4.0 安装到“C:\Program Files\JAGS\JA
当通过 rjags 在 JAGS 中运行一个相当复杂的模型时,我得到一个错误 Error: Error in node ttt3[126,509] Current value is inconsist
我对贝叶斯和 JAGS 都是新手,所以请原谅我的无知。 我收到了一个使用 JAGS 代码编写的 R 脚本(由同事发送)。 这段代码的作者定义了一组尾声样本如下: codaSamples = coda.
我通常从带有多个链的 rjags 调用 JAGS 以用于诊断目的(例如 4 个链)。之后,我经常想对后验参数估计进行一些后处理(例如,使用预测值、计算附加统计数据)。然而,此时将链存储在列表中是很麻烦
我正在研究使用组级预测变量拟合多级逻辑回归模型。我通过 R 使用 JAGS。当我用 runjags 和 R2Jags 包拟合模型时,我得到了不同的行为。 我已尝试编写一个可重现的示例来说明该问题。下面
我在加载和运行 rjags 时遇到问题。我可以安装 rjags,但它无法运行,因为它找不到 JAGS。下面列出了我的版本,并且都更新到了最新版本。我也在用Thinkpad。JAGS 版本:4.3.0R
在 Ubuntu 14.04 中从 RStudio 安装 rjags,我收到错误消息, configure: error: "JAGS module directory /usr/lib/JAGS/m
我想编写一个脚本来在新的 ubuntu 安装上安装 JAGS 和 rjags,并且它将独立于这些包的当前可用版本工作。我想知道如何在避免版本之间冲突的同时做到这一点。 我有以下 R 脚本,initia
我最近开始使用 JAGS 并在 R 中调用它。我终于用代码将 jags 链接到 R install.packages("rjags") library(rjags) 得到输出 Linked to JA
我是一名优秀的程序员,十分优秀!