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


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.  


ADD COMMENTlink modified 14 days ago by Michael Love20k • written 14 days ago by greg.deakin0
gravatar for Michael Love
14 days ago by
Michael Love20k
United States
Michael Love20k 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 14 days ago by Michael Love20k
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.2.0
Traffic: 243 users visited in the last hour