Estoy tratando de configurar un modelo de Poisson inflado a cero en R y JAGS. Soy nuevo en JAGS y necesito alguna orientación sobre cómo hacerlo.
He estado intentando con lo siguiente donde y [i] es la variable observada
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 + ....
}
}
Sin embargo, esto no funciona ya que no puede usar <- en una variable observada.
¿Alguna idea de cómo cambiar / arreglar esto? ¿Hay alguna otra forma de configurar el modelo de Poisson inflado a cero en JAGS?
Esto puede ser útil: modelos bayesianos para datos categóricos Por P. Congdon
—
user4618