Ich versuche, in R und JAGS ein Poisson-Modell ohne Inflation aufzubauen. Ich bin neu bei JAGS und brauche eine Anleitung dazu.
Ich habe Folgendes versucht, wobei y [i] die beobachtete Variable ist
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 + ....
}
}}
Dies funktioniert jedoch nicht, da Sie <- für eine beobachtete Variable nicht verwenden können.
Irgendwelche Ideen, wie man das ändert / behebt? Gibt es eine andere Möglichkeit, das Null-Inflations-Poisson-Modell in JAGS einzurichten?
r
poisson-distribution
jags
zero-inflation
George Michaelides
quelle
quelle
Antworten:
Hier ist eine einfache Lösung, bei der das Poisson Nullen liefert, wenn der Lambda-Parameter Null ist. Beachten Sie jedoch, dass JAGS dazu neigt, zu brechen, wenn Lambda genau Null ist, also "+ 0,00001".
quelle
quelle