Question: DESeq2 - Warning message
1
gravatar for madkitty
4.1 years ago by
madkitty580
Canada
madkitty580 wrote:
When comparing 13 controls and 32 samples, I get the following warning message in DESeq2. Is there a way to know which row had "a non positive estimates"? 
​
> dds <- DESeqDataSetFromMatrix(countData = as.matrix(data), colData=samples, design=~condition)
> dds <- DESeq(dds, betaPrior=FALSE)estimating size factors
estimating dispersions
gene-wise dispersion estimates
mean-dispersion relationship
final dispersion estimates
fitting model and testing
-- replacing outliers and refitting for 663 genes
-- DESeq argument 'minReplicatesForReplace' = 7 
-- original counts are preserved in counts(dds)
estimating dispersions
fitting model and testing

Warning messages:
1: In fitNbinomGLMs(objectNZ, maxit = maxit, useOptim = useOptim, useQR = useQR,  :
  1rows had non-positive estimates of variance for coefficients
2: In fitNbinomGLMs(objectNZ[fitidx, , drop = FALSE], alpha_hat = alpha_hat[fitidx],  :
  1rows had non-positive estimates of variance for coefficients

 

rna-seq deseq2 R software error • 1.5k views
ADD COMMENTlink modified 4.1 years ago by Michael Love1.8k • written 4.1 years ago by madkitty580
1
gravatar for Michael Love
4.1 years ago by
Michael Love1.8k
United States
Michael Love1.8k wrote:

hi,

Can you reproduce this bug using the release version (DESeq2 v1.8)? Typically the rows with GLM convergence issues are those with almost all zeros except for a few samples. It is valid to filter out these rows before running DESeq() as they don't have power for detection of differences anyway. For example: 

dds <- estimateSizeFactors(dds)
rmeans <- rowMeans(counts(dds, normalized=TRUE))
dds <- dds[rmeans > 1,]
ADD COMMENTlink written 4.1 years ago by Michael Love1.8k

I was unable to reproduce the error, it must have been smth out of my control. Thanks anyway

ADD REPLYlink written 4.1 years ago by madkitty580
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1613 users visited in the last hour