rm(list = ls())
library(Seurat)
pbmc<-readRDS("pbmc3k_final.rds")
pbmc$orig.ident <- colnames(pbmc)
set.seed(111) # 让每次抽样可重复
sampled.cells <- names(sample(x = pbmc@active.ident, size = 500, replace = F)) # size选择抽取细胞数
pbmc_500 <- subset(pbmc,subset = orig.ident %in% sampled.cells )
scRNA-seq | 随机抽选500个细胞
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 20、第二十计 混水摸鱼比喻趁混乱时机攫取不正当的利益。也作“浑水摸鱼”。此计是说打仗时要得于抓住敌方的可乘之隙,...
- 简介 在本教程中,我们将对PBMC的两个scATAC-seq数据集(5K和10K)和一个scRNA-seq数据集进...
- 课程学习生信技能树单细胞转录组(基础) 现在福利来了,单细胞课程免费开放了:#完整版视频上线!全!面!免!费!#全...
- 一、介绍 用于单细胞RNA-seq数据 提供严格的质量控制:将原始测序读数处理为可用于下游分析的高质量表达数据集 ...
- 本课件介绍了单细胞转录组数据分析的一个分析点:细胞类型识别。 注意区分细胞类型与细胞状态 目前大部分细胞类型的定义...