feglm function

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