【Stata】亚组分析森林图：汇总所有亚组分析的结果

小伙伴们平时在阅读文献的时候可能会见到过类似下图的亚组分析森林图，是不是觉得特别高大上？所有的亚组分析结果一目了然。其实这样的森林图制作起来并不是十分困难。今天就向大家介绍两种方法：

方法一

Excel 作图

有兴趣的小伙伴可以下载这篇文献：Meta-analyses and Forest plots using a microsoft excel spreadsheet: step-by-step guide focusing on descriptive data analysis，学习如何用Excel绘制图形。

当然，身为meta分析的搬运工，“尔云间meta分析”当然会给大家福利，关注我们的公众号，并回复“excel绘制森林图”就能获得文献和Excel模版的下载链接。

只需要在模版中修改数据，就能绘制个性化森林图，十分简单快捷（在此感谢战友“xiaoqiangge”提供方法的介绍和Excel模版）。

方法二

Stata作图

作为meta分析的爱好者，Stata是必须掌握的统计软件，它的功能十分强大，当然也可以制作高大上的森林图。

普通的森林图很多人都会做，然而遇到文中的森林图就没辙了，关键是不懂怎么录入数据。

其实，数据的录入也没有什么特别，看看下图就知道了。

关键在于“Subgroup”和“group”这两列变量，“group”是亚组的分类变量，“Subgroup”就是该亚组的名称。

然后进行HR和95%CI的对数转换，转换后运行最后的命令，即可得到亚组分析森林图的雏形，最后通过PS或其他修图软件对图形进行完善和美化，一切都OK了！

Stata的相关命令和森林图雏形如下（PS也很简单，这里就不多说了）：

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)