def testCredibleIntervalWholeNumbers(self):
simulations = np.arange(1, 101)
level = 0.9
expected = np.array([5.0, 50.0, 95.0])
obtained = utils.credible_interval(simulations, level)
np.testing.assert_array_almost_equal(expected, obtained)
def testCredibleIntervalInterpolation(self):
simulations = np.arange(1, 101)
level = 0.88
expected = np.array([6.0, 50.0, 94.0])
obtained = utils.credible_interval(simulations, level)
np.testing.assert_array_almost_equal(expected, obtained)
def testCredibleIntervalRaisesOnLargeLevel(self):
simulations = np.arange(1, 101)
level = 0.999
with self.assertRaises(ValueError):
utils.credible_interval(simulations, level)