search

aertslab / AUCell

AUCell: score single cells with gene regulatory networks
110 stars 26 forks source link

Error and warning in AUCell_exploreThresholds #30

Open rythi4 opened 1 year ago

rythi4 commented 1 year ago

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 commented 1 year ago

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 commented 1 year ago

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?