DESeq2 residuals log transformation input for WGCNA
2
1
Entering edit mode
glmazzo ▴ 10
@glmazzo-12778
Last seen 4.5 years ago

I am working with RNA seq data. I would like to compute the residuals with DESeq2 to remove confounding effect in the data before running WGCNA.

 

We computed the residuals by using  DESeq2 on our RNA Seq count matrix with this code:

fitted.common.scale = t(t(assays(dds)[["mu"]])/sizeFactors(dds))

counts(dds, normalized=TRUE) - fitted.common.scale

 

WGCNA suggests to log2 transform count data (from RNA-seq) before the analysis.

We noticed that the residuals are computed in the scale of the counts.

Is it correct to log2 transform the residuals (computed as described above) before using them in WGCNA? Or this would be a mistake?

Thank you in advance for your help.

 

deseq2 residuals wgcna • 2.2k views
ADD COMMENT
0
Entering edit mode
@ryan-c-thompson-5618
Last seen 7 months ago
Scripps Research, La Jolla, CA

For any method that can't handle count data directly, you should follow the recommendation in the DESeq2 manual and use either the rlog or variance stabilizing transformations to get log2-scale data that is approximately homoskedastic.

ADD COMMENT
0
Entering edit mode
@steve-lianoglou-2771
Last seen 13 months ago
United States

If you want to remove batch effects, first shoot your data through DESeq2's vst, then use that log2 like matrix with limma's removeBatchEffects function

ADD COMMENT

Login before adding your answer.

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