森林图简介
什么是森林图?
参考:
https://mp.weixin.qq.com/s/vwNf_sFlmhp7DeSYaQ3NxQ
森林图是以统计指标和统计分析方法为基础,用数值运算结果绘制出的图形。它在平面直角坐标系中,以一条垂直的无效线(横坐标刻度为1或0)为中心,用平行于横轴的多条线段描述了每个被纳入研究的效应量和置信区间,用一个图形描述了多个研究合并的效应量及置信区间。
标签:#微生物组数据分析 #MicrobiomeStatPlot #森林图 #R语言可视化 #Forest plot
作者:First draft(初稿):Defeng Bai(白德凤);Proofreading(校对):Ma Chuang(马闯) and Jiani Xun(荀佳妮);Text tutorial(文字教程):Defeng Bai(白德凤)
源代码及测试数据链接:
https://github.com/YongxinLiu/MicrobiomeStatPlot/项目中目录 3.Visualization_and_interpretation/ForestPlot
或公众号后台回复“MicrobiomeStatPlot”领取
这是Jakob Stokholm课题组2023年发表于Nature Medicine上的文章,第一作者为Cristina Leal Rodríguez,题目为:The infant gut virome is associated with preschool asthma risk independently of bacteria. https://doi.org/10.1038/s41591-023-02685-x
图 6 | 病毒性哮喘特征评分与早期生活暴露的关系。
森林图显示了双侧线性回归的β系数和P值,点表示点估计值,水平线表示早期暴露与1年哮喘病毒组特征评分之间相关性的95%置信区间。风险敞口按其R2贡献进行排序。对于出生季节,基线为冬季;分娩时,基线是阴道分娩。BMI,体重指数。
结果
接下来,我们研究了早期暴露是否会影响病毒组特征评分(图6)。家中有兄弟姐妹(R2=1.64%;β=−0.26(−0.41至−0.10);P=0.001),夏季出生(相对=0.83%;β=−0.27(−0.49至−0.05);P=0.017),出生体重=0.78%;β=−0.16(−0.30至−0.02);P=0.027)和妊娠期补充鱼油(R2冬季,R2(R2=0.67%;β=−0.16(−0.32至0.01);P=0.04)均与较高的病毒组评分呈负相关。出生时只养猫(R2=0.76%;β=0.21(0.02–0.41);P=0.027)与较高的病毒组评分呈正相关。在对这些暴露进行调整后,病毒组评分仍与后期哮喘相关(OR=1.32(1.09-1.61);P=0.005)。
源代码及测试数据链接:
https://github.com/YongxinLiu/MicrobiomeStatPlot/
或公众号后台回复“MicrobiomeStatPlot”领取
软件包安装
# 基于CRAN安装R包,检测没有则安装
p_list = c("forestploter", "readxl")
for(p in p_list){if (!requireNamespace(p)){install.packages(p)}
library(p, character.only = TRUE, quietly = TRUE, warn.conflicts = FALSE)}
# 加载R包 Load the package
suppressWarnings(suppressMessages(library(forestploter)))
suppressWarnings(suppressMessages(library(readxl)))
实战
使用R语言中forestploter软件包实现
# 载入数据
# Load data
dt <- read_excel("data/forest.xlsx")
# 数据处理
# Data processing
# 组别前添加空格
# Add space before group
dt$Subgroups <- ifelse(is.na(dt$HR),
dt$Subgroups,
paste0(" ", dt$Subgroups))
# 创建一列空列,用来后面存放森林图的图形部分
# Create an empty column to store the graphic part of the forest map later.
dt$` ` <- paste(rep(" ", 20), collapse = " ")
# 正常需要在图形显示数据的文本部分
# Display the text part of the data in the graph
dt$`HR (95% CI)` <- ifelse(is.na(dt$HR), "",
sprintf("%.2f (%.2f to %.2f)",
dt$HR, dt$LowerCI, dt$UpperCI))
# 设置主题
# Set theme
custom_theme <- forest_theme(
base_size = 12,
refline_col = "#EB746A",
footnote_col = "#4F4F4F",
footnote_fontface = "italic",
title_fontface = "bold",
title_size = 14
)
# 绘制森林图
# Plot
plot <- forest(
dt[, c(1, 6, 7)],
est = dt$HR,
lower = dt$LowerCI,
upper = dt$UpperCI,
ci_column = 2,
ref_line = 2,
arrow_lab = c("Placebo Better", "Treatment Better"),
xlim = c(0, 3),
ticks_at = c(0, 1, 2, 3),
theme = custom_theme
)
# 保存图形
# Save plot
pdf("results/Forest_map_Optimized.pdf", width = 7, height = 8)
print(plot)
dev.off()
#> png
#> 2
使用此脚本,请引用下文:
Yong-Xin Liu, Lei Chen, Tengfei Ma, Xiaofang Li, Maosheng Zheng, Xin Zhou, Liang Chen, Xubo Qian, Jiao Xi, Hongye Lu, Huiluo Cao, Xiaoya Ma, Bian Bian, Pengfan Zhang, Jiqiu Wu, Ren-You Gan, Baolei Jia, Linyang Sun, Zhicheng Ju, Yunyun Gao, Tao Wen, Tong Chen. 2023. EasyAmplicon: An easy-to-use, open-source, reproducible, and community-based pipeline for amplicon data analysis in microbiome research. iMeta 2: e83. https://doi.org/10.1002/imt2.83
Copyright 2016-2024 Defeng Bai baidefeng@caas.cn, Chuang Ma 22720765@stu.ahau.edu.cn, Jiani Xun 15231572937@163.com, Yong-Xin Liu liuyongxin@caas.cn
猜你喜欢
iMeta高引文章 fastp 复杂热图 ggtree 绘图imageGP 网络iNAP
iMeta网页工具 代谢组MetOrigin 美吉云乳酸化预测DeepKla
iMeta综述 肠菌菌群 植物菌群 口腔菌群 蛋白质结构预测
10000+:菌群分析 宝宝与猫狗 梅毒狂想曲 提DNA发Nature
一文读懂:宏基因组 寄生虫益处 进化树 必备技能:提问 搜索 Endnote
16S功能预测 PICRUSt FAPROTAX Bugbase Tax4Fun
生物科普: 肠道细菌 人体上的生命 生命大跃进 细胞暗战 人体奥秘
写在后面
为鼓励读者交流快速解决科研困难,我们建立了“宏基因组”讨论群,己有国内外6000+ 科研人员加入。请添加主编微信meta-genomics带你入群,务必备注“姓名-单位-研究方向-职称/年级”。高级职称请注明身份,另有海内外微生物PI群供大佬合作交流。技术问题寻求帮助,首先阅读《如何优雅的提问》学习解决问题思路,仍未解决群内讨论,问题不私聊,帮助同行。
点击阅读原文