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

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.