Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • kmcarr
    Senior Member
    • May 2008
    • 1181

    Errors in cuffdiff genes.fpkm_tracking file

    I have an mRNA-Seq experiment with a number of samples. After running them through tophat and cuffcompare I then proceeded to run the data through cuffdiff. I wasn't particularly interested in any specific differential expression testing, I just wanted to re-estimate gene level relative expression. I ran cuffdiff providing a curated combined.gtf generated by cuffcompare and four accepted_hits.sam files. According to the cufflinks manual the genes.fpkm_tracking file should contain 4 sets of qX_FPKM, qX_conf_lo and qX_conf_hi data. I was expecting q0 - q3 which is what I see in the isoforms, tss_group tracking files, but in the genes tracking file the first row (headers) lists 12 sets, q0 - q11. What's more, not all rows (loci) have the same number of data points. Many have 12 sets, consistent with the header row, some have only 4 sets which is the expected number and I've even found some rows with as many as 36 sets of FPKM and conf values.

    I am using version 0.8.2 of cufflinks. This result occurs with both the prebuilt x86_64 binaries and binaries I compiled myself.

    Has anyone else observed this behavior when providing more than two samples to cuffdiff?

    Update:

    I have now observed this behavior when only providing 2 input sam files to cuffdiff. The genes.fpkm_tracking file has q0 - q5 data sets.
    Last edited by kmcarr; 05-27-2010, 09:26 AM. Reason: Added information

Latest Articles

Collapse

  • seqadmin
    Pathogen Surveillance with Advanced Genomic Tools
    by seqadmin




    The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
    03-24-2025, 11:48 AM
  • seqadmin
    New Genomics Tools and Methods Shared at AGBT 2025
    by seqadmin


    This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

    The Headliner
    The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
    03-03-2025, 01:39 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-20-2025, 05:03 AM
0 responses
49 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-19-2025, 07:27 AM
0 responses
57 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-18-2025, 12:50 PM
0 responses
50 views
0 reactions
Last Post seqadmin  
Started by seqadmin, 03-03-2025, 01:15 PM
0 responses
201 views
0 reactions
Last Post seqadmin  
Working...