SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DotKnot issues lstbl Bioinformatics 1 09-20-2015 07:07 AM
Some issues about DE genes detection by Cuffdiff wisense RNA Sequencing 0 09-29-2013 10:01 PM
Solution found: For Cuffdiff/links 2.0.2 Make cufflinks FPKM match Cuffdiff FPKM NGSfan RNA Sequencing 4 04-16-2013 07:10 AM
cuffdiff error, issues with bam sorting rhcr56 Bioinformatics 1 10-18-2012 11:21 AM
issues found in using cufflinks/cuffcompare/cuffdiff sterding Bioinformatics 5 06-01-2011 08:04 PM

Reply
 
Thread Tools
Old 11-18-2015, 12:00 AM   #1
ClemBuntu
Member
 
Location: Lyon

Join Date: Dec 2014
Posts: 37
Default Issues with cuffdiff

Hello,

I'm performing an analysis containing 3 samples per conditions and 5 conditions (so 15 samples in total). I want to do pairwise comparison for each condition.

I'm just starting using cuffdiff.
Firsly I just tested cuffdiff with only 2 conditions, like that :

Code:
cuffdiff -L control,A -p 16 --library-type fr-firststrand --no-update-check  reference.gtf control_1.bam,control_2.bam,control_3.bam A_1.bam,A_2.bam,A_3.bam

And I obtained 3 differentially expressed genes in the gene_exp.diff file.

Then I tried to do all my comparisons in one command line :

Code:
cuffdiff --output-dir allComparisons  -L control,A,B,C,D -p 16 --library-type fr-firststrand --no-update-check  reference.gtf control_1.bam,control_2.bam,control_3.bam A_1.bam,A_2.bam,A_3.bam B_1.bam,B_2.bam,B_3.bam C_1.bam,C_2.bam,C_3.bam D_1.bam,D_2.bam,D_3.bam

And looking at the gene_exp.diff file in output, I obtained 2 differentially expressed genes between control and A which only one in common with the 3 genes I obtained before.

Which analysis I should trust ?
And what is the good way to use cuffdiff ?

Last edited by ClemBuntu; 11-18-2015 at 12:05 AM.
ClemBuntu is offline   Reply With Quote
Reply

Tags
cuffdiff

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 10:51 PM.


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