【Stata】亞組分析森林圖:匯總所有亞組分析的結(jié)果

小伙伴們平時(shí)在閱讀文獻(xiàn)的時(shí)候可能會(huì)見到過類似下圖的亞組分析森林圖,是不是覺得特別高大上?所有的亞組分析結(jié)果一目了然。其實(shí)這樣的森林圖制作起來并不是十分困難。今天就向大家介紹兩種方法:

方法一
Excel 作圖
有興趣的小伙伴可以下載這篇文獻(xiàn):Meta-analyses and Forest plots using a microsoft excel spreadsheet: step-by-step guide focusing on descriptive data analysis,學(xué)習(xí)如何用Excel繪制圖形。
當(dāng)然,身為meta分析的搬運(yùn)工,“爾云間meta分析”當(dāng)然會(huì)給大家福利,關(guān)注我們的公眾號(hào),并回復(fù)“excel繪制森林圖”就能獲得文獻(xiàn)和Excel模版的下載鏈接。
只需要在模版中修改數(shù)據(jù),就能繪制個(gè)性化森林圖,十分簡(jiǎn)單快捷(在此感謝戰(zhàn)友“xiaoqiangge”提供方法的介紹和Excel模版)。

方法二
Stata作圖
作為meta分析的愛好者,Stata是必須掌握的統(tǒng)計(jì)軟件,它的功能十分強(qiáng)大,當(dāng)然也可以制作高大上的森林圖。
普通的森林圖很多人都會(huì)做,然而遇到文中的森林圖就沒轍了,關(guān)鍵是不懂怎么錄入數(shù)據(jù)。
其實(shí),數(shù)據(jù)的錄入也沒有什么特別,看看下圖就知道了。

關(guān)鍵在于“Subgroup”和“group”這兩列變量,“group”是亞組的分類變量,“Subgroup”就是該亞組的名稱。
然后進(jìn)行HR和95%CI的對(duì)數(shù)轉(zhuǎn)換,轉(zhuǎn)換后運(yùn)行最后的命令,即可得到亞組分析森林圖的雛形,最后通過PS或其他修圖軟件對(duì)圖形進(jìn)行完善和美化,一切都OK了!
Stata的相關(guān)命令和森林圖雛形如下(PS也很簡(jiǎn)單,這里就不多說了):
gen loghr=log(hr)
gen logll=log(ll)
gen logul=log(ul)
metan loghr logll logul, lcols ( Subgroup n) by(group) nosubgroup fixed eform nooverall effect("Hazard Ratio") xlabel(0.5, 1.0, 1.5, 2.0, 3.0) texts(120)


