Saya mencoba mengatur model poisson nol-meningkat di R dan JAGS. Saya baru mengenal JAGS dan saya perlu beberapa panduan tentang cara melakukannya.
Saya sudah mencoba dengan yang berikut ini di mana y [i] adalah variabel yang diamati
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 + ....
}
}
Namun ini tidak berfungsi karena Anda tidak dapat menggunakan <- pada variabel yang diamati.
Ada ide bagaimana mengubah / memperbaikinya? Apakah ada cara lain untuk mengatur model poisson zero-inflated di JAGS?
r
poisson-distribution
jags
zero-inflation
George Michaelides
sumber
sumber
Jawaban:
Berikut adalah solusi sederhana menggunakan fakta bahwa poisson akan memberi Anda nol ketika parameter lambda nol. Perhatikan bahwa JAGS cenderung pecah jika lambda persis nol, dengan demikian "+ 0,00001".
sumber
sumber