Running DESeq2 on top variable genes?
1
0
Entering edit mode
hs.lansdell ▴ 20
@hslansdell-14246
Last seen 6.4 years ago

Since looking at the row variance and DESeq2 both act as ranking mechanisms for genes, is there any sense to taking the top 1000 or 5000 genes with the highest variance across samples from an RNA sequenced set and running the DESeq2 pipeline on that subset to look for differential genes between groups (so simple design ~condition)?

Thanks! 

deseq2 variance genes • 883 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 8 hours ago
United States

This will cause problems with DESeq2's dispersion prior, which should see counts from all the genes, or at least not subsetted by having a high or low sample variance (across all samples). Without getting into the details, it's not a problem when you subset by sample mean, but it would be a problem subsetting by sample variance.

ADD COMMENT

Login before adding your answer.

Traffic: 477 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