Winbugs中的ZINB错误代码:负二项式y [1]的顺序必须为整数

时间:2019-05-15 05:16:34

标签: regression bayesian winbugs

使用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))

0 个答案:

没有答案