Error and warning in AUCell_exploreThresholds

[rythi4]

Hi,
When I run this code I get the following error-message, but I cannot seem to fix it. Does anyone know what it could mean?

cells_assignment <- AUCell_exploreThresholds(cells_AUC, plotHist=TRUE, assignCells=TRUE)

Error in if (any(FCs > 1.5)) warning(gSetName, ":\tCheck the AUC histogram. ", :
missing value where TRUE/FALSE needed

I see the first 100 histograms, but then the error appears.

[s-aibar]

Hello,

In order to check what might be going on we would need some more information. Ideally, an example to reproduce the error.

Please, also mention which AUCell version you are using.

Thank you!

[dagarfield]

I run into this too from time to time -- never been able to replicate it and it appears to be stochastic (run it again and it vanishes).

Is there a stochastic step in the AUCell_exploreThresholds? Maybe something with mixtools? If I had to guess based on my own observations, you have a score with a sort of flat distribution and/or only very low scoring cells and that crashes things. What's a heatmap of your 101st gene set look like?