| 所在主题: | |
| 文件名: 模型代码.txt | |
| 资料下载链接地址: https://bbs.pinggu.org/a-3141381.html | |
| 附件大小: | |
|
在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、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。 (如有侵权,欢迎举报) |
|
京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明