求助贴！！！openbugs 代码部分mutiple defination，毕业论文数据处理急求大神帮助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我的代码是这样的，求助各位大神帮忙：
model
{
     for (i in 1:n) {
       for (j in 1:p) {
           Y[i , j] ~ dcat ( prob [i , j , 1: 5])
        }
        theta [i] ~ dnorm (mu[i], tau)
       e[i] ~ dnorm (0,tau)
       x[i] ~ dnorm (0,1)
       for (j in 1:p) {
           for (k in 1:4) {
             logit(P [i , j , k]) <- kappa [j , k] - alpha [j] * theta [i]
             theta[i] <- b0 + b1 * x[i] + e[i]
             mu[i] <- b0 + b1 * x[i]
           }
          P[i , j , 5] <- 1.0
         }

        for (j in 1:p) {
          prob [i , j , 1] <- P[i , j , 1]
          for (k in 2:5) {
             prob [i , j , k] <- P[i , j , k] - P[i , j , k -1]
          }
        }
     }

     for (j in 1:p) {
        alpha [j] ~ dnorm (m.alpha , pr.alpha ) I(0 , )
     }
     pr.alpha <- pow(s.alpha , -2)

     for (j in 1:p) {
        for (k in 1:4) {
          kappa.star [j , k] ~ dnorm (m.kappa , pr.kappa )
          kappa [j , k] <- ranked(kappa.star [j , 1:4 ], k)
        }
     }
     pr.kappa <- pow(s.kappa , -2)

     for (i in 1:n) {
        tau <- 1/sigma2
        log(sigma2) <- 2 * log.sigma
        log.sigma ~ dflat()
       }
   }

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：openbugs mutiple Nation 论文数据处理 求大神帮助