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

登录

该问题已被锁定！

1: 关注

5599: 浏览

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

相关问题

enhancer结合ATAC-seq进行下游靶基因的寻找策略: 2182 浏览 2 关注 1 回答 0 评论

在单细胞测序中，不同组的T细胞统计检验方法: 2620 浏览 2 关注 2 回答 0 评论

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

单细胞RNA-seq，如何计算不同cluster之间的相关性: 2248 浏览 2 关注 1 回答 0 评论

为什么单细胞测序需要对UMI去重而二代测序不需要去重？: 2730 浏览 2 关注 1 回答 0 评论

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

菌侵染植物的转录组分析中，若菌的一个基因在当前时间没有引起植物基因的变化，过了三个小时后才引起植物基因的变化，如何分析这种现象或如何将这一类基因关联起来，有没有类似的文献推荐？: 1870 浏览 2 关注 1 回答 0 评论

对突变组和对照组的某基因表达量做T检验: 1776 浏览 2 关注 1 回答 0 评论

seurat包，如何调节大热图左侧的基因字体大小: 3251 浏览 2 关注 1 回答 0 评论

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

推荐内容

问 seurat包，如何调节大热图左侧的基因字体大小: 3251 浏览 2 关注 1 回答 0 评论

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

问你好chatGPT,可以提供包含双细胞去除的单细胞分析完整流程吗？: 2605 浏览 2 关注 2 回答 0 评论

问如何使用seqtk按照比例随机提取单细胞数据？: 2824 浏览 2 关注 2 回答 0 评论

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

问 smart-seq2样本质检报告说有小片段，请问那还可以建库么？: 2537 浏览 2 关注 2 回答 0 评论

问使用seurat包，导出特定cluster的细胞-基因counts矩阵。出现问题: 5396 浏览 2 关注 1 回答 0 评论

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

问 harmony处理批次效应: 3061 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025