model{
for(i in 1:n){
y[i] ~dcat(p[i,])
p[i,1] <- phi(r1-mu[i])
p[i,2] <- phi(r1+r2-mu[i])-phi(r1-mu[i])
p[i,3] <- phi(r1+r2+r3-mu[i])-phi(r1+r2-mu[i])
p[i,4] <- phi(r1+r2+r3+r4-mu[i])-phi(r1+r2+r3-mu[i])
p[i,5] <- 1-phi(r1+r2+r3+r4-mu[i])
mu[i] <- b1*X[i,1]+b2*X[i,2]+b3*X[i,3]+b4*X[i,4]}
r1 ~ dgamma(0.1,0.1)
r2 ~ dgamma(0.1,0.1)
r3 ~ dgamma(0.1,0.1)
r4 ~ dgamma(0.1,0.1)
b1 ~ dnorm(0,0.001)
b2 ~ dnorm(0,0.001)
b3 ~ dnorm(0,0.001)
b4 ~ dnorm(0,0.001)}
