Question: If genes in one gene set are abundant in top, bottom of ranked list?
2
gravatar for ghmdsr
9 days ago by
ghmdsr20
ghmdsr20 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 • 93 views
ADD COMMENTlink modified 8 days ago by Zhilong Jia1.3k • written 9 days ago by ghmdsr20

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 9 days ago by russhh3.5k
0
gravatar for Zhilong Jia
8 days ago by
Zhilong Jia1.3k
London
Zhilong Jia1.3k 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 8 days ago by Zhilong Jia1.3k
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: 658 users visited in the last hour