def test_clipped_std_bias(self): '''Compare to www.wolframalpha.com integrals''' self.assertAlmostEqual(_clipped_std_bias(1), 0.53956, places=5) self.assertAlmostEqual(_clipped_std_bias(2), 0.879626, places=6) self.assertAlmostEqual(_clipped_std_bias(3), 0.986578, places=6) np.random.seed(1) x = np.random.normal(size=1000000) biased_std = np.std(x[np.abs(x) < 3]) self.assertAlmostEqual(biased_std, _clipped_std_bias(3), places=3)
def test_clipped_std_bias(self): '''Compare to www.wolframalpha.com integrals''' self.assertAlmostEqual(_clipped_std_bias(1), 0.53956, places=5) self.assertAlmostEqual(_clipped_std_bias(2), 0.879626, places=6) self.assertAlmostEqual(_clipped_std_bias(3), 0.986578, places=6) np.random.seed(1) x = np.random.normal(size=1000000) biased_std = np.std(x[np.abs(x)<3]) self.assertAlmostEqual(biased_std, _clipped_std_bias(3), places=3)