Hello World!
This is my first time posting here, or collaborating on a large software package, and I’m excited to be here!
About a year ago, I spent some time implementing the pdf of the Hypoexponential Distribution, which is a generalization of the Erlang distribution, which is already in Scipy. Mathematica has an implementation of the Hypoexponential, but it has numerical stability issues in certain cases (e.g. many rate parameters close to one another).
I created a Python implementation that is numerically stable. Since I couldn’t find any existing implementations of the Hypoexponential distribution online, I wanted to contribute mine to scipy, so folks don’t have to re-invent the wheel to work with it in the future.
Since I already have the code written, I am hoping it will be a light lift to add to scipy.stats, but I do need help as a first-time contributor to navigate the process. I created Issue #25639, and was pointed to ask about interest of adding this feature here.
Live Long and Prosper,
nonohuff