全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
1433 0
2018-05-14
在winbugs中,前面几步都没问题,但是update的时候,就提示 cannot bracket slice for node tt ,请问怎么能解决呢?
下面是程序:
model
{
for(i in 1:N){
r[i]~dnegbin(u[i],t)
b[i]~dnorm(0.0,tau)
logit(u[i])<-alpha0+alpha1*x1[i]+alpha2*x2[i]+alpha3*x3[i]+alpha4*x4[i]+alpha5*x5[i]+b[i]
x4[i]<-alpha41*x41[i]+alpha42*x42[i]+alpha43*x43[i]+alpha44*x44[i]
x5[i]<-alpha51*x51[i]+alpha52*x52[i]+alpha53*x53[i]+alpha54*x54[i]
}
alpha0~dnorm(0.0,1.0E-6)
alpha1~dnorm(0.0,1.0E-6)
alpha2~dnorm(0.0,1.0E-6)
alpha3~dnorm(0.0,1.0E-6)
alpha4~dnorm(0.0,1.0E-6)
alpha5~dnorm(0.0,1.0E-6)
alpha41~dnorm(0.0,1.0E-6)
alpha42~dnorm(0.0,1.0E-6)
alpha43~dnorm(0.0,1.0E-6)
alpha44~dnorm(0.0,1.0E-6)
alpha51~dnorm(0.0,1.0E-6)
alpha52~dnorm(0.0,1.0E-6)
alpha53~dnorm(0.0,1.0E-6)
alpha54~dnorm(0.0,1.0E-6)
sigma~dunif(0,100)
tau<-1/(sigma*sigma)
tt~dgamma(1000,1000)
t<-1/tt
}
list(N=240,
r=c(0,0,0,0,0,0,0,0,3,0,1,4,2,2,2,3,1,1,2,3,5,3,0,0,0,0,0,0,0,0,0,4,6,9,9,3,11,9,9,5,5,10,13,10,5,2,0,0,0,0,0,0,0,0,0,0,2,1,2,3,1,2,3数据太长了就不贴了
list(alpha0=0,alpha1=0,alpha2=0,alpha3=0,alpha4=0,alpha5=0,alpha41=0,alpha42=0,alpha43=0,alpha44=0,alpha51=0,alpha52=0,alpha53=0,alpha54=0,sigma=1,tt=1)
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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