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

Question: DESeq2 paired design with missing samples and multiple conditions
0
gravatar for greg.deakin
12 weeks ago by
greg.deakin0 wrote:

Hi,

I've got a paired (or blocks really) experiment with the setup like the below.

Five blocks and four conditions (A, B, C, D). I want to contrast all vs all (A vs B, A vs C, A vs D and etc.). The Xs in the table below indicate a sample present for each block (I think it's clearer in table format rather than colData format). Clearly there are a lot of "missing pairs", but there are at least two paired samples when comparing any of the conditions.  

Block A B C D
1   x   x
2 x x    
3 x   x x
4   x x  
5 x x x x

 

 

 

 

 

So, my question is what would be the best way to do an analysis with DESeq2?

As far as I can see, it's either keep all the data and use the design ~Block + condition then extract each contrast from the results object, or subset the data  and keep the design balanced. It's noted that block 3 and 5 (-B5) cover half the contrasts.

Reading previous forum posts the consensus seems to be subset the data or use limma-voom, but the posts I've read don't deal with multiple conditions.  

Greg

ADD COMMENTlink modified 12 weeks ago by Michael Love22k • written 12 weeks ago by greg.deakin0
Answer: DESeq2 paired design with missing samples and multiple conditions
0
gravatar for Michael Love
12 weeks ago by
Michael Love22k
United States
Michael Love22k wrote:

The case where I recommend that users use duplicateCorrelation() is a little different. That is where it is desired to estimate condition effects for different groups of subjects, but a subset of the subjects only have one level of condition,. Here there is not the same nesting of subjects within biologically-relevant groups. Here you want to estimate condition effects, which are assumed to be the same regardless of block, while controlling for block effects. You should be able to just use ~block + condition.

ADD COMMENTlink written 12 weeks ago by Michael Love22k
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: 182 users visited in the last hour