Question: DESeq2 VST question
gravatar for MarcHiggins
3.1 years ago by
MarcHiggins10 wrote:

Hello everyone,

I am using the DESeq2 command variance stabilizing transformation to normalise raw read counts across 4 RNA-Seq experiments. Once normalised I will use these VST counts in Matrix eQTL, however I am having an issue with VST.

After carrying out VST, some of my results are genes have higher VST values than others which had higher raw read counts, i.e. re-ranking has occurred for some samples. Would anyone be able to shed some light on why this might have happened?


deseq2 • 2.2k views
ADD COMMENTlink written 3.1 years ago by MarcHiggins10
gravatar for Kevin Blighe
3.1 years ago by
Kevin Blighe71k
Republic of Ireland
Kevin Blighe71k wrote:

This is most likely due to the adjustment for library size, which occurs during normalisation. For example, the samples that originally had higher raw counts may just have been sequenced to a higher overall read depth - through their adjustment for library size in relation to all other samples during normalisation, it was probably then found that these higher raw counts did not actually equate to higher expression.

For recommended reading, I would look at:

ADD COMMENTlink written 3.1 years ago by Kevin Blighe71k

Hi Kevin,

Many thanks for your reply. After contemplating this further myself I agree that this is the case.


ADD REPLYlink written 3.1 years ago by MarcHiggins10
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.3.0
Traffic: 1888 users visited in the last hour