Vglm in Rstan

You can do multinomial regression with brms. There’s an example thread here; searching the boards for “brms multinomial” will turn up more examples.