Я пытаюсь настроить модель Пуассона с нулевой раздувкой в R и JAGS. Я новичок в JAGS, и мне нужно некоторое руководство о том, как это сделать.
Я пытался со следующим, где у [я] является наблюдаемой переменной
model {
for (i in 1:I) {
y.null[i] <- 0
y.pois[i] ~ dpois(mu[i])
pro[i] <- ilogit(theta[i])
x[i] ~ dbern(pro[i])
y[i] <- step(2*x[i]-1)*y.pois[i] + (1-step(2*x[i]-1))*y.null[i]
log(mu[i]) <- bla + bla +bla + ....
theta[i] <- bla + bla + bla + ....
}
}
Однако это не работает, поскольку вы не можете использовать <- для наблюдаемой переменной.
Любые идеи, как изменить / исправить это? Есть ли другой способ настроить модель Пуассона с нулевой раздувкой в JAGS?
r
poisson-distribution
jags
zero-inflation
Джордж Михаэлидис
источник
источник
Ответы:
Вот простое решение, использующее тот факт, что пуассон даст вам нули, когда лямбда-параметр равен нулю. Однако обратите внимание, что JAGS имеет тенденцию ломаться, если лямбда равна нулю, таким образом, «+ 0,00001».
источник
источник