Driver genes

[p-gueguen]

Hello,

Thanks for this great package. I was wondering if there is a way to find which genes are most contributing to the each celltype annotation during RCTD scoring? I could not find any function to do it directly.

Thanks a lot!

[dmcable]

Good question!

It is possible to get such scores out of RCTD, but it is not necessarily something that is stored by RCTD by default.

The relevant score would be the log likelihood of each gene for each cell type. It varies across RCTD mode. For example, for doublet mode, you would want to look at the log likelihood vector for each cell type pair.

The relevant functions that compute this (for doublet mode) would be decompose_sparse and calc_log_l_vec.

I hope this helps!

Dylan