方法一: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)
森林图(与原始文献一致):
|