diff --git a/hidimstat/adaptive_permutation_threshold.py b/hidimstat/adaptive_permutation_threshold.py index 66e525d..ae40f7b 100644 --- a/hidimstat/adaptive_permutation_threshold.py +++ b/hidimstat/adaptive_permutation_threshold.py @@ -49,7 +49,7 @@ def ada_svr(X, y, rcond=1e-3): return beta_hat, scale -def _manual_inversion(X, rcond=1e-3, full_rank=False): +def _manual_inversion(X, rcond=1e-3): "Inverting taking care of low eigenvalues to increase numerical stability" X = np.asarray(X) @@ -63,9 +63,6 @@ def _manual_inversion(X, rcond=1e-3, full_rank=False): s_inv = np.zeros(np.size(s)) s_inv[:rank] = 1 / s[:rank] - if full_rank: - s_inv[rank:] = 1 / (rcond * s.max()) - X_inv = np.linalg.multi_dot([U, np.diag(s_inv), V]) return X_inv