| 查看: 2829 | 回复: 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文件出现了问题，可是不知道如何解决。谢谢！ |
回复此楼» 猜你喜欢
SEL120-34A HCl：CDK8/CDK19高选择性激酶抑制剂
已经有0人回复
-
氢氧化钾活化操作
已经有18人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有128人回复
-
散金
已经有68人回复
-
生物质超级电容器
已经有0人回复
-
ALK5 高选择性小分子 LY-3200882：可放大合成路线与质控要点解析
已经有0人回复
-
AM630（CAS 164178-33-0）：CB2 选择性配体的实验使用与供应要点
已经有0人回复
-
ADC标杆分子T-DM1全解析：从结构设计到实验室应用的技术要点
已经有0人回复
-
HILIC磁珠研究蛋白组学怎么做？
已经有0人回复
-
His标签蛋白纯化怎么才能快？
已经有0人回复