fit_multinom_model function

Fit Simple Multinomial Model