multinom_reg function

Multinomial regression