美文网首页绘图科研绘图
2024-01-16 | 个性化的富集分析图

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

作者: 千万别加香菜 | 来源:发表于2024-01-15 09:39 被阅读0次

    标新立异而已,打破常规

    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

    相关文章

      网友评论

        本文标题:2024-01-16 | 个性化的富集分析图

        本文链接:https://www.haomeiwen.com/subject/zfxhodtx.html