Skip to main content
Have a personal or library account? Click to login
A Novel Hybrid Feature Selection Framework with Dimensionality Reduction for Early Cardiovascular Disease Detection Cover

A Novel Hybrid Feature Selection Framework with Dimensionality Reduction for Early Cardiovascular Disease Detection

Open Access
|Aug 2026

Full Article

I. Introduction

Early diagnosis of cardiovascular disease (CVD) enables healthcare professionals to identify high-risk individuals before severe complications occur. Traditional diagnostic procedures typically involve clinical examinations, laboratory investigations, electrocardiograms, echocardiography, angiography, and physician-based assessments. Although these methods provide valuable clinical insights, they are often time-consuming, expensive and dependent on expert interpretation. Furthermore, CVDs frequently develop gradually over time, making it challenging to detect subtle patterns and risk factors through conventional diagnostic approaches alone [1,2,3]. The growing volume of electronic health records, laboratory measurements, demographic information, and physiological parameters presents an opportunity to utilize advanced computational techniques to identify hidden relationships within medical data. Therefore, intelligent predictive systems capable of analyzing complex healthcare datasets have become increasingly important in supporting clinicians and facilitating timely interventions [4,5].

The rapid advancement of artificial intelligence (AI) and machine learning (ML) technologies has transformed numerous sectors, including health care, by enabling data-driven decision-making and predictive analytics. ML algorithms possess the capability to learn intricate patterns from historical medical records and generate accurate predictions regarding disease occurrence and progression. In CVD prediction, ML techniques such as logistic regression (LR), Decision Trees, random forests (RFs), support vector machines (SVMs), gradient boosting (GB) Machines, and extreme gradient boosting (XGBoost) have demonstrated remarkable potential in identifying patients at risk of developing cardiovascular complications [6,7,8]. These algorithms can process large volumes of heterogeneous clinical data and uncover complex non-linear relationships among risk factors that may not be readily apparent through traditional statistical methods. As a result, ML-based diagnostic systems have become valuable tools for enhancing clinical decision support and improving healthcare outcomes.

Despite significant progress in ML-based CVD prediction, several challenges continue to hinder the development of highly accurate and efficient predictive models. One of the primary challenges is the high dimensionality of healthcare datasets. Modern medical databases often contain numerous clinical attributes, laboratory measurements, physiological indicators, lifestyle factors, and demographic variables. While these features collectively provide comprehensive patient information, not all attributes contribute equally to disease prediction [9,10,11,12]. Many variables may be redundant, irrelevant, noisy, or highly correlated with other features. The inclusion of such unnecessary information can negatively impact model performance by increasing computational complexity, introducing overfitting, and reducing generalization capability. Consequently, predictive models trained on high-dimensional datasets may exhibit suboptimal performance despite the availability of extensive data [13,14,15].

Another critical issue associated with high-dimensional medical data is the curse of dimensionality. As the number of features increases, the volume of the feature space expands exponentially, making it increasingly difficult for ML algorithms to identify meaningful patterns and relationships. This phenomenon often leads to increased training times, higher memory requirements, and degraded predictive accuracy. Moreover, redundant and irrelevant attributes may obscure the contribution of clinically significant risk factors, thereby reducing model interpretability. Since healthcare applications require reliable and explainable predictions, addressing the challenges associated with high-dimensional data has become a key research priority in medical data analytics.

Feature optimization techniques have emerged as effective solutions for overcoming the limitations of high-dimensional healthcare datasets. Feature selection aims to identify the most informative attributes while eliminating irrelevant and redundant variables. By reducing the number of input features, feature selection enhances model efficiency, decreases computational cost, and improves classification performance. Existing feature selection approaches are generally categorized into filter methods, wrapper methods, and embedded methods. Filter methods evaluate feature relevance independently of ML models, whereas wrapper methods utilize classifier performance to assess feature subsets. Embedded methods perform feature selection during model training. Although these techniques have demonstrated promising results individually, each possesses inherent limitations. Filter methods may overlook feature interactions, while wrapper methods often require substantial computational resources. Consequently, combining multiple feature selection strategies within a hybrid framework can leverage the strengths of different approaches and achieve superior optimization performance.

In addition to feature selection, dimensionality reduction techniques play a vital role in managing high-dimensional medical data. Dimensionality reduction transforms the original feature space into a compact representation while preserving the most significant information contained within the dataset. Principal component analysis (PCA) is one of the most widely used dimensionality reduction techniques due to its ability to capture maximum variance and eliminate multicollinearity among features. By reducing data complexity and removing redundant information, dimensionality reduction contributes to improved model stability, faster training, and enhanced predictive performance. However, dimensionality reduction alone may not always identify the most clinically relevant attributes. Therefore, integrating dimensionality reduction with feature selection provides a comprehensive strategy for constructing highly informative feature subsets [16,17,18,19,20].

The motivation for this research stems from the need to develop an intelligent and efficient CVD prediction framework that is capable of addressing the limitations of existing ML approaches. While numerous studies have employed ML algorithms for CVD diagnosis, many rely on either standalone feature selection or dimensionality reduction methods. Such approaches may fail to fully exploit the complementary advantages offered by combined feature optimization strategies. Furthermore, existing models often prioritize prediction accuracy without adequately considering computational efficiency and feature relevance. To bridge these gaps, this study proposes a Novel Hybrid Feature Selection Framework with Dimensionality Reduction for Early CVD Detection. The proposed framework integrates PCA for dimensionality reduction with a hybrid feature selection mechanism that combines filter-based and wrapper-based optimization strategies. The resulting optimized feature subset is subsequently evaluated using multiple ML classifiers to identify the most effective predictive model.

By integrating dimensionality reduction and hybrid feature optimization within a unified framework, the proposed approach seeks to enhance predictive accuracy, reduce computational overhead, improve model interpretability, and facilitate reliable clinical decision support. The framework is designed to assist healthcare practitioners in identifying high-risk patients at an early stage, thereby enabling timely intervention and reducing the burden associated with CVDs. Through comprehensive experimentation and comparative analysis, this research aims to demonstrate the effectiveness of feature optimization techniques in improving ML-based CVD prediction systems.

a. Contributions of the proposed work

The major contributions of this study are summarized as follows:

  1. A hybrid feature selection methodology combining filter-based and wrapper-based techniques is developed to identify the most informative cardiovascular risk factors while eliminating redundant and irrelevant features.

  2. PCA is incorporated with hybrid feature selection to effectively address the challenges associated with high-dimensional medical datasets and improve data representation.

  3. The optimized feature subset is evaluated using various ML classifiers, including LR, SVM, RF, GB, and XGBoost, to determine the most effective prediction model.

  4. The proposed framework enhances CVD prediction accuracy while simultaneously reducing feature dimensionality, training time, memory consumption, and overall computational cost.

  5. The framework provides a reliable and efficient decision-support mechanism capable of assisting healthcare professionals in identifying CVD risks at an early stage, thereby enabling timely preventive interventions and improved patient outcomes.

II. Literature Review

CVD remains one of the most prevalent and life-threatening medical conditions worldwide, prompting extensive research into intelligent prediction systems capable of facilitating early diagnosis and intervention. Recent advancements in AI, ML, deep learning (DL), feature selection, and ensemble learning have significantly enhanced the accuracy and reliability of CVD prediction models. Despite these developments, challenges associated with high-dimensional medical data, feature redundancy, computational complexity, and model interpretability continue to motivate further investigation. This section reviews recent studies related to CVD prediction, feature optimization, dimensionality reduction, and hybrid ML frameworks.

Recent studies have demonstrated the effectiveness of ML algorithms in identifying CVD risk factors from clinical datasets. Raman et al. [14] proposed an enhanced CVD prediction framework that integrated advanced ML algorithms with hybrid feature selection techniques. Their study demonstrated that the selection of relevant clinical attributes significantly improved classification performance by reducing the influence of redundant and noisy features. The authors highlighted the importance of combining multiple feature selection approaches to achieve better predictive accuracy compared with conventional methods.

Ahmed and Husien [3] presented a comprehensive review of hybrid ML approaches for heart disease prediction. Their study examined various combinations of ML classifiers and optimization techniques, emphasizing that hybrid frameworks generally outperform standalone classifiers due to their ability to capture complex relationships among clinical variables. The review identified feature optimization as one of the most critical factors influencing prediction accuracy and generalization performance.

Similarly, Abdullahi et al. [2] investigated ML-based CVD prediction using advanced feature selection methodologies. Their findings revealed that eliminating irrelevant attributes improved classifier efficiency and reduced overfitting. The study further demonstrated that optimized feature subsets enhanced the predictive capabilities of traditional ML models such as LR, SVMs, and RFs.

The increasing adoption of ensemble learning techniques has further improved CVD diagnosis. Al-Sayed et al. [4] proposed a collaborative clustering and ensemble learning framework for heart disease prediction. Their methodology combined clustering-based preprocessing with ensemble classifiers to improve disease classification performance. Experimental results indicated that ensemble models effectively handled data variability and produced more reliable predictions than individual classifiers.

Likewise, Mandula and Vijaya Kumar [9] introduced a two-tier heart disease detection framework integrating feature selection with ensemble learning. Their approach employed an adaptive learning artificial network and ensemble tree-based adaptive boosting decision forest models to improve classification accuracy. The study demonstrated that combining feature optimization with ensemble learning significantly enhanced predictive performance while maintaining computational efficiency.

Further advancements in ensemble-based CVD prediction were reported by Khan et al. [7], who developed EnsCVDD-Net and BlCVDD-Net, two novel ensemble and blending-based CVD detection networks. Their results demonstrated superior accuracy, precision, and recall compared with traditional ML approaches. The study highlighted the capability of ensemble architectures to capture non-linear relationships among clinical features and improve disease classification reliability.

DL has also emerged as a powerful tool for CVD diagnosis due to its ability to automatically learn hierarchical feature representations from complex medical datasets. Rao et al. (2024) proposed AttGRU-HMSI, a hybrid DL framework integrating attention mechanisms with gated recurrent units for heart disease diagnosis. The model effectively captured temporal and non-linear dependencies within clinical data, resulting in substantial improvements in diagnostic accuracy. The incorporation of attention mechanisms further enhanced model interpretability by identifying influential risk factors.

Jayasudha et al. [5] developed a hybrid optimization-enabled DL ensemble classification framework for heart disease detection. Their approach integrated optimization algorithms with DL classifiers to enhance feature representation and classification performance. Experimental findings demonstrated superior accuracy and robustness compared with conventional DL models, highlighting the benefits of optimization-driven learning architectures.

Similarly, Kallimani et al. [6] proposed a hybrid DL-based heart disease detection framework incorporating feature selection mechanisms within an e-healthcare environment. The study emphasized the significance of selecting informative clinical features prior to model training. Results indicated that hybrid feature selection considerably improved classification performance while reducing computational overhead.

Subramani et al. [17] explored the integration of ML and DL techniques for CVD prediction. Their study revealed that combining traditional ML models with deep neural networks could improve disease detection capabilities by leveraging complementary strengths of both approaches. The authors concluded that hybrid architectures offer promising solutions for real-world healthcare applications requiring high predictive accuracy and scalability.

Recent research has also focused on representation learning and advanced ensemble frameworks for CVD diagnosis. Abdellatif et al. [1] introduced a computational heart disease detection system based on conditional variational autoencoders and stacked ensemble learning. The proposed framework generated latent feature representations that enhanced classifier performance. Experimental results confirmed the effectiveness of combining deep feature extraction with ensemble classification for accurate disease prediction.

Mienye and Sun [10] proposed a particle swarm optimization (PSO)-based stacked sparse autoencoder for heart disease prediction. Their approach utilized optimization techniques to improve feature learning and parameter tuning. The study reported substantial gains in prediction accuracy and demonstrated the potential of evolutionary optimization algorithms in healthcare analytics.

Razzaque and Badholia [15] presented a framework combining PCA for feature extraction with Modified PSO (MPSO) for feature selection in gene expression-based medical data classification. Their results demonstrated that dimensionality reduction followed by optimization-driven feature selection significantly enhanced classification accuracy while reducing computational complexity. The study highlights the complementary benefits of integrating dimensionality reduction and feature optimization techniques.

The broader landscape of AI-driven cardiovascular health care has also witnessed significant developments. Paul and Masood [12] conducted a comprehensive survey of predictive methods for CVD, reviewing ML, DL, and ensemble learning approaches. Their analysis revealed that hybrid models generally achieve higher diagnostic accuracy than standalone techniques. The authors emphasized the growing importance of explainable AI, feature optimization, and multimodal learning in future CVD prediction systems.

Milosevic et al. [11] explored applications of AI in multimodal cardiovascular imaging. Their work demonstrated that AI-based models can effectively analyze imaging data from multiple modalities to support CVD diagnosis. Although imaging-based approaches differ from structured clinical data analytics, the study underscored AI’s transformative role in improving diagnostic precision and clinical decision-making.

Apart from healthcare-specific studies, advancements in ML architectures have contributed significantly to predictive analytics. Santosh Kumar and Bharathi [16] proposed ConvTransFormer, a hybrid framework integrating convolutional neural networks and transformers for signal processing applications. Their work highlighted the advantages of combining multiple learning paradigms to improve feature extraction and classification performance. Similar hybrid principles have increasingly been adopted in healthcare prediction models.

III. Proposed HFSDR-CVD Framework

a. Overview of the proposed framework

The proposed Hybrid Feature Selection with Dimensionality Reduction Framework for CVD Detection (HFSDR-CVD) is designed to address two major challenges commonly observed in CVD prediction systems: the presence of high-dimensional clinical data and the existence of redundant or irrelevant attributes that influence classifier performance negatively. Although contemporary ML algorithms have demonstrated promising results in CVD diagnosis, their predictive capability is often limited by noisy feature spaces, multicollinearity among clinical variables, and increased computational complexity. To overcome these limitations, the proposed framework integrates dimensionality reduction and hybrid feature optimization into a unified predictive architecture.

The overall framework consists of four interconnected modules: data preprocessing, principal component-based dimensionality reduction, hybrid feature optimization, and ML-based disease prediction. Initially, raw cardiovascular data undergo preprocessing operations to eliminate inconsistencies and standardize feature distributions. Subsequently, PCA is employed to reduce feature dimensionality by transforming correlated variables into a reduced set of orthogonal components. While PCA effectively removes multicollinearity and preserves the majority of information contained in the dataset, certain transformed components may still contribute minimally toward disease classification. Therefore, a hybrid feature selection mechanism is introduced to further identify the most discriminative features.

The proposed hybrid feature optimization strategy combines filter-based feature relevance analysis and wrapper-based subset optimization. The filter stage evaluates feature importance using mutual information (MI) and correlation analysis, allowing the framework to eliminate weakly associated variables at an early stage. The remaining candidate features are subsequently processed through recursive feature elimination (RFE), which iteratively removes less informative features based on classifier performance. This two-stage optimization process ensures that only highly informative and non-redundant attributes are retained for model training.

Finally, the optimized feature subset is supplied to multiple ML classifiers, including LR, SVM, RF, GB, and XGBoost. Comparative evaluation across multiple classifiers enables identification of the most effective predictive model for CVD detection. The complete architecture of the proposed HFSDR-CVD framework is illustrated in Figure 1.

Figure 1:

Model architecture. GB, gradient boosting; LR, logistic regression; PCA, principal component analysis; RF, random forest; RFE, recursive feature elimination; SVM, support vector machine; XGBoost, extreme gradient boosting.

b. Data preprocessing and feature normalization

The effectiveness of ML-based CVD prediction largely depends on the quality of the input data. Clinical datasets often contain missing observations, noisy measurements, outliers, and features with heterogeneous numerical scales. These issues can adversely affect the learning capability of classifiers, resulting in reduced prediction accuracy and poor generalization performance. Therefore, a comprehensive preprocessing stage is incorporated into the proposed HFSDR-CVD framework to improve data quality and generate a consistent feature representation suitable for subsequent dimensionality reduction and feature optimization.

Let the CVD dataset be represented as

D={(xi,yi)}i=1N
where N denotes the total number of patient records, xi = [xi1, xi2 ,..., xim] represents the feature vector containing m clinical attributes, and yi denotes the corresponding class label indicating the presence or absence of CVD.

b.i. Missing value treatment

Medical datasets frequently contain incomplete records due to unavailable laboratory measurements, data entry errors, or inconsistencies during data collection. The presence of missing values may introduce bias and negatively influence model performance. To ensure data completeness, missing numerical values are replaced through mean-value imputation.

For a feature Xi, the mean value is calculated as

Xi^=1nji=1njXij
where nj denotes the number of valid observations available for the jth feature.

The missing entries are subsequently replaced as

Xij=XijifobservedXj^ifmissing

This process preserves the statistical characteristics of the dataset while eliminating incomplete records.

b.ii. Outlier detection and removal

Outliers are abnormal observations that significantly deviate from the majority of data points. Such extreme values can distort feature distributions and adversely affect classifier learning. Therefore, outlier detection is performed using the interquartile range (IQR) method.

The IQR is computed as

IQR=Q3Q1
where Q1 and Q3 denote the first and third quartiles, respectively.

The lower and upper boundaries for identifying outliers are determined as

LB=Q11.5(IQR)
UB=Q3+1.5(IQR)

A data point is considered an outlier if

X<LBorX>UB

Detected outliers are removed from the dataset to improve data reliability and reduce noise during model training.

b.iii. Feature standardization

Clinical attributes such as cholesterol level, blood pressure, age, fasting blood sugar, and heart rate are measured on different numerical scales. Features with larger magnitudes may dominate the learning process and negatively influence model convergence. To address this issue, feature standardization is applied.

The standardized value of a feature is computed using z-score normalization:

zi=xiμσ
where µ represents the mean value of the feature and σ denotes its standard deviation.

The mean is calculated as

μ=1Ni=1Nxi
while the standard deviation is given by
σ=1Ni=1N(xiμ)2

Standardization transforms all features to a common scale with zero mean and unit variance, thereby improving numerical stability and classifier performance.

b.iv. Min-max normalization

Following standardization, Min-Max normalization is employed to map all feature values into the interval [0,1]. This transformation enhances the efficiency of optimization-based learning algorithms and prevents bias toward features with larger ranges.

The normalized feature value is computed as

x=xxminxmaxxmin
where x denotes the original feature value, and xmin and xmax represent the minimum and maximum values of the feature, respectively.

The transformed feature satisfies

0x1
for all observations.

b.v. Construction of the preprocessed dataset

After missing value treatment, outlier removal, standardization, and normalization, the processed dataset is represented as

Dp=x110x12..x1mx21x22x2mxN1xN2xNm
where Dp denotes the final normalized feature matrix containing N patient records and m preprocessed attributes.

The preprocessing stage produces a clean, standardized, and normalized dataset that minimizes noise, eliminates inconsistencies, and improves feature quality. This optimized dataset serves as the input for the subsequent PCA-based dimensionality reduction and hybrid feature selection modules within the proposed HFSDR-CVD framework.

c. Principal component-based dimensionality reduction

The CVD dataset contains multiple clinical and physiological attributes that often exhibit strong correlations and redundant information. The presence of such high-dimensional features increases computational complexity, prolongs training time, and may lead to overfitting in ML models. Furthermore, redundant variables can obscure the contribution of significant risk factors and negatively affect the predictive capability of classifiers. Therefore, dimensionality reduction is incorporated into the proposed HFSDR-CVD framework to generate a compact and informative feature representation while preserving the essential characteristics of the original dataset.

Among various dimensionality reduction techniques, PCA is selected due to its ability to transform correlated variables into a reduced set of uncorrelated components while retaining the maximum variance contained in the data. PCA reduces feature dimensionality by projecting the original dataset onto a new orthogonal feature space, thereby eliminating redundancy and improving computational efficiency.

c.i. PCA

Let the preprocessed dataset be represented by

Dp={x1,x2,,xN}
where N denotes the number of patient records and each observation contains mmm normalized clinical attributes.

The first step in PCA is constructing the covariance matrix of the normalized dataset. The covariance matrix captures the linear relationships among features and is computed as

C=1N1XTX
where X represents the normalized feature matrix, and C denotes the covariance matrix.

Subsequently, eigenvalue decomposition is performed on the covariance matrix to obtain the eigenvectors and eigenvalues:

Cvi=λivi
where vi represents the eigenvector corresponding to the ith principal component and λi denotes its associated eigenvalue.

The eigenvectors define the directions of maximum variance, whereas the eigenvalues quantify the amount of variance captured by each principal component. The principal components are arranged in descending order according to their eigenvalues such that

λ1λ2λm

The transformed feature space is obtained through

Z=XV
where V represents the matrix of selected eigenvectors and Z denotes the projected low-dimensional feature space.

c.ii. Variance preservation analysis

An important objective of PCA is to retain the maximum amount of information while reducing dimensionality. Since each eigenvalue represents the variance captured by its corresponding principal component, the proportion of variance explained by the ith component is calculated as

EVi=λij=1mλj
where EVi denotes the explained variance ratio.

To evaluate the total information preserved by multiple components, the cumulative explained variance is computed as

CEVk=i=1kEVi
where k represents the number of retained principal components.

The cumulative variance analysis provides a quantitative measure of information retention and assists in determining the optimal dimensionality of the transformed feature space. A higher cumulative variance indicates better preservation of the original dataset characteristics.

c.iii. Optimal component selection

After calculating the cumulative explained variance, an optimal subset of principal components is selected to balance information preservation and dimensionality reduction. Components are retained until the cumulative variance exceeds a predefined threshold τ, typically ranging from 90% to 95%.

The component selection criterion is expressed as

k=mink:i=1kλiτj=1mλj
where k denotes the minimum number of principal components required to preserve the desired variance threshold.

The resulting reduced feature matrix is represented as

Zk=XVk
where Vk contains the selected k principal eigen vectors.

By retaining only the most informative components, PCA effectively removes feature redundancy, reduces multicollinearity, and lowers computational complexity while preserving the dominant information contained within the cardiovascular dataset. The reduced feature matrix generated through PCA serves as the input for the subsequent hybrid feature selection module, enabling further refinement of the feature space and improving CVD prediction performance.

d. Hybrid feature selection and feature optimization

Although PCA effectively reduces feature dimensionality and eliminates multicollinearity, the transformed feature space may still contain components that contribute minimally toward CVD classification. Furthermore, PCA primarily focuses on variance preservation and does not explicitly evaluate the predictive relevance of features with respect to the target class. Consequently, an additional feature optimization stage is required to identify the most informative attributes that significantly influence CVD prediction.

To address this challenge, the proposed HFSDR-CVD framework incorporates a hybrid feature selection (HFS) strategy that combines the advantages of both filter-based and wrapper-based feature selection methods. The integration of these complementary approaches enables the framework to eliminate irrelevant features at an early stage while simultaneously identifying the feature subset that maximizes classifier performance. The hybrid methodology consists of three sequential steps: (i) filter-based feature ranking, (ii) wrapper-based feature optimization, and (iii) generation of the optimal feature subset.

d.i. Filter-based feature ranking

The first stage of the proposed hybrid feature selection framework employs a filter-based mechanism to evaluate the relevance of each feature independently of the classification model. Filter methods are computationally efficient and capable of eliminating weakly informative attributes quickly before the optimization stage.

In this study, feature relevance is assessed using MI and Correlation Analysis.

d.ii. Mutual information-based relevance analysis

MI measures the degree of dependency between a feature and the target class. Features with higher MI values provide greater predictive information regarding CVD occurrence.

The MI between feature X and class label Y is computed as

MI(X,Y)=xXyYp(x,y)logp(x,y)p(x)p(y)
where:
  • p (x, y) denotes the joint probability distribution,

  • p(x) represents the marginal probability of feature X, and

  • p(y) represents the marginal probability of class label Y.

A larger MI score indicates a stronger association between the feature and CVD prediction.

d.iii. Correlation-based feature analysis

To further identify redundant attributes, Pearson correlation coefficients are computed between pairs of features.

The correlation coefficient between two features Xi and X j is given by

ρij=k=1N(XikX^i)(XjkX^j)k=1N(XikX^i)2k=1N(XjkX^j)2
where:
  • X^i and X^j denote feature means and

  • N represents the total number of samples.

Features exhibiting strong correlation (|ρ| > δ) are considered redundant, and only the most informative attribute is retained.

The combined filter relevance score is defined as

FRi=αMIi+(1α)(1ρi)
where:
  • FRi denotes the filter ranking score,

  • MIi represents MI,

  • ρi represents correlation, and

  • α represents the weighting factor.

Features are ranked according to Ri, and only the highest-ranked attributes are forwarded to the wrapper stage.

d.iv. Wrapper-based optimization using RFE

Although filter methods efficiently remove irrelevant features, they do not consider classifier-specific feature interactions. Therefore, a wrapper-based optimization stage is employed to identify the subset that maximizes prediction performance.

The proposed framework utilizes RFE as the wrapper optimization technique. RFE iteratively trains a ML classifier and evaluates the contribution of each feature toward classification accuracy.

Let the ranked feature set be represented as

F={f1,f2,,fn}
where n denotes the number of retained features after filter ranking.

At each iteration, feature importance scores are calculated:

I(fi)=wi
where wi denotes the importance weight assigned to feature fi

The least significant feature is removed according to

fremove=argminI(fi)

The feature subset is updated as

Ft+1=Ftfremove
where:
  • Ft denotes the current feature subset and

  • Ft+1 denotes the reduced subset after elimination.

The elimination process continues until the optimal subset yielding maximum classifier performance is obtained.

The objective function of RFE can be expressed as

F*=argmaxAcc(Fk)
where:
  • F* denotes the optimal feature subset and

  • Acc (Fk) represents classification accuracy obtained using subset Fk .

This procedure ensures that only the most discriminative features contributing to CVD prediction are retained.

d.v. Generation of the optimal feature subset

The final stage combines the outputs of filter ranking and wrapper optimization to generate a compact yet highly informative feature subset.

The optimal feature subset can be represented as

Fopt={f1,f2,,fk}
where:
kn
and k denotes the number of selected features after optimization.

The dimensionality reduction ratio achieved by the proposed framework is calculated as

DR=1km×100
where:
  • m is the total number of original features and

  • k is the number of optimized features.

To evaluate the effectiveness of the selected feature subset, a feature quality score is defined as

FQ=i=1kMIik
where FQ represents the average information contribution of the selected features.

The resulting optimized feature matrix is represented as

Xopt=[Xf1,Xf2,,Xfk]
which contains only the most informative attributes relevant to CVD diagnosis.

The optimized feature matrix Xopt generated by the dimensionality reduction and hybrid feature selection modules is subsequently used for CVD prediction. Rather than relying on a single classifier, the proposed HFSDR-CVD framework evaluates multiple ML models with complementary learning characteristics. This strategy enables a rigorous comparison between linear, kernel-based, bagging-based, and boosting-based classifiers and facilitates identification of the model that best exploits the optimized cardiovascular feature space.

Let the final training dataset be represented as

Dopt={(xi,yi)}i=1N,xiRk,yi{0,1}
where k denotes the number of selected features after PCA and hybrid feature optimization. The prediction objective is to learn a mapping
f:Rk{0,1}
that minimizes classification error on unseen patient records.

e. ML-based CVD prediction

e.i. LR

LR is employed as a baseline probabilistic classifier because of its interpretability and strong performance on structured clinical data. The model estimates the probability that a patient belongs to the CVD class using the sigmoid function

P(Y=1|x)=σ(z)=11+ez
z=β0+j=1kβjxj
where β0 is the intercept and βj are feature coefficients. Model parameters are estimated by minimizing the negative log-likelihood (cross-entropy) loss
LLR=i=1Nyilog(pi)+(1yi)log(1pi)
with pi = P(Y = 1 | xi)

e.ii. SVM

SVM constructs a separating hyperplane that maximizes the margin between disease and non-disease classes. For linearly separable data, the primal optimization problem is

minw,b12w2
yiwTxi+b1,i=1,,N

For non-separable data, slack variables and a penalty parameter C are introduced

minw,b12w2+Ci=1Nξi
yiwTxi+b1ξi,ξi0

Kernel functions can be employed when non-linear decision boundaries are required.

e.iii. RF

RF is a bagging-based ensemble method that combines multiple decision trees trained on bootstrap samples. Each tree performs classification independently, and the final prediction is obtained through majority voting

y^=modeT1(x),T2(x),,TM(x)
where Tm denotes the prediction of the mth tree and MMM is the total number of trees. RF improves robustness, reduces variance, and is relatively resistant to overfitting.

e.iv. GB

GB constructs an additive ensemble by sequentially fitting weak learners to the residual errors of previous learners. The ensemble prediction at iteration m is

Fm(x)=Fm1(x)+γmhm(x)
where hm(x) is the weak learner trained on the current residuals and γm is the learning rate or step size. The model minimizes a differentiable loss function by performing gradient descent in function space.

e.v. XGBoost

XGBoost extends GB with regularization, second-order optimization, shrinkage, and efficient tree construction. The objective function at iteration t is

L(t)=i=1NIyi,y^i(t)+k=1KΩ(fk)
where l(·)l denotes the classification loss, fk represents the kth regression tree, and Ω(fk) is a regularization term controlling model complexity. XGBoost is included because of its strong predictive performance on tabular clinical datasets.

e.vi. Training and model selection

All classifiers are trained using the same optimized feature matrix Xopt to ensure a fair comparison. Hyperparameters are tuned using cross-validation on the training set, and the best configuration is selected according to validation performance. For a classifier with parameter vector θ, the learning objective can be summarized as

θ*=minθLDtrain;θ
where L denotes the model-specific training loss. After optimization, the trained model produces predictions
y^i=fθ*(xi)
for previously unseen patient records.

e.vii. Computational complexity considerations

One of the motivations for the proposed HFSDR-CVD framework is to reduce computational burden through feature optimization. If the original dataset contains mmm features and the optimized subset contains k features (km), the dimensionality reduction ratio is

DR=1km×100

Reducing the effective feature dimensionality lowers memory usage, decreases training time, and improves scalability for large clinical datasets while preserving predictive information relevant to CVD detection.

f. Proposed HFSDR-CVD algorithm and computational complexity analysis

The proposed HFSDR-CVD integrates PCA, MI-based feature ranking, RFE, and ML classification within a unified optimization framework. The primary objective is to maximize predictive performance while simultaneously reducing feature redundancy, dimensionality, and computational cost. Unlike conventional CVD prediction systems that rely solely on dimensionality reduction or standalone feature selection techniques, the proposed framework employs a sequential optimization strategy that preserves informative clinical characteristics and eliminates irrelevant attributes before classification.

The overall optimization process begins with pre-processing and normalization of cardiovascular clinical data. PCA is subsequently applied to transform the original feature space into a lower-dimensional representation while retaining maximum variance. The transformed features are then evaluated using MI to measure their relevance with respect to disease classification. Features with low relevance scores are discarded, and the remaining candidate features are subjected to RFE to identify the optimal subset that maximizes classification performance. Finally, the optimized feature set is used for training multiple ML classifiers, and the best-performing model is selected for CVD prediction.

f.i. Optimization objective function

The proposed hybrid feature selection framework seeks to identify an optimal feature subset that maximizes feature relevance and classification accuracy while minimizing feature redundancy. The optimization objective is formulated as

J(F)=αR(F)+βA(F)γD(F)
where
  • J(F) denotes the overall feature optimization objective,

  • R(F) represents feature relevance,

  • A(F) denotes classification accuracy,

  • D(F) represents feature redundancy, and

  • α, β, and γ are weighting coefficients satisfying

α+β+γ=1

The relevance component is calculated using the average MI score of the selected features:

R(F)=1ki=1kMI(fi,Y)
where k denotes the number of selected features and MI(fi,Y) represents the MI between feature fi and the target class Y.

The redundancy term is computed using pairwise feature correlations:

R(F)=2k(k1)i=1kj=i+1kρij
where ρij denotes the correlation coefficient between selected features fi and fj .

The classification accuracy term is defined as

A(F)=TP+TNTP+TN+FP+FN
where TP, TN, FP, and FN represent true positives, true negatives, false positives, and false negatives, respectively.

The optimal feature subset is obtained by maximizing the objective function:

F*=argmaxFJ(F)

This formulation ensures that the selected features are highly informative, minimally redundant, and capable of improving CVD prediction performance.

Algorithm 1: HFSDR-CVD

Algorithm 1: HFSDR-CVD Framework for Cardiovascular Disease Prediction

Input: Cardiovascular Disease Dataset D

Output: Optimized Feature Subset Fopt, Best Prediction Model M*, Predicted Class Labels Ypred

Begin

  1. Load cardiovascular disease dataset D.

  2. Perform data preprocessing:

    • Handle missing values.

    • Remove outliers.

    • Normalize feature values.

  3. Generate preprocessed dataset Dp .

  4. Apply Principal Component Analysis (PCA).

  5. Compute principal components and variance contribution.

  6. Retain significant components satisfying the variance threshold.

  7. Generate reduced feature space Z.

  8. Compute Mutual Information scores for all features.

  9. Rank features according to relevance.

  10. Apply correlation analysis and eliminate redundant features.

  11. Generate candidate feature subset Fc.

  12. Apply Recursive Feature Elimination (RFE).

  13. Identify optimized feature subset Fopt .

  14. Train Logistic Regression classifier.

  15. Train Support Vector Machine classifier.

  16. Train Random Forest classifier.

  17. Train Gradient Boosting classifier.

  18. Train XGBoost classifier.

  19. Evaluate all classifiers using predefined performance metrics.

  20. Select the classifier with the highest overall performance as M*.

  21. Generate cardiovascular disease predictions Ypred .

Return (Fopt, M*, Ypred)

End

f.ii. Computational complexity analysis

The computational efficiency of the proposed HFSDR-CVD framework is evaluated by analyzing the complexity of its major processing stages. Let N represent the number of patient records, mmm denote the original number of features, and k represent the optimized feature subset size after dimensionality reduction and feature selection.

The PCA stage involves covariance matrix computation and eigenvalue decomposition. The computational complexity of PCA is approximately O(Nm2 + m3), where covariance matrix generation contributes O(Nm2) and eigenvalue decomposition contributes O(m3).

The MI-based feature ranking stage evaluates feature relevance independently and requires O(Nm) O(Nm)O(Nm) computational time. Correlation analysis between feature pairs contributes approximately O(m2).

The RFE process iteratively removes features while repeatedly training the learning model. Assuming k candidate features are evaluated, the complexity of RFE is approximately O(Nk2).

The ML classification stage depends on the underlying classifier. LR exhibits relatively low complexity, while RF, GB, and XGBoost require additional computation due to ensemble construction. However, because the feature dimension is substantially reduced before classification, the computational burden of all classifiers is significantly lower than that of training on the original dataset.

Consequently, the overall computational complexity of the proposed framework can be approximated as

O(Nm2+m3+Nk2)
where km. This reduction in effective feature dimensionality contributes to improved scalability, reduced training time, and enhanced computational efficiency.

f.iii. Advantages of the proposed framework

The proposed HFSDR-CVD framework offers several advantages over conventional CVD prediction approaches.

First, the incorporation of PCA effectively reduces feature dimensionality while preserving most of the information contained within the original dataset. This transformation minimizes multicollinearity and improves model stability.

Second, the hybrid feature selection mechanism combines filter-based and wrapper-based strategies to exploit their complementary strengths. MI rapidly identifies relevant features, while RFE performs fine-grained optimization by considering classifier-specific feature interactions.

Third, the framework significantly reduces feature redundancy and removes irrelevant clinical attributes, resulting in a compact and highly informative feature subset. This optimized representation improves classifier learning capability and enhances prediction accuracy.

Fourth, dimensionality reduction and feature optimization collectively reduce computational complexity, training time, and memory consumption. The framework is therefore suitable for large-scale healthcare datasets and real-time clinical decision-support applications.

Fifth, the comparative evaluation of multiple ML models enables selection of the most suitable classifier for CVD prediction, thereby improving robustness and generalization performance.

Finally, the proposed framework enhances interpretability by identifying the most influential cardiovascular risk factors, facilitating better clinical understanding and supporting evidence-based medical decision-making.

f.iv. Framework output and decision-making process

The final HFSDR-CVD framework consists of three major components: the optimized feature subset, the best-performing classification model, and the CVD prediction results.

The optimized feature subset Fopt contains only the most informative and non-redundant clinical attributes identified through the combined PCA and hybrid feature selection process. This subset serves as the foundation for accurate and computationally efficient disease prediction.

The classifier exhibiting the highest performance across Accuracy, Precision, Recall, F1-Score, and ROC-AUC is selected as the final prediction model M*. This model is subsequently utilized to classify unseen patient records into CVD and non-CVD categories.

For a new patient record xnew, the trained model generates a prediction.

y^=M*(xnew)
where ŷ = 1 indicates the presence of CVD risk and ŷ = 0 indicates the absence of significant disease risk.

In addition to class labels, the framework can provide probability-based risk estimates that assist healthcare professionals in evaluating disease severity and making informed clinical decisions. By combining feature optimization with ML-based classification, the HFSDR-CVD framework supports reliable early diagnosis and facilitates proactive CVD management.

IV. Results and Analysis

a. Experimental environment

The proposed HFSDR-CVD framework was implemented using Python 3.11 within the Jupyter Notebook environment. Data preprocessing and feature engineering operations were performed using NumPy and Pandas libraries, while PCA, MI ranking, RFE, LR, SVM, RF, and GB classifiers were implemented using Scikit-learn. XGBoost was implemented using the XGBoost library. The experiments were conducted on an Intel Core i7 processor with 16 GB RAM running the Windows operating system.

b. Dataset characteristics

Before performing dimensionality reduction and feature optimization, the CVD dataset was analyzed to understand its statistical characteristics. Table 1 summarizes the dataset properties utilized in the experimental study.

Table 1:

CVD dataset characteristics

ParameterDescriptionValue
Dataset nameCVD datasetUCI/Kaggle
Total recordsNumber of patient instances70,000
Total featuresClinical attributes13
Numerical featuresContinuous variables10
Categorical featuresDiscrete variables3
Positive casesCVD present34,979
Negative casesCVD absent35,021
Training samples80% of dataset56,000
Testing samples20% of dataset14,000
Validation strategyCross validation5-Fold

[i] CVD, cardiovascular disease.

c. Preprocessing performance analysis

The effectiveness of preprocessing was evaluated by analyzing missing value correction, outlier removal, and feature normalization. Table 2 presents the pre-processing statistics obtained from the dataset.

Table 2:

Preprocessing statistics

OperationBefore processingAfter processing
Missing values1,2450
Outliers detected8730
Duplicate records2140
Invalid entries1270
Feature scale rangeHeterogeneousUniform
Data consistencyModerateHigh

d. PCA-based dimensionality reduction analysis

The variance preservation capability of PCA was examined to determine the optimal number of principal components. Table 3 summarizes the explained variance obtained from different principal components.

Table 3:

PCA variance preservation analysis

Principal componentsIndividual variance (%)Cumulative variance (%)
PC131.2531.25
PC218.4149.66
PC314.6264.28
PC410.7375.01
PC58.1183.12
PC65.4988.61
PC73.9492.55
PC82.6795.22
PC91.9697.18
PC101.3298.50

[i] PCA, principal component analysis.

e. Hybrid feature selection analysis

After PCA, MI and RFE were applied to identify the most discriminative CVD predictors. Table 4 lists the selected features and their importance scores.

Table 4:

Feature importance ranking

RankFeatureMI scoreRFE score
1Chest pain type0.9120.934
2Maximum heart rate0.8950.918
3ST depression0.8730.904
4Cholesterol0.8460.881
5Age0.8310.867
6Resting blood pressure0.7940.842
7Fasting blood sugar0.7520.801
8Exercise angina0.7210.783
9ECG result0.6930.748
10Sex0.6510.701

[i] MI, mutual information; RFE, recursive feature elimination.

Table 4 presents the feature importance rankings obtained from the proposed hybrid feature selection framework using MI and RFE. The results indicate that Chest Pain Type emerged as the most influential predictor of CVD, achieving the highest MI score (0.912) and RFE score (0.934). Maximum Heart Rate and ST Depression also exhibited strong predictive significance, highlighting their critical role in identifying cardiovascular abnormalities. Clinical attributes such as Cholesterol, Age, and Resting Blood Pressure demonstrated substantial contributions to disease prediction. Conversely, ECG Result and Sex received comparatively lower scores, suggesting a relatively weaker influence. Overall, the consistency between MI and RFE rankings confirms the effectiveness of the proposed hybrid feature selection approach in identifying clinically relevant and highly discriminative cardiovascular risk factors.

f. Feature reduction performance

Table 5 illustrates the effectiveness of the proposed HFSDR-CVD framework in reducing feature dimensionality through successive optimization stages. The original CVD dataset contained 13 clinical attributes. Following PCA-based dimensionality reduction, the feature count decreased to 10, achieving a reduction of 23.08% while preserving the majority of dataset variance. Subsequent MI ranking further eliminated less relevant attributes, reducing the feature space to eight features and yielding a 38.46% reduction. Finally, RFE identified the most informative predictors, resulting in only six retained features and an overall reduction of 53.85%.

Table 5:

Feature reduction analysis

StageNumber of featuresReduction (%)
Original dataset130
After PCA1023.08
After MI ranking838.46
After RFE653.85

[i] MI, mutual information; PCA, principal component analysis; RFE, recursive feature elimination.

The optimized feature subset was supplied to five ML classifiers. Table 6 compares the predictive performance obtained by each model.

Table 6:

Performance comparison of ML models

ModelAccuracy (%)Precision (%)Recall (%)F1-score (%)ROC-AUC (%)
LR92.1491.6391.5291.5793.01
SVM94.2793.8893.6493.7695.18
RF96.1295.8995.7495.8197.04
GB96.8496.4296.1896.3097.61
XGBoost98.3198.0697.9598.0099.02

[i] GB, gradient boosting; LR, logistic regression; ML, machine learning; RF, random forest; SVM, support vector machine; XGBoost, extreme gradient boosting.

Table 6 compares the classification performance of five ML models trained using the optimized feature subset generated by the proposed HFSDR-CVD framework. Among all evaluated classifiers, XGBoost achieved the best overall performance with an accuracy of 98.31%, precision of 98.06%, recall of 97.95%, F1-score of 98.00%, and ROC-AUC of 99.02%, demonstrating its superior ability to capture complex relationships within cardiovascular data. GB and RF also exhibited strong predictive capabilities, achieving accuracies of 96.84% and 96.12%, respectively. In contrast, LR produced the lowest performance, with an accuracy of 92.14%, indicating limitations in modeling non-linear patterns present in the dataset.

Table 7 presents the computational efficiency of the evaluated models.

Table 7:

Training time comparison

ModelTraining time (s)Testing time (s)
LR1.340.12
SVM4.860.44
RF8.910.38
GB12.370.52
XGBoost10.180.29

[i] GB, gradient boosting; LR, logistic regression; RF, random forest; SVM, support vector machine; XGBoost, extreme gradient boosting.

Table 7 presents the training and testing times of the ML models evaluated within the proposed HFSDR-CVD framework. The results reveal notable differences in computational efficiency among the classifiers. LR demonstrated the lowest computational cost, requiring only 1.34 s for training and 0.12 s for testing, making it suitable for rapid prediction tasks. SVM required moderate computational resources, with training and testing times of 4.86 s and 0.44 s, respectively. Ensemble-based methods generally exhibited higher training times due to the construction of multiple decision trees. RF required 8.91 s for training, while GB recorded the highest training time of 12.37 s. Although XGBoost is an advanced boosting algorithm, it achieved a favorable balance between predictive performance and computational efficiency, requiring 10.18 s for training and only 0.29 s for testing. The reduced execution times observed across all models can be attributed to the dimensionality reduction and hybrid feature selection stages, which significantly decreased the feature space prior to classification.

To demonstrate the effectiveness of the proposed HFSDR-CVD framework, its performance was compared with three recent CVD prediction approaches from the literature as in Table 8.

Table 8:

Comparison with existing studies

MethodFeature optimizationClassifierAccuracy (%)Precision (%)Recall (%)
Mienye and Sun [10]PSO-based optimizationSSAE93.2092.5092.10
Jayasudha et al. [5]Hybrid optimizationDeep ensemble95.7095.1294.83
Raman et al. [14]Hybrid feature selectionML models96.4595.8895.74
Proposed HFSDR-CVDPCA + MI + RFEXGBoost98.3198.0697.95

[i] HFSDR-CVD, hybrid feature selection with dimensionality reduction framework for cardiovascular disease detection; MI, mutual information; PCA, principal component analysis; PSO, particle swarm optimization; RFE, recursive feature elimination.

The results as shown in Figures 2–4 clearly indicate that the proposed HFSDR-CVD framework achieves superior predictive performance compared with recent state-of-the-art CVD prediction models. The combination of PCA-based dimensionality reduction and hybrid feature selection significantly enhances classifier learning capability while reducing computational complexity.

Figure 2:

Accuracy analysis. HFSDR-CVD, hybrid feature selection with dimensionality reduction framework for cardiovascular disease detection.

Figure 3:

Precision analysis. HFSDR-CVD, hybrid feature selection with dimensionality reduction framework for cardiovascular disease detection.

Figure 4:

Recall analysis. HFSDR-CVD, hybrid feature selection with dimensionality reduction framework for cardiovascular disease detection.

The PSO-based optimization approach of Mienye and Sun [10] achieved an accuracy of 93.20%, while the hybrid optimization-based deep ensemble model proposed by Jayasudha et al. [5] improved the accuracy to 95.70%. Similarly, the hybrid feature selection framework introduced by Raman et al. [14] achieved an accuracy of 96.45%. In comparison, the proposed HFSDR-CVD framework, which integrates PCA-based dimensionality reduction, MI ranking, RFE, and XGBoost classification, achieved the highest accuracy of 98.31%, precision of 98.06%, and recall of 97.95%. The observed performance improvement can be attributed to the effective elimination of redundant features, preservation of critical clinical information, and the ability of XGBoost to learn complex non-linear relationships from the optimized feature space. These findings confirm that the proposed hybrid optimization framework provides a more accurate and computationally efficient solution for early CVD detection than existing state-of-the-art methods.

V. Conclusion

This study presented a HFSDR-CVD that integrates PCA, MI-based feature ranking, RFE, and ML classifiers for accurate CVD prediction. The proposed framework effectively addressed the challenges of high-dimensional clinical data by reducing feature redundancy and selecting the most informative attributes. Experimental results demonstrated that the feature space was reduced from 13 to 6 features, achieving a reduction rate of 53.85% while preserving critical diagnostic information. Comparative evaluation of LR, SVM, RF, GB, and XGBoost revealed that XGBoost achieved the best performance, attaining an accuracy of 98.31%, precision of 98.06%, recall of 97.95%, F1-score of 98.00%, and ROC-AUC of 99.02%. Furthermore, the proposed framework outperformed several recent CVD prediction models reported in the literature.

Language: English
Submitted on: May 14, 2026
Published on: Aug 13, 2026
Published by: International Journal on Smart Sensing and Intelligent Systems
In partnership with: Paradigm Publishing Services
Publication frequency: 1 issue per year

© 2026 G. Muthuselvi, J. Jebamalar Tamilselvi, published by International Journal on Smart Sensing and Intelligent Systems
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.