Many DESeq2 P values are 0 thus preventing generation of a rank list for GSEA
1
1
Entering edit mode
4.8 years ago

Hi Everyone,

In my DESeq2 reuslts files, the smallest 100 P-values are 0 (as the smallest floating point value in R is 2.225074e-308). I am wondering if there is any way to discern the unrounded values of these P-values as they are needed in order for me to produce a rank list for gene set enrichment analysis (GSEA). Any guidance on how to approached ranking these genes would be appreciated.

Thank you

RNA-Seq gene R • 3.5k views
ADD COMMENT
0
Entering edit mode

Quick and dirty: use log-fold(change) instead. Or rank by P value, then by log-fold(change).

P.S.: see discussion here for alternatives, implications and explanations.

ADD REPLY
1
Entering edit mode
4.8 years ago
Martombo ★ 3.0k

you can use the stat column in the results object, as this (in absolute values) is monotonic when compared to the p-values. Some additional info on what metric to use for ranking can be found in this paper.

ADD COMMENT
0
Entering edit mode

Thanks for the article, I'm sure it will be very helpful!

ADD REPLY

Login before adding your answer.

Traffic: 1561 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