在医咖会既往推送的一篇文章中,向大家介绍了如何通过GraphPad Prism和Excel软件来绘制森林图,从而使得回归分析的结果能够可视化。在本期内容中,再来介绍两款常用软件——R和Stata,教大家进一步玩转森林图。
以2016年发表在JACC杂志上的这篇文章《A Prospective Natural History Studyof Coronary Atherosclerosis Using Fractional Flow Reserve》为例,利用文章中的Table 3多因素回归结果来绘制森林图。
R作为一个90后,在数据处理、统计分析、计算绘图等方面,俨然已经成为了一颗新星。尤其是R在绘图方面的突出表现,凭借强大的程序资源包,以及漂亮、高质量的图形输出,获得了一大片死忠粉的拥护。
下面我们就来介绍一下,如何通过R简单的4行代码,来绘制回归分析结果的森林图。
1. 安装forsetplot程序包(绘制森林图)和haven程序包(导入SPSS文件)
install.packages("forestplot")
install.packages("haven")
2. 加载forsetplot和haven程序包
library(forestplot)
library(haven)
3. 准备数据文件并导入R
3.1 数据录入
本文以SPSS软件数据录入为例,录入格式如下:
3.2 数据导入
#用read_sav函数导入SPSS文件并命名为ForestPlot
ForestPlot <- read_sav("C:/数据库存储路径/Forest Plot.sav")
#用attach函数将其添加到R的搜索路径中,作为当前默认数据框
确认删除