An efficient, conscise, and simple implementation of a purely on-disk B+...
Fit Gamma-Poisson Generalized Linear Models Reliably