brmultinom function

Bias reduction for multinomial response models using the Poisson trick.