model
   {
   # Calculate individual (binary) responses to each test from multinomial data
      for (j in 1 : culm[1]) {
         for (k in 1 : T) {
            r[j, k] <- response[1, k]
         }
      }
      for (i in 2 : R) {
         for (j in culm[i - 1] + 1 : culm[i]) {
            for (k in 1 : T) {
               r[j, k] <- response[i, k]
            }
         }
      }
   # Rasch model
      for (j in 1 : N) {
         for (k in 1 : T) {
            logit(p[j, k]) <- beta * theta[j] - alpha[k]
            r[j, k] ~ dbern(p[j, k])
         }
         theta[j] ~ dnorm(0, 1)
      }
   # Priors
      for (k in 1 : T) {
         alpha[k] ~ dnorm(0, 0.0001)
         a[k] <- alpha[k] - mean(alpha[])
      }
      beta ~ dflat()T(0, )
   }