def loglikeobs(self, params):

        q = 2*self.endog - 1
        X = self.exog
        return np.log(np.clip(self.cdf(q*np.dot(X,params)), FLOAT_EPS, 1))
Log-likelihood of probit model for each observation
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/discrete/discrete_model.py/Probit.loglikeobs