方法一：Stata菜单法
1. 数据沿用方法二的数据
2. 同方法二，输入命令（注：ln、log均可）
gen lnrr=ln(rr)
gen lnlci=ln(lci)
gen lnuci=ln(uci)
3. 调出菜单

4. Main选项卡，勾选如下图

5. Binary、Continuous选项卡不需设置；Effect选项卡，勾选如下图



6. Graph Opts选项卡，可设置坐标轴数值、字号等（勾选如下图），亦可不设置，系统会默认；If/in选项卡亦不需设置


7. 最后点击OK，跑出森林图，与方法二所示一致。
--------------------------------------------------------------------------------------------------------------------------------------
方法二：Stata metan命令



1.RR及其95%CI的合并
gen logrr=log(rr)
gen loglci=log(lci)
gen loguci=log(uci)
（1）固定效应模型（I-V）
metan logrr loglci loguci, eform [label(namevar=study, yearvar=year)] [by(group)] fixed [xlabel(.1, .5, 1)] effect(RR)
（2）随机效应模型（D-L, 异质性估计采用I-V固定效应模型）
metan logrr loglci loguci, eform [label(namevar=study, yearvar=year)] [by(group)] random [xlabel(.1, .5, 1)] effect(RR)

2.OR及其95%CI的合并
将上述metan命令中effect(RR)改为effect(OR)即可。

3.HR及其95%CI的合并
将上述metan命令中effect(RR)改为effect(HR)即可。

注：RR、OR、HR及其95%CI等均为precalculated effect estimates，合并时仅可选择固定效应模型（I-V）或随机效应模型（D-L, with I-V heterogeneity），metan命令菜单中的选项亦是如此，原因如下：

本例metan命令可写为：
metan logrr loglci loguci, eform label(namevar=study) by(group) random xlabel (.1, .5, 1) effect(RR)
森林图（与原始文献一致）：