The editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: non-integer counts in DESeq2
gravatar for warrena
8 months ago by
warrena0 wrote:

I have PRO-seq data with the following conditions: untreated (time 0), vehicle, and treatment (3 replicates each, 4 time points for each vehicle/treatment combination). I would like to evaluate treatment effects as a function of time. I have done the following:

1. Divided the reads in each condition by the corresponding DEseq size factor

2. Subtracted the vehicle mean from the treatment values at the corresponding time point (assuming additive effects)

3. Subtracted the min over the whole expression matrix (new min = 0)

I would like to use DESeq2 for DEG analysis, but clearly cannot enter the matrix of interest as counts. Would rounding be reasonable? Are alternative approaches advised?

ADD COMMENTlink modified 8 months ago by Michael Love22k • written 8 months ago by warrena0
Answer: non-integer counts in DESeq2
gravatar for Michael Love
8 months ago by
Michael Love22k
United States
Michael Love22k wrote:

Do not use subtracted values as input to DESeq2 or any count method.

You can account for baseline values using interaction models, and keep the data on the original count scale. See the interaction section of the DESeq2 vignette.

ADD COMMENTlink written 8 months ago by Michael Love22k
Please log in to add an answer.


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