Search
Question: DESeq: large difference in number of replicates per condition
0
gravatar for krc3004
6 weeks ago by
krc300410
krc300410 wrote:

Hi all,

I have a general question about the model estimation used by DESeq.  There have been many posts on whether or not DESeq works well with a small number/no replicates, but I'm wondering if it's appropriate to use DESeq for differential analysis across two conditions, where one condition has a small number of replicates (say, 5) and the other has a huge number (in the 100s).  The particular phenotype that we are looking at in our (clinical) data is quite rare, but we'd still like to test for differential expression.  Usually I would provide a reproducible example but these data are sensitive...

In this case, does it make sense to use DESeq?  Would it make sense to, say, randomly sample some of the replicates from the condition with 100s of replicates and run multiple tests?  I'm reading the original DESeq2 paper to try to understand how the model is built but any tips would be much appreciated.  Thank you!

ADD COMMENTlink modified 6 weeks ago by Michael Love14k • written 6 weeks ago by krc300410
1
gravatar for Michael Love
6 weeks ago by
Michael Love14k
United States
Michael Love14k wrote:

You shouldn't downsample the condition with 100s of replicates. There is nothing special to do, but you should note that the dispersion estimate will be mostly influenced from the group with more replicates. You can look at plotCounts afterward to verify that the top genes make sense, and are not affected by any artifact.

ADD COMMENTlink written 6 weeks ago by Michael Love14k

Michael, thanks very much for your help!  I will proceed without downsampling.

ADD REPLYlink written 6 weeks ago by krc300410

I would add that you should also pay attention to outlier filtering/replacement (see DESeq minReplicatesForReplace) if you expect your large group to be heterogenous. I've seen people run into trouble with that when handling large groups.

ADD REPLYlink modified 6 weeks ago • written 6 weeks ago by igor20
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 2.2.0
Traffic: 155 users visited in the last hour