带标签的火山图

先安装包，ggplot2和ggrepel，可以使用geom_text_repel()或者geom_label_repel()命令行将重叠的基因名分开。

个人认为geom_text_repel()做出来的好看一些。

volcano plot (with labels)

getwd()

###加载包

library(xlsx)

library(ggplot2)

library(ggrepel)

###读入数据

data <- read.xlsx("All_Differentially_Expressed_tRF.xlsx", stringsAsFactors = F,

                  sheetName = "AD_vs_Sham", header = T)

head(data)

###初始图

figure1 <- ggplot(data,aes(log2FC,-1*log10(p_value), label = labels)) +

   geom_point(aes(color =significant)) +

   scale_color_manual(values =c("green","black", "red")) +

geom_text_repel()

###成品图

figure <- figure1 + xlim(-2.5,2.5)+ ylim(0,2.5)+ labs(x=expression(log2(FoldChange)),y=expression(-log10(pvalue))) +

   geom_hline(yintercept=1.3, linetype = "dashed")+geom_vline(xintercept=c(-1,1), linetype = "dashed") +

   theme_bw()

###查看

figure