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

登录

该问题已被锁定！

1: 关注

6063: 浏览

单细胞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 人关注

相关问题

基因本体论分析上下调基因差异: 1866 浏览 2 关注 1 回答 0 评论

使用MAS5.0处理过的芯片数据，能不能直接取对数后用limma包找差异表达基因，为什么不能？: 2259 浏览 2 关注 1 回答 0 评论

RNA-seq，样本中存在同一基因对应不同的FPKM值: 2067 浏览 2 关注 1 回答 0 评论

bowtie2 参考基因组注释比对: 3510 浏览 2 关注 2 回答 0 评论

基因表达芯片的分类问题: 2390 浏览 3 关注 2 回答 0 评论

Bulk-RNAseq多组数据差异表达基因的筛选: 3013 浏览 2 关注 2 回答 0 评论

验证数据集基因名称: 2093 浏览 2 关注 1 回答 0 评论

EVM整合基因组注释: 2181 浏览 2 关注 1 回答 0 评论

如何根据转录组数据得到新转录本？如何验证一个基因的多个转录本？: 3639 浏览 2 关注 2 回答 0 评论

annotatePeak peak基因组注释: 3232 浏览 2 关注 1 回答 0 评论

推荐内容

问 singularity查到不到指定输入文件位置: 2835 浏览 2 关注 1 回答 0 评论

问在做RNA Velocity的时候增殖细胞群总是起点: 2572 浏览 2 关注 2 回答 0 评论

问用scanpy做单细胞分析，如何去除周期细胞的影响: 3335 浏览 2 关注 1 回答 0 评论

问 Monocle3绘制自定义轨迹错误: 3390 浏览 2 关注 1 回答 0 评论

问 celseq2转换单细胞原始数据: 2641 浏览 2 关注 2 回答 0 评论

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

问 DoubletfFinder Pk值相关计算的问题: 2804 浏览 2 关注 2 回答 0 评论

问空转bin_size: 3217 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025