This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Wilcox test using data slot (log normalized count values) of scRNA-seq between two groups

Hi,

I was wondering would one be able to do a wilcox test for a specific gene in the same cell type in two different conditions using the log normalized values in seurat (the data slot of the RNA assay)? For example, if I have the gene EGFR, and I want to see if there is a difference in expression of EGFR in epithelial cell of a control sample/samples vs the epithelial cells of a treatment sample/samples. I know I can run edgeR or DESeq2 for differential expression for full genes but would a wilcox test also work, especially if I want to test the difference for a single gene? Appreciate any insight.

scrnaseq
• 1,436 views • link
written by mropri mropri

1 answer

The best way is to generate a pseudo-bulk aggregation then apply edgeR or DESeq2. You can also directly use wilcoxon on sc-level, since many people are doing this, but be aware of possible false-positive results.

You can find details in this post: Difference between GLM and Wilcoxon rank sum test

• 0 views • link
written by yxwucq • yxwucq

Thank you, will do Wilcox for now as I only have one control and one treatment, as I get more control and treatment samples, will stick with pseudobulk and DESeq2.

• 0 views • link
written by mropri mropri

Log in to answer this question.

More posts like this