Question: To analyze the overall variance, which one to choose, vst or DESeq2 normalized counts?
gravatar for Raymond
28 days ago by
Raymond0 wrote:

Hi, All,

    I have two groups of RNAseq data, each group have 5 samples. I suspected that drug treatment would increase the overall expression variances in some pathways. So, I planned to compare all gene expression variances in that pathway or even in the top 500 expressed genes. I could calculate the variance from DEseq2 normalized counts, or variance stabilization transformed counts (vst or rld). The question is: which one should I choose? Does the VST process affect the signal I want to analyze?


Thanks and best regards,


ADD COMMENTlink modified 25 days ago by Alan0 • written 28 days ago by Raymond0

Hi Michael,

May I ask:

1) Could you specify a little bit more about what are the some of reasons for favoring VST over rlog?

2) To clarify, the VST procedure is just for purpose of clustering and visualization, right? It (the counts after VST) basically is not involved in testing for pairwise DE with the Wald test?


Thank you very much in advance. 

ADD REPLYlink written 25 days ago by Alan0

 Please check out the vignette and workflow sections on transformations. We cover both of these questions there.

ADD REPLYlink written 25 days ago by Michael Love16k

OK I see, thank you, Michael.

ADD REPLYlink written 25 days ago by Alan0
gravatar for Michael Love
28 days ago by
Michael Love16k
United States
Michael Love16k wrote:

I recommend the vst(). It will help you compare variances across genes. You could also use rlog() I just prefer vst() these days.

ADD COMMENTlink written 28 days ago by Michael Love16k

Thanks so much, Michael!  For 'across genes', do you mean the variances of genes in the same group?  I want to compare the variances of the same genes between different groups. In this case, dose vst also help?

ADD REPLYlink modified 27 days ago • written 27 days ago by Raymond0
Please log in to add an answer.


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