Feature Request: Post Estimation Score & InformationMatrix

[azev77]

It is often useful to have the score & information matrix of an estimated model (for example to do a score test).
Julia StatsBase allows StatsBase.informationmatrix(m) & StatsBase.score(m).
It would be awesome if GLM implemented these.

[kleinschmidt]

Sounds like a great first PR!

[alexfakos]

Is there a way to access the loglikelihood as a function of the parameters as each datapoint? Something like loglik_obs that maps from parameters and obs to a number.

[Tortar]

It would be also useful to have access to some other functions available in statsmodels in python: statsmodels.glmmodel.score_factor and statsmodels.glmmodel.hessian_factor