search

owkin / PyDESeq2

A Python implementation of the DESeq2 pipeline for bulk RNA-seq DEA.
https://pydeseq2.readthedocs.io/en/latest/
MIT License
573 stars 60 forks source link

BUG initialize `beta_init` when the design matrix is not full rank #281

Closed BorisMuzellec closed 3 months ago

BorisMuzellec commented 4 months ago

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.