Implements Poisson fixed effects regression and robust standard errors from Wooldridge (1999) in Julia. The implementation should be faster than GLM.jl with many fixed effects.