| 查看: 2791 | 回复: 0 | ||
hxppxh木虫 (正式写手)
|
[求助]
RNA-SEQ 过滤低表达基因遇到问题
|
|
我用trinity组装后，分析了每个样品的abundance，得到了每个样品的isoforms.results 文件，用这些文件分析差异表达基因时很顺利。 但是我想过滤掉一些低表达的基因，于是我将这些isoforms.results 文件中fpkm值小于某个范围的行删掉了，再分析差异表达基因运行rinityrnaseq_r20140413/util/abundance_estimates_to_matrix.pl ，结果出现如下错误： > library(edgeR) Loading required package: limma > > rnaseqMatrix = read.table("t_AS_fpkm10730_Trinity_genes.not_cross_norm.fpkm.tmp", header=T, row.names=1, com='') > rnaseqMatrix = round(rnaseqMatrix) > exp_study = DGEList(counts=rnaseqMatrix, group=factor(colnames(rnaseqMatrix))) > exp_study = calcNormFactors(exp_study) Error in if (any(allzero)) x <- x[!allzero, , drop = FALSE] : missing value where TRUE/FALSE needed Calls: calcNormFactors Execution halted Error, cmd: R --vanilla -q < __tmp_runTMM.R 1>&2 died with ret (256) at /diag/software/trinityrnaseq_r20140413/util/support_scripts/run_TMM_scale_matrix.pl line 98. 我觉得应该是那些isoforms.results文件出现了问题，可是不知道如何解决。谢谢！ |
回复此楼» 猜你喜欢
拉坦前列腺素(Latanoprost)：前列腺素F2α衍生物如何成为青光眼一线用药
已经有0人回复
-
卟啉家族中的明星分子：MESO-四(3,5-二溴-4-羟基苯基)卟啉的性质与制备
已经有1人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有133人回复
-
宠物寄生虫防治：氟雷拉纳阻断GABA受体，持效12周驱杀跳蚤蜱虫
已经有0人回复
-
为呼吸打开通道：依伐卡托的科学之旅
已经有1人回复
-
Dordaviprone（ONC201）：小分子药物在中线胶质瘤中的应用
已经有0人回复
-
依喜替康：新型喜树碱衍生物的研究进展
已经有1人回复
-
膀胱癌靶向治疗新选择：厄达替尼作用机制与耐药研究进展
已经有0人回复
-
从水母到实验室：腔肠素的发现历程与生物发光奥秘
已经有1人回复
-
线粒体氧化磷酸化的新靶点：S-Gboxin的发现与研究进展
已经有0人回复
-
PROTAC药物开发选择VHL配体：MDK-7526以87%出口向量占有率成首选
已经有0人回复
10