全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
3793 5
2011-06-13
各位大虾,小弟编一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))

感激不尽!!!!
二维码

扫码加我 拉你入群

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

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

全部回复
2011-6-18 01:18:44
能不能把你的likelihood写出来看一下?另外在model里,不需要加那么多括号。
按我的理解,你的log-likelihood应该是
复制代码
如果不是,那么你程序写得有误。而且为什么beta_j~dlnorm ? 需要确定一下你的参数是不是都用对了
二维码

扫码加我 拉你入群

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

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

2011-6-20 10:26:27
已经解决了,谢谢! 2# junuylia
二维码

扫码加我 拉你入群

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

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

2013-6-6 21:43:26
请问zhouyu738,您是怎么解决“can't bracket slice for node beta[*]”问题的,因为我现在也遇到了类似的问题,求解?
二维码

扫码加我 拉你入群

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

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

2015-3-19 16:39:41
zhouyu738 发表于 2011-6-20 10:26
已经解决了,谢谢! 2# junuylia
LZ是如何解决这个问题的~~?求助求助!!
二维码

扫码加我 拉你入群

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

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

2016-1-4 20:13:05
请问您的解决方案是怎样的,我也遇到相同的问题,能不能分I系那个一下经验呢
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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