395得票1回答
在JAGS中以“计数过程”形式表示参数生存模型

我正在尝试在JAGS中构建一个生存模型,允许时间变量的协变量。我希望它是参数模型 - 例如,假设生存遵循Weibull分布(但我想允许危险率变化,因此指数分布过于简单)。因此,这本质上是可以在flexsurv包中完成的内容的贝叶斯版本,该包允许参数模型中的时间变量协变量。 因此,我希望能够以...

13得票3回答
如何将mcmc.list转换为bugs对象?

我正在使用R库中的rjags。函数coda.samples生成一个mcmc.list, 例如(来自example(coda.samples)):library(rjags) data(LINE) LINE$recompile() LINE.out <- coda.samples(LINE...

11得票1回答
在WinBugs/JAGS中基于if-else条件选择不同的分布

我正在尝试编写Winbugs/Jags模型,用于建模多谷主题模型(正是这篇论文-> http://www.ryanmcd.com/papers/mg_lda.pdf)。 在这里,我想基于特定的值选择不同的分布。例如:我想做一些类似于以下的操作:`if ( X[i] > 0.5 ) { ...

10得票1回答
当响应变量为比例时的逻辑回归(使用JAGS)

我正在尝试在JAGS中拟合逻辑回归模型,但我的数据是以(# 成功数y,# 尝试数n)的形式给出,而不是二元变量。在R中,可以使用“weights”参数对这些数据进行建模,例如glm(y/n ~ ),但我不知道如何在JAGS中完成这个操作。 这里有一个简单的示例,希望能说明我想要询问的问题。请...

9得票1回答
在r中,具有iid随机效应的泊松GLM产生奇怪的输出

我正在尝试在R(通过Rstudio)中运行rjags来估计参数alpha和beta以及模型的超参数tau.nu: y_i|x_i~pois(eta_i), eta_i=exp(alpha + beta*x_i + nu_i), nu_i~N(0,tau.nu) 这是我的代码: #gen...

7得票1回答
Rjags错误信息:维度不匹配。

我想学习贝叶斯分析,基于《使用R、JAGS和Stan进行贝叶斯数据分析(2015)》这本书。 这本书中有很多例子。所以,我正在尝试在R中复制这个例子。 但是,在这个例子中,我收到了一个错误消息。 具体而言,这是这个例子的数据。 data y s 1 1 Regina...

7得票1回答
JAGS中的Cauchy先验

我正在使用rJAGS构建多级贝叶斯模型,我想为几个参数指定Cauchy prior。在JAGS中是否有办法实现这一点,还是需要切换到STAN?以下是我的JAGS模型。我想用Cauchy替换dnorm分布,但JAGS找不到标准的R Cauchy分布,例如dcauchy,pcauchy。 mod...

7得票1回答
在多个核上运行JAGS,同时运行多个链(将链子分成若干份)是否可行?

我想知道是否可以在JAGS中将3个链条分配到5或6个核心上。以下是我的代码: library(parallel) # There is no progression bar using parallel jags.parallel(data = d$data,...

7得票1回答
rjags出现语法错误

我第一次尝试使用rjags拟合一个分层模型,看起来很简单,但是我遇到了一个"Error parsing model file: syntax error on line 5 near "[""的错误提示,我完全不知道如何解释。 你能帮助我并告诉我错在哪里吗? data = list('P....