With a categorical variable t and some covariate x, I run the random-effects multinomial logit:
t
x
model <- mblogit(t ~ x, data = sa