森林图指的是对不同研究提出同一研究问题,用一系列轴来呈现转换为同一统计量的研究结果。除了meta 分析外,森林图还广泛应用于观察性研究和临床试验,例如风险分析/生存分析等结果。最大的优势是简单直观地展示单一研究和汇总研究。今天小编就用forestplot包绘制森林图,仅供参考哦~
代码如下:
1.安装并调用forestplot包
BiocManager::install("forestplot")
library(forestplot)
2.生成模拟数据,共5个模型
g1 <- rnorm(100, 8, 4)
g2 <- rnorm(100, 6, 3)
g3 <- rnorm(100, 7, 3.5)
g4 <- rnorm(100, 9, 4.5)
g5 <- rnorm(100, 5, 2.5)
Data <- data.frame(g1=g1,g2=g2,g3=g3,g4=g4,g5=g5)
head(Data)
3.计算各个模型之间的15%分位数,中位数,85%分位数,同时保留两位有效数字
Data_difference<-apply(Data(round(quantile(x,c(0.15,0.5,0.85)),digits=2))
4.设置行名