PyDESeq2
PyDESeq2 copied to clipboard
BUG initialize `beta_init` when the design matrix is not full rank
Reference Issue or PRs
Fixes #280
What does your PR implement? Be specific.
Initialize beta_init in irls even when the design matrix is not full rank.
This prevents a referenced before assignment error in the case where irls is run with a rank-deficient design matrix and irls starts diverging.