'return frozen scipy.stats distribution with mu at estimated prediction'

def predict_distribution(self, exog):

    if (not hasattr(self, result)):
        raise ValueError
    else:
        mu = np.exp(np.dot(exog, params))
        return stats.poisson(mu, loc=0)
