neurorestore/Libra

How to perform edgeR-LRT using Libra with no biological replicates?

Opened this issue · 2 comments

DE = run_de(snRNA, de_family = 'pseudobulk', de_method = 'edgeR', de_type = 'LRT',cell_type_col = "orig.ident")

This returns:
Error in mutate_at():
! Can't subset columns that don't exist.
✖ Column replicate doesn't exist.

So how can I perfrom edgeR-LRT using Libra with no biological replicates?

If you do not have biological replicates you would typically use a single-cell method (e.g., Wilcox), since pseudobulk methods require replicates. Alternatively, you could treat cells as replicates (cell_type_col = barcode) - however this will not protect you against the biases we describe in the paper.

I tried to run a single-cell method (wilcox) on a Seurat dataset (specifically ifnb, it is a PMBC dataset with 2 treatments, no biological replicates), and I also received an error "Column replicate doesn't exist". Should this be working by default?