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 02:20 PM
How DESeq2 Handles counts TheSeqGeek General 4 04-24-2014 02:32 PM
Preparing text files with counts for DESeq2 KHubbard Bioinformatics 2 10-12-2013 02:42 AM
DESeq2: Difference between condition+type vs. 3 conditions mevers Bioinformatics 3 09-02-2013 02:17 AM
cuffdiff log2(FC) when 1 condition has 0 counts crh Bioinformatics 0 06-25-2012 01:08 PM

Thread Tools
Old 01-09-2015, 08:37 AM   #1
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):

A1 A2 B1 B2
0 0 692 1350

Results table:

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?

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

Join Date: Jul 2013
Posts: 333

question answered on
Michael Love is offline   Reply With Quote

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:02 PM.

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