hierarchical multinomial marginal models