zoukankan      html  css  js  c++  java
  • cuffquant 定量的结果全部为-nan

    cuffquant 定量的过程中,当所有基因或者转录本的表达量都为0时,定量的结果就回全部是-nan  , 而不是0;

    出现这种情况有两种原因:

    1) 参考基因组搞错了,比对和定量的不是同一个参考基因组,或者gtf 文件和fasta 文件中染色体标识符对应不上,最常见的情况就是,一个以chr开头,一个不带chr;

    2) 确实是0,就是说比对灭有问题,就是没有比对上基因区间;当然这种概率非常小,如果遇到,只能说是小确幸;

    值得高兴的是,今天我就遇到了这种情况,同一批样本,其中一个样本基因的表达量全部为-nan, 找了好久的原因,以为自己哪里搞错了,原来就是这么的巧合,这种极端情况都会遇见

    真是屋漏偏逢连夜雨,还是梅雨!!!

  • 相关阅读:
    换个格式输出整数
    c++入门构造函数
    接口继承与实现
    spring mvc 之@requestmapping
    RabbitMQ
    Socket通信原理
    TCP三次握手与四次挥手
    WIN 10家庭版没有Hyper-V
    aop 切面demo
    JVM之旅------jvm内存模型
  • 原文地址:https://www.cnblogs.com/xudongliang/p/7489675.html
Copyright © 2011-2022 走看看