Bayesian inference on group differences in multivariate categorical data