Question

我是RJAGS的新手，并且正在使用rJAGS的简单模型进行练习。我的模型规格如下：

model{
#Likelihood 
for( i in 1 : N) {<br/>
 Disease[i] ~ dbern (thetacombine[i])<br/>
 thetacombine[i]<-theta1new[i]*theta0new[i]<br/>
 theta1new[i]<-pow((theta1[i]+0.0000000001),G[i])<br/>
 theta0new[i]<-pow((theta0[i]+0.0000000001),(G[i]-1))<br/>
}
  #Priors regressors <br/>
  theta1 ~ dbeta (a1,b1)<br/>
  theta0 ~ dbeta (a2,b2)<br/> 
}

我的数据如下所示：
ID G疾病
1 1 0
2 0 1
3 1 0
4 1 0
5 0 0

我指定了运行模型的代码为：

bayes_data <- list(
Disease = Disease,
G = G,
N = ntot,
a1=5,
b1=1,
a2=1,
b2=3
 )

bayes_posterior_jags <- jags.model(
 file = bayes_model_str,
 data = bayes_data,
 inits = bayes_init,
 n.chains = 3,             
 n.adapt = 100
 )<br/>

我收到一条错误消息：第9行的编译错误。索引超出范围，占用了theta0的子集。
可能是什么原因，我该如何解决？非常感谢！！！

JAGS出现错误：

0 个答案: