You can build upon a Dirichlet distribution alpha priors. I didn't test it with your data, so, my answer ill be addressing only the conception.
# K = number of parties
# T = number of periods a:e, I guess
model {
for(t in 1:T){
y[t, 1:k] ~ dmulti(alpha[t, 1:k], N[t])
# Dirichlet priors on the paramenters
alpha[t, 1:k] ~ ddirch(theta[1:k])
N[t] <- sum(y[t,1:k])
# Sample size for dmulti based on observed data
# Inference probabilities
delta[t]<-step(alpha[t,2]-alpha[t,3])
}
for(i in 1:k){ #gamma prior for the alpha vector
theta[i] ~ dgamma(0,0.01) }
}