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 ~ dnorm(0,0.0001) I(0, ) }