class TestSkewGaussian(object):
def setup(self):
self.skewGassian = SkewGaussian()
np.random.seed(seed=42)
def test_pdf(self):
x = 1
y = self.skewGassian.pdf(x, e=0., w=1., a=0.)
assert y == 0.24197072451914337
x = np.array([0, 1])
y = self.skewGassian.pdf(x, e=0., w=1., a=0.)
assert y[0] == 0.3989422804014327
assert y[1] == 0.24197072451914337
def test_pdf_skew(self):
x = 1
y = self.skewGassian.pdf_skew(x, mu=1, sigma=1, skw=0.5)
assert y == 0.39834240320473779
y = self.skewGassian.pdf_skew(x, mu=1, sigma=1, skw=-0.5)
assert y == 0.39834240320473779
def test_raise(self):
with self.assertRaises(ValueError):
skewGassian = SkewGaussian()
skewGassian.pdf_skew(x=1, mu=1, sigma=1, skw=-1)
with self.assertRaises(ValueError):
prob_density.compute_lower_upper_errors(sample=None, num_sigma=4)