全部版块 我的主页
论坛 提问 悬赏 求职 新闻 读书 功能一区 经管百科 爱问频道
1222 0
2011-08-11
我在写个程序    winbugs怎么都报错  说是empty slot not allowed in variables name 但是我检查了好多遍    变量名里面没有slot啊请教一下     到底哪里出问题了
model
{
    for(i in 1:N)
      {
         for(j in 1:(N-i+1))
           {
               z[i,j]<-y[i,j]/1000
               log(mu[i,j])<-alpha[i]+beta[j]
               zeros[i,j]<-0
               zeros[i,j]~dpois(zeros.mean[i,j])
               zeros.mean[i,j]<-(mu[i,j]-z[i,j]*log(mu[i,j])+loggam(z[i,j]+1))/scale
           }
      }
    for(i in 2:N)
      {
          for(j in (N-i+2):N)
           {
                nu[i,j]<-mu[i,j]/scale
                y[i,j]~dpois(nu[i,j])
                log(mu[i,j])<-alpha[i]+beta[j]
                z[i,j]<-scale*y[i,j]
           }
      }
    for(i in 1:N)
      {
          for(j in 1:N)
           {
                fit[i,j]<-z[i,j]*1000
           }
      }  
    scale<-3.24
    a[1]<-426.843
    alpha[1]<-log(a[1])
    for(i in 2:N)
      {   
           a[i]~dgamma(ula[i],ulb[i])
           alpha[i]<-log(a[i])
           ula[i]<-pow(ul.mean[i-1],2)/ul.var[i-1]
           ulb[i]<-ul.mean[i-1]/ul.var[i-1]
       }
    ul.mean<-c(354, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, 400)
    ul.var<-c(1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000, 1000000)
    for(j in 1:N)
       {
           p[j]~dgamma(0.0001,0.001)
       }  
    s<-sum(p[1:N])
    for(j in 1:N)
       {
           b[j]<-p[j]/s
           beta[j]<-log(b[j])
       }   
}

二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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