Good day everyone,
I am using the gllamm command in stata 15.0 for multilevel (2-leveled) multinomial logistic regression. The outcome/dependent variable is modern contraceptive use among women and it has 4 sub-categories showing the users of the different types of modern contraceptives [(1)permanent method and Long-acting reversible contraceptive users; (2)medium acting reversible contraceptive users; (3)short-acting reversible contraceptive users], each in contrast to the traditional users/those not using any method at all (this is the reference sub-category). For my study, individual women (level-1) are nested within communities (level-2), and my main independent variable; which is the economic empowerment of women, is also introduced as a random slope. In relation to this, I typed in the following command:
gllamm PLMSorNotnew CWEIcentre, i(ClusterNumber) link(mlogit) family(binomial) adapt nrf(2) eqs(inter slope) pweight(pw)

However, only one random intercept variance, one random slope variance, and thus one covariance were estimated instead of three each.

My understanding is that the single random intercept variance shows the between-community difference in the odd of generally using modern contraceptives, irrespective of the particular type, as opposed to traditional users/those not using any method at all. the single random slope variance, on the other hand, shows the between-community difference in the regression effect of women's economic empowerment on the general use of modern contraceptives, irrespective of the particular type. Please, what is the correct command to produce the three different random intercept variances, the three different random slope variances, and the three different covariances.
Thank you.