Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Cuffdiff sqrt(det(cov))==0) error

    Hi All,

    I encountered several error messages during cuffdiff run. The error message is as follows:


    [13:43:11] Loading reference and sequence.
    No fasta index found for ../ref/all.fa. Rebuilding, please wait..
    Fasta index rebuilt.
    [13:44:04] Inspecting maps and determining fragment length distributions.
    [13:46:04] Calculating initial abundance estimates.
    > Processed 23484 loci. [*************************] 100%
    [13:51:17] Learning bias parameters.
    [13:53:15] Testing for differential expression and regulation in locus.
    > Processing Locus chr1:5441872-5445261 [ ] 2%Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    > Processing Locus chr10:18647262-18648745 [**** ] 17%Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    > Processing Locus chr10:19656777-19663522 [**** ] 17%Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    > Processing Locus chr10:20340517-20356103 [**** ] 18%Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    Error: sqrt(det(cov)) == 0, 0.000000 after rounding.


    This error did not halt the program and cuffdiff actually kept running until all results were output. I am not sure what did those errors mean? Are they fatal errors that might lead to the error in my result.

    Here is what I did before running cuffdiff.

    I ranthe Tophat to get my alignment results "sample_X/accepted_hits.bam" (I have 10 samples here, so X = 1..10). I ran Cufflinks to get the "sample_X/transcripts.gtf". Then I use Cuffcompare using all transcripts.gtf files output by cufflinks to get a combined gtf file. At the end, I feed this combined gtf file into cuffdiff and two samples from my sample pool to do the diff-express analysis.

    Thanks,
    Zheng

  • #2
    Hi Zheng,

    I have got the same problem. Did you hear sth. about this problem yet?

    Regards,
    Christian

    Originally posted by endether View Post
    Hi All,

    I encountered several error messages during cuffdiff run. The error message is as follows:


    [13:43:11] Loading reference and sequence.
    No fasta index found for ../ref/all.fa. Rebuilding, please wait..
    Fasta index rebuilt.
    [13:44:04] Inspecting maps and determining fragment length distributions.
    [13:46:04] Calculating initial abundance estimates.
    > Processed 23484 loci. [*************************] 100%
    [13:51:17] Learning bias parameters.
    [13:53:15] Testing for differential expression and regulation in locus.
    > Processing Locus chr1:5441872-5445261 [ ] 2%Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    > Processing Locus chr10:18647262-18648745 [**** ] 17%Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    > Processing Locus chr10:19656777-19663522 [**** ] 17%Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    > Processing Locus chr10:20340517-20356103 [**** ] 18%Error: sqrt(det(cov)) == 0, 0.000000 after rounding.
    Error: sqrt(det(cov)) == 0, 0.000000 after rounding.


    This error did not halt the program and cuffdiff actually kept running until all results were output. I am not sure what did those errors mean? Are they fatal errors that might lead to the error in my result.

    Here is what I did before running cuffdiff.

    I ranthe Tophat to get my alignment results "sample_X/accepted_hits.bam" (I have 10 samples here, so X = 1..10). I ran Cufflinks to get the "sample_X/transcripts.gtf". Then I use Cuffcompare using all transcripts.gtf files output by cufflinks to get a combined gtf file. At the end, I feed this combined gtf file into cuffdiff and two samples from my sample pool to do the diff-express analysis.

    Thanks,
    Zheng

    Comment


    • #3
      And I have the same problem now.
      Does anyone hear something from the authors or solve the problem?

      Comment


      • #4
        I have the same error (pasted below) when running cuffcompare in a local instance of Galaxy




        An error occurred running this job: cuffcompare v1.0.3 (2401)
        cuffcompare -o cc_output -s ref.fa ./input1 ./input2
        Error running cuffcompare. You are using Cufflinks v1.0.3, which is the most recent release.
        No fasta index found for ref.fa. Rebuilding, please wait..
        Fasta index rebuilt.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Strategies for Sequencing Challenging Samples
          by seqadmin


          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
          03-22-2024, 06:39 AM
        • seqadmin
          Techniques and Challenges in Conservation Genomics
          by seqadmin



          The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

          Avian Conservation
          Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
          03-08-2024, 10:41 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 06:37 PM
        0 responses
        10 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        9 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        49 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        67 views
        0 likes
        Last Post seqadmin  
        Working...
        X