生成上图的Stata代码如下：

* Figure:Side by side horizontal bar plot (Outcomes by City & Case) using weightab

    global graph_opts ///
        title(, justification(left) color(black) span pos(11)) ///
        graphregion(color(white) lc(white) lw(med) la(center)) /// <- remove la(center) for Stata < 15
        ylab(,angle(0) nogrid) xtit(,placement(left) justification(left)) ///
        yscale(noline) xscale(noline) legend(region(lc(none) fc(none)))

    local opts lw(thin) lc(white) la(center)

    label def case 1 "Case 1" 2 "Case 2" 3 "Case 3" 4 "Case 4" , modify

    qui do "http://ds.epiman.cn/uploads/2022/06/weightab.ado"

    use "http://ds.epiman.cn/uploads/2022/06/bar-weightab.dta" , clear

    weightab ///
        correct treat_cxr re_3 re_4 treat_refer t_12 ///
        med_any med_l_any_1 med_l_any_2 med_l_any_3  med_k_any_9   ///
        if city == 2 ///
        [pweight = weight_city] ///
        , ///
        ${graph_opts} ///
        barlab ///
        barlook(1 `opts' fi(100)) ///
        title("Patna") ///
        over(case) ///
        graph ///
        legend(off) ///
        xlab(${pct}) ///
        name(Fig_1_1)
		
    weightab ///
        correct treat_cxr re_3 re_4 treat_refer t_12 ///
        med_any med_l_any_1 med_l_any_2 med_l_any_3  med_k_any_9  ///
        if city == 3 ///
        [pweight = weight_city] ///
        , ///
        ${graph_opts} ///
        barlab ///
        barlook(1 `opts'  fi(100)) ///
        title("Mumbai") ///
        over(case) ///
        graph ///
        legend(pos(5) ring(0) c(1) symxsize(small) symysize(small)) ///
        xlab(${pct}) ///
        name(Fig_1_2)
		
    graph combine Fig_1_1 Fig_1_2 ///
        , ///
        ${comb_opts} ///
        xsize(7) r(1)

* Have a lovely day!
* Source: https://worldbank.github.io/stata-visual-library/bar-weightab.html