| 查看: 2808 | 回复: 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文件出现了问题，可是不知道如何解决。谢谢！ |
回复此楼» 猜你喜欢
含氟糖皮质激素的外用抗炎利器醋酸氟轻松
已经有0人回复
-
阿莫奈韦（Amenamevir）：解旋酶-引物酶靶向分子的结构特征与合成路径全解析
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有69人回复
-
利瑞鲁单抗：靶向KIR的NK细胞检查点工具抗体，结构基础与科研应用全解析
已经有0人回复
-
Wnt-C59：PORCN靶向化学探针，Wnt信号通路阻断机制与科研应用全解析
已经有0人回复
-
Fexaramine：肠道限制性FXR化学探针，结构特征与信号通路研究全解析
已经有0人回复
-
NVS-STG2：靶向STING跨膜域的分子胶类激动剂，结构特征与寡聚化机制研究全解析
已经有0人回复
-
DOTA-GPC3肽（CAS 3066439-94-6）：双功能分子探针前体的结构特征与标记应用解析
已经有0人回复
-
科研用 TLR7/8 双激动剂怎么选？现货与质控要点
已经有0人回复
-
SEL120-34A HCl：CDK8/CDK19高选择性激酶抑制剂
已经有0人回复
