def test_qn_naive(self): assert_almost_equal(scale.qn_scale(self.normal), scale._qn_naive(self.normal), DECIMAL) assert_almost_equal(scale.qn_scale(self.range), scale._qn_naive(self.range), DECIMAL) assert_almost_equal(scale.qn_scale(self.exponential), scale._qn_naive(self.exponential), DECIMAL)
def test_qn_robustbase(self): # from R's robustbase with finite.corr = FALSE assert_almost_equal(scale.qn_scale(self.range), 13.3148, DECIMAL) assert_almost_equal(scale.qn_scale(self.stackloss), np.array([8.87656, 8.87656, 2.21914, 4.43828]), DECIMAL) # sunspot.year from datasets in R only goes up to 289 assert_almost_equal(scale.qn_scale(self.sunspot[0:289]), 33.50901, DECIMAL)
def test_qn_empty(self): empty = np.empty(0) assert np.isnan(scale.qn_scale(empty)) empty = np.empty((10, 100, 0)) assert_equal(scale.qn_scale(empty, axis=1), np.empty((10, 0))) empty = np.empty((100, 100, 0, 0)) assert_equal(scale.qn_scale(empty, axis=-1), np.empty((100, 100, 0))) empty = np.empty(shape=()) with pytest.raises(ValueError): scale.iqr(empty)
def test_qn(self): assert_almost_equal(scale.qn_scale(self.chem), 0.73231, DECIMAL)
def test_axisneg1(self): m = scale.qn_scale(self.X, axis=-1) assert_equal(m.shape, (40, 10))
def test_axis1(self): m = scale.qn_scale(self.X, axis=1) assert_equal(m.shape, (40, 30))
def test_axis0(self): m = scale.qn_scale(self.X, axis=0) assert_equal(m.shape, (10, 30))