单细胞转录组分群中,分群信息存储在Seurat数据结构的active.ident变量中。
如果要用自己表格中记录的细胞分群编号替换掉active.ident中的,可以:
cells1 <- names(xxx.Seurat@active.ident)
cells2 <- rownames(xxxTab)
先检查cells1和cells2细胞名是否都一致
替换操作:
xxx.Seurat@active.ident[cells2] <- factor(xxxTabcluster,levels = sort(unique(xxxTabcluster)))
然后使用DimPlot画图,可以看到分群信息已经变成自己表格记录中的了