Size factors should be positive

[plger]

Size factors should be positive

[plger]

Hi,
The most likely explanation is that you have cells with a very low read count. Can you check min(colSums(counts(SCE))), where SCE is your singlecellexperiment object?

[Yale73]

Hi @plger ,

Thanks for your reply. The failed ones have min(colSums(counts(SCE))) = 0. I ran the QC first and then the scDblFinder. Now they worked.

Thanks again,
Yale

[plger]

Added this to the FAQs.
Note that we recommend not filtering too much before running scDblFinder (but removing the droplets not containing cells).

[Yale73]

Hi @plger,

If you recommend not filtering too much, why not recommend us to run scDblFinder first and then filter the cells? As I found if I ran the scDblFinder first and then filter, I will keep a little more high-quality cells.

Object with 4822 cells	Singlets	Doublets
First scDblFinder	4558	264
Then filter	4493
First filter	4717
Then scDblFinder	4475	242

Thanks,
Yale

[plger]

That's actually what I recommend (see for instance in the README, or in the pipeComp paper), however your earlier issue shows that some minimal filtering is needed to at least remove empty (or nearly empty) droplets...

[Yale73]

I see. I forget the filter is to remove the empty, If with empty, I will run the earlier issue.

Thanks,
Yale