fgsea icon indicating copy to clipboard operation
fgsea copied to clipboard

fgsea hangs with nperm=1001, but not 1000

Open assaron opened this issue 3 years ago • 6 comments

# Runs fine, takes just a few seconds
fgseaRes <- fgsea(examplePathways, exampleRanks, nperm=1000, maxSize=500)

# Hangs and never completes
fgseaRes <- fgsea(examplePathways, exampleRanks, nperm=1001, maxSize=500)

It likely indicates that there are problems with BiocParallel package and it has to be reinstalled. In this case the following code hangs as well:

library(BiocParallel)
example("bplapply")

assaron avatar May 11 '21 11:05 assaron

@ChristianRohde please check if the BiocParallel example above works for you or not.

assaron avatar May 11 '21 11:05 assaron

I cannot reproduce the unexpected behavior neither on my installations on Ubuntu, MacOS or Windows any more. No need to add the nproc=1 any more as well. I could swear that I got stuck on this way too long yesterday, but now it is gone without a trace. I do not remember if I updated something on all 3 computers, but at least fgsea is now just working again as before very fast. Thank you for your support!

ChristianRohde avatar May 11 '21 13:05 ChristianRohde

Hi, I'm also having the same issue with fgsea running infinitely, seemingly out of nowhere since I had used the function successfully for at least a year before this. Setting nproc = 1 doesn't seem to help. The BiocParallel example above does work for me. I'm using the latest version of the fgsea package on MacOS. I (used to) rely heavily on this package so any advice would be really appreciated!

mvhunter1 avatar Jul 26 '21 20:07 mvhunter1

@mvhunter1 Is it the problem of some particular dataset or it hangs on different inputs?

assaron avatar Jul 27 '21 06:07 assaron

Hi, it hangs on any input I try from multiple datasets.

Edit: this seems to be specific to fgseaMultilevel as fgseaSimple seems to work fine.

mvhunter1 avatar Jul 27 '21 16:07 mvhunter1

@mvhunter1 Does it fail for the example dataset? Does setting nperm=1000 help? It looks like a different problem, not the one described in this issue, so please create another issue and describe exactly what works and what doesn't work for you providing the datasets.

assaron avatar Jul 28 '21 09:07 assaron