From 9c69a8fffec31cb201265c740ecdafd3fb90e8c7 Mon Sep 17 00:00:00 2001 From: kusch lionel Date: Tue, 7 Jan 2025 10:20:55 +0100 Subject: [PATCH] Change pvalue to evalue --- src/hidimstat/utils.py | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/hidimstat/utils.py b/src/hidimstat/utils.py index 0c1abc44..daf4e276 100644 --- a/src/hidimstat/utils.py +++ b/src/hidimstat/utils.py @@ -84,17 +84,17 @@ def _bhq_threshold(pvals, fdr=0.1): return -1.0 -def _ebh_threshold(pvals, fdr=0.1): +def _ebh_threshold(evals, fdr=0.1): """e-BH procedure for FDR control (see Wang and Ramdas 2021)""" - n_features = len(pvals) - pvals_sorted = -np.sort(-pvals) # sort in descending order + n_features = len(evals) + evals_sorted = -np.sort(-evals) # sort in descending order selected_index = 2 * n_features for i in range(n_features - 1, -1, -1): - if pvals_sorted[i] >= n_features / (fdr * (i + 1)): + if evals_sorted[i] >= n_features / (fdr * (i + 1)): selected_index = i break if selected_index <= n_features: - return pvals_sorted[selected_index] + return evals_sorted[selected_index] else: return np.inf