To analyze the overall variance, which one to choose, vst or DESeq2 normalized counts?
1
0
Entering edit mode
Raymond ▴ 20
@raymond-14020
Last seen 5.5 years ago

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,

Raymond

deseq2 • 1.9k views
ADD COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode

OK I see, thank you, Michael.

ADD REPLY
1
Entering edit mode
@mikelove
Last seen 9 hours ago
United States

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

ADD COMMENT
0
Entering edit mode

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 REPLY

Login before adding your answer.

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