Я новичок в 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 02 0 13 1 04 1 05 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. Индекс выходит за пределы диапазона, принимая подмножество тета0.в чем может быть причина и как мне это исправить?Большое спасибо !!!