Question: GSEA preranking metric for RNA Seq
gravatar for bipin
16 months ago by
bipin10 wrote:

I came across multiple posts regarding the pre-ranking metric for GSEA when using RNA seq data. However, there doesn't seem to be a consensus.

Some of the metrics I came across are:-

  • sign of log fold change * -log10(p-value[not adjusted p-val])
  • logfc shrink values from DESeq2
  • Inbuilt signal2noise from GSEA. However, this cannot be used in case of <3 replicates.

What metric do you use for ranking the genes or you know is widely used?

gsea rna-seq deseq2 • 1.5k views
ADD COMMENTlink modified 16 months ago by Kevin Blighe44k • written 16 months ago by bipin10
gravatar for Kevin Blighe
16 months ago by
Kevin Blighe44k
South America | Europe | USA
Kevin Blighe44k wrote:

It makes sense that there is no consensus, as there are countless ways to do this. My own recommendation would be to:

  1. Set an adjusted P value cut-off
  2. Rank genes based on absolute log (base 2) fold change

I believe the most widely used method is to just set an adjusted P value and log (base 2) fold change cut-off, and to then 'throw' the resulting gene list into the GSEA without any ranking.

The lack of consensus on a proper filtering strategy may in part be due to the fact that a substantial proportion of researchers do not pay much attention to the results of GSEA. GSEA results would certainly never stand as the sole evidence in a clinical test, neither would they be sufficient evidence on which conclusions could be made in most reputable journals.


ADD COMMENTlink written 16 months ago by Kevin Blighe44k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1286 users visited in the last hour