<figure><img src="image_1.png" /></figure>

HAL
open science

# Multivariate filter methods for feature selection with the $\gamma$-metric

Nicolas Ngo, Pierre Michel, Roch Giorgi

► To cite this version:

Nicolas Ngo, Pierre Michel, Roch Giorgi. Multivariate filter methods for feature selection with the $\gamma$-metric. BMC Medical Research Methodology, 2024, 24 (1), pp.307. 10.1186/s12874-024-02426-9.
hal-04848056

HAL Id: hal-04848056

https://hal.science/hal-04848056v1

Submitted on 19 Dec 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# Multivariate filter methods for feature selection with the γ-metric

Nicolas Ngo<fnref n="1" />*, Pierre Michel<fnref n="2" />† and Roch Giorgi<fnref n="3" />†

**Abstract**

**Background** The γ-metric value is generally used as the importance score of a feature (or a set of features) in a classification context. This study aimed to go further by creating a new methodology for multivariate feature selection for classification, whereby the γ-metric is associated with a specific search direction (and therefore a specific stopping criterion). As three search directions are used, we effectively created three distinct methods.

**Methods** We assessed the performance of our new methodology through a simulation study, comparing them against more conventional methods. Classification performance indicators, number of selected features, stability and execution time were used to evaluate the performance of the methods. We also evaluated how well the proposed methodology selected relevant features for the detection of atrial fibrillation, which is a cardiac arrhythmia.

**Results** We found that in the simulation study as well as the detection of AF task, our methods were able to select informative features and maintain a good level of predictive performance; however in a case of strong correlation and large datasets, the γ-metric based methods were less efficient to exclude non-informative features.

**Conclusions** Results highlighted a good combination of both the forward search direction and the γ-metric as an evaluation function. However, using the backward search direction, the feature selection algorithm could fall into a local optima and can be improved.

**Keywords** Atrial fibrillation, Classification, Feature selection, γ-metric

**Introduction**

Adding features to classification models is not always beneficial to the problem at hand. In fact, some features may degrade the model's predictive performance. This study aims to propose a new methodology for

selecting relevant features in classification problems. The motivation for this work stems from the need to improve the detection of atrial fibrillation (AF) from electrocardiogram (ECG) data. AF is the most common cardiac arrhythmia, characterized by an irregular cardiac rhythm and an often rapid heart rate. Its prevalence increases with age [1], and it is associated with a significantly higher risk of stroke - up to five times greater [2, 3].

AF is typically detected through a 24-hour ECG recording, which is performed in a hospital or specialized facility. If confirmed, patients can be prescribed anticoagulants to reduce the risk of blood clot formation, and subsequently, the risk of stroke. However, AF is often asymptomatic and can be intermittent, with episodes lasting from a few minutes to a few days before the heart

<footnote>†Pierre Michel and Roch Giorgi contributed equally to this work.</footnote>

<footnote>*Correspondence:</footnote>

<footnote>Nicolas Ngo
nicolas.NGO@univ-amu.fr</footnote>

<footnote n="1">¹ Aix Marseille Univ, Inserm, IRD, SESSTIM, Sciences Économiques & Sociales de la Santé & Traitement de l'Information Médicale, ISSPAM, Marseille, France</footnote>

<footnote n="2">² Aix Marseille Univ, CNRS, AMSE, Aix-Marseille School of Economics, Marseille, France</footnote>

<footnote n="3">³ Aix Marseille Univ, APHM, Inserm, IRD, SESSTIM, ISSPAM, Hop Timone, BioSTIC, Marseille, France</footnote>
returns to normal sinus rhythm (NSR), sometimes for
longer periods than 24–48 hours. This irregularity makes
it challenging to detect AF unless the patient is continu-
ously monitored with ECG.

Emerging technologies, particularly mobile health (mHealth), are being developed to monitor patients non-invasively and continuously, facilitating long-term data collection for suspected AF cases [4–6]. These tools allow for the collection of physiological data, accurately reflecting cardiac activity. From this data, various features (e.g., the average time between two consecutive heartbeats) can be derived and used as inputs for predictive classification models to detect abnormal heart activity.

However, in real-world data analysis, it is often unclear
whether each extracted feature is relevant for classifying
heart rhythms (i.e., NSR vs AF). As the dimensionality of
a dataset increases (i.e., the number of features grows),
the risk of including irrelevant or redundant features also
increases. These extraneous features can reduce the pre-
dictive capacity of classification models. Feature selection
addresses this issue by reducing the dataset's dimension-
ality, which can enhance model performance, reduce
computation time, and simplify model complexity.
Moreover, by eliminating irrelevant features, the result-
ing models are better equipped to generalize to new data,
avoiding overfitting. In turn, such a model can facilitate
the interpretation of results by researchers, clinicians,
and decision-makers.

Feature selection is a statistical process that can be
applied to both regression and classification models. In
this study, we focused on distinguishing AF from NSR
using Holter-based ECG data and applied feature selec-
tion in the context of supervised classification models.
Feature selection methods are commonly grouped into
three categories [7]. The first is ‘filter methods’; these
methods are applied before model construction. They
evaluate the importance of individual features (or set
of features) based on predefined metrics. The second
category is ‘wrapper methods’; these methods rely on a
classification model to assess features and by evaluating
subsets of features using model performance criteria.
The last category is ‘embedded methods’; In these meth-
ods, feature selection is integrated into model-building
process. The feature selection process can be univari-
ate, where each feature is evaluated independently and
assigned an importance score. However, univariate meth-
ods do not account for potential correlations between
features. To address this, multivariate methods evaluate
subsets of features, considering the relationships among
them.

In 2017, a new metric called the γ-metric was proposed
by Pons et al. [8] to evaluate the discriminatory power
of numeric features in classification tasks. The authors

demonstrated that features with the highest γ-metric
values yielded the highest accuracy in univariate logistic regression models for AF detection. In a more recent
study, [9] examined the performance of the γ-metric as
a univariate filter for feature selection, comparing it with
three other univariate filters in the context of AF detection. The results showed that γ-metric produced comparable results to existing ranking methods, though it was univariate and only considered the relevance of individual features.

In the present study, we developed a new methodol-
ogy for multivariate filter feature selection in classifica-
tion, using the γ-metric as an evaluation function. The
γ-metric is associated with a specific search direction
(and therefore a specific stopping criterion). As three
search directions are used, we effectively created three
distinct methods. We assess the performance of these
methods through a simulation study and compare them
to seven conventional feature selection methods. Finally,
we applied all ten methods (the three novel and the seven
conventional) to the detection of AF using ECG data.

The paper is organized as follows: **Methods** section explains the computation of the γ-metric and its integration into three new feature selection methods, along with a description of the conventional feature selection methods. **Simulation study** section details the simulation study used to validate the three novel methods and compares their performance with existing methods. **Applica-tion** section presents the application of these methods to AF detection and their related results. Finally, **Discus-sions** section concludes the paper with a discussion of the simulation results and potential perspectives for using the γ-metric in classification models.

Methods

To select the best subset of features in a classification
model, one approach is to evaluate all possible combi-
nations of candidate features and then choose the sub-
set that optimizes an evaluation function. However, for
a dataset with *p* features, this requires evaluating 2^p - 1
combinations, which can be extremely time-consuming
and computationally expensive [10]. For example, with
only 20 features, over a million combinations must be
evaluated. Therefore, attempting to explore the entire
feature space is impractical. Instead, more efficient strat-
egies should be employed to explore the feature space in
a way that ensures a reasonable solution can be found
without evaluating every possible feature combination.
Various feature selection methods have been developed
based on different mathematical concepts and explora-
tion strategies, but the objective remains the same: to
identify the best subset of relevant features as closely as
possible [11]. Feature space exploration generally follows
four key steps [10]: (i) Defining a subset generation procedure, often referred to as search direction, (ii) selecting an evaluation function to assess feature relevance (in this study, the γ-metric is used), (iii) establishing a stopping criterion, and (iv) validating the feature selection process, (i.e., validating steps (i), (ii), and (iii)).

In step (i), the subset generation procedure defines how the feature space will be explored, specifically how features are added or removed from candidate subsets. This process, known as the search direction, determines how features are selected for evaluation. For example, in an exhaustive search, all possible subsets of features are evaluated.

Step (ii) involves the evaluation function, which measures the relevance of the feature (or subset of features) generated in the previous step. This function determines how important each feature or subset is for the model in question.

Step (iii) defines the stopping criterion, which indicates when the feature selection process should halt. Two methods with identical generation procedures and evaluation functions can yield different results if they employ different stopping criterion.

Finally, step (iv) is the validation of the feature selection process itself. This occurs outside the selection process and ensures that the chosen subset of features is indeed relevant and performs well in the selected model.

In the following sections, we describe the γ-metric, which is used as the evaluation function in step (ii). We then present three proposed multivariate filter feature selection methods, each employing a different search direction (step (i)) and using the γ-metric as an evaluation function (step (ii)) with a strict stopping criterion (step (iii)). Finally, we compare the performance of these methods with more conventional feature selection methods through a simulation study and by applying them to the task of discriminating between AF and NSR using Holter-based ECG data (step (iv)).

## γ-metric as an evaluation function

When developing the γ-metric for classification tasks, Pons et al. [8] main idea was to represent the K classes by p-dimensional ellipsoids, with p being the number of features. Each ellipsoid is characterised by the position of its center and the length and direction of its axes. The γ-metric computes the distance between the centers of each ellipsoid by taking into account the overlap of the ellipsoids. If an overlap exists, then the γ-metric will be negative; otherwise, it will be positive. In order to compute the γ-metric, we consider a set of n observations {Xᵢ}ᵢ=1,...,n with Xᵢ ∈ ℝp. These observations belong to one of the K classes, i.e., Yᵢ = 1, ..., K ∀i = 1, ..., n. The

first step is to group all observations from each class k and compute the associated class covariance matrix:

$$W_{k,p} = \text{Cov}(X_i|Y_i = k), \qquad (1)$$

Wₖ,ₚ is a diagonalizable, symmetrical positive semi-definite p × p matrix, with eigenvalues {λₖ,ₚ}ᵢ=1,...,p, and eigenvectors {uₖ,ₚ}ᵢ=1,...,p representing, respectively, the length and direction of the p axes of the ellipsoid. This ellipsoid is centered at μₖ, the mean vector of observations for class k. Hence, each class is represented by an ellipsoid, and the γ-metric represents the sum of the distances between each pair of ellipsoids. More specifically, for each pair of ellipsoids, the γ-metric represents the distance between the centroids minus the distance between the centroids and the borders of each ellipsoid; this ensures that any overlap of the ellipsoid is taken into account. For K = 2, there are only two ellipsoids, and the γ-metric is the distance between the two centroids minus the centroid-border distance of each. For K > 2, the distance between all pairs of ellipsoids is computed and then summed. More details on the computation of the γ-metric can be found in Appendix A.

In Eq. 1, the covariance matrix can be estimated either with the empirical estimator or the shrinkage estimator, depending on the ratio between features and observations. In reality, in the context of a large-scale dataset (i.e., p ≫ n) the empirical covariance matrix is a poor estimation of the covariance matrix, as it can lead to an ever-increasing number of null eigenvalues and to singularity of the covariance matrix [12]. This could pose a problem for the computation of the γ-metric, since the eigenvalues are used in the computation of the distances involved (see Appendix A for details on the computation of the γ-metric). Therefore, for the present study, we used the shrinkage estimation [13] of the covariance matrix when the number of observations for at least one class was smaller than the number of features.

## Multivariate feature selection method using the γ-metric

We used the γ-metric as an evaluation function to develop a novel multivariate feature selection methodology. Among all the subsets of candidate features, we looked for the one which maximized the γ-metric value. Specifically, we used three search directions (see below). Each search direction used a specific means to generate and explore the feature space. Therefore, each was separately associated with the γ-metric as an evaluation function. The three corresponding algorithms had the same strict stopping criterion. Hence there were three distinct feature selection methods using the γ-metric as an evaluation function. All three are described below:

Backward search (GAMMA_BACK). For this method, the value of the γ-metric is first computed for the whole
set of features. Next, the γ-metric value of all the possible subsets with one less feature is evaluated. If the value of the γ-metric of any given subset is strictly greater than that of the first candidate subset, then the former subset is retained; otherwise the algorithm stops and returns the first candidate subset. This process is repeated with the highest value being compared to the γ-metric value of the subset retained in the previous step. If the γ-metric value is strictly greater, the iterations continue either until there are no more features to remove or the γ-metric value cannot be increased by removing a feature of the candidate subset.

Forward search (GAMMA_FORW). In this method, the value of the γ-metric is first computed for each feature individually; the first candidate subset will be the feature with the highest value. Next, the γ-metric of all the possible subsets with one more feature is evaluated. If the γ-metric value of the given subset is strictly greater than that of the first candidate subset, then the former subset is retained; if not, the algorithm stops and returns the first set. The same process is then repeated: the γ-metric value of all the possible subsets with one more feature is evaluated; then the highest value is compared to the γ-metric value of the subset retained in the previous step. If the γ-metric value is strictly greater, the iterations continue either until there are no more features to add or the γ-metric value cannot be improved by adding a feature to the candidate subset.

Best first search (GAMMA_BF). With this method, the value of the γ-metric is first computed for each feature; the first candidate subset is formed by the feature with the highest γ-metric value. Next, the γ-metric of all the possible subsets with one more feature is evaluated. If the γ-metric value of a given subset is strictly greater than that of the first candidate subset, then the former subset is retained; if not the algorithm returns the first set. In the next step, the γ-metric value of all possible neighbour subsets, with one more feature, is evaluated and the highest value is then compared to the γ-metric value of the subset retained in the previous step. If the γ-metric value is strictly greater, the iterations continue; if not, the candidate subset is not directly returned. The best first search provides the possibility to go back to the second-best subset of features of the previous step and to continue the iterations with these poorer candidate features. In this way, this search direction is less likely to return a local maximum for the evaluation function. The iterations continue either until there are no more features to add to the candidate subset or the γ-metric value cannot be improved or the number of 'go back to a less optimal candidate subset' iterations reaches a maximum.

The simulation study and healthcare application was conducted using R to compute the γ-metric. The

shrinkage estimation was performed with the R package *corpcor* [14] and the *cov.shrink* function.

## Conventional feature selection methods

Several other feature selection methods have been proposed for supervised classification purposes which are more conventional in nature. In order to compare the performance of our novel γ-metric multivariate filter feature selection methods with these existing methods, we considered seven feature selection methods (4 filter methods, 2 wrapper methods and one embedded method). The general principle behind each of these methods is described below. Related technical details for some methods are provided in Appendix B.

Chi-squared filter (CHI2), this univariate filter method uses the Chi-squared statistic [15] to measure the dependence between the feature and the class. To do this, continuous features are discretized, the Chi-squared statistic is estimated, and Cramer's V is used as a ranking score for each feature. Cramer's V values close to 1 indicate a strong association between the feature and the class (see Appendix B: Chi-squared section).

Correlation-based feature selection (CFS) [16] is a multivariate filter feature selection method. It evaluates subsets of features on the basis of the hypothesis that "good feature subsets contain features highly correlated with the class, yet uncorrelated to each other" [16]. CFS ranks features according to an evaluation function based on correlations. This method assumes that irrelevant features have a low correlation with the class, and should be ignored. For a classification problem, CFS first discretizes numeric features and then entropy measure is used to estimate the degree of association between discrete features [16, 17]. The best-first search direction is used to generate the subsets to evaluate. The mathematical computation of the criterion used in CFS can be found in Appendix B: Correlation-based feature selection section.

Least absolute shrinkage and selection operator regression (LASSO) [18] is a multivariate embedded feature selection method. In a linear regression equation, the LASSO method adds a penalty term that discourages the model from assigning too much importance to any single feature. The penalty applied here is the L1 norm, which is the sum of the absolute values of the regression coefficients. This method allows some coefficients to be shrunk exactly to zero, effectively performing feature selection. A penalty parameter, which controls the strength of the regularization, is calibrated using a cross-validation.

Random forest importance (RFI) is a multivariate filter method. It is embedded in the random forest algorithm. More specifically, it computes the mean decrease accuracy (MDA) score for each feature in order to rank them [19]. This score describes how much accuracy the
model loses by permuting values of the feature. The idea is that for each feature, its score is computed by comparing the accuracy of the full prediction with the accuracy of the prediction when the feature values are randomly permuted. A high MDA value means that the permutations of the feature greatly impacted the accuracy, hence that the feature is important for a good accuracy. A low MDA value means that the values of the feature does not impact the prediction of the model. Once the score is computed for each feature, they can be ranked by order of importance. The features with the best ranks are then selected.

Stepwise AIC selection (STEP) is a multivariate wrapper method. The idea is to compare the AIC (a measure of the goodness-of-fit and complexity of a model) of the models by removing or adding features to the model. Start with the model with no features (only the intercept) and add the feature that reduces the AIC the most. Then repeat by adding or removing the feature that reduce the most the AIC. The methods is stopped when we cannot decrease the value of the AIC by removing or adding a feature.

Symmetrical uncertainty (SU) [20] is an univariate filter method based on the entropy and the information gain. This method is a variant of the mutual information [21] where we compute the entropy of all the data and the entropy of the candidate subset of features. In the mutual information, the entropy tends to be biased toward features with a large number of different values, in the SU method a normalization of the mutual information is applied to lower this bias (see Appendix B: Symmetrical uncertainty section).

Support Vector Machine Recursive Feature Elimination (SVM-RFE) is a multivariate wrapper method [22]. The method uses a recursive feature elimination searching procedure associated with a support vector machine classification model [23, 24]. The idea is to train an SVM model with all the features at first step. In the process of building an SVM model, a weights vector $v$ of the features is estimated. This weights vector is a linear combination of the training sample and the weight of each feature can be used as a ranking criteria of the features. The features with the smallest values are then eliminated and we train again the SVM with the remaining features.

To apply the feature selection methods CFS, CHI2, RFI and SU we used the implementation of the R package `FSelector` [17]. In addition, since CHI2, RFI and SU implementations return only a ranking of the features, with regard to their importance, we used a cutoff based on the biggest difference of importance score for each method. This way, the number of features selected was not fixed beforehand. For the γ-based feature selection methods (`GAMMA_BACK`, `GAMMA_FORW`, and

`GAMMA_BF`), the `FSelector` package offer the possibility to use a custom evaluation function with their implementation of the backward, forward and best first search direction functions. So we plugged in the γ-metric value function. For the LASSO feature selection method, we used the `glmnet` package [25, 26]. The step function from `stats` package was used to apply the stepwise AIC feature selection method. Last the SVM-RFE method was applied with the package `mlr3` [27].

## Simulation study

### Design of the simulation study and assessment of the feature selection methods' performances

To validate each of our three feature selection methods using the γ-metric as an evaluation function, we considered three distinct scenarios of binary classification problems. In these scenarios we explored the trade-off between number of observations and the number of features; class balancing and separability; and the effect of multicollinearity. Accordingly, we incorporated informative features with a fixed non-null effect and non-informative features with a null effect. In scenario 1, to assess whether each method selected informative features (and did not select non-informative features) in classical context, we considered situations where the number of observations was much higher than the number of features, with strong positive and negative effects. We also included a feature with a much lower effect. In scenario 2, we considered (i) more features than observations, (ii) both balanced and unbalanced classes, and (iii) situations where the two classes were strongly separated or not by the features. In scenario 3, we explored the efficiency of the methods in complex contexts with various pattern of correlations between the features, we considered (i) constant and non-constant levels of correlation, and (ii) three different levels of correlation. In all scenarios, we aimed to test whether the feature selection methods can guarantee good classification performance, select all the features that are truly informative and disregard all the non-informative features.

The classification performance was assessed using three criteria: (i) the area under the curve (AUC) which is an overall indicator of the model predictive performance, (ii) the sensitivity (iii) and specificity of the model. The feature selection process for each method was assessed by computing the average number of informative and non-informative features selected. The stability of the feature selection methods was assessed using the Jaccard index [28]. Additionally, we measured the execution time of the feature selection process, in seconds, of each method.
We applied three γ-metric based feature selection methods, plus the seven conventional methods presented above. The performance of the ‘FULL’ model, which uses all available features without any selection, is included as a reference for comparison. This model serves as a baseline to illustrate the impact of feature selection on predictive performance.

In each of the three scenarios for each of the eleven methods (the three γ-metric based methods, the seven conventional methods and the FULL method which returns all available features), we generated two datasets with specific parameters: a training dataset and a validation dataset. After applying a given feature selection method, the selected features were included in a logistic regression model developed from the training dataset. Using the validation dataset, we assessed the classification performance indicators. For each scenario, we repeated the generation of the training and validation datasets (both generated with the same sample size), the feature selection step, and the assessment of the classification performance 50 times. At each repetition of the simulation, we ensured that all the feature selection methods are applied on the same data (also assessment of the classification performances was done on the same validation dataset for all models). Stability was assessed using the Jaccard index [28] with pairwise comparison of the selected features over all the repetitions.

## Data generation

To generate data, we considered the following logistic equation, with $ρ_i$ being the probability that the observation i is a class 1 type:

$$ρ_i = \frac{\exp(X_iβ)}{1 + \exp(X_iβ)}, \quad \forall i \in 1, ..., n.$$

With $β$ the vector of effects, including an intercept $β_0$. We considered $p$ informative features associated with a non-null effect (i.e., $β_j \neq 0, \forall j \in \{1, ..., p\}$). We also added $p'$ features having a null-effect (i.e., $β_j = 0, \forall j \in \{p + 1, ..., p + p'\}$). X is thus the matrix of

$R^{n \times (p+p'+1)}$ where each column is a vector of $n$ realisations of a Gaussian distribution. The class of the observation i is then defined by $Y_i \sim Bernoulli(ρ_i)$. The different values of $β$ for both scenario 1 and 2 are reported in Table 1 (see Appendix D for details on the choice of the $β$ values). The generation of X described above applies for both scenario 1 and 2. For scenario 3, in order to generate data with multicollinearity, we used a multivariate Gaussian distribution with a variance-covariance matrix $Σ$. There was, in the generation process, groups of features, that were correlated together, but features from different groups were not correlated together. One group of features was generated with no correlation (within and with the other groups). All features had a variance of 1. Hence, the covariance matrix $Σ$ was a block diagonal matrix, with each block being the variance-covariance matrix, noted $Σ_g$, of a group of features. The construction of the covariance matrix was similar as the work of [29] with $Σ_g$ defined as:

$$Σ_g = \begin{pmatrix} 1 & α_{ij} & \cdots & α_{ij} \\ α_{ij} & 1 & \cdots & α_{ij} \\ \vdots & \vdots & \ddots & \vdots \\ α_{ij} & α_{ij} & \cdots & 1 \end{pmatrix}, \forall i, j = 1, ..., s_g.$$

With $s_g$ the size of each group, and $α_{ij}$ the correlation between feature i and feature j. We considered constant and non-constant levels of correlation within a group. For a constant level of correlation, $Σ_g$ is a matrix with 1 on its diagonal and $α_{max}$ as off-diagonal elements:

$$α_{ij} = \begin{cases} α_{max}, & \forall i \neq j \\ 1, & \text{otherwise.} \end{cases}$$

For a non-constant level of correlation, the level of correlation within a group will depend on the index of features. The correlation between feature i and feature j within a group will be affected by the value of $|i - j|$:

$$α_{ij} = \begin{cases} α_{max} \exp\{-w(|i - j| - 1)\}, & \forall i \neq j \\ 1, & \text{otherwise.} \end{cases} \quad (2)$$

<table><caption>**Table 1** Values of $β$ for the data generation process in scenario 1 and scenario 2</caption><thead><tr><th></th><th>β<sub>0</sub></th><th>β<sub>1</sub></th><th>β<sub>2</sub></th><th>β<sub>3</sub></th><th>p'</th><th>n</th></tr></thead><tbody><tr><td>Scenario 1</td><td>0.00</td><td>3.00</td><td>-2.00</td><td>0.50</td><td>22</td><td>2000</td></tr><tr><td>Scenario 2</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Unbalanced/strong</td><td>-2.65</td><td>3.60</td><td>-2.20</td><td>-1.00</td><td>197</td><td>100</td></tr><tr><td>Balanced/strong</td><td>0.00</td><td>3.60</td><td>-4.00</td><td>-1.00</td><td>197</td><td>100</td></tr><tr><td>Unbalanced/weak</td><td>-2.75</td><td>0.60</td><td>-2.50</td><td>-1.00</td><td>197</td><td>100</td></tr><tr><td>Balanced/weak</td><td>0.50</td><td>0.60</td><td>-2.50</td><td>-1.00</td><td>197</td><td>100</td></tr></tbody></table>

Values of $β$ associated with non-informative features are not reported since they were all null. Last two columns are $p'$ number of non-informative features and $n$ number of observations
The correlation level $\alpha_{ij}$ will decrease in (2), if the value $|i - j|$ increase. For instance, in a group, two consecutive features will have the maximum correlation level ($|i - j| - 1 = 0$ and $\alpha_{ij} = \alpha_{max}$). The first and the last feature of a group will have the minimum correlation level ($|i - j| - 1 = s_g - 2$ and $\alpha_{ij} = \alpha_{max} \exp\{-w(s_g - 2)\}$). Values of the parameter $w$ allowed us to control this minimum level of correlation $c$ between two features within a group:

$$ \alpha_{\max} \exp\{-w(|i - j| - 1)\} > c, $$

if we wish to achieve at least a correlation of $c$ between the first and last feature of a group, we set $w$ as

$$ w < \frac{\log(\alpha_{\max}/c)}{s_g - 2}. $$

The group of features independent from all other features had the identity matrix as a variance-covariance matrix. $\Sigma$ had the following format:

$$ \Sigma = \begin{pmatrix} \Sigma_g & 0 & \cdots & 0 & 0 & 0 \\ 0 & \Sigma_g & \cdots & 0 & 0 & 0 \\ \vdots & \vdots & \ddots & \vdots & \vdots & \vdots \\ 0 & 0 & \cdots & \Sigma_g & 0 & 0 \\ 0 & 0 & \cdots & 0 & \Sigma_g & 0 \\ 0 & 0 & \cdots & 0 & 0 & I_{sg \times s_g} \end{pmatrix}. \quad (3) $$

In scenario 1, we considered $n = 2000$ observations, and $p' = 22$ non-informative features. In scenario 2, we considered $n = 100$ observations, and $p' = 197$ non-informative features. We also considered different fixed effects $\beta$ for the informative features in order to generate classes with a strong or weak separation. Furthermore, we generated balanced and unbalanced classes. In scenario 1 and scenario 2, $x_1, x_2$, and $x_3$ are the only informative features. Table 1 summarises the previous information for scenario 1 and scenario 2, with the values of $\beta$ used in the generation process.

In scenario 3, we considered $n = 2000$ observations and the same value of $\beta$ fixed at 1.5, for all informative features. The intercept $\beta_0$ was set to 0. The datasets were generated with 10 groups of 10 features. In the first 5 groups, the first feature was informative and the others were non-informative. Groups 6, 7, 8, and 9 had only non-informative features and group 10 was the independent group with the first feature being informative and the others non-informative. Hence, in scenario 3, $p' = 94$ and $p = 6$. We considered $\Sigma_g$ with constant and non-constant level of dependence with $\alpha_{max} \in \{0.9, 0.6, 0.3\}$ to test high, medium, and low correlation within a group of features. For non-constant level of dependence we set $c \in \{0.35, 0.25, 0.1\}$. Visualisation of $\Sigma$, $\Sigma_g$ and the $\beta$ vector for this scenario are provided in Appendix E.

## Performance indicators

We evaluated feature selection methods primarily based on their predictive power using classification models (AUC, specificity and sensitivity at maximum Youden's index), the quantity and importance of the selected features (number of informative and non-informative features selected), and the stability of the selection (Jaccard index). Additionally, we measured the execution time for the selection process (runtime).

(1) **AUC**: Is a well-known indicator of a classification model's predictive power. It is computed from the ROC curve, which plots sensitivity against 1-specificity at various probability cutoffs. The AUC is simply the area under this ROC curve and can be interpreted as the probability that the model will rank a randomly chosen positive example higher than a negative one. A higher AUC indicates better model performance.

(2) **Sensitivity and specificity**: Using the ROC curve, we also determined the optimal sensitivity and specificity based on Youden's index. This index is calculated as the point on the ROC curve that maximizes the value of specificity + sensitivity - 1.

(3) **Feature selection**: In each scenario of the simulation study, we have prior knowledge of which features are truly informative and which are non-informative. For each simulation repetition, we assessed how many of the selected features were informative or non-informative. At the conclusion of the simulation, we calculated the average number of selected informative and non-informative features for each method.

(4) **Stability**: The stability of feature selection was evaluated using the average over all pairwise similarity comparisons between the different set of selected features [30]. For a feature selection method we compute its stability over all the repetition using the formula:

$$ \text{Stability}(S) = \frac{2}{R(R-1)} \sum_{i=1}^{R-1} \sum_{j=i+1}^{R} J(S_i, S_j), $$

where $S$ is the set of all selected features at each iteration by a given method, and $R$ is the number of repetitions. The Jaccard index [28], $J(S_i, S_j) = \frac{|S_i \cap S_j|}{|S_i \cup S_j|}$ measures the similarity between the sets of selected features at repetition $i$ and $j$. A Jaccard index of 1 means the two sets are identical, while an index of 0 indicates completely different sets.

(5) **Runtime**: It is the execution time of the feature selection process in seconds. The simulation study
was performed on a device with 16.0Go of RAM, AMD Ryzen 7 5700X 8-Core Processor 3.40 GHz.

## Feature selection methods ranking

Unless a single method outperforms all others across every indicators, it can be challenging to determine which method performs better overall, as each indicator measures a different aspect of performance. To address this, we applied a multiple criteria decision-making (MCDM) method to rank the methods based on their performance across all criteria (in our case, a criteria is a performance indicator). We applied the Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) method [31, 32]. TOPSIS ranks feature selection methods based on their geometrical distance to an ideal solution, where higher-ranked methods are closer to this ideal and farthest from the negative ideal (the worst solution). Before calculating these distances, the indicators are normalized, and weights are applied based on the relative importance of each indicator as defined by the user. A step-by-step explanation of the TOPSIS score is detailed in Appendix C. In our study, we prioritize the performance of the classification models, specifically their predictive power. This means that the indicators

AUC, sensitivity and specificity are given the highest weight. Following that, we assign importance to the number of selected features - preferring methods that select fewer features, provided they achieve similar predictive performance. Finally, the stability of the feature selection methods and their running time are considered of lower importance, and thus, they receive lower weights in the ranking process. Table 2 summarizes the ranking and corresponding weights for each indicator.

## Results

### Scenario 1

The results for feature selection methods are presented in Table 3 and TOPSIS scores are presented in Fig. 1. The GAMMA_BF and GAMMA_FORW methods performed similarly, ranking as the top two methods in this scenario. Both demonstrated strong predictive performance, with AUC values of 85.92%, sensitivity of 85.98% and specificity of 86.59%. These methods also selected nearly all of the informative features (an average of 2.92 out of 3), while keeping the number of non-informative features low (1.62 out of 22 on average). The distinction between these methods lies in their execution time that

<table><caption>**Table 2** Ranks, weights ($w_i$), and categories of the performance indicators used for evaluating feature selection methods</caption><thead><tr><th></th><th>AUC</th><th>Sensitivity</th><th>Specificity</th><th>NSIF</th><th>NSNIF</th><th>Stability</th><th>Runtime (s)</th></tr></thead><tbody><tr><th>Rank</th><td>1</td><td>2</td><td>3</td><td>4</td><td>5</td><td>6</td><td>7</td></tr><tr><th>Weight w<sub>i</sub></th><td>0.2500</td><td>0.2143</td><td>0.1786</td><td>0.1429</td><td>0.1071</td><td>0.0714</td><td>0.0357</td></tr><tr><th>Category</th><td>Benefit</td><td>Benefit</td><td>Benefit</td><td>Benefit</td><td>Cost</td><td>Benefit</td><td>Cost</td></tr></tbody></table>

The Category column specifies whether the indicator is a benefit (higher values preferred) or a cost (lower values preferred)

Abbreviations: NSIF Number of Selected Informative Features, NSNIF Number of Selected Non-Informative Features

<table><caption>**Table 3** Results for scenario 1, displaying all performance indicators averaged over the 50 repetitions</caption><thead><tr><th>FSM</th><th>AUC</th><th>Se</th><th>Spe</th><th>NSIF</th><th>NSNIF</th><th>Stability</th><th>Runtime (s)</th><th>Rank</th></tr></thead><tbody><tr><th>FULL</th><td>85.77</td><td>86.11</td><td>86.18</td><td>3.00</td><td>22.00</td><td>1.00</td><td>0.04</td><td>11</td></tr><tr><th>CFS</th><td>83.76</td><td>84.30</td><td>84.02</td><td>1.82</td><td><strong>0.00</strong></td><td>0.849</td><td>0.23</td><td>6</td></tr><tr><th>CHI2</th><td>85.09</td><td>86.09</td><td>84.88</td><td>1.98</td><td><strong>0.00</strong></td><td><strong>0.980</strong></td><td><strong>0.16</strong></td><td>4</td></tr><tr><th>GAMMA_BACK</th><td><strong>85.92</strong></td><td>85.98</td><td><strong>86.59</strong></td><td><strong>2.92</strong></td><td><strong>1.62</strong></td><td>0.805</td><td>5.68</td><td>3</td></tr><tr><th><strong>GAMMA_BF</strong></th><td><strong>85.92</strong></td><td>85.98</td><td><strong>86.59</strong></td><td><strong>2.92</strong></td><td><strong>1.62</strong></td><td>0.805</td><td>0.45</td><td><strong>2</strong></td></tr><tr><th><strong>GAMMA_FORW</strong></th><td><strong>85.92</strong></td><td>85.98</td><td><strong>86.59</strong></td><td><strong>2.92</strong></td><td><strong>1.62</strong></td><td>0.805</td><td>0.44</td><td><strong>1</strong></td></tr><tr><th>LASSO</th><td>84.02</td><td>84.48</td><td>84.21</td><td>2.90</td><td>2.26</td><td>0.504</td><td>2.38</td><td>5</td></tr><tr><th>RFI</th><td>85.25</td><td><strong>86.21</strong></td><td>85.05</td><td>2.00</td><td><strong>0.00</strong></td><td><strong>1.00</strong></td><td>9.22</td><td>7</td></tr><tr><th>STEP</th><td>85.87</td><td><strong>86.41</strong></td><td>86.05</td><td><strong>3.00</strong></td><td>3.16</td><td>0.404</td><td>1.75</td><td>8</td></tr><tr><th>SU</th><td>77.28</td><td>78.18</td><td>77.27</td><td>1.08</td><td><strong>0.00</strong></td><td>0.925</td><td><strong>0.14</strong></td><td>9</td></tr><tr><th>SVM-RFE</th><td>85.91</td><td>86.16</td><td><strong>86.42</strong></td><td><strong>3.00</strong></td><td>2.44</td><td>0.489</td><td>23.36</td><td>10</td></tr></tbody></table>

The best value of each indicator is highlighted with bold and underline text. The second-best value is highlighted in bold. The indicators for the FULL model are not highlighted in the table since the focus of the study is to compare the feature selection methods

Abbreviations: FSM Feature Selection Methods, Se Sensitivity, Spe Specificity, NSIF Number of Selected Informative Features, NSNIF Number of Selected Non-Informative Features
<figure><img src="image_4.png" /><figcaption>**Fig. 1** TOPSIS score of the feature selection methods on scenario 1</figcaption></figure>

was slightly longer for GAMMA_BF (0.45s against 0.44s for GAMMA_FORW).

Other methods also showed good performance. CFS, CHI2, RFI, and SU were able to select only informative features, but fewer than the γ-metric based methods, resulting in lower predictive power and rankings (6th, 4th, 7th, and 9th for CFS, CHI2, RFI, and SU, respectively). Notably, RFI exhibited excellent stability, with a Jaccard index of 1.00, indicating that it consistently selected the same two informative features. However, RFI was unable to select all the informative features, limiting its predictive performance.

STEP and SVM-RFE were able to select all informative features in each iteration but also included more non-informative features. SU performed the worst, with an average of only 1.08 informative features selected and significantly lower predictive power (AUC = 77.28%, sensitivity = 78.18%, and specificity = 77.27%).

Figure 1 illustrates how closely the methods actually performed in this scenario.

### Scenario 2

Results for all feature selection methods and indicators (for 200 features) are presented in Table 4, with TOPSIS scores shown in Fig. 2. In this scenario, LASSO consistently emerged as the best feature selection method across all cases (as reflected in its TOPSIS score in Fig. 2). In the strong separation case, GAMMA_BF and GAMMA_FORW also performed well, ranking 2nd and 3rd with unbalanced classes, and 4th and 5th with balanced classes, respectively. Both methods achieved the second-best AUC values (81.46% for unbalanced and 87.46% for balanced classes) and the highest sensitivity (88.70% for unbalanced classes).

For balanced classes, CFS and CHI2 were ranked 3rd and 2nd, respectively. These methods excelled at ignoring non-informative features (averaging 0.40 for CFS and 0.46 for CHI2), achieving the best and second-best AUC values (87.86% for CFS and 88.41% for CHI2) and sensitivity (90.52% for CFS and 91.12% for CHI2).

In the case of weak separation between classes, all methods demonstrated lower predictive performance, with the highest AUC at 72.15% (LASSO). For unbalanced classes, GAMMA_BF and GAMMA_FORW selected more non-informative features (17.52 on average out of 197) and fewer informative ones (1.82 on average out of 3). GAMMA_BACK managed to select the most informative features (2.60 on average) but also included a high number of non-informative features (93.62 on average). SVM-RFE exhibited a similar pattern, selecting many non-informative features (19.32 on average for unbalanced classes, 24.86 for balanced classes) and a moderate number of informative features (2.12 for unbalanced classes and 2.48 for balanced classes).

As shown in Fig. 2, LASSO remains the top-performing method, especially in cases with strong classes separation, with GAMMA_BF and GAMMA_FORW consistently close behind in strong separation scenarios.

### Scenario 3

Results for all feature selection methods and indicators are presented in Table 5, with TOPSIS scores shown in Fig. 3.

**Constant correlation:** First, let's consider the case where the correlation is constant (the left side of Table 5). In the case of a strong correlation $\alpha_{max} = 0.9$, LASSO
<table><caption>**Table 4** Results for scenario 2, displaying all performance indicators averaged over 50 repetitions</caption><thead><tr><th rowspan="2">FSM</th><th colspan="6">Unbalanced classes</th><th colspan="6">Balanced classes</th></tr><tr><th>AUC</th><th>Se</th><th>Spe</th><th>NSIF</th><th>NSNIF</th><th>Stability</th><th>Rank</th><th>FSM</th><th>AUC</th><th>Se</th><th>Spe</th><th>NSIF</th><th>NSNIF</th><th>Stability</th><th>Runtime</th><th>Rank</th></tr></thead><tbody><tr><td colspan="17">Weak separation</td></tr><tr><td>FULL</td><td>53.23</td><td>62.19</td><td>49.75</td><td>3.00</td><td>197.00</td><td>1.00</td><td>0.01</td><td>11</td><td>FULL</td><td>53.53</td><td>57.75</td><td>52.52</td><td>3.00</td><td>197.00</td><td>1.00</td><td>0.01</td><td>11</td></tr><tr><td>CFS</td><td>71.96</td><td><strong>88.23</strong></td><td>78.80</td><td>1.18</td><td>0.36</td><td>0.638</td><td>1.39</td><td>4</td><td><strong>CFS</strong></td><td>77.59</td><td>80.80</td><td>81.52</td><td>1.24</td><td>0.28</td><td>0.681</td><td>1.77</td><td><strong>2</strong></td></tr><tr><td><strong>CHI2</strong></td><td><strong>72.06</strong></td><td><strong>88.38</strong></td><td>78.94</td><td>1.22</td><td>0.38</td><td>0.616</td><td><strong>0.36</strong></td><td><strong>2</strong></td><td><strong>CHI2</strong></td><td><strong>77.61</strong></td><td><strong>81.31</strong></td><td>81.11</td><td>1.32</td><td>0.44</td><td>0.552</td><td><strong>0.48</strong></td><td>5</td></tr><tr><td>GAMMA_BACK</td><td>54.71</td><td>61.86</td><td>56.23</td><td><strong>2.60</strong></td><td>93.62</td><td>0.317</td><td>271.90</td><td>10</td><td>GAMMA_BACK</td><td>53.16</td><td>58.03</td><td>54.18</td><td><strong>2.82</strong></td><td>115.22</td><td>0.422</td><td>299.86</td><td>10</td></tr><tr><td>GAMMA_BF</td><td>68.27</td><td>78.59</td><td>74.40</td><td>1.82</td><td>17.52</td><td>0.119</td><td>60.52</td><td>8</td><td>GAMMA_BF</td><td>70.33</td><td>73.12</td><td>73.91</td><td><strong>2.50</strong></td><td>28.82</td><td>0.130</td><td>136.31</td><td>9</td></tr><tr><td>GAMMA_FORW</td><td>68.27</td><td>78.59</td><td>74.40</td><td>1.82</td><td>17.52</td><td>0.119</td><td>26.41</td><td>7</td><td>GAMMA_FORW</td><td>70.33</td><td>73.12</td><td>73.91</td><td><strong>2.50</strong></td><td>28.82</td><td>0.130</td><td>60.18</td><td>7</td></tr><tr><td><strong>LASSO</strong></td><td><strong>72.15</strong></td><td>87.08</td><td><strong>79.84</strong></td><td>1.54</td><td>0.92</td><td>0.619</td><td>1.85</td><td><strong>1</strong></td><td><strong>LASSO</strong></td><td><strong>78.34</strong></td><td><strong>81.17</strong></td><td><strong>81.81</strong></td><td>2.00</td><td>1.08</td><td>0.510</td><td>2.41</td><td><strong>1</strong></td></tr><tr><td>RFI</td><td>70.98</td><td>86.33</td><td>77.21</td><td>1.02</td><td><strong>0.14</strong></td><td><strong>0.865</strong></td><td><strong>0.45</strong></td><td>5</td><td>RFI</td><td>77.23</td><td>80.27</td><td>80.86</td><td>1.04</td><td><strong>0.00</strong></td><td><strong>0.961</strong></td><td>0.68</td><td>6</td></tr><tr><td>STEP</td><td>70.82</td><td>58.12</td><td><strong>86.06</strong></td><td>1.94</td><td>5.10</td><td>0.136</td><td>3.60</td><td>6</td><td>STEP</td><td>73.72</td><td>76.71</td><td>72.72</td><td>2.34</td><td>7.36</td><td>0.137</td><td>5.70</td><td>3</td></tr><tr><td>SU</td><td>71.90</td><td>88.01</td><td>78.95</td><td>1.18</td><td><strong>0.34</strong></td><td><strong>0.659</strong></td><td><strong>0.45</strong></td><td>3</td><td>SU</td><td>77.46</td><td>80.21</td><td><strong>81.61</strong></td><td>1.18</td><td><strong>0.26</strong></td><td><strong>0.733</strong></td><td><strong>0.61</strong></td><td>4</td></tr><tr><td>SVM-RFE</td><td>64.34</td><td>65.03</td><td>75.19</td><td><strong>2.12</strong></td><td>19.32</td><td>0.090</td><td>10.07</td><td>9</td><td>SVM-RFE</td><td>67.30</td><td>71.72</td><td>69.38</td><td>2.48</td><td>24.86</td><td>0.104</td><td>14.63</td><td>8</td></tr><tr><td colspan="17">Strong separation</td></tr><tr><td>FULL</td><td>54.36</td><td>60.03</td><td>53.04</td><td>3.00</td><td>197.00</td><td>1.00</td><td>0.01</td><td>11</td><td>FULL</td><td>54.44</td><td>61.21</td><td>53.30</td><td>3.00</td><td>197.00</td><td>1.00</td><td>0.01</td><td>11</td></tr><tr><td>CFS</td><td>78.82</td><td>86.43</td><td>84.10</td><td>1.54</td><td>0.36</td><td>0.586</td><td>1.77</td><td>6</td><td>CFS</td><td><strong>87.86</strong></td><td><strong>90.52</strong></td><td>90.37</td><td>2.04</td><td><strong>0.40</strong></td><td>0.711</td><td>2.49</td><td>3</td></tr><tr><td><strong>CHI2</strong></td><td><strong>79.81</strong></td><td>87.23</td><td>85.17</td><td>1.70</td><td>0.44</td><td>0.561</td><td><strong>0.38</strong></td><td>4</td><td><strong>CHI2</strong></td><td><strong>88.41</strong></td><td><strong>91.12</strong></td><td>90.74</td><td>2.10</td><td>0.46</td><td><strong>0.713</strong></td><td><strong>0.43</strong></td><td><strong>2</strong></td></tr><tr><td>GAMMA_BACK</td><td>56.49</td><td>61.30</td><td>59.48</td><td><strong>2.88</strong></td><td>92.36</td><td>0.316</td><td>274.23</td><td>10</td><td>GAMMA_BACK</td><td>55.45</td><td>61.06</td><td>56.06</td><td><strong>2.80</strong></td><td>107.92</td><td>0.382</td><td>268.72</td><td>10</td></tr><tr><td>GAMMA_BF</td><td><strong>81.46</strong></td><td><strong>88.70</strong></td><td>84.96</td><td>2.30</td><td>4.50</td><td>0.357</td><td>17.85</td><td>3</td><td>GAMMA_BF</td><td><strong>87.86</strong></td><td>89.33</td><td><strong>91.19</strong></td><td>2.24</td><td>1.04</td><td>0.525</td><td>7.36</td><td>5</td></tr><tr><td><strong>GAMMA_FORW</strong></td><td><strong>81.46</strong></td><td><strong>88.70</strong></td><td>84.96</td><td>2.30</td><td>4.50</td><td>0.357</td><td>10.09</td><td><strong>2</strong></td><td>GAMMA_FORW</td><td><strong>87.86</strong></td><td>89.33</td><td><strong>91.19</strong></td><td>2.24</td><td>1.04</td><td>0.525</td><td>6.05</td><td>4</td></tr><tr><td><strong>LASSO</strong></td><td><strong>82.84</strong></td><td><strong>87.94</strong></td><td><strong>87.41</strong></td><td>2.36</td><td>2.02</td><td>0.536</td><td>1.95</td><td><strong>1</strong></td><td><strong>LASSO</strong></td><td><strong>87.74</strong></td><td>89.04</td><td><strong>91.27</strong></td><td>2.26</td><td>1.30</td><td>0.698</td><td>2.17</td><td><strong>1</strong></td></tr><tr><td>RFI</td><td>76.13</td><td>85.21</td><td>81.39</td><td>1.22</td><td><strong>0.00</strong></td><td><strong>0.819</strong></td><td>0.51</td><td>8</td><td>RFI</td><td>84.94</td><td>88.19</td><td>86.98</td><td>1.78</td><td><strong>0.00</strong></td><td><strong>0.780</strong></td><td>0.61</td><td>7</td></tr><tr><td>STEP</td><td>81.07</td><td>76.48</td><td><strong>88.12</strong></td><td>2.48</td><td>3.58</td><td>0.243</td><td>3.53</td><td>5</td><td>STEP</td><td>83.76</td><td>84.88</td><td>84.33</td><td><strong>2.56</strong></td><td>3.90</td><td>0.235</td><td>4.12</td><td>8</td></tr><tr><td>SU</td><td>77.44</td><td>85.75</td><td>82.50</td><td>1.36</td><td><strong>0.22</strong></td><td><strong>0.668</strong></td><td><strong>0.48</strong></td><td>7</td><td>SU</td><td>86.58</td><td>89.97</td><td>88.29</td><td>1.92</td><td>0.42</td><td>0.65</td><td><strong>0.53</strong></td><td>6</td></tr><tr><td>SVM-RFE</td><td>72.16</td><td>76.14</td><td>78.41</td><td><strong>2.68</strong></td><td>16.86</td><td>0.116</td><td>10.48</td><td>9</td><td>SVM-RFE</td><td>75.52</td><td>76.71</td><td>79.98</td><td>2.52</td><td>18.26</td><td>0.110</td><td>11.54</td><td>9</td></tr></tbody></table>

The best value of each indicator is highlighted with bold and underline. The second-best value is highlighted in bold. The indicators for the FULL model are not highlighted in the table since the focus of the study is to compare the feature selection methods

Abbreviations: FSM Feature Selection Methods, Se Sensitivity, Spe Specificity, NSIF Number of Selected Informative Features, NSNIF Number of Selected Non-Informative Features
<figure><img src="image_5.png" /><figcaption>**Fig. 2** TOPSIS scores of the feature selection methods on scenario 2</figcaption></figure>

was ranked 1st, followed by CFS in 2nd. Both methods achieved an AUC above 85%, with CFS and LASSO reaching 85.68% and 85.62%, respectively. They also had high sensitivity values: 86.56% for CFS and 86.04% for LASSO, as well as among the best specificity (85.54% for CFS and 85.79% for LASSO). In this scenario, most feature selection methods were able to select all the informative features, but none could fully exclude the non-informative ones. GAMMA_BF and GAMMA_FORW were ranked 9th and 8th, respectively. Both selected all the informative features, but they also selected a significant number of non-informative features (an average of 57.56 out of 94 features). RFI was less efficient at selecting informative features under strong correlation (2.70 on average out of 6) and exhibited the worst predictive performance (AUC = 69.82%, sensitivity = 69.90%, and specificity = 70.74%).

When the correlation level $\alpha_{max}$ was lower ($\alpha_{max} = 0.6$), RFI performed much better compared to the strong correlation case. It was able to select all the informative features (6.00 on average) and excluded all the non-informative ones. RFI also achieved the best performance in terms of AUC (86.03%), specificity (86.78%), and the second-best sensitivity (85.94%). Additionally, RFI demonstrated perfect stability (1.00), consistently selecting only informative features. LASSO and CFS ranked 2nd and 3rd, respectively. LASSO selected only a few non-informative features (2.76 on average), although it also selected slightly fewer informative features. CFS had the second-best AUC (85.94%) and specificity (86.68%). GAMMA_BF and GAMMA_FORW performed well in terms of predictive metrics (AUC = 85.31%, sensitivity = 85.40%, and specificity = 85.98%), but they selected some non-informative features (36.48 on average for both).

At a correlation level $\alpha_{max}$ of 0.3, CFS ranked 1st with an AUC of 86.23% (tied with RFI) and a sensitivity of 87.02%. CFS, SU, and RFI were the only methods that selected only informative features, ranking 1st, 2nd, and 3rd, respectively. In this case, GAMMA_BF and GAMMA_FORW improved rankings and achieved the second-best specificity (86.16%).

**Non-constant correlation:** Next, we look at the case of non-constant correlation (the right side of Table 5). For a strong correlation level ($\alpha_{max} = 0.9$), CFS and LASSO once again ranked 1st and 2nd, respectively, with CFS achieving the best AUC (85.81%) and sensitivity (86.09%). In this setting, the feature selection methods generally selected fewer non-informative features. For instance, GAMMA_BF and GAMMA_FORW selected an average of 35.48 non-informative features compared to 57.56 in the constant correlation case. Both GAMMA_BF and GAMMA_FORW achieved over 85% in AUC, sensitivity, and specificity. RFI selected very few non-informative features (0.32 on average) but was only able to select 3.00 informative features, resulting in the lowest predictive performance (AUC = 71.63%, sensitivity = 74.33%, and specificity = 69.95%).

For a lower correlation level ($\alpha_{max} = 0.6$), CFS clearly ranked 1st, demonstrating the highest predictive performance, with all indicators exceeding 86%. RFI ranked 2nd with strong AUC (85.91%), sensitivity (85.99%), and specificity (86.53%). Moreover, RFI did not select any non-informative features and maintained a high stability score (0.993). GAMMA_BF and GAMMA_FORW ranked 6th and 5th, with a high specificity (86.37% for both) and an AUC of 85.61%. However, both methods selected a similar number of non-informative features (16.92 on average).
<table><caption>**Table 5** Results for scenario 3, displaying all indicators average over the 50 repetitions</caption><thead><tr><th rowspan="2">FSM</th><th colspan="10">Constant correlation</th><th colspan="10">Non constant correlation</th></tr><tr><th>AUC</th><th>Se</th><th>Spe</th><th>NSIF</th><th>NSNIF</th><th>Stability</th><th>Runtime</th><th>Rank</th><th>FSM</th><th>AUC</th><th>Se</th><th>Spe</th><th>NSIF</th><th>NSNIF</th><th>Stability</th><th>Runtime</th><th>Rank</th></tr></thead><tbody><tr><td colspan="18" style="text-align: left; font-weight: bold;">α<sub>max</sub> = 0.9</td></tr><tr><td>FULL</td><td>84.66</td><td>84.80</td><td>85.18</td><td>6.00</td><td>94.00</td><td>1.00</td><td>0.03</td><td>11</td><td>FULL</td><td>85.64</td><td>84.46</td><td>84.57</td><td>6.00</td><td>94.00</td><td>1.00</td><td>0.00</td><td>11</td></tr><tr><td>CFS</td><td><strong>85.68</strong></td><td><strong>86.56</strong></td><td>85.54</td><td><strong>5.94</strong></td><td>19.82</td><td>0.393</td><td>9.46</td><td><strong>2</strong></td><td>CFS</td><td><strong>85.81</strong></td><td><strong>86.09</strong></td><td>86.29</td><td><strong>6.00</strong></td><td>8.08</td><td><strong>0.625</strong></td><td>4.20</td><td><strong>1</strong></td></tr><tr><td>CHI2</td><td>85.40</td><td>85.95</td><td>85.47</td><td><strong>6.00</strong></td><td>45.00</td><td><strong>1.00</strong></td><td><strong>0.58</strong></td><td>5</td><td>CHI2</td><td>85.40</td><td>85.36</td><td><strong>86.34</strong></td><td><strong>6.00</strong></td><td>41.42</td><td><strong>0.897</strong></td><td><strong>0.52</strong></td><td>7</td></tr><tr><td>GAMMA_BACK</td><td>84.90</td><td>85.26</td><td>85.16</td><td><strong>6.00</strong></td><td>64.24</td><td>0.712</td><td>68.64</td><td>10</td><td>GAMMA_BACK</td><td>84.93</td><td>85.41</td><td>85.21</td><td><strong>6.00</strong></td><td>57.16</td><td>0.570</td><td>73.57</td><td>10</td></tr><tr><td>GAMMA_BF</td><td>84.95</td><td>85.45</td><td>85.23</td><td><strong>6.00</strong></td><td>57.56</td><td><strong>0.725</strong></td><td>42.92</td><td>9</td><td>GAMMA_BF</td><td>85.29</td><td>85.69</td><td>85.64</td><td><strong>6.00</strong></td><td>35.48</td><td>0.408</td><td>25.94</td><td>9</td></tr><tr><td>GAMMA_FORW</td><td>84.95</td><td>85.45</td><td>85.23</td><td><strong>6.00</strong></td><td>57.56</td><td><strong>0.725</strong></td><td>31.95</td><td>8</td><td>GAMMA_FORW</td><td>85.29</td><td>85.69</td><td>85.64</td><td><strong>6.00</strong></td><td>35.48</td><td>0.408</td><td>18.89</td><td>8</td></tr><tr><td><strong>LASSO</strong></td><td><strong>85.62</strong></td><td>86.04</td><td><strong>85.79</strong></td><td>5.92</td><td><strong>7.24</strong></td><td>0.373</td><td>2.99</td><td><strong>1</strong></td><td>LASSO</td><td>85.48</td><td><strong>85.90</strong></td><td>85.89</td><td><strong>5.94</strong></td><td><strong>3.48</strong></td><td>0.562</td><td>2.91</td><td><strong>2</strong></td></tr><tr><td>RFI</td><td>69.82</td><td>69.90</td><td>70.74</td><td>2.70</td><td>15.30</td><td>0.551</td><td>18.36</td><td>7</td><td>RFI</td><td>71.63</td><td>74.33</td><td>69.95</td><td>3.00</td><td><strong>0.32</strong></td><td>0.489</td><td>17.72</td><td>5</td></tr><tr><td>STEP</td><td>85.46</td><td><strong>86.40</strong></td><td>85.16</td><td><strong>6.00</strong></td><td>11.08</td><td>0.283</td><td>18.71</td><td>3</td><td>STEP</td><td>85.33</td><td>84.94</td><td><strong>86.57</strong></td><td><strong>6.00</strong></td><td>9.54</td><td>0.299</td><td>15.39</td><td>3</td></tr><tr><td>SU</td><td>85.40</td><td>85.95</td><td>85.47</td><td><strong>6.00</strong></td><td>45.00</td><td><strong>1.00</strong></td><td><strong>0.68</strong></td><td>6</td><td>SU</td><td>80.92</td><td>80.24</td><td>82.47</td><td>5.02</td><td>33.10</td><td>0.583</td><td><strong>0.62</strong></td><td>6</td></tr><tr><td>SVM-RFE</td><td>85.50</td><td>85.75</td><td><strong>86.07</strong></td><td><strong>6.00</strong></td><td><strong>10.02</strong></td><td>0.310</td><td>55.89</td><td>4</td><td>SVM-RFE</td><td><strong>85.55</strong></td><td>85.73</td><td>86.24</td><td><strong>6.00</strong></td><td>8.28</td><td>0.348</td><td>55.58</td><td>4</td></tr><tr><td colspan="18" style="text-align: left; font-weight: bold;">α<sub>max</sub> = 0.6</td></tr><tr><td>FULL</td><td>84.60</td><td>85.56</td><td>84.47</td><td>6.00</td><td>94.00</td><td>1.00</td><td>0.00</td><td>11</td><td>FULL</td><td>84.72</td><td>84.72</td><td>85.50</td><td>6.00</td><td>94.00</td><td>1.00</td><td>0.00</td><td>11</td></tr><tr><td>CFS</td><td><strong>85.94</strong></td><td>85.89</td><td><strong>86.68</strong></td><td><strong>6.00</strong></td><td>5.96</td><td>0.559</td><td>3.94</td><td>3</td><td>CFS</td><td><strong>86.01</strong></td><td><strong>86.06</strong></td><td><strong>86.61</strong></td><td><strong>6.00</strong></td><td><strong>0.24</strong></td><td><strong>0.935</strong></td><td>1.88</td><td><strong>1</strong></td></tr><tr><td>CHI2</td><td>85.45</td><td><strong>86.20</strong></td><td>85.39</td><td><strong>6.00</strong></td><td>44.98</td><td><strong>0.999</strong></td><td><strong>0.52</strong></td><td>7</td><td>CHI2</td><td>85.72</td><td>85.84</td><td>86.19</td><td><strong>6.00</strong></td><td>27.76</td><td>0.655</td><td><strong>0.51</strong></td><td>9</td></tr><tr><td>GAMMA_BACK</td><td>84.83</td><td>85.77</td><td>84.68</td><td><strong>6.00</strong></td><td>62.52</td><td>0.670</td><td>68.40</td><td>10</td><td>GAMMA_BACK</td><td>85.24</td><td>85.16</td><td>85.98</td><td><strong>6.00</strong></td><td>40.44</td><td>0.383</td><td>94.79</td><td>10</td></tr><tr><td>GAMMA_BF</td><td>85.31</td><td>85.40</td><td>85.98</td><td><strong>6.00</strong></td><td>36.48</td><td>0.412</td><td>27.75</td><td>9</td><td>GAMMA_BF</td><td>85.61</td><td>85.60</td><td>86.37</td><td><strong>6.00</strong></td><td>16.92</td><td>0.473</td><td>13.93</td><td>6</td></tr><tr><td>GAMMA_FORW</td><td>85.31</td><td>85.40</td><td>85.98</td><td><strong>6.00</strong></td><td>36.48</td><td>0.412</td><td>20.63</td><td>8</td><td>GAMMA_FORW</td><td>85.61</td><td>85.60</td><td>86.37</td><td><strong>6.00</strong></td><td>16.92</td><td>0.473</td><td>10.29</td><td>5</td></tr><tr><td><strong>LASSO</strong></td><td>85.22</td><td>85.29</td><td>85.84</td><td><strong>5.88</strong></td><td><strong>2.76</strong></td><td>0.607</td><td>2.42</td><td><strong>2</strong></td><td>LASSO</td><td>85.45</td><td>85.51</td><td>86.08</td><td>5.92</td><td>2.06</td><td>0.669</td><td>2.42</td><td>4</td></tr><tr><td>RFI</td><td><strong>86.03</strong></td><td><strong>85.94</strong></td><td><strong>86.78</strong></td><td><strong>6.00</strong></td><td><strong>0.00</strong></td><td><strong>1.00</strong></td><td>17.55</td><td><strong>1</strong></td><td>RFI</td><td><strong>85.91</strong></td><td><strong>85.99</strong></td><td><strong>86.53</strong></td><td><strong>5.98</strong></td><td><strong>0.00</strong></td><td><strong>0.993</strong></td><td>17.03</td><td><strong>2</strong></td></tr><tr><td>STEP</td><td>85.41</td><td>85.42</td><td>86.08</td><td><strong>6.00</strong></td><td>12.56</td><td>0.267</td><td>18.60</td><td>4</td><td>STEP</td><td>85.28</td><td>85.60</td><td>85.75</td><td><strong>6.00</strong></td><td>13.26</td><td>0.261</td><td>18.25</td><td>7</td></tr><tr><td>SU</td><td>84.28</td><td>83.98</td><td>85.25</td><td>5.72</td><td>25.18</td><td>0.511</td><td><strong>0.62</strong></td><td>5</td><td>SU</td><td>84.20</td><td>83.99</td><td>85.16</td><td>5.64</td><td>0.56</td><td>0.852</td><td><strong>0.61</strong></td><td>3</td></tr><tr><td>SVM-RFE</td><td>85.59</td><td>85.91</td><td>85.96</td><td><strong>6.00</strong></td><td>8.68</td><td>0.325</td><td>59.76</td><td>6</td><td>SVM-RFE</td><td>85.52</td><td>85.54</td><td>86.28</td><td><strong>6.00</strong></td><td>9.42</td><td>0.315</td><td>58.75</td><td>8</td></tr><tr><td colspan="18" style="text-align: left; font-weight: bold;">α<sub>max</sub> = 0.3</td></tr><tr><td>FULL</td><td>84.73</td><td>84.38</td><td>85.84</td><td>6.00</td><td>94.00</td><td>1.00</td><td>0.00</td><td>11</td><td>FULL</td><td>84.64</td><td>84.95</td><td>85.07</td><td>6.00</td><td>94.00</td><td>1.000</td><td>0.00</td><td>11</td></tr><tr><td>CFS</td><td><strong>86.23</strong></td><td><strong>87.02</strong></td><td>86.07</td><td><strong>6.00</strong></td><td><strong>0.00</strong></td><td><strong>1.000</strong></td><td>1.81</td><td><strong>1</strong></td><td>CFS</td><td><strong>86.05</strong></td><td><strong>86.18</strong></td><td><strong>86.67</strong></td><td><strong>6.00</strong></td><td><strong>0.00</strong></td><td><strong>1.000</strong></td><td>1.80</td><td><strong>3</strong></td></tr><tr><td>CHI2</td><td><strong>86.19</strong></td><td><strong>86.87</strong></td><td><strong>86.17</strong></td><td><strong>6.00</strong></td><td>1.92</td><td>0.841</td><td><strong>0.50</strong></td><td>4</td><td>CHI2</td><td><strong>86.04</strong></td><td><strong>86.23</strong></td><td><strong>86.61</strong></td><td><strong>6.00</strong></td><td>0.02</td><td><strong>0.994</strong></td><td><strong>0.51</strong></td><td><strong>2</strong></td></tr><tr><td>GAMMA_BACK</td><td>85.44</td><td>85.69</td><td>85.93</td><td><strong>6.00</strong></td><td>28.70</td><td>0.310</td><td>109.10</td><td>10</td><td>GAMMA_BACK</td><td>85.45</td><td>85.68</td><td>85.91</td><td><strong>6.00</strong></td><td>23.14</td><td>0.325</td><td>118.05</td><td>10</td></tr><tr><td>GAMMA_BF</td><td>85.71</td><td>86.01</td><td><strong>86.16</strong></td><td><strong>6.00</strong></td><td>17.10</td><td>0.440</td><td>13.97</td><td>7</td><td>GAMMA_BF</td><td>85.56</td><td>85.69</td><td>86.15</td><td><strong>6.00</strong></td><td>15.88</td><td>0.371</td><td>12.49</td><td>7</td></tr><tr><td>GAMMA_FORW</td><td>85.71</td><td>86.01</td><td><strong>86.16</strong></td><td><strong>6.00</strong></td><td>17.10</td><td>0.440</td><td>10.33</td><td>6</td><td>GAMMA_FORW</td><td>85.56</td><td>85.69</td><td>86.15</td><td><strong>6.00</strong></td><td>15.88</td><td>0.371</td><td>9.10</td><td>6</td></tr><tr><td>LASSO</td><td>85.85</td><td>86.63</td><td>85.71</td><td><strong>5.98</strong></td><td>4.56</td><td>0.587</td><td>2.38</td><td>5</td><td>LASSO</td><td>85.36</td><td>85.84</td><td>85.52</td><td><strong>5.90</strong></td><td>2.40</td><td>0.646</td><td>2.40</td><td>5</td></tr><tr><td>RFI</td><td><strong>86.23</strong></td><td><strong>87.02</strong></td><td>86.07</td><td><strong>6.00</strong></td><td><strong>0.00</strong></td><td><strong>1.000</strong></td><td>17.43</td><td>3</td><td>RFI</td><td><strong>86.05</strong></td><td><strong>86.18</strong></td><td><strong>86.67</strong></td><td><strong>6.00</strong></td><td><strong>0.00</strong></td><td><strong>1.000</strong></td><td>17.50</td><td>4</td></tr><tr><td>STEP</td><td>85.28</td><td>85.96</td><td>85.33</td><td><strong>6.00</strong></td><td>15.56</td><td>0.251</td><td>21.37</td><td>8</td><td>STEP</td><td>85.27</td><td>85.51</td><td>85.85</td><td><strong>6.00</strong></td><td>14.14</td><td>0.255</td><td>19.87</td><td>8</td></tr></tbody></table>
<table><caption>**Table 5 (continued)**</caption><thead><tr><th rowspan="2">FSM</th><th colspan="5">Constant correlation</th><th colspan="5">Non constant correlation</th></tr><tr><th>AUC</th><th>Se</th><th>Spe</th><th>NSIF</th><th>NSNIF</th><th>Stability</th><th>Runtime</th><th>Rank</th><th>FSM</th><th>AUC</th><th>Se</th><th>Spe</th><th>NSIF</th><th>NSNIF</th><th>Stability</th><th>Runtime</th><th>Rank</th></tr></thead><tbody><tr><td>SU</td><td>85.73</td><td>86.39</td><td>85.77</td><td>5.90</td><td><u>0.00</u></td><td><u>0.967</u></td><td><u>0.62</u></td><td>2</td><td><u>SU</u></td><td><u>86.05</u></td><td><u>86.18</u></td><td><u>86.67</u></td><td><u>6.00</u></td><td><u>0.00</u></td><td><u>1.000</u></td><td><u>0.59</u></td><td><u>1</u></td></tr><tr><td>SVM-RFE</td><td>85.54</td><td>85.90</td><td>85.93</td><td><u>6.00</u></td><td>10.64</td><td>0.291</td><td>60.70</td><td>9</td><td>SVM-RFE</td><td>85.47</td><td>86.01</td><td>85.74</td><td><u>6.00</u></td><td>9.52</td><td>0.305</td><td>60.91</td><td>9</td></tr></tbody></table>

The best value of each indicator is highlighted with bold and underline. The second-best value is highlighted in bold. The indicators for the FULL model are not highlighted in the table since the focus of the study is to compare the feature selection methods

Abbreviations: FSM Feature Selection Methods, Se Sensitivity, Spe Specificity, NSIF Number of Selected Informative Features, NSNIF Number of Selected Non-Informative Features
<figure><img src="image_6.png" /><figcaption>**Fig. 3** TOPSIS scores of the feature selection methods on scenario 3</figcaption></figure>

At the lowest correlation level ($\alpha_{max} = 0.3$), SU, RFI, and CFS achieved similar results, except for differences in running time. They ranked 1st, 4th, and 3rd, respectively, with an AUC of 86.05%, sensitivity of 86.18%, and specificity of 86.67%. These three methods successfully selected all and only the informative features. CHI2 had comparable results and ranked 2nd, mainly thanks to its short execution time (0.51s on average). GAMMA_BF and GAMMA_FORW ranked 7th and 6th, respectively, with an average of 15.88 non-informative features selected, while still selecting all informative features.

## Application

To evaluate the performance of the three proposed novel $\gamma$-metric feature selection methods, we compared each method's ability to select the most relevant features for distinguishing between AF and NSR using Holter-based ECG data. This evaluation was based on the classification performance of models built with the features selected by each method.

The Holter ECG recordings used in this study were from patients who exhibited only one type of cardiac activity - either AF or NSR. These data, which were previously collected as part of routine medical care, were obtained from the Department of Cardiology and Rhythmology at Marseille University Hospital Center (Timone Hospital) in France. The sample consisted of 34 files, each corresponding to a 24-hour Holter recording for and individual patient. Recordings were conducted between November 2016 and February 2017. The sample included 18 men and 16 women, with a median age of 62 years. Among these patients, 11 (32.4%) had episodes of AF.

## Data

The recordings consisted of RR-intervals, or beat-to-beat time intervals. Each patient's recording was segmented into a total 41 661 segments, each lasting 60 seconds. For each segment, we calculated various heart rhythm variability indicators, as outlined by [33]. These included: (i) the standard deviation of all RR intervals (SDNN); (ii) the standard deviation of the averages of 5-second RR intervals (SDANN); (iii) the mean of the standard deviations of 5-second RR intervals (SDNNidx); (iv) the percentage of differences between successive RR intervals greater than 50 ms (pNN50); (v) the standard deviation of successive differences (SDSD); (vi) the root-mean-square of successive differences (RMSSD); (vii) the interquartile range of the differences between successive RR intervals (IRRR); (viii) the median of the absolute differences between adjacent RR intervals (MADRR); (ix) the triangular interpolation of the RR interval histogram (TINN); (x) the integral of the RR interval histogram density divided by its height (HRV.index); (xi) the means (denoted mn.0df to mn.2df); and (xii) the standard deviations (denoted sd.0df to sd.2df) of the RR-intervals derivatives (up to the 2nd order). Altogether, the dataset included 16 features, which are described in Table 6 for the training sample.

## Feature selection and classification

In this healthcare-focused application, the dataset was divided into a training sample and a validation sample, using 60% of the dataset for the training sample. Sampling was done so segments from the same patient was exclusively in the training sample or exclusively in the validation sample. Feature selection was performed on the training sample for each feature selection method, and the selected features were then used to build a
<table><caption>**Table 6** Description of the training sample for AF and NSR. Mean and standard (in brackets) deviations are displayed for each feature</caption><thead><tr><th></th><th>AF (N = 6 087)</th><th>NSR (N = 16 712)</th><th>p-value †</th></tr></thead><tbody><tr><td>mn.0df (mean (SD))</td><td>0.73 (0.17)</td><td>0.89 (0.14)</td><td>&lt; 0.001</td></tr><tr><td>sd.0df (mean (SD))</td><td>0.14 (0.04)</td><td>0.05 (0.03)</td><td>&lt; 0.001</td></tr><tr><td>mn.1df (mean (SD))</td><td>-0.04 (0.02)</td><td>0.00 (0.01)</td><td>&lt; 0.001</td></tr><tr><td>sd.1df (mean (SD))</td><td>0.29 (0.08)</td><td>0.04 (0.06)</td><td>&lt; 0.001</td></tr><tr><td>mn.2df (mean (SD))</td><td>0.13 (0.08)</td><td>0.01 (0.03)</td><td>&lt; 0.001</td></tr><tr><td>sd.2df (mean (SD))</td><td>0.79 (0.35)</td><td>0.10 (0.24)</td><td>&lt; 0.001</td></tr><tr><td>SDNN (mean (SD))</td><td>144.06 (35.58)</td><td>49.37 (34.62)</td><td>&lt; 0.001</td></tr><tr><td>SDANN (mean (SD))</td><td>57.84 (25.01)</td><td>38.16 (28.03)</td><td>&lt; 0.001</td></tr><tr><td>SDNNIDX (mean (SD))</td><td>136.12 (32.62)</td><td>27.01 (20.52)</td><td>&lt; 0.001</td></tr><tr><td>pNN50 (mean (SD))</td><td>76.79 (7.23)</td><td>11.68 (16.09)</td><td>&lt; 0.001</td></tr><tr><td>SDSD (mean (SD))</td><td>201.73 (49.18)</td><td>36.81 (38.55)</td><td>&lt; 0.001</td></tr><tr><td>RMSSD (mean (SD))</td><td>200.48 (48.76)</td><td>36.54 (38.27)</td><td>&lt; 0.001</td></tr><tr><td>IRRR (mean (SD))</td><td>186.24 (48.08)</td><td>60.91 (42.72)</td><td>&lt; 0.001</td></tr><tr><td>MADRR (mean (SD))</td><td>129.02 (33.84)</td><td>19.58 (15.68)</td><td>&lt; 0.001</td></tr><tr><td>TINN (mean (SD))</td><td>256.39 (54.54)</td><td>121.61 (46.07)</td><td>&lt; 0.001</td></tr><tr><td>HRV.index (mean (SD))</td><td>16.41 (3.49)</td><td>7.78 (2.95)</td><td>&lt; 0.001</td></tr></tbody></table>

logistic regression classifier. The performance of each logistic regression classifier was subsequently evaluated on both the training and validation samples.

Since the truly informative features were not known beforehand, it was not feasible to define if the selected features were truly informative or not. Hence the evaluation of the performances rely solely on the AUC, the sensitivity and specificity (at maximum Youden index), the number of selected features in total, and the running time of the feature selection process. Feature selection was run only once, so no stability index could be computed. For all the above reasons, TOPSIS score was not computed on the healthcare-case application. For this study, AF was treated as the positive case. Sensitivity, in this context, represents the percentage of correctly predicted AF cases, while specificity represents the percentage of correctly predicted NSR cases.

## Results

The results for discriminating AF from NSR are presented in Tables 7 and 8 for the training and validation samples, respectively. The results from the FULL method highlight the importance of including a feature selection step. On the training sample, the model built with all features showed near-perfect sensitivity and specificity (99.61% and 99.83%, respectively), but its performance dropped significantly when applied to new data. On the validation sample, the FULL model's sensitivity decreased to 95.08%, while its specificity remained high at 99.7%.

<table><caption>**Table 7** Training sample results for AF detection</caption><thead><tr><th>FSM</th><th>AUC</th><th>Sensitivity</th><th>Specificity</th><th>NSF</th><th>Runtime (s)</th></tr></thead><tbody><tr><td>FULL</td><td>99.72</td><td>99.61</td><td>99.83</td><td>16</td><td>0.00</td></tr><tr><td>CFS</td><td>99.95</td><td>99.46</td><td>99.41</td><td>3</td><td>0.11</td></tr><tr><td>CHI2</td><td>99.87</td><td>99.77</td><td>99.96</td><td>15</td><td>0.83</td></tr><tr><td>GAMMA_BACK</td><td>99.99</td><td>99.95</td><td>99.90</td><td>7</td><td>2.12</td></tr><tr><td>GAMMA_BF</td><td>99.99</td><td>99.95</td><td>99.90</td><td>7</td><td>0.91</td></tr><tr><td>GAMMA_FORW</td><td>99.99</td><td>99.95</td><td>99.90</td><td>7</td><td>0.80</td></tr><tr><td>LASSO</td><td>99.73</td><td>99.51</td><td>99.95</td><td>7</td><td>20.84</td></tr><tr><td>RFI</td><td>78.42</td><td>71.35</td><td>86.93</td><td>1</td><td>178.29</td></tr><tr><td>STEP</td><td>99.99</td><td>99.98</td><td>99.95</td><td>6</td><td>52.84</td></tr><tr><td>SU</td><td>99.87</td><td>99.77</td><td>99.96</td><td>15</td><td>0.87</td></tr><tr><td>SVM-RFE</td><td>99.81</td><td>99.66</td><td>99.96</td><td>14</td><td>43.97</td></tr></tbody></table>

NSF Number of Selected Features

The CHI2, SU, and SVM-RFE methods selected nearly all features and achieved very high specificity (99.67%, 99.67%, and 99.75%, respectively), but their sensitivity was considerably lower (90.90%, 90.90%, and 87.17%). In contrast, the RFI method, which selected only one feature, performed the worst on the validation sample, with an AUC of 66.22%, sensitivity of 56.52%, and specificity of 79.36%. The γ-metric based methods - GAMMA_BACK, GAMMA_BF, and GAMMA_FORW - had similar results. They achieved the second-highest AUC (98.74%) and sensitivity (95.08%) and maintained a high specificity of 97.4%. LASSO selected seven features but showed lower sensitivity (82.87%) and AUC (91.27%). Finally, the CFS method delivered the best overall results, selecting three features and achieving an AUC of 99.51%, sensitivity of 98.61%, and specificity of 96.33%.

<table><caption>**Table 8** Validation sample results for AF detection</caption><thead><tr><th>FSM</th><th>AUC</th><th>Sensitivity</th><th>Specificity</th><th>NSF</th><th>Runtime (s)</th></tr></thead><tbody><tr><td>FULL</td><td>95.48</td><td>91.26</td><td>99.7</td><td>16</td><td>0.00</td></tr><tr><td>CFS</td><td><strong>99.51</strong></td><td><strong>98.61</strong></td><td>96.33</td><td><strong>3</strong></td><td><strong>0.11</strong></td></tr><tr><td>CHI2</td><td>95.28</td><td>90.90</td><td><strong>99.67</strong></td><td>15</td><td>0.83</td></tr><tr><td>GAMMA_BACK</td><td><strong>98.74</strong></td><td><strong>95.08</strong></td><td>97.40</td><td>7</td><td>2.12</td></tr><tr><td>GAMMA_BF</td><td><strong>98.74</strong></td><td><strong>95.08</strong></td><td>97.40</td><td>7</td><td>0.91</td></tr><tr><td>GAMMA_FORW</td><td><strong>98.74</strong></td><td><strong>95.08</strong></td><td>97.40</td><td>7</td><td><strong>0.80</strong></td></tr><tr><td>LASSO</td><td>91.27</td><td>82.87</td><td>99.66</td><td>7</td><td>20.84</td></tr><tr><td>RFI</td><td>66.22</td><td>56.52</td><td>79.36</td><td><strong>1</strong></td><td>178.29</td></tr><tr><td>STEP</td><td>98.27</td><td>94.75</td><td>98.07</td><td>6</td><td>52.84</td></tr><tr><td>SU</td><td>95.28</td><td>90.90</td><td>99.67</td><td>15</td><td>0.87</td></tr><tr><td>SVM-RFE</td><td>93.46</td><td>87.17</td><td><strong>99.75</strong></td><td>14</td><td>43.97</td></tr></tbody></table>

The best value of each indicator is highlighted with bold and underline. The second-best value is highlighted in bold. The indicators for the FULL model are not highlighted in the table since the focus of the study is to compare the feature selection methods

Abbreviation: NSF Number of Selected Features

<footnote>† The *p*-value is based on the Mann-Whitney U test</footnote>
## Discussions

In this article, we proposed a novel multivariate filter methodology for feature selection, based on three distinct methods that utilize the γ-metric as an evaluation function for classification. These methods - GAMMA_BACK, GAMMA_BF, and GAMMA_FORW - differ in their search direction. Unlike the univariate methods presented by [9], the methods explored here are multi-variate. We also incorporated shrinkage estimation of the covariance matrix into the γ-metric calculation and compared these methods to seven traditional feature selection methods: CFS, CHI2, LASSO, RFI, STEP, SU, and SVM-RFE. Both simulation studies and real-world data on AF detection were used for comparison. The three γ-metric based methods effectively identified features with non-null effects, although they were less efficient at excluding non-informative features. GAMMA_BACK tended to select more features compared to GAMMA_BF and GAMMA_FORW. Additionally, the γ-metric could only be computed for numerical features, limiting the feature selection to numerical data.

The capacity of the γ-metric based methods to detect informative features was illustrated in all scenario and also in the application, Table 9 summarizes the best results and key conclusions from each scenario. When the informative features had varying effect size (Scenario 1), the γ-metric based methods consistently selected the three informative features, including $x_3$, which had the smallest β coefficient. Only STEP and SVM-RFE selected more informative features in this scenario. GAMMA_BACK, GAMMA_BF, and GAMMA_FORW outperformed the other methods and were ranked 1st, 2nd and 3rd according to the TOPSIS score. We tested the capacity of the methods to perform feature selection in the cases where the number of features was greater than the number of observations (Scenario 2), with both balanced and unbalanced classes, as well as weak and strong class separation. In strong separation cases, GAMMA_BF and GAMMA_FORW successfully selected most of

the informative features and only a few non-informative ones, even achieving the 3rd and 2nd rank respectively for unbalanced data and 5th and 4th rank with balanced classes. We tested if the methods were able to disregard non-informative features that were correlated to informative features (Scenario 3) with different levels of correlation. The γ-metric methods consistently selected all informative features. GAMMA_BF and GAMMA_FORW were more successful at excluding non-informative features when the correlation level was low. In the healthcare application, GAMMA_BACK, GAMMA_BF and GAMMA_FORW demonstrated strong predictive performance, selecting seven features, striking a balance between methods that selected almost all features but had lower predictive performance (e.g., SVM-RFE with 87.17% sensitivity vs. 95.08% for γ-metric based methods) and RFI, which selected only one feature but had a very low sensitivity of 56.52%, compared to 95.08% for the proposed methods. Only CFS outperformed them, selecting three features and achieving a high AUC of 99.51%, sensitivity of 98.61%, and specificity of 96.33%. For AF detection, [34] reported the accuracy of general practitioners (92% specificity and 80% sensitivity) and practice nurses (85% specificity and 77% sensitivity) in distinguishing AF from NSR cases (99% specificity and 83% sensitivity). Similarly, [35] described the performance of a computer-based algorithm for diagnosing primary cardiac rhythms, including AF, with a specificity of 98.9% and sensitivity of 90.8% for AF diagnosis.

Despite their effectiveness, the three methods based on the proposed γ-metric methods can still be improved, as they have certain limitations. In all scenarios, they selected some non-informative features, with GAMMA_BACK being particularly prone to this. In Scenario 1, the number of selected non-informative features was very low and did not impact the predictive performance of the models. This was more pronounced in Scenario 2, where GAMMA_BACK selected nearly half the features, compared to GAMMA_BF and

<table><caption>**Table 9** Summary of the best results and key conclusions for each scenario in the simulation study</caption><thead><tr><td></td><td>Best results</td><td>Conclusions</td></tr></thead><tbody><tr><td>Scenario 1: Informative features with different effect size.</td><td>γ-metric methods selected almost every time each informative features and very few non-informative features.</td><td>In classical feature selection tasks, the proposed methods outperformed the other methods.</td></tr><tr><td>Scenario 2: Large dataset with balanced/unbalanced classes and weak/strong separation.</td><td>GAMMA_FORW and GAMMA_BF could select very few features among which mostly the informative ones for the strong separation cases.</td><td>GAMMA_FORW and GAMMA_BF were able to perform well with few observations w.r.t the number of features.</td></tr><tr><td>Scenario 3: Different levels of correlation between informative and non-informative features.</td><td>GAMMA_FORW and GAMMA_BF could select all informative features and less non-informative features with low correlation levels.</td><td>GAMMA_BF and GAMMA_FORW were able to disregard non-informative features that were less correlated to the informative features.</td></tr></tbody></table>
GAMMA_FORW, which performed better in cases of strong class separation. This is likely because the γ-metric value is predominantly influenced by informative features, and adding or removing non-informative features does not significantly affect its value. GAMMA_BACK, a backward search method, may terminate prematurely, retaining non-informative features if removing them doesn't drastically change the γ-metric. In contrast, GAMMA_FORW's forward search direction stops adding features once all informative ones are selected, which helps avoid non-informative features. In both methods, the γ-metric is primarily impacted by the inclusion or exclusion of informatives features. Another limitation is the restriction to numerical features, as the γ-metric requires covariance matrix calculations.

Future work could address these limitations. One potential improvement is to introduce a penalty in the γ-metric calculation, based on the number of features in the set. This would favor smaller feature sets and help the algorithm ignore non-informative features more effectively. The γ-metric would then account for both a distance criterion and the dimensionality of the feature subsets. Additionally, exploring alternative search directions may improve performance. As demonstrated, search strategies have a significant impact on results, even when using the same evaluation function. Future search could explore method like genetic algorithms [36] combined with the γ-metric to avoid local optima by introducing randomness during the search.

Lastly, expanding the γ-metric to handle qualitative features would be a valuable direction for future research. For example, existing methods for calculating covariance for categorical features could be adapted [37] proposed a variance definition for categorical features, while [38] discussed the use of polychoric correlation [39] for ordinal features in principal component analysis.

## Conclusions

Based on the results from both the simulation study and the healthcare application, the proposed feature selection methods utilizing the γ-metric as an evaluation function were effective in identifying informative features. Combining the γ-metric with a forward search strategy, such as in GAMMA_FORW, produced excellent results in traditional feature selection tasks and reasonable performance in scenarios with high feature correlation and large datasets. However, the backward search method, GAMMA_BACK, was more prone to getting stuck in local optima, resulting in the selection of more features than necessary.

## Appendix A: Computation of the γ-metric

Let S be a set of n observations noted {$X_i$}$_{i=1,...,n}$, characterised by p features, where $X_i \in S \subset \mathbb{R}^p$. S is divided into K classes such that we have an integer vector Y where $Y_i = 1, ..., K$ $\forall i = 1, ..., n$. For each k ∈ {1, ..., K}, $W_{k,p}$ is the covariance matrix of the corresponding sub-sample of observations belonging to class k:

$$W_{k,p} = \mathrm{Cov}(X_i|Y_i = k) \quad (4)$$

where $W_{k,p}$ is a diagonizable $p \times p$ symmetrical positive semi-definite matrix in which all eigenvalues {$\lambda_{k,j}$}$_{j=1,...,p}$ are positive ($\forall j = 1, ..., p$, $\lambda_{k,j} \ge 0$). Let {$u_{k,j}$}$_{j=1,...,p}$ be the normalized eigenvectors associated with eigenvalues {$\lambda_{k,j}$}$_{j=1,...,p}$. These eigenvectors represent the directions of the p axes of length $\sqrt{\lambda_{k,j}}$ in a p-dimensional ellipsoid centred in $\mu_k$, which is the mean vector of observations in class k. Each class k (with k ∈ {1, ..., K}) is thus represented by an ellipsoid in $\mathbb{R}^p$. Let $k_1, k_2 \in \{1, ..., K\}$ such that $k_1 \ne k_2$ along the mean-mean axis given by $\mu_{k_1} \mu_{k_2} = \mu_{k_1} - \mu_{k_2}$ can be defined as follows:

$$d_{k_1,k_2} = \frac{1}{\alpha_{k_1,k_2}} (||\mu_{k_1} \mu_{k_2}|| - (d_{k_1,k_1k_2} + d_{k_2,k_1k_2})),$$

where $\alpha_{k_1,k_2}$ is a normalisation factor defined as:

$$\alpha_{k_1,k_2} = \sqrt{\sum_{j=1}^{p} \lambda_{k_1,j}} + \sqrt{\sum_{j=1}^{p} \lambda_{k_2,j}},$$

and $d_{k_1,k_1k_2}$ and $d_{k_2,k_1k_2}$ are defined as:

$$d_{k_1,k_1k_2} = \frac{1}{\sqrt{\sum_{j=1}^{p} \frac{\tilde{\mu}_{k_1,j}^2}{\lambda_{k_1,j}}}} \quad \text{and} \quad d_{k_2,k_1k_2} = \frac{1}{\sqrt{\sum_{j=1}^{p} \frac{\tilde{\mu}_{k_2,j}^2}{\lambda_{k_2,j}}}},$$

where $\tilde{\mu}_{k_1,j}$ (respectively $\tilde{\mu}_{k_2,j}$) represents the coordinates of the normalised vector $\mu_{k_1}$, $\mu_{k_2}$ expressed in the orthogonal basis formed by the eigenvectors of ellipsoid $k_1$ (respectively $k_2$). If $U_{k_1}$ (respectively $U_{k_2}$) is the matrix whose columns correspond to the eigenvectors of ellipsoid $k_1$ (respectively $k_2$), then the normalized mean-mean vector $\tilde{\mu}_{k_1}$ (respectively $\tilde{\mu}_{k_2}$) can be written as:

$$\tilde{\mu}_{k_1} = U_{k_1}^{-1} \frac{\mu_{k_1} \mu_{k_2}}{||\mu_{k_1} \mu_{k_2}||_2} \quad \text{and} \quad \tilde{\mu}_{k_2} = U_{k_2}^{-1} \frac{\mu_{k_1} \mu_{k_2}}{||\mu_{k_1} \mu_{k_2}||_2}.$$

$d_{k_1,k_1k_2}$ represents the distance between $\mu_{k_1}$ and the border of the ellipsoid. Any point on this border is determined by drawing a segment between $\mu_{k_1}$ and the border of the ellipsoid in the same direction as $\mu_{k_1} \mu_{k_2}$. Similarly, $d_{k_2,k_1k_2}$ is the distance between $\mu_{k_2}$ and the border of the ellipsoid, and any point on this border is determined by drawing a segment between $\mu_{k_2}$ and the border of the ellipsoid in the same direction as vector $-\mu_{k_1} \mu_{k_2}$. Finally,
the γ-metric for a set of K classes of observations in S ⊂ ℝP is defined as follows:

$$\gamma_p = \sum_{k_1=1}^{K} \sum_{k_2<k}^{K} d_{k_1,k_2}.$$

# Appendix B: Conventional feature selection methods' computation

## Chi-squared

The Chi-squared statistic is computed for each discretized feature and class as follows:

$$\chi^2 = \sum_{i=1}^{2} \sum_{j=1}^{k} \frac{(A_{ij} - E_{ij})^2}{E_{ij}},$$

With k being the number of classes, Aij the number of observations in interval i and class j, and Eij the expected frequency of Aij. Cramer's V is used to obtain values between 0 and 1, with 1 being interpreted as a very strong correlation between the two features.

$$E_{ij} = \frac{R_i \times C_j}{N},$$

with R_i = Σ_{j=1}^{k} A_{ij} the number of observations in interval i, C_j = Σ_{i=1}^{k} A_{ij} the number of observations in class j, and N = Σ_{i=1}^{2} R_i the total number of observations.

## Correlation-based feature selection

To compute the importance score of the subsets we compute:

$$M = \frac{p_S \bar{r}_{cf}}{\sqrt{p_S + p_S(p_S - 1)\bar{r}_{ff}}},$$

with p_S being the number of features in subset S, $\bar{r}_{cf}$ the average feature-class correlation, and $\bar{r}_{ff}$ is the average feature-feature correlation.

## Symmetrical uncertainty

Symmetrical uncertainty method is a variant of mutual information, or information gain, to lower the bias of the features with a large number of different values. To do so we operate a normalization of the mutual information. The mutual information measures the dependence between a feature $f_i$ and class C and being computed as follows:

$$\text{IG}(f_i, C) = H(f_i) - H(f_i|C),$$

with $H(f_i)$ is the entropy of feature $f_i$ and $H(f_i|C)$ the conditional entropy of $f_i$ given C, which is defined as follows:

$$H(f_i) = - \sum_{x \in \text{supp}(f_i)} P(x) \log_2(P(x)),$$

$$H(f_i|C) = - \sum_{k \in \{1, \dots, K\}} P(c_k) \sum_{x \in \text{supp}(f_i)} P(x|c_k) \log_2(P(x|c_k)),$$

with $P(x)$ being the probability that $f_i$ takes value $x$, $\text{supp}(f_i)$ the support of $f_i$, $P(c_k)$ the probability that C takes value $c_k$, and $P(x|c_k)$ the probability that $f_i$ takes the value $x$ when C takes the value $c_k$. The symmetrical uncertainty is a normalize variant of $\text{IG}(f_i, C)$:

$$\text{SU}(f_i, C) = 2 \frac{\text{IG}(f_i, C)}{H(f_i) + H(C)}.$$

# Appendix C: technique for order preference by similarity to ideal solution

The TOPSIS method can be implemented using the following procedure:

**Step 1:** Normalize the indicator values

$$r_{ij} = \frac{x_{ij}}{\sqrt{\sum_{k=1}^{J} x_{ik}^2}}, \qquad (5)$$

with $x_{ij}$ represents the performance of feature selection method j on indicator i. The total number of methods is J, the total number of indicators is m.

**Step 2:** Calculate the weighted indicator values:

$$v_{ij} = w_i r_{ij}, \qquad (6)$$

where $w_i$ is the weight assigned to indicator i.

**Step 3:** Define the ideal and worst solutions

Each indicator is classified as either benefit indicator (high values indicate better performance) or a cost indicator (lower values indicate better performance). The ideal solution $S^+$ maximizes all benefit indicators and minimizes all cost indicators, while the worst solution $S^-$ minimizes all benefit indicators and maximizes all cost indicators. The ideal and worst solutions can be defined as:

$$S^+ = \{ \langle \max_j v_{ij} | i \in I' \rangle, \langle \min_j v_{ij} | i \in I'' \rangle \} = \{ v_i^+ | i \in \{1, \dots, m\} \} \qquad (7)$$

$$S^- = \{ \langle \min_j v_{ij} | i \in I' \rangle, \langle \max_j v_{ij} | i \in I'' \rangle \} = \{ v_i^- | i \in \{1, \dots, m\} \}, \qquad (8)$$

where $I'$ and $I''$ represent the sets of benefit and cost indicators, respectively.

**Step 4:** Compute the distances
For each feature selection method j, we calculate the Euclidean distance to both the ideal solution $D_j^+$ and the worst solution $D_j^-$:

$$D_j^+ = \sqrt{\sum_{i=1}^{m} (v_{ij} - v_i^+)^2} \quad \text{and} \quad D_j^- = \sqrt{\sum_{i=1}^{m} (v_{ij} - v_i^-)^2}. \quad (9)$$

Step 5: Compute the relative closeness to the ideal solution. The relative closeness $R_j^+$ to the ideal solution is calculated as:

$$R_j^+ = \frac{D_j^-}{D_j^+ + D_j^-} \quad (10)$$

A high value of $R_j^+$ indicates that the method is closer to the ideal solution, with $R_j^+ = 1$ representing the best possible performance and $R_j^+ = 0$ representing the worst.

Step 6: Rank the feature selection methods

The feature selection methods are ranked in decreasing order based on their $R_j^+$ values. The method with the highest $R_j^+$ offers the best compromise across all indicators.

To compute the $R_j^+$ values, we must first define the weights $w_i$ for each indicator. These weights reflect the emphasis placed on each indicator when assessing the performance of feature selection methods. Each weight $w_i$ lies within the range [0, 1] and must satisfy the condition $\sum_{i=1}^{m} w_i = 1$.

## Appendix D: Justification of the choice of the $\beta$ values

In scenario 2, we generated classes based on two levels of balance (balanced/unbalanced) and separation (strong/ weak), using three features with non-null effects and an intercept, $\beta_0$. As outlined in the data generation section, we computed a probability $\rho_i$ for each observation. Each observation's class was determined through a Bernoulli process with parameter $\rho_i$. This section explains the choice of $\beta = (\beta_0, \beta_1, \beta_2, \beta_3)$ values according to the scenario we aimed to generate.

To recap, the probability $\rho_i$ for each observation $i \in \{1, \dots, n\}$ is given by:

$$\rho_i = \frac{\exp(X_i \beta)}{1 + \exp(X_i \beta)} = \frac{1}{1 + \exp(-X_i \beta)}.$$

For simplicity, we focus on two features $X_1$ and $X_2$ with non-null effects. We express the equation of the decision boundary where $\rho_i = \frac{1}{2}$, which represents observations with a 50% probability of being classified as 1. Thus $\forall i \in \{1, \dots, n\}$, we have:

$$\frac{1}{1 + \exp\left\{-\left(\beta_0 + x_{i1}\beta_1 + x_{i2}\beta_2\right)\right\}} = \frac{1}{2}$$

This simplifies to:

$$\forall \beta_2 \neq 0, x_2 = -\frac{\beta_0}{\beta_2} - \frac{\beta_1}{\beta_2}x_1,$$

Observations near this line have an approximately 50% probability of being in class 1. The values of $\beta = (\beta_0, \beta_1, \beta_2)$ determine the slope (i.e., $-\frac{\beta_1}{\beta_2}$) and intercept (i.e., $-\frac{\beta_0}{\beta_2}$) of this line, thereby influencing the distribution of the two classes.

In the case of a balanced dataset, where observations are centered around 0, the decision line should pass through the origin (0, 0). To achieve this, the intercept must be zero, i.e., $-\frac{\beta_0}{\beta_2} = 0$. By varying values of $\beta_0$, we control the degree of class imbalance.

For example, as shown in Fig. 4, when only the values of $\beta_0$ change (top panels of Fig. 4), the decision line shifts, reducing the balance between the two classes. Conversely, when $\beta_0 = 0$, altering $\beta_1$ and $\beta_2$ does not affect class balance, as all observations are drawn from a standard normal distribution.

<figure><img src="image_7.png" /><figcaption>**Fig. 4** Illustration of data generation with different values of $\beta$ and two features. The straight black line shows the coordinates of points where $\rho_i = 1/2$. The percentage of observations labeled as 1 at the end of the generation procedure is displayed in the bottom-right corner. The top panels depict the class of each observation for varying values of $\beta_0$, while $\beta_1$ and $\beta_2$ are held constant. The bottom panels show the class of each observation with $\beta_0$ fixed at zero and, varying values of $\beta_1$ and $\beta_2$</figcaption></figure>

For class separation, we used the Bernoulli distribution to generate random class labels based on the probability $\rho_i$ calculated for each observation. Rather than using a strict decision rule such as $Y_i = 1$ when $\rho_i \ge 0.5$, we allowed observations with $\rho_i \ge 0.5$ to still be classified as 0, introducing overlap between the classes. When the probabilities are concentrated around 0 or 1, the overlap is minimal, indicating strong separation.
<figure><img src="image_8.png" /><figcaption>**Fig. 5** Distribution of the probability $\rho_i$ based on the chosen $\beta$ values, illustrating strong (left panel) and weak (right panel) separation between classes</figcaption></figure>

To demonstrate this, we used the same matrix $X$ to generate the results shown in Fig. 5, which compares strong separation (left panel) and weak separation (right panel). In the strong separation case, the distribution of $\rho_i$ is concentrated near 0 and 1, while in the weak separation case, the distribution is flatter. In both scenarios, approximately 50% of the observations are labeled as 1. Thus, we control the degree of class overlap (i.e., separability) by adjusting the distribution of $\rho_i$.

## Appendix E: Visualisation of the correlation matrices and $\beta$ vector chosen in scenario 3 of the simulation study

In scenario 3, we generated data with multicollinearity by using a multivariate Gaussian distribution with a specific covariance matrix, $\Sigma$, which was structured as a block diagonal matrix. The following figures illustrate the matrix $\Sigma$ and the block matrix $\Sigma_g$ for a group of features under conditions of constant (Fig. 6) and non-constant (Fig. 7) correlation levels.

<figure><img src="image_9.png" /><figcaption>**Fig. 6** The $\Sigma$ matrix (top panels) used for data generation in scenario 3, illustrating constant correlation at different levels. The $\Sigma_g$ matrix (bottom panels) represents the correlation structure within a group of features</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 7** The $\Sigma$ matrix (top panels) used for data generation in scenario 3, showing non-constant correlations at different levels. The $\Sigma_g$ matrix (bottom panels) represents the correlation structure within a group of features</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 8** Illustration of the $\beta$ vector used in the generation process for scenario 3. Blue cells indicate $\beta$ coefficients corresponding to informative features, while red cells indicate coefficients corresponding to non-informative features</figcaption></figure>

Additionally, we employed a specific $\beta$ vector with values of 1.5 or 0, distinguishing between informative and non-informative features. As shown in Fig. 8, each group contained ten features. In the first five groups, only the first feature was informative, while the rest were non-informative. Groups 6 through 9 consisted entirely of non-informative features, and group 10 was independent, with the first feature being informative and the others non-informative. As a result, only $\beta_1$, $\beta_{11}$, $\beta_{21}$, $\beta_{31}$, $\beta_{41}$ and $\beta_{91}$ were non-zero. The intercept term, $\beta_0$, was set to 0.

### Acknowledgements

Not applicable.

### Authors' contributions

NN, PM, and RG designed the study. NN and PM designed the computational framework. NN wrote the first version of the manuscript. RG and PM contributed to the analysis of the data, the interpretations of the findings, and drafting the manuscript. All authors read and approved the final manuscript.
## Funding

The project leading to this publication has received funding from the French government under the "France 2030" investment plan managed by the French National Research Agency (reference: ANR-17-EURE-0020) and from Excellence Initiative Aix-Marseille University - A*MIDEX.

## Data availability

The dataset analyzed for the healthcare application is available from the corresponding author upon reasonable request. The R code used to generate the simulation study datasets, as well as to perform the analyses, and produce the tables and figures in this manuscript, is available in the repository https://github.com/NicolasNgo/Multivariate-filter-methods-for-feature-selection-with-the-gamma-metric.

## Declarations

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors declare no competing interests.

Received: 27 September 2024 Accepted: 26 November 2024

Published online: 19 December 2024

## References

1. Feinberg WM, Blackshear JL, Laupacis A, Kronmal R, Hart RG. Prevalence, age distribution, and gender of patients with atrial fibrillation: analysis and implications. Arch Intern Med. 1995;155(5):469–73.

2. Benjamin EJ, Wolf PA, D'Agostino RB, Silbershatz H, Kannel WB, Levy D. Impact of atrial fibrillation on the risk of death: the Framingham Heart Study. Circulation. 1998;98(10):946–52.

3. Wolf PA, Abbott RD, Kannel WB. Atrial fibrillation as an independent risk factor for stroke: the Framingham Study. Stroke. 1991;22(8):983–8.

4. Guo Y, Lane DA, Wang L, Zhang H, Wang H, Zhang W, et al. Mobile health technology to improve care for patients with atrial fibrillation. J Am Coll Cardiol. 2020;75(13):1523–34.

5. Turakhia MP, Kaiser DW. Transforming the care of atrial fibrillation with mobile health. J Interv Card Electrophysiol. 2016;47:45–50.

6. McConnell MV, Turakhia MP, Harrington RA, King AC, Ashley EA. Mobile health advances in physical activity, fitness, and atrial fibrillation: moving hearts. J Am Coll Cardiol. 2018;71(23):2691–701.

7. Kumar V, Minz S. Feature selection. SmartCR. 2014;4(3):211–29.

8. Pons JF, Haddi Z, Deharo JC, Charaï A, Bouchakour R, Ouladsine M, et al. Heart rhythm characterization through induced physiological variables. Sci Rep. 2017;7(1):5059.

9. Michel P, Ngo N, Pons JF, Delliaux S, Giorgi R. A filter approach for feature selection in classification: application to automatic atrial fibrillation detection in electrocardiogram recordings. BMC Med Inform Decis Making. 2021;21(Suppl 4):130.

10. Dash M, Liu H. Feature selection for classification. Intell Data Anal. 1997;1(1–4):131–56.

11. Chandrashekar G, Sahin F. A survey on feature selection methods. Comput Electr Eng. 2014;40(1):16–28.

12. Ledoit O, Wolf M. A well-conditioned estimator for large dimensional covariance matrices. J Multivar Anal. 2004;88(2):365–411.

13. Schäfer J, Strimmer K. A shrinkage approach to large-scale covariance matrix estimation and implications for functional genomics. Stat Appl Genet Mol Biol. 2005;4(1):Article32.

14. Schäfer J, Opgen-Rhein R, Zuber V, Ahdesmaki M, Silva APD, Strimmer K. corpcor: Efficient Estimation of Covariance and (Partial) Correlation.

2021. R package version 1.6.10. https://CRAN.R-project.org/package=corpcor. Accessed 26 Sept 2024.

15. Jin X, Xu A, Bie R, Guo P. Machine learning techniques and chi-square feature selection for cancer classification using SAGE gene expression profiles. In: Data Mining for Biomedical Applications: PAKDD 2006 Workshop, BioDM 2006, Singapore, April 9, 2006. Proceedings. Berlin, Heidelberg: Springer; 2006. p. 106–15.

16. Hall MA. Correlation-based feature selection for machine learning [PhD Thesis]. The University of Waikato; 1999.

17. Romanski P, Kotthoff L, Schratz P. FSelector: Selecting Attributes. 2023. R package version 0.34. https://CRAN.R-project.org/package=FSelector. Accessed 26 Sept 2024.

18. Tibshirani R. Regression shrinkage and selection via the lasso. J R Stat Soc Ser B Stat Methodol. 1996;58(1):267–88.

19. Han H, Guo X, Yu H. Variable selection using mean decrease accuracy and mean decrease gini based on random forest. In: 2016 7th IEEE International Conference on Software Engineering and Service Science (ICSESS). Beijing, CN: IEEE; 2016. p. 219–24.

20. Witten IH, Frank E. Data mining: practical machine learning tools and techniques with Java implementations. Acm Sigmod Record. 2002;31(1):76–7.

21. Quinlan JR. Induction of decision trees. Mach Learn. 1986;1:81–106.

22. Guyon I, Weston J, Barnhill S, Vapnik V. Gene selection for cancer classification using support vector machines. Mach Learn. 2002;46:389–422.

23. Boser BE, Guyon IM, Vapnik VN. A training algorithm for optimal margin classifiers. In: Proceedings of the fifth annual workshop on Computational learning theory. New York: Association for Computing Machinery; 1992. p. 144–52.

24. Vapnik V. Statistical Learning Theory. New-York: Wiley-Interscience; 1998.

25. Friedman J, Tibshirani R, Hastie T. Regularization Paths for Generalized Linear Models via Coordinate Descent. J Stat Softw. 2010;33(1):1–22. https://doi.org/10.18637/jss.v033.i01.

26. Tay JK, Narasimhan B, Hastie T. Elastic Net Regularization Paths for All Generalized Linear Models. J Stat Softw. 2023;106(1):1–31. https://doi.org/10.18637/jss.v106.i01.

27. Lang M, Binder M, Richter J, Schratz P, Pfisterer F, Coors S, et al. mlr: A modern object-oriented machine learning framework in R. J Open Source Softw. 2019. https://doi.org/10.21105/joss.01903.

28. Jaccard P. Distribution de la flore alpine dans le bassin des Dranses et dans quelques régions voisines. Bull Soc Vaudoise Sci Nat. 1901;37:241–72.

29. Halsteinslid E. Addressing collinearity and class imbalance in logistic regression for statistical fraud detection [Master's Thesis]. University of Oslo; 2019.

30. Saeks Y, Abeel T, Van de Peer Y. Robust feature selection using ensemble feature selection techniques. In: Machine Learning and Knowledge Discovery in Databases: European Conference, ECML PKDD 2008, Antwerp, Belgium, September 15-19, 2008, Proceedings, Part II 19. Berlin, Heidelberg: Springer; 2008. p. 313–25.

31. Hwang CL, Yoon KP. Multiple decision making: Methods and applications. New York: Springer-Verlag; 1981.

32. Kou G, Yang P, Peng Y, Xiao F, Chen Y, Alsaadi FE. Evaluation of feature selection methods for text classification with small datasets using multiple criteria decision-making methods. Appl Soft Comput. 2020;86:105836.

33. Rodríguez-Liñares L, Méndez AJ, Lado MJ, Olivieri DN, Vila X, Gómez-Conde I. An open source tool for heart rate variability spectral analysis. Comput Methods Prog Biomed. 2011;103(1):39–50.

34. Mant J, Fitzmaurice DA, Hobbs FR, Jowett S, Murray ET, Holder R, et al. Accuracy of diagnosing atrial fibrillation on electrocardiogram by primary care practitioners and interpretative diagnostic software: analysis of data from screening for atrial fibrillation in the elderly (SAFE) trial. Brmj. 2007;335(7616):380.

35. Poon K, Okin PM, Kligfield P. Dianostic performance of a computer-based ECG rhythm algorithm. J Electrocardiol. 2005;38(3):235–8.

36. Leardi R, Boggia R, Terrible M. Genetic algorithms as a strategy for feature selection. J Chemom. 1992;6(5):267–81.

37. Gini C. Variabilità e mutabilità: contributo allo studio delle distribuzioni e delle relazioni statistiche. [Fasc. I.]. Bologna: Tipogr. di P. Cuppini; 1912.
38. Kolenikov S, Angeles G. Socioeconomic status measurement with discrete proxy variables: Is principal component analysis a reliable answer? Rev Income Wealth. 2009;55(1):128–65.

39. Olsson U. Maximum likelihood estimation of the polychoric correlation coefficient. Psychometrika. 1979;44(4):443–60.

**Publisher’s Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.