Error in py_call_impl about scicsr HOT 3 CLOSED

Hi,

This is not a technical (i.e. software dependency) issue - but a data-related one. Your error/warning messages suggest the neighbour graph constructed to describe the data is disconnected:

WARNING: Biased k-NN graph is disconnected
Error in py_call_impl(callable, dots$args, dots$keywords) :
ValueError: Input matrix is not weakly connected. Therefore it has no unique stationary distribution. Separate disconnected components and handle them separately

The warning/error messages are directly from CellRank (hence the py_call...). My guess is the dataset you selected has disparate sets of cells where (transcriptomically) they just are too distinct and therefore end up as separate 'clusters' and therefore the program struggles with predicting transitions between them. We did encounter this when we tested these on datasets and usually it is an indication of us needing to subset data ...

from scicsr.

Hello, I have solved this problem so far. I guess this may still be a technical problem. I speculate that some information may be lost during the conversion process using convertSeuratToH5ad. This problem does not occur when using other software packages that convert rds to h5ad (for example: scdior). Thank you very much for answering my question.

Hi,

This is not a technical (i.e. software dependency) issue - but a data-related one. Your error/warning messages suggest the neighbour graph constructed to describe the data is disconnected:

WARNING: Biased k-NN graph is disconnected
Error in py_call_impl(callable, dots$args, dots$keywords) :
ValueError: Input matrix is not weakly connected. Therefore it has no unique stationary distribution. Separate disconnected components and handle them separately

The warning/error messages are directly from CellRank (hence the py_call...). My guess is the dataset you selected has disparate sets of cells where (transcriptomically) they just are too distinct and therefore end up as separate 'clusters' and therefore the program struggles with predicting transitions between them. We did encounter this when we tested these on datasets and usually it is an indication of us needing to subset data ...

from scicsr.

Glad you have solved your problem - I have been looking into changing the engine behind conversion of the Seurat and AnnData objects - it seems there are better alternatives to SeuratDisk in doing this. Will add this as a To-Do. Thanks again for trying sciCSR!

from scicsr.