美文网首页seurat
跟着Cell学单细胞转录组分析(十四):细胞比例柱状图---连线

跟着Cell学单细胞转录组分析(十四):细胞比例柱状图---连线

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

    之前做单细胞比例图的时候做过堆叠柱状图,跟着Cell学单细胞转录组分析(六):细胞比例计算及可视化),看那篇Cell文章封面也给我一个启发,我们做的柱状图每个样本是一个单独的柱子,如果有很多细胞类型,从左到右看起来很费劲,那么如果把相同的细胞类型用线连起来的话,看起来就方便多了,包括不同组中细胞比例变化一目了然。

    一般的做法如下:

    table(scedata$orig.ident)#查看各组细胞数
    prop.table(table(Idents(scedata)))
    table(Idents(scedata), scedata$orig.ident)#各组不同细胞群细胞数
    Cellratio <- prop.table(table(Idents(scedata), scedata$orig.ident), margin = 2)#计算各组样本不同细胞群比例
    Cellratio <- as.data.frame(Cellratio)
    
    allcolour=c("#DC143C","#0000FF","#20B2AA","#FFA500","#9370DB","#98FB98","#F08080","#1E90FF","#7CFC00","#FFFF00",
                "#808000","#FF00FF","#FA8072","#7B68EE","#9400D3","#800080","#A0522D","#D2B48C","#D2691E","#87CEEB","#40E0D0","#5F9EA0",
                "#FF1493","#0000CD","#008B8B","#FFE4B5","#8A2BE2","#228B22","#E9967A","#4682B4","#32CD32","#F0E68C","#FFFFE0","#EE82EE",
                "#FF6347","#6A5ACD","#9932CC","#8B008B","#8B4513","#DEB887")
    library(ggplot2)
    ggplot(Cellratio) + 
      geom_bar(aes(x =Var2, y= Freq, fill = Var1),stat = "identity",width = 0.7,size = 0.5,colour = '#222222')+ 
      theme_classic() +
      labs(x='Sample',y = 'Ratio')+
      scale_fill_manual(values = allcolour)+
      theme(panel.border = element_rect(fill=NA,color="black", size=0.5, linetype="solid"))
    
    

    使用连线连接:



    详细内容请至我的公众号《KS科研分享与服务》,感谢支持

    相关文章

      网友评论

        本文标题:跟着Cell学单细胞转录组分析(十四):细胞比例柱状图---连线

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