GO富集分析bar图

最新推荐文章于 2024-03-29 15:10:38 发布

cujiangcao2

最新推荐文章于 2024-03-29 15:10:38 发布

阅读量406

点赞数

文章标签： python 开发语言

本文链接：https://blog.csdn.net/cujiangcao2/article/details/126246979

版权

rm(list = ls())
library(tidyverse)
df<-read.delim('CHN_LD.USA_LD.CV_LD.same.deg_upgene_goterm',sep='\t',header = T)
head(df)
#GO_term_order=factor(as.integer(rownames(df)),labels=df$GO_Name)
GO_term_order=factor(as.integer(rownames(df)),labels=df$GO_Name)
ggplot(df,aes(GO_term_order,-log10(df$corrected.p.value.BH.method.),fill=Class))+geom_bar(stat="identity", width=0.8)+
  theme_classic()+
  scale_fill_manual(values=c('Biological process'="#8DA1CB", 'Cellular component'="#66C3A5",'Molecular function'= "#FD8D62"))+
  theme(axis.text.x=element_text(size=20, color="gray50",angle = 70,vjust = 1, hjust = 1 ),
        axis.text.y=element_text(size=24),
        axis.title.y = element_text(size=24),
        axis.title.x = element_blank(),
        legend.title = element_blank(),
        legend.text = element_text(size = 15),
        plot.margin=unit(c(0.2,0.2,0.2,1),'inches')  #上右下左
        )+
  ylab("-log(10P-value)")