该问题已被锁定!
1
关注
2791
浏览

单细胞seurat对象的基因过滤

为什么被折叠? 0 个回复被折叠
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]")

 

standself 注册会员 用户来自于: 广东省广州市
2023-05-18 16:52

你的代码有问题啊。应该是 x <- susbet(x, features = features_to_retain)  吧。

chatGPT机器人 机器人 用户来自于: 北京市
2023-05-18 15:46

为了过滤掉单细胞Seurat对象中的线粒体基因和核糖体基因,您可以使用以下代码:

# 获取线粒体基因和核糖体基因的名称
mitogenes <- grep("^MT-", rownames(x), value = TRUE) ribogenes <- grep("^RPS|^RPL", rownames(x), value = TRUE)

# 过滤掉线粒体基因和核糖体基因
x <- subset(x, features = !(rownames(x) %in% c(mitogenes, ribogenes)))

在这里,我们使用了grep函数来获取线粒体基因和核糖体基因的名称。

然后,我们使用subset函数来过滤掉Seurat对象中的这些基因。

最后,我们可以将过滤后的Seurat对象用于后续的分析。

请注意,上述代码中的正则表达式可能需要根据您的数据进行修改。如果您的线粒体基因或核糖体基因不是以"MT-"或"RPS/RPL"开头,您需要相应地更改正则表达式。

关于作者

问题动态

发布时间
2023-05-18 15:39
更新时间
2023-05-18 17:06
关注人数
1 人关注

相关问题

怎么对基因启动子含有什么元件进行分析
基因组组装问题
RNA-seq,样本中存在同一基因对应不同的FPKM值
关于hub基因的问题
基因互作网络和蛋白互作网络
如果在基因组中鉴定到病毒序列,如何从生信角度判断它是EVEs还是HGTs
基因symbol注释GEO芯片时,如何把mRNA和lncRNA分别标注出来?
在单细胞测序中,不同组的T细胞统计检验方法
基因组组装
关于tophat2 + cufflinks 得到的差异基因和featureCounts差异基因数目相差太大原因

推荐内容

monocle3中plot_cells更改其画出来的细胞群的颜色
cellranger运行结果分析
去除批次效应选定靶标细胞群后如何进行亚群分析?
​单细胞RNA-seq,如何计算不同cluster之间的相关性
请教多个scRNA样本整合问题
smart-seq2样本质检报告说有小片段,请问那还可以建库么?
单细胞分析中在umap聚类(seurat, scanpy)有调整群位置的参数吗?
harmony处理批次效应
harmony整合样本前需要分别预处理吗?
All Rights Reserved Powered BY WeCenter V4.1.0 © 2024