Question: If genes in one gene set are abundant in top, bottom of ranked list?
2
gravatar for ghmdsr
2.4 years ago by
ghmdsr40
ghmdsr40 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 • 783 views
ADD COMMENTlink modified 2.4 years ago by Zhilong Jia1.6k • written 2.4 years ago by ghmdsr40

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 2.4 years ago by russhh5.5k
0
gravatar for Zhilong Jia
2.4 years ago by
Zhilong Jia1.6k
London
Zhilong Jia1.6k 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 2.4 years ago by Zhilong Jia1.6k
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: 838 users visited in the last hour