69-R可视化12-用easylabel轻松手动添加标记

参考：

easylabel (r-project.org)[1]

前言

先前我介绍过ggrepel 这个包：[[67-R可视化11-用ggrepel更加美观的添加标记（火山图的实现）]]

其实现的文本标记更加美观，且个性化程度远远优于：[[66-R可视化10-自由的在ggplot上添加文本（柱状图加计数）]]

那么我们可不可以直接指哪打哪，连代码都不用呢？

基于shiny 的easylabel 包或许是你的选择。

这里我以火山图为例。

1-准备工作

加载包、制作假数据：

代码语言：javascript复制

my_packages<- c("RColorBrewer", "paletteer", "ggplot2",
                "ggpubr", "tidyverse", "ggrepel", "devout",
                "ggforce", "easylabel")
tmp <- sapply(my_packages, function(x) library(x, character.only = T)); rm(tmp, my_packages)

# make fake data
DEG <- data.frame(
  genes = paste0("gene", 1:100),
  fold_change = round(runif(100, -10, 10), 2),
  P_value = runif(100, 0.001, 0.1)
)
DEG$group <-  ifelse((DEG$P_value > 0.05) | (abs(DEG$fold_change) < 2) , "no-Significant",
                     ifelse(DEG$fold_change > 2, "increase", "decrease"))
DEG$group <- factor(DEG$group, levels = c("decrease", "no-Significant", 
                                          "increase"))
sig_DEG <- DEG[!DEG$group %in% "no-Significant",]

2-直接操作easylabel

开始我以为，这个包会像ggthemeassistant 一样，直接操作绘图对象p，直接修改代码，更加方便且个性化。

没想到，它是需要直接以数据框操作的。

这里因为以火山图为例，先使用火山图绘图函数easyVolcano ：需要注意这里的数据框的行名需要是基因名

代码语言：javascript复制

# 开始之前修改一下数据框
DEG$new_P <- -log10(DEG$P_value)
easylabel(DEG)
rownames(DEG) <- DEG$genes

easyVolcano(DEG, x = "fold_change", 
            y = "P_value",
            fccut = 1, fdrcutoff = 0.05,
            ylim = c(0, 6), xlim = c(-5, 5),
            colScheme = colScheme,
            vline = c(-1, 1))

这里参数分别表示：