2. ggplot2画细胞亚群的比例和条形图

2.1 按照细胞celltype分群进行绘制细胞亚群ratio的百分比图

prop.table(table(sce$celltype))
table(sce$celltype, sce$group)#各组不同细胞群细胞数
Cellratio <- as.data.frame(prop.table(table(Idents(sce), sce$group), margin = 2))#计算各组样本不同细胞群比例
Cellratio
colnames(Cellratio) <- c("Celltype","Group","Ratio")
colourCount = length(unique(Cellratio$Celltype))
Cellratio$Group <-factor(Cellratio$Group,ordered=TRUE,levels=c("Sham","LHT")) #修改因子水平 
library(ggplot2)
p1 <- ggplot(Cellratio) + 
      geom_bar(aes(x =Ratio, y= Group, fill = Celltype),stat = "identity",width = 0.7,size = 0.3,colour = '#222222')+ 
      theme_classic() +
      labs(x='Ratio',y = 'Group')+
      coord_flip()+
      theme(panel.border = element_rect(fill=NA,color="black", size=0.5, linetype="solid"),
            text = element_text(size = 16), #修改所有的字体大小
            axis.line = element_line(size = 0), #修改x,y轴的线的宽度
            axis.text.x=element_text(size=16,colour = "black",angle = 45,hjust = 1,vjust = 1),#修改x轴的字体大小和颜色位置 
            axis.text.y=element_text(size=16,colour = "black"),#修改y轴的字体颜色和大小
            legend.title = element_text(size = 16), #修改图例的title的字体大小
            legend.text = element_text(size = 16)) #修改图例的字体的大学
                                    
p1
# ggsave(filename = "cellcluster_bar_percent.pdf",width = 3, height =6)

2.2 按照细胞Celltype分群进行ratio绘制条形图

p2 <- ggplot(Cellratio,aes(x=Celltype,y=Ratio,fill=Group))+
  geom_bar(stat='identity',
           width=0.5,
           colour='black',
           position='dodge') + 
  theme(panel.background = element_rect('white'), #改变背景色
        # panel.border = element_rect(fill = NA,colour = 'black',size = 1,linetype = 'solid'), #改变边框颜色
        axis.line = element_line(size = 1), #改变x,y轴线的宽度
        text = element_text(size = 16), #改变图中所有字体的大小
        axis.text.x = element_text(size = 16,angle = 45, hjust = 1, vjust = 1,colour = 'black'), #改变坐标轴上标注的字体大小
        axis.text.y = element_text(size = 16,colour = 'black'),
        legend.title = element_text(size = 16),
        legend.text = element_text(size = 16))
p2

2.3 按照细胞Celltype分群进行绘绘制细胞数量的条形图

Cellnumbers <- as.data.frame(table(Idents(sce), sce$group))
colnames(Cellnumbers) <- c("Celltype","Group","Numbers")
Cellnumbers$Group <-factor(Cellnumbers$Group,ordered=TRUE,levels=c("Sham","LHT")) #修改因子水平
p3 <- ggplot(Cellnumbers,aes(x=Celltype,y=Numbers,fill=Group))+
  geom_bar(stat='identity',
           width=0.5,
           colour='black',
           position='dodge') + 
  theme(panel.background = element_rect('white'), #改变背景色
        # panel.border = element_rect(fill = NA,colour = 'black',size = 1,linetype = 'solid'), #改变边框颜色
        axis.line = element_line(size = 1), #改变x,y轴线的宽度
        text = element_text(size = 16), #改变图中所有字体的大小
        axis.text.x = element_text(size = 16,angle = 45, hjust = 1, vjust = 1,colour = 'black'), #改变坐标轴上标注的字体大小
        axis.text.y = element_text(size = 16,colour = 'black'),
        legend.title = element_text(size = 16),
        legend.text = element_text(size = 16))
p3 
# ggsave(filename = "cellcluster_numbers_bar.pdf",width = 10, height =5,units = "cm")
plot <- p1+ p2 + p3
plot
library(ggpubr)
ggarrange(p1,p2,p3,nrow = 3)
ggsave(filename = "celltype_numbers_ratio_bar.pdf",width = 8, height =12)
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 218,284评论 6 506
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,115评论 3 395
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,614评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,671评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,699评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,562评论 1 305
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,309评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,223评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,668评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,859评论 3 336
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,981评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,705评论 5 347
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,310评论 3 330
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,904评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,023评论 1 270
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,146评论 3 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,933评论 2 355

推荐阅读更多精彩内容