我是用seraut自帶函數(shù)作圖的，當然很多人用ggplot做的更好看，再次提醒有人會用ggplot調(diào)好色的代碼，記得簡書發(fā)給我，我就用很蠢的自帶函數(shù)作圖。

VlnPlot(experiment.aggregate, features = "SDC1",group.by = "celltype3")

圖片.png

選取大于0的

VlnPlot(subset(experiment.aggregate,SDC1>0), features = "SDC1",group.by = "celltype3")

圖片.png

疊加

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()

圖片.png

換位置

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()

圖片.png