该问题已被锁定!
2
关注
948
浏览

关于tophat2 + cufflinks 得到的差异基因和featureCounts差异基因数目相差太大原因

为什么被折叠? 0 个回复被折叠
孟浩巍 超级管理员 用户来自于: 北京市
2018-09-08 21:01
我建议: cuffdiff结果里面满足下面3个条件的算显著基因: 1. sample1 FPKM 或者 sample2 FPKM至少有1个大于1; 2. log2foldchange 的绝对值大于1; 3. qvalue < 0.05;   DESeq2 或者 edgeR的结果中满足下面3个条件 算显著基因; 1. sample1 read count 或者 sample2 read count至少有1个大于10条; 2. log2foldchange 的绝对值大于1; 3. qvalue < 0.05;   然后把这两个gene list做overlap再做比较。
qliu 初级会员 用户来自于: 华中农业大学
2018-09-08 20:57
图中mut开头的表示我用sujunc+featureCounts+DESeq2得到的差异基因,zsl表示我用cuffdiff得到的差异基因, 请问为啥和cufflinks的结果差别这么大呢,谢谢了 [attach]167[/attach] 我之前比较过不用的比对软件得到bam然后再用reads得到差异基因(不是本组数据结果),他们几乎没什么差别,详细见, [attach]168[/attach] [attach]169[/attach]   后者比较我之前记的笔记 [url=https://www.jianshu.com/p/aa7e019ca8a3]不同比对软件和HTseq+featureCounts之间的比较[/url]  

关于作者

qliu 初级会员

这家伙很懒,还没有设置简介

问题动态

发布时间
2018-09-08 20:49
更新时间
2018-09-08 21:01
关注人数
2 人关注

相关问题

怎么寻找某个基因的不同可变剪切形式?
一个cohort中,某基因count大于多少才有意义?
使用seurat包,导出特定cluster的细胞-基因counts矩阵。初始数据命名为pbmc
绵羊转录组测序后差异表达基因太多了
小麦660K芯片里,参考碱基那一列的A/C,哪个才是参考基因组的碱基?
根据GFF和fasta等文件提取某一基因的ATG位点信息,比如具体的位置?
新人求教关于测序以及免疫的问题
GEO数据芯片数据基因名转换
关于蛋白质文件的疑问
关于几个数据库对GO注释的疑问

推荐内容

trimmomatric使用报错
耐药细菌的比较基因组分析流程和机制探索
All Rights Reserved Powered BY WeCenter V4.1.0 © 2024