Two independent multinomial variates have ordered response categories.

Using independent Dirichlet priors, explain how to simulate to find the posterior probability that those two distributions are stochastically ordered.

(Altham (1969) found an exact expression for this probability.)