搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  模型代码.txt
资料下载链接地址: https://bbs.pinggu.org/a-3141381.html
附件大小:
331.76 KB   举报本内容
在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
}
整个的代码在附件中

求求大家帮帮我QAQ


    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

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

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

GMT+8, 2026-1-22 06:59