def score_obs(self, params):

        y = self.endog
        X = self.exog
        L = self.cdf(np.dot(X, params))
        return (y - L)[:,None] * X
Logit model Jacobian of the log-likelihood for each observation
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/discrete/discrete_model.py/Logit.score_obs