[问答] winbugs程序help [推广有奖]

3关注
42粉丝

已卖：431份资源

学科带头人

41%

还不是VIP/贵宾

威望: 0 级
论坛币: 2302 个
通用积分: 908.3324
学术水平: 114 点
热心指数: 120 点
信用等级: 83 点
经验: 52009 点
帖子: 1552
精华: 1
在线时间: 2357 小时
注册时间: 2005-1-13
最后登录: 2024-5-21

楼主

zhangtao 发表于 2011-6-20 13:04:32 |AI写论文

是否 +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))
thanks!