# 安装包
if (!requireNamespace("ggalluvial", quietly = TRUE)) {
install.packages("ggalluvial")
}if (!requireNamespace("ggplot2", quietly = TRUE)) {
install.packages("ggplot2")
}
# 加载包
library(ggalluvial)
library(ggplot2)
矩阵气泡图
注记
Hiplot 网站
本页面为 Hiplot Matrix Bubble
插件的源码版本教程,您也可以使用 Hiplot 网站实现无代码绘图,更多信息请查看以下链接:
彩色矩阵气泡用于可视化不同细胞(列)中多个基因(行)的表达矩阵数据。
环境配置
系统: Cross-platform (Linux/MacOS/Windows)
编程语言: R
依赖包:
ggalluvial
;ggplot2
数据准备
数据结构:
- <1st-col>: (String) 单元格样本名称作为 X 轴
- <2nd-col>: (String) 基因名作为Y轴
- <3rd-col>: (Numeric) 基因表达
- <4rd-col>: (Numeric) 组
# 加载数据
<- read.delim("files/Hiplot/118-matrix-bubble-data.txt", header = T)
data
# 整理数据格式
1] <- factor(data[, 1], levels = unique(data[, 1]))
data[, 2] <- factor(data[, 2], levels = unique(data[, 2]))
data[,
# 查看数据
head(data)
x y value group
1 cell1 gene1 1 groupC
2 cell1 gene2 5 groupC
3 cell1 gene3 10 groupC
4 cell1 gene4 15 groupC
5 cell1 gene5 20 groupC
6 cell1 gene6 25 groupC
可视化
# 矩阵气泡图
<- ggplot(data = data, aes(x = x, y = y, size = value, color = y)) +
p geom_point(alpha = 1) +
labs(title = "Matrix Bubble") +
guides(color = FALSE) +
theme(panel.background = element_blank(),
panel.grid.major = element_line(colour = "gray"),
strip.background = element_blank(),
panel.border = element_rect(colour = "black", fill = NA),
panel.spacing = unit(0, "lines"),
plot.title = element_text(size = 12, hjust = 0.5),
text = element_text(family = "Arial"),
legend.title = element_text(size = 10),
axis.text.x = element_text(angle=0, hjust=0.5, vjust=1)) +
facet_grid(~group, scales = 'fixed', margins = F) +
scale_color_manual(values = c(
"#3B4992FF","#EE0000FF","#008B45FF","#631879FF","#008280FF","#BB0021FF",
"#5F559BFF","#A20056FF","#808180FF","#1B1919FF"))
p
