求助winbugs报错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我在写个程序    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])
       }   
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 程序