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)