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)