def test_scaling(): """Test the scaling function""" shape = (400, 10) u = np.random.randn(*shape) mean = 100 * np.random.rand(shape[1]) + 1 Y = u + mean Y_, mean_ = mean_scaling(Y) assert_almost_equal(Y_.mean(0), 0, 5) assert_almost_equal(mean_, mean, 0) assert Y.std() > 1
def test_scaling(): """Test the scaling function""" rng = np.random.RandomState(42) shape = (400, 10) u = rng.standard_normal(size=shape) mean = 100 * rng.uniform(size=shape[1]) + 1 Y = u + mean Y_, mean_ = mean_scaling(Y) assert_almost_equal(Y_.mean(0), 0, 5) assert_almost_equal(mean_, mean, 0) assert Y.std() > 1