View Single Post
Old 12-10-2016, 03:34 PM   #1
Rebekah_h4545
Junior Member
 
Location: Missouri

Join Date: Dec 2016
Posts: 1
Unhappy Help with Cuffdiff output and CummeRbund

Dear community,

I have RNA-seq data of 3 different samples, mock, condition 1 and condition 2. I have used the tuxedo suite for my analysis. I bypassed cuffmerge and used only cuffcompare. I finished with cuffdiff. The cuffdiff output looks perfect, all the files are there. The problem lies with cummeRbund. After making my cuffData.db file I immediately receive 22 warnings. I went back to my cuffdiff results (gene_exp.diff and cds_exp.diff file) and saw that they had "-inf" and "=-nan" values in log2(fold change) and test_stat column.

How can I have so many errors when I ran everything in the correct way? I ran cufflinks and tophat with an annotation file and I have proper files (seemingly proper files) for everything up until cuffdiff. I don't know how to circumvent these Nan or inf values for visualization in cummeRbund. I can get an adobe pdf of various plots, but they are damaged files and can't be repaired and therefore not viewed. My organism is Arabidopsis thaliana.

Any help would be appreciated.
Attached Images
File Type: png Screen Shot 2016-12-10 at 6.23.15 PM.png (35.5 KB, 5 views)
File Type: png Screen Shot 2016-12-10 at 6.25.44 PM.png (125.2 KB, 4 views)
Rebekah_h4545 is offline   Reply With Quote