R语言：统计结果汇总森林图画法

Posted 2021-04-18 荟医点

图表是文章的灵魂，一张高颜值的图表可以迅速引起编辑、审稿人和读者的兴趣。在meta分析中，森林图是一种重要又很常见的元素，用来展示每个研究的结果，而本文主要为大家介绍统计结果汇总森林图。什么是统计结果汇总森林图？大家看一张已经发表在论文中（百度随便找的图不知道是哪篇文献中的）的插图（图片来源：Medsci）就明白了。

上图容纳了不同的亚组中（同一变量的不同水平）的统计结果，并以森林图形式呈现，十分直观地展示了不同亚组之间的区别。话不多说，直接开画！

──────§──────

一、数据准备

如上图是我随便“编造”出的一组结果：第一列为亚组名称，第二列为样本量，第三列为两亚组结果OR值的比较（或交互作用）的p值，后三列为OR及95%CI。

二、加载所需要的包

library(haven) #读入SPSS数据所需library(forestplot) #画森林图所需

有些包第一次使用需要先安装，安装方法为：

install.packages("包的名字")

三、读入数据

---

a <- read_sav("d:\森林图.sav")attach(a)

四、画图

---

forestplot(as.matrix(a[,1:3]), #设置用于文本展示的列，此处我们用数据的前三列作为文本，在图中展示mean=a$OA,  #设置均值lower=a$LA,  #设置均值的上限upper=a$UA,  #设置均值的下限is.summary=c(T,T,F,T,F,F,F,T,F,F,F,T,F,F),  #用于定义数据中的每一行是否是汇总值或者是否需要突出显示zero = 1,  #设置参考值boxsize = 0.7,  #设置点估计的方形的大小graphwidth=unit(70,"mm"),  #置信区间图的宽度lineheight="auto",  #设置图形中行间距，此处设置为自动colgap = unit(2,'mm'),  #设置图形中列间距lwd.zero = 2, #设置参考线的粗细lwd.ci = 2, #设置区间估计线的粗细xticks.digits=2,  #设置X轴允许存在的位数xticks=(c(0.5,1,1.5,2,2.5,3,3.5,4)), txt_gp=fpTxtGp(ticks=gpar(cex=0.9),xlab=gpar(cex=0.9),title=gpar(cex=1),cex=0.9), #设置字体col=fpColors(box='#458B00',summary="#8B008B",lines='black',zero = '#7AC5CD'),  #设置颜色，从左到右分别为点估计方形、汇总值、区间估计线、参考线xlab="OR (95%CI)",  #设置森林图横坐标坐标轴名称lwd.xaxis=2,  #设置X轴粗细lty.ci = "solid",  #设置点估计形式graph.pos =3  #设置森林图的位置，此处设置为3，则出现在第三列)

五、结果展示

如图，一张森林图就画好啦 ！代码中有些红色的是微信自动报错，可能是微信中“添加代码”这一功能基于的语言（具体什么语言我也不知道 ）与R语言某些规则不同，可以无视。大家画出来有不适合的地方多调调里面的参数就可以了，还有不明白的大家可以评论或者后台提问 。

────── § ──────

声明：本公众号推文的目的在于传播医学知识，文章仅代表作者观点。素材来源于网络，如有侵权请与我们获得联系。建议、投稿、合作相关事宜欢迎联系[e-mail]