def hessian(self, params):
        X = self.exog
        L = self.cdf(np.dot(X,params))
        return -np.dot(L*(1-L)*X.T,X)
Logit model Hessian matrix of the log-likelihood
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/discrete/discrete_model.py/Logit.hessian