Question: If genes in one gene set are abundant in top, bottom of ranked list?
2
gravatar for ghmdsr
5 months ago by
ghmdsr30
ghmdsr30 wrote:

Let's suppose that there are 100 genes in one specific gene set and 20000 genes are used in making ranked list.

50 genes in the gene set are located in top of ranked list(pvalue very low,sign of fold-change= +), and the other 50 genes are located in the bottom of ranked list.(pvalue very low, sign of fold-change= -)

If then, I think Enrichment plot is similar below. How can I get ES between two peaks if two peaks in this picture are same? I used ranked list which uses inverse of pvalue and sign of fold change

enter image description here

rna-seq gene • 226 views
ADD COMMENTlink modified 5 months ago by Zhilong Jia1.4k • written 5 months ago by ghmdsr30

which tool are you using for enrichment? I believe you can modify the power to which your normalised ranks are raised, try squaring them

ADD REPLYlink written 5 months ago by russhh3.9k
0
gravatar for Zhilong Jia
5 months ago by
Zhilong Jia1.4k
London
Zhilong Jia1.4k wrote:

You mean GSEA. The ES is 0 if it is unweighted or related to the stats used (inverse of pvalue and sign of fold change) if weighted (stats)your stats. Details are as below. The GSEA implementation used is fgsea or from Github. and parameter gseaParam: GSEA weight parameter (0 is unweighted, suggested value is 1).

library(fgsea)
data(exampleRanks)
pathway_test <- list(top50_bot50=names(c(head(exampleRanks, 50), tail(exampleRanks, 50) )) )

# unweighted
fgseaRes <- fgsea(pathways = pathway_test, stats = exampleRanks, minSize=15, maxSize=500, nperm=10000, gseaParam=0)

enter image description here

plotEnrichment(pathway_test$top50_bot50, exampleRanks, gseaParam=0)

enter image description here

# weighted
fgseaRes1 <- fgsea(pathways = pathway_test, stats = exampleRanks, minSize=15, maxSize=500, nperm=10000, gseaParam=1)

enter image description here

 plotEnrichment(pathway_test$top50_bot50, exampleRanks, gseaParam=1)

enter image description here

ADD COMMENTlink written 5 months ago by Zhilong Jia1.4k
Please log in to add an answer.

Help
Access

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