Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Gene filtering for FGSEA #143

Open
SirKuikka opened this issue Nov 15, 2023 · 0 comments
Open

Gene filtering for FGSEA #143

SirKuikka opened this issue Nov 15, 2023 · 0 comments

Comments

@SirKuikka
Copy link

Hi,

I have scRNA-seq data and I am using log2 fold changes as the gene-level stats.

How should I filter the genes prior to FGSEA?

I have all 36 000 genes in a vector like this:

> head(fcs1)
MIR1302-2HG     FAM138A       OR4F5  AL627309.1  AL627309.3  AL627309.2 
 0.00000000  0.00000000  0.00000000  0.05287563  0.00000000  0.00000000 
> length(fcs1)
[1] 36601
> length(unique(fcs1))
[1] 14121

And the analysis warns of 3.8% of ties in the gene stats

> fgsea_results1 <- RunFGSEA(fcs1)
  |===========================================================================| 100%

Warning message:
In preparePathwaysAndStats(pathways, stats, minSize, maxSize, gseaParam,  :
  There are ties in the preranked stats (3.87% of the list).
The order of those tied genes will be arbitrary, which may produce unexpected results.
> 

Are those ties causing any issues? How would you filter the genes prior to GSEA?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant