Question: DESeq2 blocking by condition over time
0
gravatar for akp
3.2 years ago by
akp0
akp0 wrote:

I am using deseq2 to find differential significant taxa (count data). The data is longitudinal and three groups of mice and I am using sth like


diagdds = phyloseq_to_deseq2(physeq.fecal, ~ condition + week_experiment +
                             team + condition:week_experiment)

since I have three groups, I wonder how DEseq2 is calculating the fold change and p-value !? There is neither any warning nor any error ! should I slice the data 3 different combination of each group ?!

The second question I have is about condition:week_experiment term. I wonder how I can set it up in a way that two groups have the same initial value at the beginning (same intercept).

phyloseq deseq2 time course • 744 views
ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by akp0
Answer: DESeq2 blocking by condition over time
0
gravatar for Michael Love
3.2 years ago by
Michael Love26k
United States
Michael Love26k wrote:

I'd recommend that you partner with a statistician who can help set up a design that works for you and interpret results. There are in fact many different designs you could use which "work" in that they don't produce an error about rank-deficiency or linear dependency between columns, but that doesn't mean these are the right design. And furthermore, you will need to know how to properly interpret the estimated coefficients or build contrasts.

ADD COMMENTlink written 3.2 years ago by Michael Love26k
Answer: DESeq2 blocking by condition over time
0
gravatar for akp
3.2 years ago by
akp0
akp0 wrote:

I think, I haven't explained myself well, let me try one more time. Bacterial abundance of three different groups across 10 weeks is measure. Therefore, I chose this formulation/design to answer what are the taxa( in RNA world genes) that are distinguishing these three groups over the time.

diagdds = phyloseq_to_deseq2(physeq.fecal, ~ condition + week_experiment +
                             team + condition:week_experiment)

 In this setup, I wonder how does the fold change is calculated when there are three conditions, which I think I should look into the results function and it depends on the contrast

res = results(diagdds, contrast=c("condition","A","B"), cooksCutoff = FALSE)

 

I think I got my answer now

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by akp0

condition specifies the difference across three conditions at time=0,

week gives the difference across time for the reference condition

team accounts for differences in the team variable

and condition:week creates terms which provide for condition-specific differences at each time point after time=0.

You can use a likelihood ratio test, removing the interaction term, to find taxa where there are any condition-specific differences over time. We have a time series example in the workflow.

http://www.bioconductor.org/help/workflows/rnaseqGene/

ADD REPLYlink written 3.2 years ago by Michael Love26k
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: 194 users visited in the last hour