fit_multinomial function

Fit Multinomial Model on Covariate