model {
for (i in 1:N) {
Score ~ dcat(p[i,])
p[i,1] <- 1 - Q[i,1]
for (r in 2:5) {
p[i,r] <- Q[i,r-1] - Q[i,r]
}
p[i,6] <- Q[i,5]
for (r in 1:5) {
logit(Q[i,r]) <- b.apd*lAPD - c[r]
}
}
for (i in 1:5) {
dc ~ dunif(0, 20)
}
c[1] <- dc[1]
for (i in 2:5) {
c <- c[i-1] + dc
}
b.apd ~ dnorm(0, 1.0E-03)
or.apd <- exp(b.apd)
}