estimate_multinomial_rsp function

Estimate proportions of each level of a variable