SEQanswers

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

Reply
 
Thread Tools
Old 08-05-2014, 11:28 AM   #1
rutgers2015
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
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

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
friducha
Junior Member
 
Location: UK

Join Date: Jan 2015
Posts: 9
Default

Hi,

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
Reply

Tags
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