def probitloglike(params, endog, exog):
 q = ((2 * endog) - 1)
 X = exog
 return np.add.reduce(stats.norm.logcdf((q * np.dot(X, params))))

def probitloglike ( params , endog , exog ) : q = ( ( 2 * endog ) - 1 ) X = exog return np . add . reduce ( stats . norm . logcdf ( ( q * np . dot ( X , params ) ) ) )
log likelihood for the probit .