winbugs程序问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位大虾，小弟编一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))

感激不尽！！！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 程序 winbugs

能不能把你的likelihood写出来看一下？另外在model里，不需要加那么多括号。
按我的理解，你的log-likelihood应该是

1. l_j=log(alpha)-log(beta_j)-beta_j*log(y_j)+alpha*beta_j*(y_n^beta_n)/n

复制代码

如果不是，那么你程序写得有误。而且为什么beta_j~dlnorm ? 需要确定一下你的参数是不是都用对了

已经解决了，谢谢！ 2# junuylia

请问zhouyu738,您是怎么解决“can't bracket slice for node beta[*]”问题的，因为我现在也遇到了类似的问题，求解？

zhouyu738 发表于 2011-6-20 10:26
已经解决了，谢谢！ 2# junuylia

LZ是如何解决这个问题的~~？求助求助！！

请问您的解决方案是怎样的，我也遇到相同的问题，能不能分I系那个一下经验呢