When using the glm function: which is the best distribution for multinomial regression --- is it "binomial"? If the glm function doesn't have an appropriate distribution for multinomial regression, which r package has a better function? Thank you!
You may try the nnet::multinom()
You can find the documentation at here.