alexdobin / STAR

RNA-seq aligner
MIT License
1.77k stars 497 forks source link

STARsolo knee-like cell filtering #2086

Closed mkazanov closed 3 months ago

mkazanov commented 3 months ago

Hi Alex,

According to the STARsolo documentation, there are two options for cell filtering: knee-like and EmptyDrop-like filtering. While the EmptyDrop method references Lun et al.'s paper, the knee-like approach is not described.

Could you please provide a reference for the knee-like method, if available, or offer a more detailed explanation?

Thank you in advance, Marat

alexdobin commented 3 months ago

Hi Marat,

This is a very simplistic calculation that replicates the old versions (<=2.2) of CellRanger. 1). Calculate robust maximum of UMI counts: default 99 percentile in the expected number of cells (3000 by default). 2). The UMI cutoff for cell calling is robust maximum divided by default 10.

These three parameters can be changed, the default values are: --soloCellFilter CellRanger2.2 3000 0.99 10

mkazanov commented 3 months ago

Thank you!