def pdf(self, X):
        y = self.endog
        return np.exp(stats.poisson.logpmf(y, np.exp(X)))
Poisson model probability mass function
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/discrete/discrete_model.py/Poisson.pdf