def pdf(self, X):
        X = np.asarray(X)
        return np.exp(-X)/(1+np.exp(-X))**2
The logistic probability density function
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/discrete/discrete_model.py/Logit.pdf