Question: DESeq2 residuals log transformation input for WGCNA
0
gravatar for glmazzo
2.5 years ago by
glmazzo0
glmazzo0 wrote:

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 wgcna residuals • 1.1k views
ADD COMMENTlink modified 2.5 years ago by Steve Lianoglou12k • written 2.5 years ago by glmazzo0
Answer: DESeq2 residuals log transformation input for WGCNA
0
gravatar for Ryan C. Thompson
2.5 years ago by
The Scripps Research Institute, La Jolla, CA
Ryan C. Thompson7.4k wrote:

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 COMMENTlink modified 2.5 years ago • written 2.5 years ago by Ryan C. Thompson7.4k
Answer: DESeq2 residuals log transformation input for WGCNA
0
gravatar for Steve Lianoglou
2.5 years ago by
Denali
Steve Lianoglou12k wrote:

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 COMMENTlink written 2.5 years ago by Steve Lianoglou12k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 418 users visited in the last hour