Próbuję skonfigurować model Poissona z napompowaniem zerowym w R i JAGS. Jestem nowy w JAGS i potrzebuję wskazówek, jak to zrobić.
Próbowałem z następującymi, gdzie y [i] jest obserwowaną zmienną
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 + ....
}
}
Nie działa to jednak, ponieważ nie można użyć <- dla obserwowanej zmiennej.
Wszelkie pomysły, jak to zmienić / naprawić? Czy istnieje inny sposób skonfigurowania modelu Poissona z napompowaniem zerowym w JAGS?
r
poisson-distribution
jags
zero-inflation
George Michaelides
źródło
źródło
Odpowiedzi:
Oto proste rozwiązanie wykorzystujące fakt, że poissona da ci zera, gdy parametr lambda wynosi zero. Zauważ jednak, że JAGS ma tendencję do łamania się, jeśli lambda ma dokładnie zero, więc „+ 0,00001”.
źródło
źródło