在winbugs里面建好了贝叶斯框架下的条件自回归NB模型initial采用的是gen inits

但是在update的时候遇到了undefined real result

请问大神们是哪里出现了问题？
才开始学着用winbugs，虽然看了用户手册，但还是不明白怎么改QAQ
model代码如下：
#Model
model
{
    for (i in 1:N) {
        thi~dgamma(1,beta)
        thita<-log(thi)
        log(mu) <- alpha0+ alpha1*log(DVMT)  + alpha2 *pop/1000 +alpha3*sixt+alpha4*income+alpha5*job+alpha6*pub+alpha7*home+alpha8*park+alpha9*air+alpha10*street+alpha11*bike+alpha12*yiji+alpha13*erji+alpha14*path+alpha15*jiaochakou+alpha16*fwy+fai+thita
        sev ~ dpois(mu)
    }
    # CAR prior distribution for random effects:
    fai[1:N] ~ car.normal(adj[], weights[], num[], tau)
    for (k in 1:sumNumNeigh) {
        weights[k] <- 1
    }
    # Other priors:
    be~dnorm(0,0.01)
    beta<-exp(be)
    alpha0 ~ dnorm(0, 1000)
        alpha1 ~ dnorm(0, 1000)
        alpha2 ~ dnorm(0, 1000)
    alpha3 ~ dnorm(0, 1000)
    alpha4 ~ dnorm(0, 1000)
        alpha5 ~ dnorm(0, 1000)
        alpha6 ~ dnorm(0, 1000)
        alpha7 ~ dnorm(0, 1000)
        alpha8 ~ dnorm(0, 1000)
        alpha9 ~ dnorm(0, 1000)
        alpha10 ~ dnorm(0, 1000)
        alpha11 ~ dnorm(0, 1000)
        alpha12 ~ dnorm(0, 1000)
        alpha13 ~ dnorm(0, 1000)
        alpha14 ~ dnorm(0, 1000)
        alpha15 ~ dnorm(0, 1000)
        alpha16 ~ dnorm(0, 1000)
    tau ~ dgamma(0.5, 0.0005)
    # prior on precision
    # standard deviation
}
整个的代码在附件中

模型代码.txt
大小:(331.76 KB)

 马上下载

求求大家帮帮我QAQ