单细胞 10X 和seurat对象学习

单细胞seurat数据的基础知识

rm(list = ls()) 
library(Seurat)
#注意这个报错
#Warning: Feature names cannot have underscores ('_'), replacing with dashes ('-')
folders=list.files('./',pattern='[123]$')
folders
scList = lapply(folders,function(folder){ 
  CreateSeuratObject(counts = Read10X(folder), 
                     project = folder,
                     min.cells = 3, min.features = 200)
})

BM <- merge(scList[[1]], 
            y = c(scList[[2]],scList[[3]]), 
            add.cell.ids = c("BM1","BM2","BM3"), 
            project = "BM")
BM
GM <- merge(scList[[4]], 
            y = c(scList[[5]],scList[[6]]), 
            add.cell.ids = c("GM1","GM2","GM3"), 
            project = "GM")