Question: DESeq2 - Effect of DDS Import of All vs Select Conditions
0
12 months ago by
pvd21070
pvd21070 wrote:

Hi Everyone,

We've recently completed an RNASeq experiment with 35 samples and 10 different conditions. Almost all of our biological questions are about pairwise comparisons between two of these conditions.

As of now, my workflow has been to import all of the samples as follows:

dds <- DESeqDataSetFromTximport(txi.salmon, sampleTable, ~condition)
dds <- dds[ rowSums(counts(dds)) > 10, ]

nrow(dds) -> 33,270

After I've done this, I've ran DESeq on the enitre dds dataset:

dds <- DESeq(dds)

I've then proceeded to run results on the pairwise comparison of interest:

YCNT <- results(dds, contrast = c("condition", "YB6CNT", "YBJCNT"), alpha = 0.05, lfcThreshold = 1)

In this instance, I find very few significant differences using the log2 and padj thresholds (which is not that surprising since this is a control condition comparison)

> table(YCNT.05\$padj < 0.05)

FALSE  TRUE
33072    24 

My question is, if I'm interested in just this comparison, is there any difference in how the model runs if I were to only import the samples corresponding to these two conditions to begin with, as opposed to importing all 35 samples and then running the pairwise results contrast?

For example, I could imagine that the initial DDS filtering would return less rows with rowsum >=10 and that might make the FDR less strict. I'm not sure if that is accounted for in the model, but was just curious if one approach is more appropriate than the other.

Thanks a lot.

Paul

deseq2 • 190 views
ADD COMMENTlink
modified 12 months ago by Michael Love22k • written 12 months ago by pvd21070
Answer: DESeq2 - Effect of DDS Import of All vs Select Conditions
1
12 months ago by
Michael Love22k
United States
Michael Love22k wrote:

This is one of the frequently asked questions in the DESeq2 vignette.

ADD COMMENTlink written 12 months ago by Michael Love22k

sorry about that, thanks!

ADD REPLYlink written 12 months ago by pvd21070
Please log in to add an answer.

Content
Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 129 users visited in the last hour