Wilcox test using data slot (log normalized count values) of scRNA-seq between two groups
1
0
Entering edit mode
8 weeks ago
mropri ▴ 150

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 • 326 views
ADD COMMENT
0
Entering edit mode
8 weeks ago
yxwucq • 0

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

ADD COMMENT
0
Entering edit mode

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.

ADD REPLY

Login before adding your answer.

Traffic: 2287 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6