Club STATA 19-04-13
FOREST PLOT con metan
clear
input str12 strata RR LIC LSC
Men 0.69 0.51 0.94
Women 0.73 0.5 1.07
"< 70 yr" 0.73 0.52 1.05
"≥70 yr" 0.71 0.51 0.98
No 0.67 0.45 1.01
Yes 0.71 0.53 0.96
No 1.25 0.64 2.45
Yes 0.65 0.5 0.84
No 0.95 0.64 1.42
Yes 0.6 0.44 0.8
Never 0.67 0.47 0.94
Ever 0.75 0.54 1.03
"< median" 0.76 0.53 1.08
"≥ median" 0.67 0.48 0.93
end
list
gen logRR=ln(RR)
gen SE =ln(LSC/LIC)/3.92
foreach v of varlist RR LIC LSC{
tostring `v', generate(`v's) force format(%3.2f)
}
generate str title = "("+RRs+ "; "+LICs+ "-"+LSCs+ ")"
lab var title "RR (95% CI)"
l strata RR LIC LSC title
ge g=_n
gen inter=""
replace inter="0.62" if g==1
replace inter="0.84" if g==3
replace inter="0.63" if g==5
replace inter="0.06" if g==7
replace inter="0.06" if g==9
replace inter="0.75" if g==11
replace inter="0.72" if g==13
lab var inter "P for interaction"
lab var strata "Subgroups"
gen i=round(floor((g+1)/2))
list g i strata
lab def i 1 "Sex" 2 "Age" 3 "Diabetes" 4 "Hypertension" ///
5 "Dyslipidemia" 6 "Smoking" 7 "BMI (kg/m2)"
lab val i i
sort i
metan logRR SE, by(i) eform
metan logRR SE, by(i) eform ///
lcols(strata) rcols(title inter) nostats nosubgroup nooverall notable ///
textsize (95) astext(45) boxsca(190) ///
boxopt(mcolor(purple)) xlabel(.5, 1, 2) xtick(.05, .75, 1.33, 20) ///
favours ("Mediet better" # "Control better")
+--------------------------------------------------+
| strata RR LIC LSC title |
|--------------------------------------------------|
1. | Men .69 .51 .94 (0.69; 0.51-0.94) |
2. | Women .73 .5 1.07 (0.73; 0.50-1.07) |
3. | < 70 yr .73 .52 1.05 (0.73; 0.52-1.05) |
4. | ≥70 yr .71 .51 .98 (0.71; 0.51-0.98) |
5. | No .67 .45 1.01 (0.67; 0.45-1.01) |
|--------------------------------------------------|
6. | Yes .71 .53 .96 (0.71; 0.53-0.96) |
7. | No 1.25 .64 2.45 (1.25; 0.64-2.45) |
8. | Yes .65 .5 .84 (0.65; 0.50-0.84) |
9. | No .95 .64 1.42 (0.95; 0.64-1.42) |
10. | Yes .6 .44 .8 (0.60; 0.44-0.80) |
|--------------------------------------------------|
11. | Never .67 .47 .94 (0.67; 0.47-0.94) |
12. | Ever .75 .54 1.03 (0.75; 0.54-1.03) |
13. | < median .76 .53 1.08 (0.76; 0.53-1.08) |
14. | ≥ median .67 .48 .93 (0.67; 0.48-0.93) |
+--------------------------------------------------+
|