示例#1
0
def compute_univ_cdf():
    eta, lam = 20, 1.5
    skst = SkStJR(ndim=1, eta=eta, lam=lam)
    print(skst.cdf(np.zeros(1)))

    eta, lam = 20, 1.5
    skst = SkStBL(ndim=1, eta=eta, lam=lam)
    print(skst.cdf(np.zeros(1)))

    eta, lam = 100, 1.5
    skst = SkStDM(ndim=1, eta=eta, lam=lam)
    print(skst.cdf(np.zeros(1) + 10))
示例#2
0
def compute_cdf():
    eta, lam = [20, 5], [1.5, .5]
    skst = SkStJR(ndim=len(lam), eta=eta, lam=lam)
    print(skst.cdf(np.zeros(2)))

    eta, lam = 20, [1.5, .5]
    skst = SkStBL(ndim=len(lam), eta=eta, lam=lam)
    print(skst.cdf(np.zeros(2)))

    eta, lam = 100, [1.5, -2]
    skst = SkStDM(ndim=len(lam), eta=eta, lam=lam)
    print(skst.cdf(np.zeros(2)))
示例#3
0
def compute_univ_cdf():
    eta, lam = 20, 1.5
    skst = SkStJR(ndim=1, eta=eta, lam=lam)
    print(skst.cdf(np.zeros(1)))

    eta, lam = 20, 1.5
    skst = SkStBL(ndim=1, eta=eta, lam=lam)
    print(skst.cdf(np.zeros(1)))

    eta, lam = 100, 1.5
    skst = SkStDM(ndim=1, eta=eta, lam=lam)
    print(skst.cdf(np.zeros(1)+10))
示例#4
0
def compute_cdf():
    eta, lam = [20, 5], [1.5, .5]
    skst = SkStJR(ndim=len(lam), eta=eta, lam=lam)
    print(skst.cdf(np.zeros(2)))

    eta, lam = 20, [1.5, .5]
    skst = SkStBL(ndim=len(lam), eta=eta, lam=lam)
    print(skst.cdf(np.zeros(2)))

    eta, lam = 100, [1.5, -2]
    skst = SkStDM(ndim=len(lam), eta=eta, lam=lam)
    print(skst.cdf(np.zeros(2)))
示例#5
0
    def test_cdf(self):
        """Test cdf."""

        eta, lam = 20, 1.5
        skst = SkStJR(ndim=1, eta=eta, lam=lam)
        cdf = skst.cdf(np.zeros(1))

        self.assertIsInstance(cdf, float)

        eta, lam = [20, 5], [1.5, .5]
        skst = SkStJR(ndim=len(lam), eta=eta, lam=lam)
        cdf = skst.cdf(np.zeros(2) - 10)

        self.assertIsInstance(cdf, float)
示例#6
0
def compute_quantile():
    eta, lam = 20, 1.5
    skst = SkStJR(ndim=1, eta=eta, lam=lam)
    cdf = skst.cdf(np.zeros(1) - 2)
    print(skst.ppf(cdf))

    eta, lam = 20, 1.5
    skst = SkStBL(ndim=1, eta=eta, lam=lam)
    cdf = skst.cdf(np.zeros(1) - 2)
    print(skst.ppf(cdf))

    eta, lam = 100, 1.5
    skst = SkStDM(ndim=1, eta=eta, lam=lam)
    cdf = skst.cdf(np.zeros(1) - 2)
    print(skst.ppf(cdf))
示例#7
0
def compute_quantile():
    eta, lam = 20, 1.5
    skst = SkStJR(ndim=1, eta=eta, lam=lam)
    cdf = skst.cdf(np.zeros(1) - 2)
    print(skst.ppf(cdf))

    eta, lam = 20, 1.5
    skst = SkStBL(ndim=1, eta=eta, lam=lam)
    cdf = skst.cdf(np.zeros(1) - 2)
    print(skst.ppf(cdf))

    eta, lam = 100, 1.5
    skst = SkStDM(ndim=1, eta=eta, lam=lam)
    cdf = skst.cdf(np.zeros(1) - 2)
    print(skst.ppf(cdf))
示例#8
0
    def test_quantile(self):
        """Test quantile."""

        eta, lam = 20, 1.5
        skst = SkStJR(ndim=1, eta=eta, lam=lam)
        arg = -2.
        cdf = skst.cdf(arg)
        ppf = skst.ppf(cdf)

        self.assertAlmostEqual(ppf, arg)

        arg = -.1 * np.ones(3)
        cdf = skst.cdf_vec(arg)
        quantiles = skst.ppf_vec(cdf)

        npt.assert_array_almost_equal(arg, quantiles)