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)