SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
DESEq2 seems to under report counts for SE reads Starr_Hazard Bioinformatics 3 11-19-2014 01:20 PM
How DESeq2 Handles counts TheSeqGeek General 4 04-24-2014 01:32 PM
Preparing text files with counts for DESeq2 KHubbard Bioinformatics 2 10-12-2013 01:42 AM
DESeq2: Difference between condition+type vs. 3 conditions mevers Bioinformatics 3 09-02-2013 01:17 AM
cuffdiff log2(FC) when 1 condition has 0 counts crh Bioinformatics 0 06-25-2012 12:08 PM

Reply
 
Thread Tools
Old 01-09-2015, 07:37 AM   #1
gmatt
Junior Member
 
Location: Saint Louis

Join Date: May 2014
Posts: 2
Default DESeq2 and zero counts in one condition

Does anyone know how DESeq2 handles genes that have zero counts in one condition and >0 counts in another?

My results output shows that these genes have a positive log2 fold-change value, but I do not understand how DESeq2 arrives at this number if it is taking the log of a ratio, in which the numerator is divided by zero.

Count data (letters are conditions; numbers are replicates):

Gene1:
A1 A2 B1 B2
0 0 692 1350

Results table:

Gene1:
baseMean log2FoldChange lfcSE stat pvalue padj
585.1394 11.78521 1.458913 8.078076 6.579644e-16 3.884571e-15


I know that DESeq1 gave an 'Inf' value in these cases, but how does DESeq2 arrive at a real number value?

Thanks
gmatt is offline   Reply With Quote
Old 01-09-2015, 10:32 AM   #2
Michael Love
Senior Member
 
Location: Boston

Join Date: Jul 2013
Posts: 333
Default

question answered on support.bioconductor.org
Michael Love is offline   Reply With Quote
Reply

Tags
deseq, deseq rnaseq, deseq2, zero counts

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 01:06 AM.


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