feglm.nb function

Efficiently fit negative binomial glm's with high-dimensional $k$-way fixed effects