J'essaie de mettre en place un modèle de poisson zéro gonflé dans R et JAGS. Je suis nouveau chez JAGS et j'ai besoin de conseils sur la façon de le faire.
J'ai essayé avec ce qui suit où y [i] est la variable observée
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 + ....
}
}
Cependant, cela ne fonctionne pas car vous ne pouvez pas utiliser <- sur une variable observée.
Avez-vous des idées pour changer / corriger cela? Existe-t-il une autre façon de configurer le modèle de poisson zéro gonflé dans JAGS?
r
poisson-distribution
jags
zero-inflation
George Michaelides
la source
la source
Réponses:
Voici une solution simple utilisant le fait que le poisson vous donnera des zéros lorsque le paramètre lambda est nul. Notez cependant que JAGS a tendance à se casser si lambda est exactement nul, donc le "+ 0,00001".
la source
la source