用seurat自带函数画
https://www.jianshu.com/p/e879cb6dfd26
Idents(experiment.aggregate) <- "celltype3"markers <- c('IGHM','FCER2','IL4R','CD27','SUGCT','LRMP','SDC1','CD38','JCHAIN')markers <- CaseMatch(markers, rownames(experiment.aggregate))B2 <- VlnPlot(experiment.aggregate, features = markers, pt.size = 0, stack = T)+ggsci::scale_fill_jco()
用factor控制group_list的顺序
Idents(experiment.aggregate) <- "celltype3"
markers <- c('IGHM','FCER2','IL4R','CD27','SUGCT','LRMP','CD38','SDC1','JCHAIN')
markers <- CaseMatch(markers, rownames(experiment.aggregate))
group_list <- factor(experiment.aggregate@meta.data$celltype3,levels = c('IGG+ B cells','IGA+ B cells','GC B cells','Memory B cells','Naive B cells'))
experiment.aggregate@meta.data$celltype4 <- group_list
B2 <- VlnPlot(experiment.aggregate, features = markers, pt.size = 0, stack = T,group.by = 'celltype4')+
ggsci::scale_fill_jco()
