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!
附件列表