我刚从R2jags转自R2OpenBUGS，注意到了一些我不明白的事情。在用jags()进行仿真和用as.mcmc()转换输出后，第一个样本总是有很大的偏差，并且通常离收敛参数估计很远。使用bugs()运行相同的数据时，不会出现此示例。这几乎就像第一个样本是真正的第一个样本，从烧伤阶段。可重复的代码，包括一个糟糕的初始估计，以在jags()的第一个示例中显示糟糕<
浏览 0提问于2014-10-31得票数 1
回答已采纳
1回答
我正在通过R包runjags运行JAGS模型。我刚刚从Jags3.4更新到Jags4.0.0，并注意到了一些与更新相关的意外行为。它似乎不影响结果(但这是非常令人费解的)。然而，我在玩这个的时候注意到了几次，对于一些变量来说，后验和前驱几乎是一样的(表示没有更新)。我现在似乎无法重新创建更新失败，但下面是一个可重复的代码示例，演示了奇怪的警告信息。run.
浏览 5提问于2015-10-09得票数 1
回答已采纳
1回答
在JAGS/BUGS中指定模型时，"is as“符号~非常有用。当使用要求我指定可能性的MCMC方法时，如何在R中做到这一点？
比方说，我想估计三个多元正态分布的参数。在JAGS中，我会通过指定pars[1:n] ~ dmnorm(mu[1:3], sigma[1:3, 1:3])来做到这一点。如果所有内容都被正确指定，JAGS将继续在给定分布下估计这些参数。在R中，有类似的函数，如mvtnorm
浏览 2提问于2020-04-10得票数 1
回答已采纳
我有一个函数调用(对jags.parallel)，当给定一个数值参数(如n.iter = 100 )时，它可以工作，但当参数使用变量值n.iter = n.iter时，它就失败了。这看起来可能是jags.parallel中的一个错误 library(R2jags)
model.file <- system.file(package="R2jags=jags.da
浏览 2提问于2013-05-24得票数 4
回答已采纳
1回答
我用以下命令成功地获得了最终样本： data = data, n.adapt = 100)
c('theta', &
浏览 4提问于2016-04-25得票数 0
1回答
我可以安装rjags，但它不能运行，因为它找不到JAGS。下面列出了我的版本，并全部更新到最新版本。我也在用Thinkpad。is set tobut no JAGS installationcan be found there 如果转到安装JAGS的</
浏览 1680提问于2021-10-13得票数 2
我想使用JAGS到R在Tweedie分布变量上运行一个模型。我知道JAGS没有Tweedie分布作为标准，但可以指定一个作为复合Gamma/Poisson。不幸的是，我不知道如何在JAGS中编写它。我基于从各种来源收集的代码编写了以下代码，以简单地尝试并从Tweedie随机变量中恢复平均值、幂和phi参数。它目前不会运行，因为y上的父值无效，大概是因为yi出现在表达式的右侧和左侧。这是在源
浏览 16提问于2020-08-19得票数 0
回答已采纳
1回答
为了进行贝叶斯分析，我使用了这两个包，但是有一些我不明白的区别：其次，在rjags中，我可以说这两段代码是相似的吗？=3)
jsample <- coda.samples(jmod, n.
浏览 1提问于2018-06-08得票数 4
回答已采纳
这应该很快就能回答：
我想把字符串传递给这个参数。解决此问题的一种方法是将字符串写入临时文件并传递文件名。但我希望有一种直接的方法将
浏览 2提问于2013-02-20得票数 7
回答已采纳
我试图使用以下代码(改编自Gelman和Hill‘model中的代码)来估计Jags中的可变系数/拦截有序概率模型。然而，它给我一个“与未观察到的父母在initialization.Try上设置适当的初始值不一致的观察节点”。我哪里出问题了？有人能帮帮我吗？提前谢谢!！rm(list=ls(all=TRUE));library(mvtnorm)library(foreign)
l
浏览 3提问于2015-02-25得票数 0
我试图按照"Lahoz-Monfort JJ，Guillera-Arroita G，Tingley R (2015)统计方法“中提供的指令在R2jags中运行这个脚本。似乎一切正常，但我想不出要看结果的命令.有人能帮忙吗？data = jags.data, inits = jags.inits, parameters.to.save= jags.params, che
浏览 5提问于2021-10-02得票数 0
目前正在讨论当前的np.random.dirichlet函数，因为它不适用于小参数：
mtrand.pyx in mtrand.RandomSta
浏览 2提问于2015-11-11得票数 0
回答已采纳
2回答
R的runjags软件包太棒了。并行功能和使用extend.jags函数的能力使我的生活变得更好了。然而，有时，在我运行一个模型后，我意识到烧成阶段应该更长。如何从run.jags输出中削减额外的样本，以便重新估计参数分布并检查收敛性？jags.object <- run.jags(model, n.chains=3, data=data, monitor =c('a','b&#x
浏览 3提问于2016-06-23得票数 5
回答已采纳
4回答
我想升级一个R包到更新的版本，它已经可供使用。我试过了但它什么也做不了！控制台上没有输出，没有错误，什么也没有。我使用了与install.packages相同的语法，但可能我做错了什么。我一直在看?update.packages，但是我无法弄清楚它是如何工作的，在哪里指定包等等。没有例子。我还尝试使用install.packages更新包，使其再次“安装”，但上面写着“警告:软件包‘R2jags</e
浏览 4提问于2014-01-30得票数 47
回答已采纳
1回答
我最近更新了我的R、Rstudio和rjags包。message:我正在运行一个jags模型，其中包含作为'dgenpareto‘的广义Pareto我使用以下命令行来运行我的jags模型：
> jags.parallel(model.file = bayes_gpd,data = jags.data,parameter
浏览 10提问于2022-08-01得票数 0
2回答
我正在用群组水平预测来拟合一个多层次的logistic回归模型。我通过R使用JAGS，当我将模型与runjags和R2Jags包相匹配时，我得到了不同的行为。，，为什么这
浏览 4提问于2020-02-19得票数 3
回答已采纳
jags.model需要一个包含BUGS模型的文件名作为它的第一个参数。为了将所有内容包含在一个脚本中，我有时使用writeLines命令将BUGS模型写入文件。Y1[i] ~ dnorm(Beta0, Beta1) Beta0 ~ dunif(1, 5)}, "m1.jags")
mod1 <- jags.model("m1.ja
浏览 2提问于2012-05-21得票数 4
回答已采纳
似乎runjags突然(在更新到2.0.3-2版之后)在查找JAGS二进制文件时遇到了问题，并发出一个错误：我通过将这一行放到我的Rprofile中修复了这个问题
.runjags.options <- list(jagspath ="c:/Program
浏览 4提问于2016-11-01得票数 0
2回答
是否有可能误用JAGS作为从具有已知参数的模型生成数据的工具？我需要从预定义的模型中采样数据点，以便进行模拟研究，并测试我在R中开发的模型的威力。不幸的是，该模型有些棘手(带有AR和VAR组件的层次结构)，而且我无法在R中直接模拟数据，而在互联网上搜索时，我发现了一个，其中数据是在JAGS中使用data{}块生成的。在这篇文章中，作者直接估算了<
浏览 5提问于2016-07-10得票数 1
回答已采纳
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号