model
   {
       for( i in 1 : N ) {
          beta[i , 1 : 2] ~ dmnorm(mu.beta[], R[ , ])
          for( j in 1 : T ) {
             Y[i, j] ~ dnorm(mu[i , j], tauC)
             mu[i, j] <- beta[i, 1] + beta[i, 2] * x[j]
          }
       }

       mu.beta[1 : 2] ~ dmnorm(mean[], prec[ , ])
       R[1 : 2 , 1 : 2] ~ dwish(Omega[ , ], 2)
       tauC ~ dgamma(0.001, 0.001)
       sigma <- 1 / sqrt(tauC)
   }