SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq2 multiple contrasts no "control" condition aprilw Bioinformatics 2 04-30-2015 10:21 AM
DESeq2 paired plus condition bioBob Bioinformatics 3 02-24-2015 08:23 AM
DESeq2 and zero counts in one condition gmatt RNA Sequencing 1 01-09-2015 11:32 AM
DESeq2: Difference between condition+type vs. 3 conditions mevers Bioinformatics 3 09-02-2013 02:17 AM

Reply
 
Thread Tools
Old 07-16-2015, 04:03 PM   #1
friducha
Junior Member
 
Location: UK

Join Date: Jan 2015
Posts: 9
Default What happened to per-condition and gene-est-only in DESeq2?

In DESeq, there was an option to estimate the gene dispersions within each condition instead of across all conditions, and per gene instead of sharing dispersion estimates across genes.

However, in DESeq2, I can't find these options anywhere -- how to specify these tweaks to the dispersion estimates now?

Thank you!
friducha is offline   Reply With Quote
Old 07-17-2015, 07:26 AM   #2
Michael Love
Senior Member
 
Location: Boston

Join Date: Jul 2013
Posts: 333
Default

In DESeq2 we only estimate a single dispersion value per gene. See the methods description in the paper:

http://genomebiology.com/2014/15/12/550
Michael Love is offline   Reply With Quote
Old 07-17-2015, 11:46 AM   #3
friducha
Junior Member
 
Location: UK

Join Date: Jan 2015
Posts: 9
Default

Thanks, Michael.

So it is clear to me now that the dispersion values are now per-gene, so gene-est-only no longer applies.

Quoting from the paper:

Quote:
Many methods for differential expression analysis of RNA-seq data perform such information sharing across genes for variance (or, equivalently, dispersion) estimation. edgeR[2],[3] moderates the dispersion estimate for each gene toward a common estimate across all genes, or toward a local estimate from genes with similar expression strength, using a weighted conditional likelihood. Our DESeq method [4] detects and corrects dispersion estimates that are too low through modeling of the dependence of the dispersion on the average expression strength over all samples.

What about estimating the dispersions within each condition instead of across all samples? (per-condition)?

Thanks.
friducha is offline   Reply With Quote
Reply

Tags
bioconductor, deseq, deseq2, rna-seq, rnaseq

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 04:32 PM.


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