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
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.
Please check out the vignette and workflow sections on transformations. We cover both of these questions there.
OK I see, thank you, Michael.