Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 error in data.frame (multiple treatments and multiple replicates) KYR Bioinformatics 3 03-19-2019 10:05 AM
DESeq2 Contrasts cacti Bioinformatics 4 08-04-2015 12:20 PM
[DESeq2] Multiple replicates, multiple treatments KYR Bioinformatics 4 08-02-2014 07:01 AM
Contrasts with DESeq2 mistrm RNA Sequencing 3 07-24-2014 10:03 AM
Differential expression analysis on RPKMs - Contrasts and contrasts of contrasts giorgifm Bioinformatics 6 08-16-2013 11:50 AM

Thread Tools
Old 08-05-2014, 11:28 AM   #1
Junior Member
Location: New Jersey

Join Date: Aug 2014
Posts: 2
Default DESeq2 - using contrasts vs. multiple runs for DE analysis

Hello all,

I have a time series dataset with 4 time points (no other factors besides time). What is the better, more statistically sound strategy to test genes in each time point against one another for DE:

1) Run DESeq on the entire data set (4 time points * 3 replicates each) and extract each contrast using results(dds, contrast = c("time",x,y) for a total of 6 results tables ( 1 v 2, 1 v 3, 1 v 4, 2 v 3, 2 v 4, 3 v 4)

2) Run DESeq six times on six different data sets, with 2 time points in each, and use the results table from each run.

I plan on extracting all genes that have at least LFC of 1 and padj < 0.05 . Applying this filter on a result table from (2) gives me less genes than on the corresponding table in (1). Does this mean that running the analysis individually for each comparison is more stringent?

Any help would be appreciated.
rutgers2015 is offline   Reply With Quote
Old 08-05-2014, 11:42 AM   #2
Devon Ryan
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480

Option 1 is generally the better idea. Option 2 isn't more stringent, you just have fewer total samples so that'll affect your variance estimates and degrees of freedom.
dpryan is offline   Reply With Quote
Old 07-14-2015, 11:34 AM   #3
Junior Member
Location: UK

Join Date: Jan 2015
Posts: 9


I am in a similar situation with my data.

I have three conditions, with 4 replicates for CondA, 2 replicates for CondB, and 2 replicates for CondC.

When I run DESeq using all 8 samples, I get 197 DE genes with a padj < 0.1 for contrast CondB/CondA.

However, when I take out the samples from CondC from the counts matrix, for the same contrast CondB/CondA, I get only 170 DE genes.

So in my case, the second method is more stringent... although I would like to understand why this happens.

Thanks for your help.
friducha is offline   Reply With Quote

deseq2, differential anals

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 09:40 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO