ENH: Glm add score_obs

so far just adds analytic (glm generic) score_factor and score_obs

score_factor is the same as score residuals in Stata.

I checked some examples against Stata, but unittests are against discrete and OLS/explicit

TODO: Hessian, observed and expected (information matrix) should be relatively easy to add, but not yet done.