Club stata 19-04-13 forest plot con

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

gen SE =ln(LSC/LIC)/3.92

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

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) |

+--------------------------------------------------+