使用Winbugs运行ZINB回归时遇到问题,当我在“规范工具”选项卡中单击“ gen inits”时,它一直显示“负二项式y [1]的阶数必须为整数”。这是我的ZINB winbugs代码。请检查一下。谢谢。
model{
for(i in 1:300){
y[i] ~ dnegbin(p.ind[i], r.ind[i])
p.ind[i] <- r.ind[i]/(r.ind[i]+lambda.ind[i]*(1-u[i]))
u[i] ~ dbern(p[i])
r.ind[i] ~ dgamma(2,1)
logit(p[i])<-a[1]+a[2]*x1[i]+a[3]*x2[i]+a[4]*x3[i]
log(lambda.ind[i]) <- b[1]+b[2]*x1[i]+b[3]*x2[i]+b[4]*x3[i]}
for(k in 1:4){
a[k] ~ dnorm(0,0.001)
b[k] ~ dnorm(0,0.001)}
}
list(a=c(0,0,0,0),b=c(0,0,0,0))