SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
multiple lines per single gene in gene_exp.diff output of cuffdiff reut Bioinformatics 2 03-09-2015 04:16 AM
Cuffdiff splicing.diff Greg Bioinformatics 1 11-04-2013 04:04 PM
Why can't cuffdiff identify splicing difference Jolin Bioinformatics 0 10-14-2011 08:30 AM
multiple FPKM problem for single gene in gene_exp.diff after running cuffdiff ngs RNA Sequencing 4 03-30-2011 01:55 PM
Strange cuffdiff result (disagreement between tracking and diff files) jebe Bioinformatics 0 07-01-2010 10:42 AM

Reply
 
Thread Tools
Old 06-28-2010, 02:35 PM   #1
Greg
Member
 
Location: British Columbia

Join Date: Oct 2009
Posts: 31
Default Cuffdiff NOTEST for splicing.diff

Using the tophat/cufflinks combo to find differential alternative splicing.

The cuffdiff manual says that if there "not enough alignments for testing" then you will get NOTEST and the treatments will not be compared. Does anyone have an idea as to what the cutoff for this is? Many of the genes that get NOTEST in my data set (single end) have a large number of reads (several hundred). Could it be because of insufficient coverage? or maybe operator error?

thanks!
Greg is offline   Reply With Quote
Old 06-30-2010, 05:14 PM   #2
Greg
Member
 
Location: British Columbia

Join Date: Oct 2009
Posts: 31
Default

So I just found the -c option to set the number of required reads. It does change the number of tests that happen but there are still less then I would think. Any ideas? What are the other criteria for doing the test or not? Would more info on my workflow be helpful?
Greg is offline   Reply With Quote
Old 07-14-2010, 02:49 PM   #3
griffon42
Member
 
Location: New York

Join Date: Jan 2009
Posts: 23
Default

I'm experiencing the same issue. I'm using cufflinks 0.8.3 (the 7/2 release) and the -c option doesn't seem to work as described in the manual. As extreme examples, I have loci that are "OK" for stat testing and "yes" for significant but have only 2 reads vs. 0 reads across samples.

As mentioned, changing -c DOES affect the number of loci tested, but not by much and it does not fix situations like that described above.

Does anyone know how to manage this situation? The differential expression output files i'm generating are rather unmanageable.

Thanks.
griffon42 is offline   Reply With Quote
Old 07-14-2010, 09:43 PM   #4
jetspeeder
Member
 
Location: U.S.

Join Date: Jun 2010
Posts: 12
Default

I have a lot of NOTEST and FAIL as well,

Also do you guys know the significance of this and why the percentages are high for missed and wrong exons and what the signifies? Thanks.

#> Genomic sequence: chr1
# Query mRNAs : 20972 in 20055 loci (3356 multi-exon transcripts)
# (633 multi-transcript loci, ~1.0 transcripts per locus)
# Reference mRNAs : 1024 in 912 loci (966 multi-exon)
# Corresponding super-loci: 695
#--------------------| Sn | Sp | fSn | fSp
Base level: 63.4 16.4 - -
Exon level: 37.6 9.6 44.0 11.3
Intron level: 52.6 31.5 53.4 31.9
Intron chain level: 18.6 5.4 24.7 7.1
Transcript level: 0.0 0.0 0.1 0.0
Locus level: 19.5 0.9 22.9 1.0
Missed exons: 2554/8302 ( 30.8%)
Wrong exons: 26036/32327 ( 80.5%)
Missed introns: 3161/7289 ( 43.4%)
Wrong introns: 7825/12186 ( 64.2%)
Missed loci: 208/912 ( 22.8%)
Wrong loci: 17823/20055 ( 88.9%)
jetspeeder is offline   Reply With Quote
Old 07-15-2010, 09:41 AM   #5
Greg
Member
 
Location: British Columbia

Join Date: Oct 2009
Posts: 31
Default

I've had the same experience. Could only guess...
Greg is offline   Reply With Quote
Old 11-14-2011, 05:22 AM   #6
SEQond
Member
 
Location: Italy

Join Date: Jul 2010
Posts: 27
Default

So anyone has an idea on how to find the ideal -c "minimum number of alignments in a locus for needed to conduct significance testing" (--min-alignment-count) ?
SEQond is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:46 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO