runmlwin causing MLwiN to load all observations
Posted: Mon Oct 17, 2022 4:13 pm
Hi There - I am trying to run a 5-level model in which individuals are nested in households, villages, districts, and states in India. However, when I execute my runmlwin code in Stata, it opens MLwiN and the MLwiN begins to read in all the data file observations as though I am loading a new sheet directly in to MLwiN. This does not happen when I run a 4-level model without the individuals. Any thoughts on why this might be occurring? My code is below. Thanks!
*4-LEVEL CODE (works fine)
xi: runmlwin santype1 constant,level4(hv024: constant, residuals(res_state)) level3(shdist: constant, residuals(res_dist)) level2(hv001: constant, residuals(res_clust)) ///
level1(case:, ) discrete(distribution(binomial) link(logit) denominator(constant)) mlwinsettings(optimat)nopause
xi: runmlwin santype1 constant,level4(hv024: constant, residuals(res_state1)) level3(shdist: constant, residuals(res_dist1)) level2(hv001: constant, residuals(res_clust1)) ///
level1(case:, ) discrete(distribution(binomial) link(logit) denominator(constant)) mlwinsettings(optimat) mcmc(on) initsprevious nopause
ge beta0 = _b[constant]
*5-LEVEL CODE (not working)
xi: runmlwin santype1 constant,level5(hv024: constant, residuals(res_state)) level4(NFHS5_Code: constant, residuals(res_dist)) level3(hv001: constant, residuals(res_clust)) ///
level2(case: constant, residuals(res_case)) level1(id:, ) discrete(distribution(binomial) link(logit) denominator(constant)) mlwinsettings(optimat)nopause
xi: runmlwin santype1 constant,level5(hv024: constant, residuals(res_state1)) level4(NFHS5_Code: constant, residuals(res_dist1)) level3(hv001: constant, residuals(res_clust1)) ///
level2(case: constant, residuals(res_case1)) level1(id:, ) discrete(distribution(binomial) link(logit) denominator(constant)) mlwinsettings(optimat) mcmc(on) initsprevious nopause
ge beta0 = _b[constant]
*4-LEVEL CODE (works fine)
xi: runmlwin santype1 constant,level4(hv024: constant, residuals(res_state)) level3(shdist: constant, residuals(res_dist)) level2(hv001: constant, residuals(res_clust)) ///
level1(case:, ) discrete(distribution(binomial) link(logit) denominator(constant)) mlwinsettings(optimat)nopause
xi: runmlwin santype1 constant,level4(hv024: constant, residuals(res_state1)) level3(shdist: constant, residuals(res_dist1)) level2(hv001: constant, residuals(res_clust1)) ///
level1(case:, ) discrete(distribution(binomial) link(logit) denominator(constant)) mlwinsettings(optimat) mcmc(on) initsprevious nopause
ge beta0 = _b[constant]
*5-LEVEL CODE (not working)
xi: runmlwin santype1 constant,level5(hv024: constant, residuals(res_state)) level4(NFHS5_Code: constant, residuals(res_dist)) level3(hv001: constant, residuals(res_clust)) ///
level2(case: constant, residuals(res_case)) level1(id:, ) discrete(distribution(binomial) link(logit) denominator(constant)) mlwinsettings(optimat)nopause
xi: runmlwin santype1 constant,level5(hv024: constant, residuals(res_state1)) level4(NFHS5_Code: constant, residuals(res_dist1)) level3(hv001: constant, residuals(res_clust1)) ///
level2(case: constant, residuals(res_case1)) level1(id:, ) discrete(distribution(binomial) link(logit) denominator(constant)) mlwinsettings(optimat) mcmc(on) initsprevious nopause
ge beta0 = _b[constant]