Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Cufflinks, differentially expressed genes statsteam Bioinformatics 5 11-15-2013 11:28 AM
Pathway analysis of Differentially expressed genes himanshu04 RNA Sequencing 12 09-05-2013 08:04 AM
too many Differential expressed genes using cuffdiff papori Bioinformatics 4 04-30-2012 01:42 PM
Finding diferentially expressed genes Amative Bioinformatics 8 04-09-2012 09:44 PM
DESeq and EdgeR: too many differentially expressed genes!?!? cutcopy11 Bioinformatics 5 12-08-2011 12:14 AM

Thread Tools
Old 07-13-2012, 09:59 AM   #1
Location: belgium

Join Date: Jul 2011
Posts: 11
Default cuffdiff2 fails to find diff expressed genes!!??


I processed 8 breast cancer samples, 4ER+ and 4ER- using tophat2/cufflinks2/cuffmerge2/cuffdiff2/cummeRbund.

The expression of the gene GREB1 is nearly equal to 0 for the 4 ER-
and ~100 for ER+:

cuff <- readCufflinks()
Neg_0 Neg_1 Neg_2 Neg_3 Pos_0 Pos_1 Pos_2 Pos_3
XLOC_011466 0.182529 0.791468 0.721075 0.15861 35.9609 231.247 64.1414 39.4893

However, the result of cuffdiff is:
diffData(getGene(cuff, 'GREB1'))
gene_id sample_1 sample_2 status value_1 value_2
log2_fold_change test_stat p_value q_value significant
1 XLOC_011466 Neg Pos OK 0.460534 92.726
7.65352 -1.1267 0.25987 0.999954 no

How can you explain this. It seems not consistent for me. I would have
expected a nice q value, a significant test, etc.

thanks for feedbacks,

ps. data were downloaded from SRA/GEO : GSE27003 and are available at inSilico:
colinmolter is offline   Reply With Quote

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 06:08 AM.

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