老师，图1是我用图2的代码实现的，现在我想仿照提取细胞亚类的方式，使用图3的代码单独提取TIP这个maker基因，不提取整个cluster3，这个代码可以实现吗，或者是否有其他代码可以使用

默认排序时间排序

1 个回答

omicsgene - 生物信息 3天前

擅长：重测序,遗传进化,转录组,GWAS

基因是在细胞中表达，不明白你：提取TIP这个maker基因是什么意思？
是想提取TIP高表达的细胞吗？你可以这样：

Rscript $scripts/subset_seurat_obj.r --rds  cluster3.qs --subset 'TIP >1' -p cluster3.TIP