| 所在主题: | |
| 文件名: 333.rar | |
| 资料下载链接地址: https://bbs.pinggu.org/a-926655.html | |
本附件包括:
|
|
| 附件大小: | |
|
winbugs程序,在specification时均已经完成,且无错误。而在“update”环节出现如下提示:“can't bracket slice for node beta[*]”,点击一次update会出现不同的错误提示中beta[*]不同。查资料讲是变量的先验分布参数太离散,但修改后依然有此问题,故求高手指点。程序如下:
model{ for(i in 1:n){zeros[i]<-0 zeros[i]~dpois(l[i])} for(j in 1:n){l[j]<--log(alpha)-(log(beta[j]))-((beta[j])*log(y[j]))+(alpha*(beta[j])*((pow(y[n],beta[n])))/n)} alpha~dgamma(0.1,0.001) for(j in 1:n){beta[j]~dlnorm(mu,tau)} mu~dnorm(0,0.01) tau~dgamma(0.1,0.1) s<-sqrt(1/tau) } #Data list(n=19,y=c(0.3945,0.7435,0.9055,1.1075,1.4631,1.8113,1.9136,2.1390,2.3514,2.7121,2.7612,2.7807,2.8761,2.9606,2.9802,3.0674,3.1310,3.1719,3.2531)) #Initial values list(alpha=44,mu=0.01,tau=0.01,beta=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)) thanks! |
|
熟悉论坛请点击新手指南
|
|
| 下载说明 | |
|
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。 2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。 3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。 (如有侵权,欢迎举报) |
|
京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明