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.
Related Posts with multilevel multinomial modelling using the gllamm command in stata
Adding 2-Way Plot Options to User Written CommandsLet's say there's a command we're using to do some causal effect estimation. We run the command, and…
poisson regression with instrumental variables and high dimensional fixed effectsI have a poisson regression model with large number of fixed effects and an endogenous independent v…
Converting R date into STATA dateHi, I am trying to convert R dates into STATA dates but unable to do so. For example, the R date i…
Filling a number of missing variables sharing prefix with "replace"Hi Statalisters, I've merged two datasets which are cross-sectional(year of 2004) and time series(y…
Calculating difference between observations in a panel with missing valuesHello! I am working with online auction data (see -dataex- below), where id is the listing identifi…
Subscribe to:
Post Comments (Atom)
0 Response to multilevel multinomial modelling using the gllamm command in stata
Post a Comment