单细胞转录组画小提琴VlnPlot只显示需要类型细胞

时间:2025-02-21 07:35:38

探序基因肿瘤研究院  整理

在R语言中,单细胞转录组数据画小提琴图中,我们如果想只展示出需要的特定细胞类型,可以把这需要展示的细胞类型单独构建成一个Seurat对象,再画图。

例如我们在Seurat对象中的meta.data表格中,有一列,记录了细胞类型。needgrp变量传入要显示的细胞类型名字,为一个向量。

GetSubSeurat <- function(temp.Seurat,needgrp)
{
  a <- unique(temp.Seurat@meta.data$group)
  ChoClus <- a[-which(a %in% needgrp)]
  temp.Seurat.sub <- subset(x = temp.Seurat, subset = group %in% ChoClus)
  temp.Seurat.sub@meta.data$anno1 <- temp.Seurat.sub@meta.data$group
  return(temp.Seurat.sub)
}

注意,meta.data表格中,要有group这一列,记录细胞的分群信息。