Question: Should I use RUVseq on raw counts or on normalized counts coming from EDASeq package
0
gravatar for VeloSand
5 weeks ago by
VeloSand0
France
VeloSand0 wrote:

Hi, I'm working on an RNAseq experiment with spikes added and I would like to use RUVseq package. I had a look in the documentation and started with the very nice example on zebrafish. In this example RUVg is run on the object obtained after betweenLaneNormalization of the EDA package which is equivalent to runing RUVg on the normalized counts obtained by this normalisation method (I checked). Without this I would have run RUVg directly on the counts (either via a matrix or via an object). Is there any guidelines on whether we should run RUVg on counts or on "pre-normalized" counts? Thanks in advance for any help.

ruvseq • 63 views
ADD COMMENTlink modified 5 weeks ago by davide risso830 • written 5 weeks ago by VeloSand0
Answer: Should I use RUVseq on raw counts or on normalized counts coming from EDASeq pac
0
gravatar for davide risso
5 weeks ago by
davide risso830
University of Padova
davide risso830 wrote:

Hi,

in our experience, it is usually preferable to apply RUV to the normalized counts. However, that means that the assumption of the normalization method hold in your data.

In particular, usually normalizations assume that the majority of genes are not DE and/or that there is a roughly equal number of up- and down-regulated genes. If you think that these are reasonable assumptions in your experiment, we found that RUV on normalized data leads to slightly more robust results.

The fact that you are specifically mentioning spike-ins makes me think that perhaps you expect a very different amount of RNA among the conditions, or a lot of gene deregulation. If that's the case, you might be better off not normalizing the data and apply RUV directly on the raw counts (using the spike-ins as negative controls).

I hope this answers your question.

ADD COMMENTlink written 5 weeks ago by davide risso830

Thanks a lot for your answer. You are right I'm expecting a lot of gene deregulations for some of the tested conditions. I tested both approach (RUVg on raw counts and RUVg after EDA UQ normalization. Even in this case, RUVg after EDA seems to work better (even though it is not perfect).

ADD REPLYlink written 5 weeks ago by VeloSand0
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: 448 users visited in the last hour