def getDistribution(data):
    kernel = stats.gaussian_kde(data)


class rv(stats.rv_continuous):

    def _cdf(self, x):
        return kernel.integrate_box_1d(-numpy.Inf, x)
Creating new distributions in scipy