Note: this post has been firstly posted on Biostars, Moderators suggested to post it here being this more appropriate for my purpose.
Hi, I am not going to ask how to solve issue related to the maxit argument with nbinomialWladTest within DESeq2 but for me, far from this scientific field, it's really hard to go through the vignette, other posts and lessons and trying to end up with a clear explanation of what is the convergence of beta and why using a larger maxit argument with nbinomWaldTest may solve the issue.
- I assume that first of all this may come up when the design of the experiment is not well balanced;
- Referring to this link Dr Love explains that when there is a single count in a row of 0s (so I guess when a certain gene in a sample is ≥1 where the same gene has 0 counts in every other samples) the GLM code may have some problems converging the betas. Now, I am a bit rusty with my stats knowledge but if I got this right: convergence of betas is required to fit a line in our experiment so that is possible to estimate dispersions of our parameters of interest. right? effect of X on Y in Y=βX+ε
Is that correct?
- Why using a larger maxit argument may solve this convergence problem? I presume that this step removes lowly expressed genes with low power so that makes the fit more reliable? Thanks