全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
2476 2
2011-06-20
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!
附件列表

333.rar

大小:846 Bytes

 马上下载

本附件包括:

  • 333.odc

二维码

扫码加我 拉你入群

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

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

全部回复
2011-9-6 14:25:34
对各个参数分别用不同的方法抽样,在OpenBUGS里
二维码

扫码加我 拉你入群

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

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

2014-1-27 19:32:18
你好,我看您几年前有篇关于winbugs的问题cannot bracket slice for node。。。我也遇到这个问题,不知道您解决了没有呢
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群