harmony处理批次效应 - 做专业的生信问答网站

登录

该问题已被锁定！

2: 关注

1285: 浏览

harmony处理批次效应

单细胞RNA_seq

假如现在有sampleA，sampleB两个样本，在合并之前，是否需要单独对sampleA 和sampleB做FindVariableFeatures，ScaleData，RunPCA呢？还是直接合并sampleA 和sampleB的seurat对象，合并之后统一去做？然后再进行harmony

好问题 0 评论收藏举报

查看全部 2 个回答

前台管理员用户来自于: 北京市
2023-05-25 13:53

个人建议是，在做样本合并之前，先对每个样本单独进行FindVariableFeatures，得到每个样本的高变异基因（hvf），而后取样本hvf的并集，作为后续合并完成的SeuratObject的hvf，这样可以保证每个数据集的差异特征尽可能的保留。而后，再对合并完成的SeuratObject进行ScaleData和RunPCA的操作，然后再RunHarmony进行批次效应矫正。

赞同 1 0评论

问题动态

发布时间: 2023-05-23 15:44

更新时间: 2023-05-25 13:53

关注人数: 2 人关注

相关问题

测了实验组与对照组的mRNA与miRNA的转录组，miRNA-mRNA互作分析网络怎么处理呢？: 851 浏览 2 关注 2 回答 0 评论

奇异值分解在处理芯片数据时的意义是什么？: 651 浏览 1 关注 0 回答 0 评论

如何在差异表达前消除批次效应: 800 浏览 2 关注 1 回答 0 评论

stringtie和gffcompare处理转录本的问题: 698 浏览 2 关注 1 回答 0 评论

请问RNA-seq采用poly A(+)策略建库，处理数据时若不去除rRNA会对后续分析有何影响？: 770 浏览 2 关注 1 回答 0 评论

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

xpclr的XPCLR_score的inf值应该怎么处理: 841 浏览 2 关注 1 回答 0 评论

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

缺失值处理: 820 浏览 2 关注 1 回答 0 评论

Harmony 和Seurat的FindIntegrationAnchors算法各有何优劣: 825 浏览 2 关注 1 回答 0 评论

推荐内容

问单细胞RNA-seq分析流程: 861 浏览 2 关注 2 回答 0 评论

问整合scRNA和scATAC相关问题请教: 1043 浏览 2 关注 3 回答 0 评论

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

问 scRNA-seq有什么推荐的监督聚类方式吗？: 923 浏览 2 关注 1 回答 0 评论

问 cellranger运行结果分析: 695 浏览 2 关注 1 回答 0 评论

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

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

问单细胞seurat对象的基因过滤: 2569 浏览 1 关注 3 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2024