rho_mlogit function

Fitting disease models via multinomial logistic models