# 安装包
if (!requireNamespace("data.table", quietly = TRUE)) {
install.packages("data.table")
}
if (!requireNamespace("jsonlite", quietly = TRUE)) {
install.packages("jsonlite")
}
if (!requireNamespace("GOplot", quietly = TRUE)) {
install.packages("GOplot")
}
# 加载包
library(data.table)
library(jsonlite)
library(GOplot)GO 条形图
注记
Hiplot 网站
本页面为 Hiplot GOBar Plot 插件的源码版本教程,您也可以使用 Hiplot 网站实现无代码绘图,更多信息请查看以下链接:
GO 条形图用于显示按 z-score 或 adjusted p-value 的负对数排序的 Z-score 彩色条形图。
环境配置
系统: Cross-platform (Linux/MacOS/Windows)
编程语言: R
依赖包:
data.table;jsonlite;GOplot
sessioninfo::session_info("attached")─ Session info ───────────────────────────────────────────────────────────────
setting value
version R version 4.5.2 (2025-10-31)
os Ubuntu 24.04.3 LTS
system x86_64, linux-gnu
ui X11
language (EN)
collate C.UTF-8
ctype C.UTF-8
tz UTC
date 2026-01-28
pandoc 3.1.3 @ /usr/bin/ (via rmarkdown)
quarto 1.8.27 @ /usr/local/bin/quarto
─ Packages ───────────────────────────────────────────────────────────────────
package * version date (UTC) lib source
data.table * 1.18.0 2025-12-24 [1] RSPM
ggdendro * 0.2.0 2024-02-23 [1] RSPM
ggplot2 * 4.0.1 2025-11-14 [1] RSPM
GOplot * 1.0.2 2016-03-30 [1] RSPM
gridExtra * 2.3 2017-09-09 [1] RSPM
jsonlite * 2.0.0 2025-03-27 [1] RSPM
RColorBrewer * 1.1-3 2022-04-03 [1] RSPM
[1] /home/runner/work/_temp/Library
[2] /opt/R/4.5.2/lib/R/site-library
[3] /opt/R/4.5.2/lib/R/library
* ── Packages attached to the search path.
──────────────────────────────────────────────────────────────────────────────
数据准备
加载的数据是具有七个列的GO富集的结果:category, GO id, GO term, gene count, gene name, logFC, adjust pvalue and zscore。
# 加载数据
data <- data.table::fread(jsonlite::read_json("https://hiplot.cn/ui/basic/gobar/data.json")$exampleData$textarea[[1]])
data <- as.data.frame(data)
# 整理数据格式
colnames(data) <- c("category","ID","term","count","genes","logFC","adj_pval","zscore")
data <- data[data$category %in% c("BP","CC","MF"),]
data <- data[!is.na(data$adj_pval),]
data$adj_pval <- as.numeric(data$adj_pval)
data$zscore <- as.numeric(data$zscore)
# 查看数据
head(data) category ID term count genes logFC adj_pval
1 BP GO:0007507 heart development 54 DLC1 -0.9707875 2.17e-06
2 BP GO:0007507 heart development 54 NRP2 -1.5153173 2.17e-06
3 BP GO:0007507 heart development 54 NRP1 -1.1412315 2.17e-06
4 BP GO:0007507 heart development 54 EDN1 1.3813006 2.17e-06
5 BP GO:0007507 heart development 54 PDLIM3 -0.8876939 2.17e-06
6 BP GO:0007507 heart development 54 GJA1 -0.8179480 2.17e-06
zscore
1 -0.8164966
2 -0.8164966
3 -0.8164966
4 -0.8164966
5 -0.8164966
6 -0.8164966
可视化
# GO 条形图
p <- GOBar(data, display = "multiple", order.by.zscore = T,
title = "GO Enrichment Barplot ",
zsc.col = c("#EF8A62","#F7F7F7","#67A9CF")) +
theme(plot.title = element_text(hjust = 0.5),
axis.text.x = element_text(size = 8))
p
如示例图所示,x 轴表示每个 GO term,y 轴表示 -log(adj_pvalue),每个条按 z-score 着色。如果显示分面,则面板的宽度将与 x 刻度的长度成正比。
