我正在尝试在R和JAGS中建立零膨胀泊松模型。我是JAGS的新手,我需要一些有关该操作的指南。
我一直在尝试以下方法,其中y [i]是观察到的变量
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 + ....
}
}
但是,这不起作用,因为您不能在观察到的变量上使用<-。
任何想法如何更改/解决此问题?还有其他方法可以在JAGS中建立零膨胀泊松模型吗?
这可能会有所帮助:分类数据的贝叶斯模型作者
—
P。Congdon