model {
  for (i in 1:T) {
    y[Event[i],Nominee[i]] <- Choice[i]
    for (j in 1:D) {Z[j,Event[i],Nominee[i]] <- X[i,j]}
    }
  for (j in 1:D) {
    delta[j] ~ dnorm(0,0.1)
#    delta[j] ~ dnorm(0,0.001)
    }
  for (i in 1:N) {
    y[i,1:(P[i])] ~ dmulti(p[i,1:(P[i])], 1)
    for (k in 1:(P[i])) {
      log(phi[i,k]) <- inprod2(delta[],Z[,i,k])
      p[i,k]  <- phi[i,k] / sum(phi[i,1:(P[i])])
      }
    }
}