Suppose I have RNAseq data for three conditions, A, B, and C. I can take all the data I have, build a sample table, feed it to DESeq, and use results(dds, contrast=) to compare each set of conditions.
Alternatively, I can take only the data for each pairwise comparison and run DESeq for each pair. .
I've tried this, and the results are fairly different between the two approaches. Only about 60% of significant genes overlap with my data.
How do I interpret this effect? Obviously, the data I'm not contrasting is still going to affect the linear model being created, which will affect the results table. But is it making the model better or worse? Which do I trust?