cummeRbund sigMatrix wrong counts
1
0
Entering edit mode
Till ▴ 30
@till-9696
Last seen 4.2 years ago

Hi,

I'm using cummRbund to analyze some data. However the sigMatrix function seems to produce a graph with incorrect counts of DEGs, in my case they are much higher than they should be (as reported by getSig with same alpha).This is the same error as already reported at Seqanswers (http://seqanswers.com/forums/showthread.php?t=51230) it seems.

Is that actually a bug or am I missing something?

Cheers

 

cummeRbund • 870 views
ADD COMMENT
0
Entering edit mode
r.chereji • 0
@rchereji-8534
Last seen 5.2 years ago
United States

I noticed the same problem: the number of DEGs reported by sigMatrix is higher than the numbers reported by getSig using the same alpha. I would appreciate if someone could clarify why these 2 functions report different numbers and how one could use these 2 functions to get identical results.

Here is a short example:

library(cummeRbund)
cuff = readCufflinks(dir=system.file("extdata", package="cummeRbund"))

mySigMat = sigMatrix(cuff,level='genes',alpha=0.05)
mySigMat

setOfGenes = getSig(cuff,x='iPS',y='hESC',alpha=0.05,level='genes')
length(setOfGenes)

In the previous example, sigMatrix reports 120 DEGs when comparing hESC vs iPS, while getSig reports 117 DEGs. Why are these numbers different? 

 

Thanks!

ADD COMMENT

Login before adding your answer.

Traffic: 259 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6