代码library(SingleR) counts<-pbmc@assays$RNA@counts clusters<-pbmc@meta.data$seurat_clusters ann=pbmc@meta.data$orig.ident singler = CreateSinglerObject(counts, annot = ann, \”pbmc\”, min.genes = 0, species = \”Human\”, citation = \”\”, ref.list = list(), normalize.gene.length = F, variable.genes = \”de\”, fine.tune = F, do.signatures = T, clusters = clusters, do.main.types = T, reduce.file.size = T, numCores = 1) singler$seurat = pbmc singler$meta.data$xy = pbmc@reductions$tsne@cell.embeddings clusterAnn=singler$singler[[2]]$SingleR.clusters.main$labels write.table(clusterAnn,file=\”07.clusterAnn.txt\”,quote=F,sep=\”\\t\”,col.names=F) write.table(singler$other,file=\”07.cellAnn.txt\”,quote=F,sep=\”\\t\”,col.names=F)
评论前必须登录!
注册