multinomial model code snippets in r