GSEApy icon indicating copy to clipboard operation
GSEApy copied to clipboard

Getting precise FDR values and plotting if values are 0

Open ThomasThaewel opened this issue 1 year ago • 1 comments

Thank you so much for the great implementation, I just have a couple of questions here. If I use the standard gp.gsea function on my data (merged single cell dataset), I do get a lot of terms where FDR, NOM p-val and FWER are 0. I use permutation_num=1000 like suggested in the snRNA-seq example, I read in the closed issues that this seems to be an issue if used with ssGSEA but is this the same case here for multiple samples? The problem I run into then is that I can't use the suggested dot plot (gp.plot.dotplot(df_fin, column='FDR q-val', top_term=15, figsize=(4,8), cutoff=0.25)) since the FDR values of all top terms are 0. Is there a way to get more precise FDR values or enable plotting without manually changing the values? Thanks!

ThomasThaewel avatar Jun 04 '23 16:06 ThomasThaewel