2024-01-16 | 个性化的富集分析图

标新立异而已,打破常规

library(ggplot2)
library(openxlsx)
library(RColorBrewer)
library(ggthemes)

kegginput <- read.xlsx("kegg.xlsx")
x=kegginput$logp
y=factor(kegginput$Term, levels=kegginput$Term)

p1 = ggplot(kegginput,aes(x=x,y=y))+ 
     geom_point(aes(size=Count,color=-0.5*log(pvalue)))+
     geom_text(aes(x=0, y=y,label=y, color=-0.5*log(pvalue)), 
               hjust=(0))+
     scale_color_gradient(low = "BLUE", high = "OrangeRed")  + 
     theme_few()+
     labs(size="Count", color=expression(-log[10](pvalue)), title="KEGG pathway", x=expression(-log[10](pvalue)), y="")+
     theme(axis.text.x=element_text(face="bold"),
           axis.ticks = element_blank(),
           axis.text.y=element_blank())
p1
图片2.png
library(ggplot2)
library(openxlsx)
library(RColorBrewer)
library(ggthemes)

goinput <- read.xlsx("go.xlsx")
x1=goinput$logp
y1=factor(goinput$Term, levels=goinput$Term)

p2 = ggplot(goinput, aes(x=x1,y=y1))+ 
     geom_bar(stat="identity",width = 0.4,position = position_dodge(0.7),
              aes(fill=-0.5*log(pvalue)))+
     geom_point(aes(size=Count,color=-0.5*log(pvalue)), color="gray")+
     geom_text(aes(x=0, y=y1, label=y1), color="black",
               hjust=(0))+
     scale_color_gradient(low="SpringGreen", high="OrangeRed")+ 
     scale_fill_gradient(low="SpringGreen", high="OrangeRed")+ 
     theme_few()+ 
     labs(size="Count", fill=expression(-log[10](pvalue)), title="GO terms", x=expression(-log[10](pvalue)), y="")+
     theme(axis.text.x = element_text(face="bold"),
           axis.ticks = element_blank(),
           axis.text.y = element_blank())
p2
图片1.png
library(ggplot2)
library(openxlsx)
library(RColorBrewer)
library(ggthemes)
library(patchwork)
goinput <- read.xlsx("KEGG.xlsx")
x3=goinput$logp
y3=factor(goinput$Term, levels=goinput$Term)

p3 = ggplot(goinput, aes(x=x3,y=y3))+ 
  geom_bar(stat="identity",width = 0.5,position = position_dodge(0.7),
           aes(fill=-0.5*log(pvalue)), color="black")+
  geom_point(aes(size=Count,color=-0.5*log(pvalue)), color="gray")+
  coord_flip()+
  scale_color_gradient(low =  "#00BFC4", high = "#F8766D")+ 
  scale_fill_gradient(low =  "#00BFC4", high = "#FCA311")+ 
  theme_bw()+ 
  labs(size="Count", fill=expression(-log[10](pvalue)), 
       title="KEGG pathway", x=expression(-log[10](pvalue)), y="")+
  theme(axis.text.y = element_text(size=10, face="bold", color="black"),
        axis.ticks = element_blank(),
        axis.text.x = element_text(size=10, face="bold", color="black", 
                                   angle=90,  hjust = 1, vjust = 1))
p3
2.3-KEGG富集.png
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容