Question: pre-ranked GSEA within R: options?
1
gravatar for Daniel Schmolze
5.2 years ago by
Daniel Schmolze20 wrote:
I want to do a GSEA entirely from within R, using genes ranked by my own metric. At the moment I'm saving my ranked genes to a .rnk file, then calling Broad's GSEA java program, then reading the resulting output back into R (all I care about are the p-values in gsea_report_for_na_pos_####.xls and gsea_report_for_na_neg_####.xls). Cumbersome to say the least. As far as I can tell, the Broad GSEA R script won't accept pre-ranked genes, but maybe I'm wrong? If not, I'm interested in other options. I'd like to specifically stick with the Broad GSEA algorithm if possible.
• 969 views
ADD COMMENTlink modified 8 months ago by atakanekiz20 • written 5.2 years ago by Daniel Schmolze20
Answer: pre-ranked GSEA within R: options?
0
gravatar for predeus
2.0 years ago by
predeus0
United States
predeus0 wrote:

use fgsea package on Bioconductor

ADD COMMENTlink written 2.0 years ago by predeus0
Answer: pre-ranked GSEA within R: options?
0
gravatar for atakanekiz
8 months ago by
atakanekiz20
atakanekiz20 wrote:

I've been playing with GSEA analyses in the last couple of days by using my single cell rnaseq data. I automated sample subsetting, gene ranking (S2N method, but can be replaced by other metrics), performing GSEA (using fgsea package) and reporting the top hits. If you are still interested in this type of analyses, I'm happy to share my code. Let me know.

ADD COMMENTlink written 8 months ago by atakanekiz20
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 16.09
Traffic: 426 users visited in the last hour