Question: DESeq2 rlog as input for GAGE
0
gravatar for Tom
3.3 years ago by
Tom10
Tom10 wrote:

Dear all,

In the vignette of GAGE pathway analysis package, a joint workflow with DESeq2 is shown, where only the column of log2 fold change from DESeq2 output is used as an input for GAGE.

Is rlog matrix calculated by DESeq2 applicable as an input for GAGE pathway analysis?

Any comments would be appreciated.

Best,

Tom

gage deseq2 • 593 views
ADD COMMENTlink modified 3.3 years ago • written 3.3 years ago by Tom10
Answer: DESeq2 rlog as input for GAGE
0
gravatar for Luo Weijun
3.3 years ago by
Luo Weijun1.5k
United States
Luo Weijun1.5k wrote:
Not sure exact what rlog does in stabilizing variance. It may be safee to use just the normalized raw counts followed by log2 transformation, something like: cnts=count(dds, normalized=T) cnts.log=log2(cnts+1) cnts.log can then be used as gage input. Remember to specify ref and samp (column number vectors) though.
ADD COMMENTlink written 3.3 years ago by Luo Weijun1.5k
1

The idea of the vst() and rlog() functions in DESeq2 is to stabilize the variance of transformed counts across the range of values. These two functions are a motivated, statistical approach to transformation rather than picking either an arbitrary psuedocount and/or filter threshold for each dataset which will reduce the large noise component of log(count + x). 

vst() and rlog() return data which is log2-like, but the variance has been stabilized for low counts, and should be appropriate for downstream usage.

ADD REPLYlink modified 3.3 years ago • written 3.3 years ago by Michael Love26k
Answer: DESeq2 rlog as input for GAGE
0
gravatar for Tom
3.3 years ago by
Tom10
Tom10 wrote:

Luo and Michael,

Thank you for your helpful comments. I will use rlog transformed data as an input for downstream GAGE analysis.

Best,

Tom

ADD COMMENTlink written 3.3 years ago by Tom10
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: 410 users visited in the last hour