Я пытаюсь использовать регрессию для обоих параметров биномиальной пуассоновской смеси. Количество событий распределено по Пуассону, но я могу обнаружить эти события с некоторой биномиальной вероятностью. Когда я делаю только пуассоновскую регрессию, она работает, но регрессия для биномиальной части, похоже, вызывает ошибку «узел несовместим с родителями».
Мой код указан ниже.
mod_string = " model{
## Likelihood
for(i in 1:length(Count)){
Count[i] ~ dbin(p[i],n[i])
n[i]~dpois(lambda[i])
lambda[i]=exp(b[1]+b[2]*Time[i])
logit(p[i])=a[1]+a[2]*Time[i]
}
for (i in 1:2){
b[i]~ dnorm(0.0,1.0/1e6)
a[i]~ dnorm(0.0,1.0/1e6)
}
}"
Это говорит так:
Error in jags.model(textConnection(mod_string), data = data_jags, n.chains = 3) :
Error in node Count[2]
Node inconsistent with parents
В чем здесь проблема?