急！update出现undefined real result如何解决？

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

在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)

2020-4-25 09:05:57 上传

整个的代码

求求大家帮帮我QAQ

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖