#multilevel model

WhatsApp: +91-8754446690

cap log close log using "$log_regress/xtmelogit_all outcomes.txt", text replace display "File = $S_FN" timer clear 2 timer on 2 ** Outcomes foreach z in $outcomes { ** Model for Reliability adjustment using empirical bayesian methods xtmelogit `z' xbp_1`z' || hospid4merge:, var est save "$estimates/xtmelogit_`z'", replace estat icc esti stat cap drop re_1`z' predict re_1`z', reffect lab var re_1`z' "random effect" cap drop se_1`z' predict se_1`z', reses lab var se_1`z' "stan. error" cap drop mlp_1`z' predict mlp_1`z' lab var mlp_1`z' "prediction, multilevel" ** Generate adjusted rates, confidence intervals cap drop mean_`z' egen mean_`z' = mean(`z') cap drop xb_mean_`z' gen xb_mean_`z' = logit(mean_`z') cap drop re_xb_1`z' gen re_xb_1`z' = xb_mean_`z' + re_1`z' cap drop rradj_1`z' gen rradj_1`z' = invlogit(re_xb_1`z') lab var rradj_1`z' "rr-adj'd rate" cap drop xb_ciub_1`z' cap drop xb_cilb_1`z' gen xb_ciub_1`z' = re_xb_1`z' +1.96*se_1`z' gen xb_cilb_1`z' = re_xb_1`z' -1.96*se_1`z' cap drop ciub_1`z' cap drop cilb_1`z' gen ciub_1`z' = invlogit(xb_ciub_1`z') gen cilb_1`z' = invlogit(xb_cilb_1`z') } // closes outcomes loop log close timer off 6 save "$data_pt/temp2_pt_mega_0509.dta", replace ** timer cap log close log using "$timers/regression_xtmelogit_all.txt", text replace display "timers for multi-level model" qui:timer list 6 display r(t6)/60 " minutes" log close