Search
Question: Issues with DEXSeq: rows did not converge in beta
0
gravatar for Konstantin Okonechnikov
6 days ago by

HI! 

I am using DEXSeq to compare exons within two conditions. The applied strategy was the same as in tutorial.

However, after running testForDEU(...) function the following messages were reported:

using supplied model matrix

5 rows did not converge in beta, labelled in mcols(object)$fullBetaConv. Use larger maxit argument with nbinomLRT

I tried to set maxit value, however this did not work:

Error in testForDEU(dxdFull, maxit = 500) : unused argument (maxit = 500)

Afterwards I tried to run next step and got many warning messages:

Fit for gene/exon ENSG00000241769.3+ENSG00000241489.3+ENSG00000010404.13+ENSG00000238039.1 threw the next warning(s): the matrix is either rank-deficient or indefinite

Why is this occurring? How is it possible to fix?

Also, I did not clean up the exon counts before running DEXSeq. Should I  perform some filtering in advance (i.e. exclude lowly expressed genes)? 

 

 

 

 

 

 

 

 

 

 

ADD COMMENTlink written 6 days ago by Konstantin Okonechnikov0

Hi Konstantin. The first one should be a warning and not an error, right? For those 5 exons you will get NAs instead of p-values. It is a problem of convergence of the GLMs. It is likely that these exons have very low counts or very dispersed counts. As for the estimateExonFoldChanges, its the same problem, do you get many NAs in the fold changes column?

ADD REPLYlink written 6 days ago by Alejandro Reyes1.4k
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.2.0
Traffic: 249 users visited in the last hour