之前做單細(xì)胞比例圖的時(shí)候做過堆疊柱狀圖，（跟著Cell學(xué)單細(xì)胞轉(zhuǎn)錄組分析(六):細(xì)胞比例計(jì)算及可視化），看那篇Cell文章封面也給我一個(gè)啟發(fā)，我們做的柱狀圖每個(gè)樣本是一個(gè)單獨(dú)的柱子，如果有很多細(xì)胞類型，從左到右看起來(lái)很費(fèi)勁，那么如果把相同的細(xì)胞類型用線連起來(lái)的話，看起來(lái)就方便多了，包括不同組中細(xì)胞比例變化一目了然。

一般的做法如下：

table(scedata$orig.ident)#查看各組細(xì)胞數(shù)
prop.table(table(Idents(scedata)))
table(Idents(scedata), scedata$orig.ident)#各組不同細(xì)胞群細(xì)胞數(shù)
Cellratio <- prop.table(table(Idents(scedata), scedata$orig.ident), margin = 2)#計(jì)算各組樣本不同細(xì)胞群比例
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"))

使用連線連接：

詳細(xì)內(nèi)容請(qǐng)至我的公眾號(hào)《KS科研分享與服務(wù)》，感謝支持