数据和代码获取:请查看主页个人信息!!!
大家好,今天我将介绍高分杂志组合图绘制,图形来源于Nature的一篇文章:
该图展示了转录组多样性的主成分(PCs)与基因组(黑色标签)和临床(蓝色标签)变量之间的关系。原文使用混合效应线性模型进行显著性的确定,感兴趣的同学可以进行深入学习。
本文主要对可视化方法进行展示,话不多说,我们直接开始:
Step1:数据载入
rm(list=ls())
pacman::p_load(tidyverse,aplot,ggtext)
data <- read.csv('data.csv', row.names = 1)
大家可以看到:原文图中有三组Y文本非黑色字体,这个并非P图实现,而是通过【ggtext包】实现,因此我们在绘图前需要对数据分组颜色进行赋值。
Step2:定义Y轴字体颜色
d <-
data %>%
filter(variable_type %in% c("genomic", "clinical")) %>%
mutate(axis_cols = if_else(variable_type == "clinical",
true = "#3399FF",
false = "