This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem with Deseq2

Hello every one I'm trying to obtain differentially expressed genes in galaxy using Deseq2. But unfortunately, I get the following error. I need to mention that my data is pair end and the stringtie results are in a list.

Fatal error: An undefined error occurred, please check your input carefully and contact your administrator.
Tool generated the following standard error:

estimating size factors
Error in estimateSizeFactorsForMatrix(counts(object), locfunc = locfunc,  : 
  every gene contains at least one zero, cannot compute log geometric means
Calls: DESeq ... estimateSizeFactors -> .local -> estimateSizeFactorsForMatrix
galaxy deseq2 sequencing rna-seq
• 3,916 views • link
updated by Lila M Lila M • written by Yoosef • Yoosef

every gene contains at least one zero, cannot compute log geometric means

How can this be, so what kind of data is this?

• 0 views • link
written by ATpoint ATpoint

How do you generate your countData? Maybe the problem is in the upstream steps, because the matrix seems to have no counts for your genes.

• 0 views • link
written by Lila M Lila M

1 answer

It seems that you have data with zero counts. Filter to get rid of them. DESeq2 best practices recommend a filtering of low counts reads before the actual analysis

• 0 views • link
written by Antonio R. Franco Antonio R. Franco

I do not think filtering is a general best practice. From the vignette:

While it is not necessary to pre-filter low count genes before running the DESeq2 functions, there are two reasons which make pre-filtering useful: by removing rows in which there are very few reads, we reduce the memory size of the dds data object, and we increase the speed of the transformation and testing functions within DESeq2.

https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#pre-filtering

In any case it will not help here I guess as every gene contains a zero so removal of lowly-expressed genes will not change that. So either one or many corrupted samples or unusual data. @OP, please elaborate.

• 0 views • link
written by ATpoint ATpoint

Log in to answer this question.

More posts like this