pbmc.counts <- Read10X(data.dir ="pbmc3k/filtered_gene_bc_matrices/hg19/")#.数据导入,导入示例数据
过滤检测少于200个基因的细胞和少于3个细胞检测出的基因（min.cells = 3）#1个数据集，总共得到2700个细胞和13714个基因
pbmc <- CreateSeuratObject(counts = pbmc.counts) #2创建Seurat对象
pbmc <- NormalizeData(object = pbmc)
pbmc <- FindVariableFeatures(object = pbmc)
pbmc <- ScaleData(object = pbmc)
pbmc <- RunPCA(object = pbmc)
pbmc <- FindNeighbors(object = pbmc)
pbmc <- FindClusters(object = pbmc)
pbmc <- RunTSNE(object = pbmc)
pbmc <- RunUMAP(object = pbmc)
DimPlot(object = pbmc, reduction = "tsne")
网友评论