我正在尝试在JAGS中构建一个生存模型,允许时间变量的协变量。我希望它是参数模型 - 例如,假设生存遵循Weibull分布(但我想允许危险率变化,因此指数分布过于简单)。因此,这本质上是可以在flexsurv包中完成的内容的贝叶斯版本,该包允许参数模型中的时间变量协变量。 因此,我希望能够以...
我正在尝试编写Winbugs/Jags模型,用于建模多谷主题模型(正是这篇论文-> http://www.ryanmcd.com/papers/mg_lda.pdf)。 在这里,我想基于特定的值选择不同的分布。例如:我想做一些类似于以下的操作:`if ( X[i] > 0.5 ) { ...
我正在尝试在JAGS中拟合逻辑回归模型,但我的数据是以(# 成功数y,# 尝试数n)的形式给出,而不是二元变量。在R中,可以使用“weights”参数对这些数据进行建模,例如glm(y/n ~ ),但我不知道如何在JAGS中完成这个操作。 这里有一个简单的示例,希望能说明我想要询问的问题。请...
我正在尝试在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...
我想学习贝叶斯分析,基于《使用R、JAGS和Stan进行贝叶斯数据分析(2015)》这本书。 这本书中有很多例子。所以,我正在尝试在R中复制这个例子。 但是,在这个例子中,我收到了一个错误消息。 具体而言,这是这个例子的数据。 data y s 1 1 Regina...
我正在使用rJAGS构建多级贝叶斯模型,我想为几个参数指定Cauchy prior。在JAGS中是否有办法实现这一点,还是需要切换到STAN?以下是我的JAGS模型。我想用Cauchy替换dnorm分布,但JAGS找不到标准的R Cauchy分布,例如dcauchy,pcauchy。 mod...
我想知道是否可以在JAGS中将3个链条分配到5或6个核心上。以下是我的代码: library(parallel) # There is no progression bar using parallel jags.parallel(data = d$data,...
我第一次尝试使用rjags拟合一个分层模型,看起来很简单,但是我遇到了一个"Error parsing model file: syntax error on line 5 near "[""的错误提示,我完全不知道如何解释。 你能帮助我并告诉我错在哪里吗? data = list('P....