Jags中的Rank函数

文章/答案/技术大牛

发布

1回答

rank、jags、rjags

我在让Jags中的'rank‘函数工作时遇到问题。下面是改编自Winbugs的代码、模型和数据。不起作用的一点是： myorder[i] <- rank(aux.u[],i) ，这将抛出以下错误： RUNTIME ERROR: Incorrect number of parametersin function rank 我知道Jags中的函数<

浏览 30提问于2019-04-11得票数 1

1回答

用于在Jags中建模均值的函数

r、jags

我正在使用jags来模拟贝叶斯框架中的工程反问题。我想知道我是否可以在jags模型中包含一个函数来定义µ参数。] ,1/sigma^2 )} fi ~ dnorm( 24 , 1/4 )当我包含函数时有什么方法可以在模型中包含一个函数吗？谢谢

浏览 0提问于2018-09-04得票数 1

1回答

从R中的runjags对象中删除链

r、bayesian、jags、runjags、rjags

我有一个runjags对象，它有两个混合非常好的链(链1和链3)，而另一个链没有混合(链2)。如何将runjags对象修剪为只包含链1和链3？这是一个使用runjags生成JAGS模型的可重复示例(尽管这里的链混合得很好)。=3,实现此目的的一种方法是将runjags对象转换为mcmc.list，然后使用以下代码删除链： trim.jags <-as.mcmc.list(jags

浏览 2提问于2017-08-22得票数 1

回答已采纳

3回答

在函数中使用jags.parallel ( get(name，envir = envir)中的R语言错误:找不到对象'y‘)

r、parallel-processing、r2jags、jags.parallel

从命令行或脚本中使用jags.parallel很好。我可以在上运行这个修改过的示例。(data=jags.data, inits=jags.inits, jags.params, 但是，如果我将它包装在一个函数中(jagsfit)然后我得到了错误: get中的错误(name，envir = envir)：对象'y‘根据我找到的</

浏览 2提问于2014-05-21得票数 3

回答已采纳

1回答

rjags与r2jags的差异

r、jags、r2jags、rjags

为了进行贝叶斯分析，我使用了这两个包，但是有一些我不明白的区别：其次，在rjags中，我可以说这两段代码是相似的吗？=3) jsample

浏览 1提问于2018-06-08得票数 4

回答已采纳

1回答

奇怪的jags.parallel错误/避免函数调用中的延迟计算

r、lazy-evaluation、jags、r2jags、jags.parallel

我有一个函数调用(对jags.parallel)，当给定一个数值参数(如n.iter = 100 )时，它可以工作，但当参数使用变量值n.iter = n.iter时，它就失败了。这看起来可能是jags.parallel中的一个错误 library(R2jags)但这不是：

浏览 2提问于2013-05-24得票数 4

回答已采纳

1回答

我在R2jags中看不到模型的结果--怎么了？

r、r2jags

我试图按照"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

1回答

在JAGS/BUGS中指定模型时，"is as“符号~非常有用。当使用要求我指定可能性的MCMC方法时，如何在R中做到这一点？比方说，我想估计三个多元正态分布的参数。在JAGS中，我会通过指定pars[1:n] ~ dmnorm(mu[1:3], sigma[1:3, 1:3])来做到这一点。如果所有内容都被正确指定，JAGS将继续在给定分布下估计这些参数。在R中，有类似的函数，如mvtnorm包中</e

浏览 2提问于2020-04-10得票数 1

回答已采纳

1回答

在R包中捆绑.JAGS模型文件

r、package、exec、bundle、jags

我正在编写一个使用作为依赖项的R包。导出的函数需要在内部调用rjags::jags.model("myModel.JAGS")。我觉得我应该将myModel.JAGS文件捆绑在exec文件夹中，即使它不是一个严格意义上的“脚本”。那么我应该如何访问它呢？path <- path.package('myPackage') file <- file.path(path, 'exec'

浏览 4提问于2017-02-12得票数 0

回答已采纳

1回答

runjags模型中的未使用变量警告

jags、r2jags、runjags

我正在通过R包runjags运行JAGS模型。我刚刚从Jags3.4更新到Jags4.0.0，并注意到了一些与更新相关的意外行为。首先，当我运行模型时，我现在得到一个警告消息WARNING: Unused variable(s) in data table:，后面是模型中引用并作为数据提供的数据对象列表。我现在似乎无法重新创建更新失败，但下面是一个可重复的代码示例，演示了奇怪的警告信息。run.jags帮助页上的代码

浏览 5提问于2015-10-09得票数 1

回答已采纳

1回答

附加() ->对象未找到错误(R2Jags)

r、jags、r2jags

我正在使用R2Jags包。我首先得到的模型如下：# The parameters that we are monitoring (mustjags(jags.data, parameters.to.save=jags.params, model.file=model.loc, n.chains = 3, n.burnin=50

浏览 3提问于2015-04-11得票数 0

回答已采纳

2回答

为什么与rjags和R2Jags匹配的模型输出是不同的？

r、bayesian、jags、rjags、r2jags

我正在用群组水平预测来拟合一个多层次的logistic回归模型。我通过R使用JAGS，当我将模型与runjags和R2Jags包相匹配时，我得到了不同的行为。，，为什么这

浏览 4提问于2020-02-19得票数 3

回答已采纳

1回答

如何在R中保存JAGS模型对象？

r、jags

我使用包rjags在R中执行MCMC，并且我希望保存函数jags.model的输出，以便以后在另一个R会话中使用。下面是正态分布平均值的一个简单例子：N <- 1000model.str <- 'model {for (i in 1:N) {list(x = x, N = N))我可以生成这样的mu示例： coda.samples

浏览 1提问于2014-08-21得票数 5

回答已采纳

1回答

具有收敛性测试的并行RJAGS

r、mcmc、jags

问题是，如果我需要根据诊断值恢复运行，那么重新编译的链将从第一个初始化的先验值重新启动，而不是在链停止的参数空间中的位置。如果我不重新编译模型，RJAGS会抱怨。当链停止时，是否有一种方法来存储链的位置，这样我就可以从我停止的地方重新初始化？这里我将给出一个非常简单的例子。[[i]]$recompile() jags.samples <- coda.samples(j

浏览 0提问于2015-04-06得票数 6

回答已采纳

2回答

在runjags中运行JAGS模型后延长刻录期

r、bayesian、jags、runjags

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

回答已采纳

1回答

如何使字符向量的内容在R中显示为文件？

r、jags

jags.model需要一个包含BUGS模型的文件名作为它的第一个参数。为了将所有内容包含在一个脚本中，我有时使用writeLines命令将BUGS模型写入文件。Y1[i] ~ dnorm(Beta0, Beta1) Beta0 ~ dunif(1, 5)}, "m1.jags") mod1 <- jags.model("m1.jags&q

浏览 2提问于2012-05-21得票数 4

回答已采纳

1回答

jags线性回归，当先验相互依赖时，如何设置先验？

r、linear-regression、bayesian、mcmc、jags

library(rjags)for (i in 1:n){mu[i] <- beta0 + beta1*x[= list(x = c(5,1,2,3,4),update(

浏览 2提问于2020-12-23得票数 0

回答已采纳

1回答

贾格斯(只是另一个吉布斯·采样器)能处理常微分方程吗？

jags、winbugs

我搜索了JAGS的手册和一篇关于常微分方程( ode )的帖子()。我认为JAGS可以，因为它类似于WinBUGS (通过WBdiff接口完成)。然而，如果我让JAGS在我的ode代码中阅读，它甚至无法识别D(y[...], t)表达式。 JAGS能处理ode吗？也许我错过了JAGS的插件，比如WBdiff。

浏览 10提问于2022-10-15得票数 1

回答已采纳

1回答

将R2jags对象转换为Stanreg (rstanarm)对象

r、r2jags、rstanarm

我用R2jags做了一个模型。我喜欢jags语法，但是我发现R2jags产生的输出并不容易使用。我最近读到了关于rstanarm包的内容。它有许多有用的功能，并且得到了tidybayes和bayesplot包的良好支持，以便轻松地进行模型诊断和可视化。但是，我不喜欢用rstanarm编写模型的语法。理想情况下，我希望在这两个领域中取得最佳效果，即用R2jags编写模型，并将输出转换为Stanreg对象以使用rstanarm函数

浏览 4提问于2020-09-09得票数 0

1回答

jags.parallel - get(name，envir = envir)错误:第一个参数无效

r、winbugs、jags、r2jags、jags.parallel

在使用jags.parallel时，我得到以下错误：+nc, ni, nb, nt);使用jags函数的相同调用运行正常。我只找到了，但只有一个推测性的建议不适用，也不适用于这里。可重现<e

浏览 5提问于2013-07-23得票数 5

回答已采纳

点击加载更多