Estou tentando configurar um modelo de poisson inflado com zero no R e JAGS. Eu sou novo no JAGS e preciso de algumas orientações sobre como fazer isso.
Eu tenho tentado com o seguinte onde y [i] é a variável 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 + ....
}
}
No entanto, isso não funciona, pois você não pode usar <- em uma variável observada.
Alguma idéia de como alterar / corrigir isso? Existe outra maneira de configurar o modelo de poisson inflado com zero no JAGS?
r
poisson-distribution
jags
zero-inflation
George Michaelides
fonte
fonte
Respostas:
Aqui está uma solução simples, usando o fato de que o poisson fornecerá zeros quando o parâmetro lambda for zero. Observe, no entanto, que o JAGS tende a quebrar se lambda for exatamente zero, portanto, o "+ 0,00001".
fonte
fonte