大家看看这个森林图有何不同?这个森林图不是meta分析的,而是logistic回归的,高大上啊!若是咱们用此森林图展示logistic回归的结果,妙哉啊。但是,怎么画呢?今天咱们用R软件画一画logistic回归的森林图。
首先,按下图将logistic回归的结果整理好。v1为各变量名称,v2为OR(95%CI),v3为OR点估计,v4为OR的95%CI下限,v5为OR的95%CI上限。存储为csv格式,待会有大用哦。
接下来就要画森林图了,大家准备好了吗?
1. 下载、安装R和RStudio,百度搜索即可。
2. 安装forestplot程序包,install.packages("forestplot")。
3. 加载forestplot程序包,library(forestplot)。
4. 读入数据,dat<-read.csv(file.choose())。键入该命令后,会弹出浏览、选择文件的对话框,选择整理好的数据。
5. 绘制森林图
6. 运行程序后,RStudio右下角就出现我们想要的森林图了。点击Zoom可在新窗口查看森林图,通过Export可将森林图存储成适合的格式。
完整版程序如下:
install.packages("forestplot")
library(forestplot)
dat<-read.csv(file.choose()) forestplot(as.matrix(dat[,1:2]), dat$v3, dat$v4, dat$v5, zero = 1, clip = c(0.4,15), graph.pos = 2, xticks = c(0.1,1,10,100), boxsize=0.05, xlog=TRUE)