单细胞seurat对象的基因过滤 - 做专业的生信问答网站

登录

该问题已被锁定！

1: 关注

5729: 浏览

单细胞seurat对象的基因过滤

单细胞seurat对象x，在质控(去掉线粒体基因表达比例>10%的细胞等)之后，想去除线粒体基因和核糖体基因等，应该如何操作呢？尝试了如下代码：mito_genes <- grep("^mt-", rownames(x), value = TRUE)
x <- subset(x, features = !rownames(x) %in% mito_genes)没有成功

好问题 0 评论收藏举报

查看全部 3 个回答

Ezio_2023 注册会员用户来自于: 云南省昆明市呈贡县
2023-05-18 16:01

如果你的不是人的话，上面的应该不行，建议从gtf里面提取线粒体的基因来做，我这个是猴子的

mt.genes <- c("ND1","ND2","COX1","COX2","ATP8","ATP6","COX3","ND3","ND4L","ND4","ND5","ND6","CYTB")
obj[["percent.mt"]] <- PercentageFeatureSet(obj, features = mt.genes[mt.genes %in% rownames(RTT4)])
obj[["percent.ribo"]] <- PercentageFeatureSet(obj, pattern = "^RP[L|S]")

赞同 0 4评论

关于作者

: SeaWestKing 注册会员
这家伙很懒，还没有设置简介

0: 回答

0: 文章

3: 问题

问题动态

发布时间: 2023-05-18 15:39

更新时间: 2023-05-18 17:06

关注人数: 1 人关注

相关问题

单细胞数据整合: 2640 浏览 2 关注 1 回答 1 评论

批量获得基因间区及内含子的突变距离最近的基因外显子的距离: 1791 浏览 2 关注 1 回答 0 评论

riboseq中frame指的是什么，是如何确定的，通过riboseq如何知道一个基因的开放阅读框: 2159 浏览 2 关注 1 回答 0 评论

在很多分析表达谱芯片的教程中，使用limma 包寻找差异表达基因，都是先找到差异表达探针，再把探针注释为基因id，为何不先进行注释再计算差异表达呢？这两种方法有什么区别呢？: 1942 浏览 2 关注 1 回答 0 评论

3D基因组里，如何分清compartment、TAD和chromatin loop: 3440 浏览 2 关注 2 回答 0 评论

单细胞亚类聚类和相关细胞分析: 1983 浏览 2 关注 1 回答 0 评论

如何按一个列表对基因型文件进行过滤，剔除不需要的样本？: 2774 浏览 2 关注 1 回答 0 评论

如何理解单细胞测序中用fastQC对低质量细胞进行过滤: 2396 浏览 2 关注 1 回答 0 评论

转录本坐标转换成基因组坐标: 3282 浏览 4 关注 2 回答 0 评论

3D基因组里compartment里一般是包含好多TAD的，但这图为什么compartment数量比TAD多这么多呢: 2510 浏览 2 关注 2 回答 0 评论

推荐内容

问单细胞——小鼠T细胞注释的markers: 3469 浏览 2 关注 2 回答 0 评论

问运行harmony后counts数变成小数: 1917 浏览 1 关注 1 回答 0 评论

问使用seurat包，导出特定cluster的细胞-基因counts矩阵。初始数据命名为pbmc: 2824 浏览 2 关注 1 回答 0 评论

问去除批次效应选定靶标细胞群后如何进行亚群分析？: 3224 浏览 2 关注 2 回答 0 评论

问如何理解单细胞测序中用fastQC对低质量细胞进行过滤: 2396 浏览 2 关注 1 回答 0 评论

问单细胞转录因子: 2097 浏览 2 关注 1 回答 0 评论

问单细胞monocle3，路径很不对，有什么参数是值得调的吗: 3163 浏览 2 关注 1 回答 0 评论

问 GEO数据读入: 1759 浏览 2 关注 1 回答 0 评论

问在用seurat做多样本整合的时候的问题: 2804 浏览 2 关注 1 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025