I am running a model in which I am looking at differences in odds of binge drinking (1=yes, 0=no) across 47 groups. I would like to obtain the predicted probabilities for these groups (not individual level probabilities). Here is my model:
Code: Select all
runmlwin binge cons, ///
level2(group: cons) ///
level1(idcode:) ///
discrete(distribution(binomial) link(logit) denominator(denominator) pql2) ///
rigls maxiterations(100) ///
nopause
runmlwin binge cons, ///
level2(group: cons, residuals(u, savechains("binge_null_u.dta", replace))) ///
level1(idcode:) ///
discrete(distribution(binomial) link(logit) denominator(denominator)) ///
mcmc(burnin(5000) chain(50000) thinning(50)) initsprevious ///
nopause