美文网首页R语言ggplot2绘图作图
复现Nature medicine图表---堆叠柱状图显示每个样

复现Nature medicine图表---堆叠柱状图显示每个样

作者: KS科研分享与服务 | 来源:发表于2022-03-30 11:11 被阅读0次

    学习单细胞文章的时候,看到NM上一篇文章差异基因的显示方法。

    image.png

    图形如下:

    image.png

    先解释下这样展示的意义(举例可能并不是很恰当):假设我们需要研究某个细胞群分为两个亚群,想要看这两个亚群之间的差异基因变化情况,就可以使用这种方法。先计算每个样本中这两群细胞的差异基因,再将所有样本合并,找出需要关注的基因,就可以作图了。

    因为我们演示数据并没有涉及这样的意义,所以偷懒用了不同群细胞的差异基因构建了用来演示的示例数据。图形本身是柱状图,而且是堆叠柱状图。那么作图需要以下因素,首先是纵坐标的LogFC,横坐标的基因名,还有分组。具体的数据形式如下(差异基因筛选不再赘述):

    image.png

    接下来用ggplot作图即可,为了让基因排列按照我们拍好的顺序,将其转化为因子排序,这样作图显示基因名顺序就不会变了。

    
    diff <- read.csv("diff.csv",header = T,row.names = 1)
    library(ggplot2)
    library(forcats)
    diff$gene  <- as.factor(diff$gene )
    diff$gene  <- fct_inorder(diff$gene 
    

    先做一个基本的柱状图。

    
    ggplot()+
      geom_bar(data=diff,
               aes(x=avg_log2FC,y=gene,fill=cluster),
               stat="identity",
               colour = "grey50")#柱子边框颜色
    
    image.png

    接下来,改变主题,图形方向,填充颜色即可。

    
    ggplot()+
      geom_bar(data=diff,
               aes(x=avg_log2FC,y=gene,fill=cluster),
               stat="identity",
               colour = "grey50")+
      scale_fill_manual(values = c("#2271b6","#6bafd6",
                                   "#9ecbe2","#fd9272","#fee1d3"))+
      theme_bw()+
      labs(x="Cumulative avg. logFC",y="Gene expression")+
      coord_flip()+
      theme(axis.text.x=element_text(angle=90,hjust = 1,vjust=0.5))#X轴
    
    image.png

    这样就完成了,就是一摸一样有没有?解锁新的数据展示,快尝试下吧。想要示例数据的小伙伴可以到我的公众号《KS科研分享与服务》,感谢支持!

    相关文章

      网友评论

        本文标题:复现Nature medicine图表---堆叠柱状图显示每个样

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