model{
  y ~ dmulti(theta,n)
  theta ~ ddirch(alpha)
  p1 <- step(theta[1]-0.5)
  p2 <- step(theta[2]-theta[3])
}