美文网首页
KEGG分析柱形图做法

KEGG分析柱形图做法

作者: RaoZC | 来源:发表于2020-06-07 22:21 被阅读0次

    该柱形图纵坐标是路径名,横坐标是基因数量,颜色梯度为-log2(q-value)。代码如下:

    setwd(".../06. WGCNA/turquoise")
    a = read.table("05. KEGG-R.txt",header = T,sep="\t")
    head(a)
    library(ggplot2)
    ggplot(a,aes(y=reorder(Pathways,Order),x=GeneNumbers,fill=X.LOG2.q.value.))+
      geom_bar(stat = "identity")+
      scale_fill_gradient(a$X.LOG2.q.value.,low = "blue",high = "red")+
      labs(title = "Pathways Enrichment",
           x = "Gene numbers", 
           y = "Pathways",
           fill = "-Log2(q-value)")+
      theme(axis.title.x = element_text(face = "bold",size = 16),
            axis.title.y = element_text(face = "bold",size = 14),
            legend.title = element_text(face = "bold",size = 14))+
      theme_get()
    

    得到的结果如下:


    KEGG柱形图

    我就没搞清楚那个legend的title怎么回事。待搞清楚了再补上。
    得到的结果在AI上面修改,添加分类,可以得到如下结果:


    AI修改

    相关文章

      网友评论

          本文标题:KEGG分析柱形图做法

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