R-forestplot包HR结果绘制森林图
本⽂⾸发于“⽣信补给站”微信公众号
更多关于R语⾔，ggplot2绘图，⽣信分析的内容，敬请关注⼩号，给您⼲货。
上⼀篇简单的介绍了COX⽣存分析结果绘制森林图Forest plot（森林图） | Cox⽣存分析可视化，本⽂将介绍根据数据集合的基本信息以及点估计值（置信区间区间）的结果直接绘制森林图的⽅法。
其中点估计值（置信区间）的结果可以是COX也可以是logistic回归等其他⽅法的结果，适⽤范围更⼴。
⼀准备数据
#载⼊R包
library(forestplot)
#数据来源：/forest-plot-with-horizontal-bands/
data <- read.csv("ForestPlotData.csv", stringsAsFactors=FALSE)
#查看数据
head(data)
⼆绘制森林图
广州华尔街英语学费2.1 简单森林图
对数据进⾏部分修改，⽅便⾏名和列名字输出
## 构建tabletext，更改列名称，展⽰更多信息
bee是什么意思np <- ifel(!is.na(data$Count), paste(data$Count," (",data$Percent,")",p=""), NA)
## The rest of the columns in the table.
tabletext <- cbind(c("Subgroup","\n",data$Variable),
c("No. of Patients (%)","\n",np),
新东方刘韧c("4-Yr Cum. Event Rate\n PCI","\n",data$PCI.Group),
北京十二中初中一部c("4-Yr Cum. Event Rate\n Medical Therapy","\n",data$Medical.Therapy.Group),
c("P Value","\n",data$P.Value))
##绘制森林图
forestplot(labeltext=tabletext, graph.pos=3,
备忘录英文mean=c(NA,NA,data$Point.Estimate),
lower=c(NA,NA,data$Low), upper=c(NA,NA,data$High),
boxsize=0.5)
如上图所⽰基本信息OK了，但是可以在以下⼏个⽅⾯进⾏优化：
district
添加线条，区分Subgroup期待你
更改箱线图的宽度，颜⾊和⼤⼩
更改字体⼤⼩，更易区分
济南会计培训
添加标题和横坐标轴标⽰
2.2 优化森林图
## 定义亚组，⽅便后⾯线条区分
subgps <- c(4,5,8,9,12,13,16,17,20,21,24,25,28,29,32,33)
data$Variable[subgps] <- paste(" ",data$Variable[subgps])
forestplot(labeltext=tabletext,
graph.pos=3, #为Pvalue箱线图所在的位置
mean=c(NA,NA,data$Point.Estimate),
lower=c(NA,NA,data$Low), upper=c(NA,NA,data$High),
#定义标题
title="Hazard Ratio Plot",
##定义x轴
xlab=" <---PCI Better--- ---Medical Therapy Better--->",
##根据亚组的位置，设置线型，宽度造成“区块感”
hrzl_lines=list("3" = gpar(lwd=1, col="#99999922"),
"7" = gpar(lwd=60, lineend="butt", columns=c(2:6), col="#99999922"),
"15" = gpar(lwd=60, lineend="butt", columns=c(2:6), col="#99999922"),
"23" = gpar(lwd=60, lineend="butt", columns=c(2:6), col="#99999922"),
"31" = gpar(lwd=60, lineend="butt", columns=c(2:6), col="#99999922")),
#fpTxtGp函数中的cex参数设置各个组件的⼤⼩
txt_gp=fpTxtGp(label=gpar(cex=1.25),
ticks=gpar(cex=1.1),
xlab=gpar(cex = 1.2),
title=gpar(cex = 1.2)),
##fpColors函数设置颜⾊
zolocol=fpColors(box="#1c61b6", lines="#1c61b6", zero = "gray50"),
#箱线图中基准线的位置
zero=1,
cex=0.9, lineheight = "auto",
colgap=unit(8,"mm"),
#箱⼦⼤⼩，线的宽度
lwd.ci=2, boxsize=0.5,
#箱线图两端添加⼩竖线，⾼度
ci.vertices=TRUE, ci.vertices.height = 0.4)
如此即绘制完成了，颜⾊，⼤⼩，间隔等需要根据实际情况进⾏调整。更多：
R|⽣存分析（1）：⽣存分析介绍以及绘制KM曲线
Nomogram（诺莫图） | Logistic、Cox⽣存分析结果可视化
Forest plot（森林图） | Cox⽣存分析可视化
maftools|TCGA肿瘤突变数据的汇总，分析和可视化
maftools | 从头开始绘制发表级oncoplot（瀑布图）
ggalluvial|炫酷桑基图（Sankey），你也可以秀
ggplot2|详解⼋⼤基本绘图要素
ggplot2|ggpubr进⾏“paper”组图合并
pheatmap|暴⾬暂歇，“热图”来袭
ggplot2-plotly|让你的⽕⼭图“活”过来
ggplot2| 绘制KEGG⽓泡图
ggplot2|绘制GO富集柱形图
mylf是什么意思
绘图系列|R-corrplot相关图
绘图系列|R-VennDiagram包绘制韦恩图
R|clusterProfiler-富集分析
【觉得不错，右下⾓点个“在看”，期待您的转发，谢谢！】
