Sto cercando di impostare un modello poisson gonfiato a zero in R e JAGS. Sono nuovo di JAGS e ho bisogno di alcune indicazioni su come farlo.
Ho provato con quanto segue in cui y [i] è la variabile osservata
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 + ....
}
}
Tuttavia, ciò non funziona poiché non è possibile utilizzare <- su una variabile osservata.
Qualche idea su come cambiare / riparare questo? Esiste un altro modo per impostare il modello poisson gonfiato a zero in JAGS?