R实战 | 森林图绘制

forestplot.jpg

R实战 | 森林图绘制

回答一下R绘图实战|GSEA富集分析图 - 知乎 (zhihu.com)里会员的问题。简单的森林绘制。

Multivariate logistic regression determined independent candidate diagnostic biomarkers.

Ref：He S, Deng Z, Li Z, et al. Signatures of 4 autophagy-related genes as diagnostic markers of MDD and their correlation with immune infiltration. J Affect Disord. 2021;295:11-20. doi:10.1016/j.jad.2021.08.005

22

森林图

森林图是一种简单直观地展示单一研究和汇总研究的图形，广泛应用于荟萃分析、观察研究和临床试验中。

Odds ratio(OR)

在病例对照研究（case-control study）中，分析暴露风险因素与疾病（或者用药）的关联程度；主要是反映暴露与疾病之间关联强度的指标，OR常适用于病例对照研究，也可以运用于前瞻性的研究（当观察时间相等时）

OR值的统计学意义：

• OR>1，暴露与疾病的危险度增加，两者呈正相关
• OR<1，暴露与疾病的危险度减少，两者呈负相关
• OR=1，暴露与疾病的危险度无关，两者呈不相关

绘制

library(survival)
library(survminer)
library(eoffice) #导出ppt用

# 构建模型
# 以colon数据为示例数据
# 时间变量是time，结局变量是status，自变量选择 sex，rx和 adhere
model <- coxph( Surv(time, status) ~ sex + rx + adhere,
                data = colon )
# 森林图
ggforest(model,  #coxph得到的Cox回归结果
         data = colon,  #数据集
         main = 'Hazard ratio of colon',  #标题
         cpositions = c(0.05, 0.15, 0.35),  #前三列距离
         fontsize = 1, #字体大小
         refLabel = 'reference', #相对变量的数值标签，也可改为1
         noDigits = 3 #保留HR值以及95%CI的小数位数
)

Snipaste_2022-03-04_22-16-36

参考

GitHub - cardiomoon/autoReg: R package autoReg

木舟笔记2022年度VIP企划

关于木舟笔记2022年度VIP会员企划
https://mp.weixin.qq.com/s/cTf38LqiV07ILX3PeVEQFA

往期内容

1. 跟着Nature学作图 | 配对哑铃图+分组拟合曲线+分类变量热图
2. (免费教程+代码领取)|跟着Cell学作图系列合集
3. 跟着Nat Commun学作图 | 1.批量箱线图+散点+差异分析
4. 跟着Nat Commun学作图 | 2.时间线图
5. 跟着Nat Commun学作图 | 3.物种丰度堆积柱状图
6. 跟着Nat Commun学作图 | 4.配对箱线图+差异分析
7. 跟着 Nat Med. 学作图 | GSVA+limma差异通路分析+发散条形图

---