Next Article in Journal
Evaluating an Immersive Virtual Classroom as an Augmented Reality Platform in Synchronous Remote Learning
Next Article in Special Issue
Computer Vision Tasks for Ambient Intelligence in Children’s Health
Previous Article in Journal
Healthcare Trust Evolution with Explainable Artificial Intelligence: Bibliometric Analysis
Previous Article in Special Issue
Exploring the State of Machine Learning and Deep Learning in Medicine: A Survey of the Italian Research Community
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cost-Sensitive Models to Predict Risk of Cardiovascular Events in Patients with Chronic Heart Failure

by
Maria Carmela Groccia
1,†,
Rosita Guido
1,†,
Domenico Conforti
1,†,
Corrado Pelaia
2,
Giuseppe Armentaro
2,
Alfredo Francesco Toscani
2,
Sofia Miceli
2,
Elena Succurro
2,
Marta Letizia Hribal
2 and
Angela Sciacqua
2,*,†
1
Department of Mechanical, Energy and Management Engineering, University of Calabria, Ponte Pietro Bucci 41C, 87036 Arcavacata di Rende, Italy
2
Department of Medical and Surgical Sciences, University Magna Graecia, 88100 Catanzaro, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 19 August 2023 / Revised: 30 September 2023 / Accepted: 1 October 2023 / Published: 3 October 2023
(This article belongs to the Special Issue Computer Vision, Pattern Recognition and Machine Learning in Italy)

Abstract

:
Chronic heart failure (CHF) is a clinical syndrome characterised by symptoms and signs due to structural and/or functional abnormalities of the heart. CHF confers risk for cardiovascular deterioration events which cause recurrent hospitalisations and high mortality rates. The early prediction of these events is very important to limit serious consequences, improve the quality of care, and reduce its burden. CHF is a progressive condition in which patients may remain asymptomatic before the onset of symptoms, as observed in heart failure with a preserved ejection fraction. The early detection of underlying causes is critical for treatment optimisation and prognosis improvement. To develop models to predict cardiovascular deterioration events in patients with chronic heart failure, a real dataset was constructed and a knowledge discovery task was implemented in this study. The dataset is imbalanced, as it is common in real-world applications. It thus posed a challenge because imbalanced datasets tend to be overwhelmed by the abundance of majority-class instances during the learning process. To address the issue, a pipeline was developed specifically to handle imbalanced data. Different predictive models were developed and compared. To enhance sensitivity and other performance metrics, we employed multiple approaches, including data resampling, cost-sensitive methods, and a hybrid method that combines both techniques. These methods were utilised to assess the predictive capabilities of the models and their effectiveness in handling imbalanced data. By using these metrics, we aimed to identify the most effective strategies for achieving improved model performance in real scenarios with imbalanced datasets. The best model for predicting cardiovascular events achieved mean a sensitivity 65%, a mean specificity 55%, and a mean area under the curve of 0.71. The results show that cost-sensitive models combined with over/under sampling approaches are effective for the meaningful prediction of cardiovascular events in CHF patients.

1. Introduction

In the world, more than 64 million people suffer from heart failure (HF) [1]. The term chronic heart failure (CHF) refers to a clinical syndrome characterised by symptoms and signs due to structural and/or functional abnormalities of the heart, resulting in a reduction in cardiac output and/or an increase in intracardiac pressure at rest or under stress. CHF is defined as a condition characterised by asthenia and dyspnoea, as well as signs of pulmonary and/or systemic venous congestion. Although these symptoms manifest over a medium-long period, patients seem to be asymptomatic during that time, while the causes that lead to HF often exist for a long time before symptoms appear. The early identification of the organic causes that drive HF is very important to ameliorate the prognosis, but also to allow the precocious detection of decompensation, in case of an already developed disease [2,3]. CHF can be classified according to different criteria. The updated European Society of Cardiology Guidelines recommend considering the classification based on the measurement of the left ventricle ejection fraction (LVEF). In addition, the New York Heart Association (NYHA) classification is used to define the gravity of symptoms. The higher the NYHA class, the worse the clinical condition, and the poorer the overall outcome. However, it is possible that patients with mild symptoms may have an increased risk of mortality and/or hospitalisation [4,5].
According to the guidelines of the European Society of Cardiology, we recognise three phenotypes of HF, one with reduced ejection fraction (HFrEF) with an LVEF ≤ 40%; one HF with mildly reduced ejection fraction (HFmrEF) with an LVEF between 41 and 49%; and one HF with preserved ejection fraction (HFpEF) with LVEF ≥ 50%. This characterisation is not only echocardiographic but also aetiological and clinical. In fact, HFrEF is predominantly due to ischaemic heart disease, affects mostly young male patients and responds well to pharmacological treatment with inhibitors of the renin–angiotensin–aldosterone system, beta-blockers and glyflozines. In contrast, the HFmrEF and HFpEF forms are prevalent in the elderly and in women, and in this phenotype, an important influence on the aetiology and symptomatology is determined by comorbidities, especially chronic kidney disease, anaemia, arterial hypertension, atrial fibrillation and chronic obstructive pulmonary disease. In these two phenotypes, the correct treatment of comorbidities and treatment with diuretics and glyflozines is crucial [6]. In addition to the NYHA functional class, the INTERMACS (Interagency Registry for Mechanically Assisted Circulatory Support) classification plays an important role when it comes to terminal forms of HF, which is used to assess patients who could benefit from mechanical supportive circulation [7].
HF is widespread in Europe. Indeed, about 15 million people suffer from this disease, whose estimated incidence is 2–3% per year in a population older than 40 years of age. Current epidemiological projections forecast that such numbers will further increase in parallel with the progressive increment in life expectancy [8]. In particular, there is a CHF prevalence of 6.4% in subjects older than 65 years. These considerations highlight the importance of a better and earlier disease identification, hopefully coupled with greater public awareness of the relevance of CHF detection, associated with a more rational use of health resources [9]. Therefore, it is noteworthy to consider the problem of recurrent hospitalisations for CHF patients. An Italian study showed that the one-year hospitalisation rate of CHF patients was about 22%. The outcome after one year demonstrated, in the same study, that the mortality rate was greater in NYHA III–IV patients than in NYHA I–II patients (14.5% and 4%, respectively) [10]. Epidemiological data based on hospital discharge showed that CHF was the second cause of hospitalisation, with an incidence rate of 4–5 cases for every one-thousand inhabitants. Moreover, one patient out of four had a hospital readmission after one month. Fifty percent of patients were re-hospitalised during the six months following the first hospital admission, and their prognosis worsened. According to these reports, health costs were around 2% of the total healthcare expenditures, and hospitalisation could account for 60–70% of the global CHF cure burden [11,12,13]. Obviously, frailty and comorbidities may have a negative impact on health costs and prognosis [14].
The importance of early diagnosis is well recognised. Predictive models based on common clinical variables are potentially useful for the early detection of decompensation risk. In this regard, several studies have been published, based on the evaluation of tele-monitored patients undergoing investigation related to common clinical parameters, such as weight, systolic blood pressure (SBP), diastolic blood pressure (DBP), heart rate (HR), and arterial oxygen saturation [15,16]. A relevant advantage of these systems relies on the implementation of easy home checking, affordable by patients themselves or their caregivers. This approach can be associated with a better compliance of most patients, who prefer non-invasive procedures which do not require going to hospital.
Machine learning (ML) is an interesting area of research within heart failure. ML techniques have been applied to many aspects of heart failure such as diagnosis, classification, and prediction [17]. Some limitations of ML in this field were discussed in [18]. Knowledge discovery techniques can be properly taken into account due to their effective impact in the cardiovascular domain for prediction tasks [19,20]. ML and KD techniques can significantly contribute to disease identification and make a real-time effective clinical decision. Many ML methodologies have already been investigated in predicting the presence of adverse events in CHF patients, such as destabilisation, re-hospitalisation, and mortality [21,22].
It is important to point out that in many applications of ML, such as medical diagnosis, datasets are often imbalanced. Typically, the number of patients is far less than that of healthy individuals. In order to solve this problem, several methods are proposed in the literature. In order to generate a balanced dataset, over-sampling methods add more data to the smaller class, making it the same size as the larger class [23,24]; under-sampling methods sample the larger class in order to have the same size as the smaller class [25]. Cost-sensitive learning approaches take the costs of misclassification errors into account [26,27,28,29,30]. They assign a higher misclassification cost for objects belonging to the minority class with respect to the misclassification cost for objects belonging to the majority class. For instance, the approach presented in [28] employs support vector machines (SVMs) as a classification method and assigns penalty coefficients to both positive and negative instances. Combinations of sampling techniques with cost-sensitive learning methods were shown to be effective in addressing the class imbalance problem [31].
In this paper, we designed a knowledge discovery (KD) task and implemented it with a two-fold purpose. Firstly, the predictive capability of clinical variable (CV) events are investigated analysing real collected data spanning five years from the ambulatory of the Geriatrics Division at the “Mater Domini” University Hospital in Catanzaro, Italy. Secondly, different ML models for predicting cardiovascular deterioration events in CHF patients were developed. The analysis focused on clinical decompensation events and major CV events were selected. The KDD analysis was defined as a predictive task stated as a supervised binary classification problem [32]. The real data were analysed and a dataset was constructed from it. The dataset exhibited an imbalanced distribution due to the under-representation of event cases. Dealing with this imbalanced data was one of the most challenging aspects of applied ML. To address the issue, a pipeline was specifically developed to handle imbalanced data. Subsequently, various ML models were trained, and their hyper-parameters were optimised using a grid search approach. The performance of these learned models was then evaluated using appropriate metrics. To further tackle the class imbalance problem, three distinct approaches were implemented and tested: cost-sensitive learning methods, data resampling methods, and a combination of cost-sensitive learning methods with data resampling methods. The results demonstrated that combining sampling methods with cost-sensitive learning models yielded promising values for sensitivity and balanced accuracy. Moreover, several computational experiments were carried out to optimise the hyper-parameters of the ML models to improve the performance on the real-world and imbalanced dataset. The ML approach adopted in this study can be broken down into four main steps: (1) Data preprocessing—this step involved operations such as data cleaning, handling missing data, data transformation, and reducing data imbalance; (2) Features selection—this step aimed to reduce overfitting, training times, and improve accuracy by selecting the most relevant subset of variables to build the predictive model; (3) Model building—in this stage, parameter and hyper parameter values for the ML model are chosen to optimise its performance; (4) Cross-validation approach—the dataset was divided into two separate groups—a training set and a test set—for validating the ML model. The model is trained on the training set and then tested on the test set. By following these steps, the study successfully developed predictive models for identifying cardiovascular deterioration events in CHF patients based on the real-world data collected from the Geriatrics Division at the “Mater Domini” University Hospital.
Figure 1 illustrates the knowledge discovery process that we designed and implemented. It is detailed in the following sections.
The remainder of this paper is structured as follows. Section 2 offers a comprehensive description of our data and outlines the data processing method employed. Section 3 presents an overview of the ML models selected for this research paper (i.e., support vector machine, artificial neural network, naïve Bayes, decision tree, and random forest), the three methods that we adopt to balance the class distribution of the unbalanced dataset, and the parameter tuning approach for the ML models. Moving on, Section 4 presents the experimental results and discusses the model performance metrics. Lastly, Section 5 provides the concluding remarks for this paper.

2. Real Data Collection and Dataset Construction

The data were collected over five years, as part of a pilot study conducted at the CHF outpatient clinic of the Geriatrics Division at the “Mater Domini” University Hospital in Catanzaro, Italy. The key steps undertaken are described below, and the statistical analysis was performed using R software, version 4.0.1 [33].

Study Population

A total of 154 patients suffering from CHF participated, comprising 119 men (77.3%) and 35 women (22.7%). However, only 50 patients (i.e., 32.5% of the total) who voluntarily provided their consent were enrolled in the pilot study. During the baseline assessment, all patients underwent medical history examination requiring a full physical examination. The medical history evaluation mainly focused on CV events, respiratory, and metabolic comorbidities, while key haemodynamic and anthropometric parameters were measured. The set of outpatients was meticulously monitored with a five-year follow-up conducted every 3 months. During each visit, six parameters were measured: weight, heart rate (HR), respiratory rate (RR), body temperature (BT), systolic blood pressure (SBP), and diastolic blood pressure (DBP). Clinical decompensation, with or without hospitalisation, and major CV events such as acute coronary syndrome, myocardial infarction, percutaneous transluminal coronary angioplasty (PTCA), surgical coronary artery revascularisation, stroke, death, and hospitalisation for any reason were reported. All clinical events had to be validated by source data (hospital records, death certificates or other original documents). Throughout the five-year follow-up, some patients missed a few medical visits, while others did not attend any appointments from the second-year onwards. Only eight patients completed the full five-year follow-up period.

Patient Characteristics

Among the 50 enrolled patients, 14 were women (28%) and 36 were men (72%). Table 1 summarises the demographic and clinical characteristics of the patients. The average age of the population was 72.5 ± 14.2 years. Among the patients, 3 (6%) showed a NYHA I class, 38 (76%) were in NYHA II class, and 9 (18%) in NYHA III class. The aetiology of CHF was represented as follows: ischaemic heart disease for 23 patients (46%); idiopathic dilated cardiomyopathy for 9 patients (18%); arterial hypertension for 4 patients (8%); valvular disease for 8 patients (16%); valvular disease that coexists with hypertension for 4 patients (8%); and alcoholic habit for 2 patients (4%). With regard to CV history, 1 patient (2%) had been revascularised by PTCA; 7 patients (14%) had undergone coronary artery bypass grafting; 13 patients (23%) had atrial flutter; 3 patients (6%) had undergone pacemaker implantation; 2 (4%) had implantable cardioverter defibrillator (ICD); 1 (2%) had undergone cardiac re-synchronisation therapy; 21 patients (42%) had mitral insufficiency; 4 patients (8%) had aortic insufficiency; 28 patients (56%) were hypertensive; and 2 patients (4%) suffered from a transitory ischaemic attack (TIA). In regard to the other comorbidities, 11 patients (22%) were diabetics, 1 patient (2%) had hypothyroidism, 4 patients (8%) suffered from chronic renal failure, 5 patients (10%) had chronic obstructive pulmonary disease (COPD), 1 patient (2%) suffered from bronchial asthma, 4 patients (8%) suffered from obstructive sleep apnoea, 4 patients (8%) suffered from gastrointestinal diseases and 3 patients (6%) suffered from liver disease. With regard to pharmacological treatment, 47 patients (94%) were on ACE-I/ARB, 29 patients (58%) were taking diuretic therapy, and 44 patients (88%) were receiving beta-blockers. No patients used corticosteroids or NSAIDs. Finally, 25 patients (50%) were taking triple therapy, 18 patients (36%) were receiving dual therapy with ACE-I/ARB - beta-blockers, 2 patients (4%) were on dual therapy with ACE-I/ARB and diuretic, whilst 5 patients (10%) were on single therapy.

Events during Follow-Up

During an average follow-up of 60 months, 19 patients presented a CV event. Among those patients who experienced an adverse event, 8 patients developed a second episode of decompensation. Notably, patients who manifested episodes of clinical instability were found to be older ( p < 0.05 ). Figure 2 displays the distribution of events based on the nine subtypes of CHF. Among the patients with ischaemic aetiology, 9 individuals (39.13%) had an event, while 2 patients (22.2%) were affected by idiopathic dilated cardiomyopathy; 2 patients (50%) had hypertensive aetiology; 3 patients (37.5%) had valvular disease; and 2 patients (50%) presented both valvular disease and hypertension. Interestingly, patients with alcoholic aetiology remained in stable clinical condition without any complications, in contrast to other types of CHF. Approximately half of the patients who experienced CV deterioration had a history of hypertension, and just under half had mitral insufficiency, while no patient with transient ischemic attack (TIA) presented further complications.

2.1. Data Preprocessing

A dataset consisting of 187 instances and suitable for the prediction task was created based on the collected data. Originally, the dataset was in a wide format with 50 rows, each containing personal data and a medical history of patients, along with visit dates, vital signs recorded at the visit, events occurring between the current and previous visits, and the corresponding event dates.
To perform the classification task, the data were converted into a long format where each row represents a patient’s visit. Input errors were corrected, outliers were discarded, and missing values were statistically imputed. Categorical variables with n values (e.g., aetiology, CV history, and other diseases) were converted into n numerical binary variables, while numerical data were expressed in their correct units of measurement. The resulting dataset comprises 794 instances and 37 features. Each instance was labelled as positive if there were CV deterioration events between two consecutive visits, and negative otherwise. Instances representing CV deterioration events were assigned to Class 1, while instances without any events were assigned to Class 2.
Imbalanced datasets are common in real-world applications and often become the focus of significant research efforts in knowledge discovery and data engineering. A dataset is considered imbalanced when one class, known as the minority class, is under-represented compared to the other class, which is the majority class. Relative imbalances frequently occur in practical scenarios, prompting extensive research in knowledge discovery and data engineering. For ML algorithms, imbalanced datasets pose a challenge because they tend to be overwhelmed by the abundance of majority class instances during the learning process. Therefore, methods are needed to improve recognition rates and address the issue of imbalanced data. In our dataset, the significant disparity between the number of negative instances (majority class) and positive instances (minority class) makes it imbalanced. Specifically, there are only 31 positive instances, accounting for a mere 4.6% of the entire dataset.

2.2. Feature Selection

Our aim was to develop a predictive model for the early detection of cardiovascular events in patients with CHF, utilising a limited set of basic clinical parameters. Through the feature selection process, vital signs were identified as the most informative factors. Vital signs, i.e., HR, RR, DPB, and SBP, were identified as the key factors for monitoring CHF and assessing the patient’s overall condition. These parameters play a crucial role in the predictive model’s design due to their significance in both CHF monitoring and patient assessment. Following the feature selection process, duplicate instances were removed from the dataset.

3. Machine Learning Process

In this section, we illustrate an overview of the entire machine learning process. We constructed, tested, and compared five ML predictive models, which are briefly introduced below. The details of this process are described in the subsequent section.
Supervised learning models were developed to accurately predict the risk of major events in CHF patients. The ML models implemented to develop the related prediction models are support vector machines (SVMs), artificial neural network, naïve Bayes, decision tree, and random forest.
SVM is based on the statistical learning theory [34,35] and is the most widely used ML technique available nowadays. It searches for an optimal hyperplane that separates the patterns of two classes by maximising the margin. Let X be a dataset with N instances X = ( x 1 , , x N ) , where x i , i = 1 , , N , denotes an instance with m features, and y i { ± 1 } its label. Finding the optimal hyperplane means solving the quadratic programming model (1)–(3)
m i n 1 2 | | w | | 2 + C 1 N ϵ i
y i ( w T ϕ ( x i ) + b ) 1 + ϵ i 0 i = 1 , , N
ϵ i 0 i = 1 , , N
where C, named penalty parameter, is a trade-off between the size of the margin and the slack variable penalty. In a non-linearly separable dataset, the SVM basically maps inputs into high-dimensional feature spaces by the so-called kernel functions. A kernel function denotes an inner product in a feature space, measures similarity between any pair of inputs x i and x j , and is usually denoted by K ( x i , x j ) = ϕ ( x i ) , ϕ ( x j ) [36]. Here, we used three kernel functions: linear kernel K ( x i , x j ) = x i , x j , polynomial kernel K ( x i , x j ) = ( x i , x j + 1 ) d , and the RBF kernel K ( x i , x j ) = e x p ( γ | | x i x j | | 2 ) . The linear kernels are a special case of polynomial kernels as the degree d is set to 1 and they compute similarity in the input space, whereas the other kernel functions compute similarity in the feature space.
Artificial neural networks are computational models, consisting of a number of artificial neural units. They emulate biological neural networks [37]. In this study, we used a feed-forward artificial neural network named multilayer perceptron (MLP) with a three-layer structure of neurons: an input layer, one or more hidden layers with a variable number of neurons, and an output classification layer. The neurons in the MLP are trained with the back propagation learning algorithm [38].
Naïve Bayes [39] is a probabilistic ML algorithm based on the Bayes Theorem. It assumes that a particular feature in a class is unrelated to the presence of any other feature.
Decision trees are a non-parametric supervised learning method [40]. One of their main advantages is that they are simple to understand and interpret, and they can be visualised.
Random forest [41] consists of individual decision trees that operate as an ensemble. Each tree is built by applying bagging, which is the general technique of bootstrap aggregation. A simple majority vote of all trees gives the final result. RF had good accuracy results in medical diagnosis problems [42,43].

3.1. Dealing with Imbalance Data: Cost-Sensitive Learning and Methods for Model Assessment

The problem addressed here is one of the most challenging issues in applied ML, as the event cases are under-represented. It is strongly important to correctly identify instances from the minority class compared to the majority class. To handle this problem, specific methods need to be employed. Usually, misclassification errors are treated equally but they are different depending on the class. In this work, we use and test cost-sensitive algorithms which involve the use of different misclassification costs.
A general cost matrix C o s t denotes the cost of each class misclassification [44]. A cost matrix for datasets with two classes is illustrated in Table 2.
In general, c i j is the cost of predicting an instance belonging to class i as belonging to class j. The goal of this type of learning is to minimise the total misclassification cost of a model on the training set. Formally, given a cost matrix C o s t and an instance x, the cost R ( i | x ) of classifying x into class i is R ( i | x ) = j p ( j | x ) c i j , where p ( j | x ) is the probability estimation of classifying an instance into class j. In the above cost matrix, c 12 represents the cost of a false positive misclassification and c 21 is the cost of a false negative misclassification. Usually, as we also assume here, there is no cost for correct classifications, that is, c 11 = c 22 = 0 . We tested different cost matrices, as detailed in the next section.

3.2. Hybrid Method for Imbalanced Dataset and Hyper-Parameter Optimisation Approach

Our dataset suffers from class imbalance, a common issue wherein classifiers developed on such data tend to be biased towards negative predictions due to the majority class having no-event cases. To address this problem and enable the better generalisation to new data, it is essential to use appropriate methods for handling imbalanced classes. The main methods for sampling-based imbalance correction are based on over-sampling and under-sampling. Over-sampling methods involves adding more data to the smaller class equalising its size with the larger class. On the other hand, under-sampling methods involve randomly selecting data from the larger class, matching its size with the smaller class. To tackle the class imbalance in our study, we adopted a hybrid method that balances the class distribution by combining over-sampling and under-sampling approaches. This approach adds data in the minority class while simultaneously removing data from the majority class, achieving a more balanced representation of both classes in the dataset.
To assess the models, we conducted a k-fold cross-validation process, ensuring the use of k independent sets to test the model, effectively simulating unseen data. The procedure basically consists in randomly partitioning the dataset into k equal-sized folds. During each of the k rounds, the k -th fold is the test set while the remaining folds are used as the training set. The test set is never used during the training of the model, preventing overfitting. Each fold is used exactly once as a test set, ensuring that each instance is used for testing exactly once. Figure 3 provides a schematic representation of the k-fold cross-validation process. The performance metrics are averaged across the k estimates from each test fold. In a well-designed k-fold cross-validation procedure, it is crucial to determine the training and test partitions before applying any oversampling technique. As thoroughly discussed in [45], conducting k-fold cross-validation on imbalanced data could have overly optimistic estimates if oversampling methods are applied to the entire dataset. To ensure a reliable evaluation of the model’s ability to generalise to real-world data, the oversampling should only be applied to the folds designated as training data. This approach maintains the integrity of the test sets, providing a more accurate assessment of the model’s performance.
The impact of hyper-parameters on the performance of an ML model is widely recognised. Generally, the hyper-parameters of each model are adjusted to find a hyper-parameter setting that maximises the model performances and enables accurate predictions on unseen data. Recent studies proved that hyper-parameter tuning with class weight optimisation are efficient in handling imbalanced data [46,47,48,49,50]. In this study, we employed the grid search optimisation strategy [51] to determine the optimal hyper-parameters for each ML model on the entire dataset. This strategy involves exploring all specified hyper-parameter combinations within a multi-dimensional grid. Each combination is evaluated using a performance metric to assess its effectiveness in enhancing the model’s performance. Multiple combinations of hyper-parameters are evaluated during the Grid Search optimisation process. Among these combinations, the one that yields the best performance is selected. Subsequently, this optimal set of hyper-parameters is utilised to train the ML model on the entire dataset.

3.3. Performance Metrics for Imbalanced Dataset

The predictive performance of the constructed ML models is evaluated and compared using various metrics, including the receiver operating characteristic (ROC—area under the curve (AUC), sensitivity, specificity, balanced accuracy, and G-mean values.
Let P and N be the numbers of positive and negative instances in the dataset, respectively. Let T P and T N be the numbers of instances correctly predicted as positive and negative, respectively, and F P and F N be the number of instances predicted as positive and negative while actually belonging to the opposite class, respectively.
  • AUC measures the classifier’s ability to avoid false classification [52]. It is the area under the curve of the true positive ratio vs. the false positive ratio that indicates the probability that the model will rank a positive case more highly than a negative case. A model whose predictions are 100% correct has an AUC of 1.0.
  • Sensitivity, also referred to as true positive rate or recall, measures the proportion of positive instances that are correctly identified, i.e., it is the ability to predict a CV event:
    Sens = T P T P + F N
  • Specificity (also known as true negative rate) is used to determine the ability to correctly classify. It measures the proportion of negatives that are correctly identified, and is defined as
    Spec = T N T N + F P
  • The accuracy metric can be misleading for our imbalanced dataset. As it is equally important to accurately predict the events of the positive and negative classes for the addressed problem, we used the balanced accuracy metric [53,54], which is defined as the arithmetic mean of sensitivity and specificity:
    B-acc = Sens + Spec 2 = T P R + T N R 2
  • Another useful metric is the so-called geometric mean or G-Mean that balances both sensitivity and specificity by combining them. It is defined as
    G-Mean = s q r t ( S e n s S p e c )
The sensitivity and specificity are also known as quality parameters and used to define the quality of the predicted class.

4. Results and Discussions

For predicting cardiovascular deterioration events in CHF patients, we only used the vital signs as predictive features for a deterioration event. With the aim of finding an ML model that is able to predict CV events, we carried out computational experiments with three methods, named as follows:
  • Method1: cost-sensitive learning methods;
  • Method2: data resampling methods;
  • Method3: cost-sensitive learning methods combined with data resampling methods.
In cost-sensitive learning, we considered the cost of misclassifying a positive instance and the cost of misclassifying a negative instance as hyper-parameters to be tuned during the model training process. Referring to the notation of Table 2, we fixed the misclassification cost c 12 = 1 and explored different values for
c 21 { 1 , 1.5 , 2 , 5 , 10 , 20 , 30 , 40 , 50 , 100 , 200 , 300 , 400 } , while setting c 11 = c 22 = 0 .
The Waikato Environment for Knowledge Analysis (Weka, version 3.8.2, [55]) was utilised for constructing and evaluating the classification models. For each ML model, we used the grid search algorithm to select the optimal values of the cost matrix.
Table 3 presents the hyper-parameter tuning performed with a data resampling approach for SVM and MLP. The final column displays the optimal hyper-parameter values, as follows:
 SVM: 
We tested SVM with three kernel functions, i.e., the linear kernel, polynomial kernel, and RBF kernel. The hyper-parameter C and those related to the polynomial kernel and RBF kernel were optimised by searching for the best values within the specified range, as reported in Table 3. The incremental value used for optimisation is denoted in the fourth column as “Step.”
 MLP: 
The model parameter optimisation involves choosing the number of hidden layers, the number of neurons in each layer, the number of epochs, the learning rate, and the momentum. We tested different MLP models consisting of one input layer, one output layer, and one hidden layer. The number of neurons in the hidden layer was set according to the formula
n r . n e u r o n s = n r . o f i n p u t f e a t u r e s + n r . o f c l a s s e s 2
We optimised the learning rate, the momentum, and the number of epochs in a defined range of values, as reported in Table 3. The incremented value is denoted in the fourth column as “Step”.
Table 3 presents the values tested for optimising the hyper-parameters of SVM and MLP models through the grid search algorithm. The results of the hyper-parameter tuning for both SVM and MLP models using a cost matrix with equal misclassification costs (i.e., false negative cost and false positive cost) are shown as optimal hyper-parameter values in the last column.
Table 4 displays the hyper-parameter values set for each model, alongside the corresponding cost matrix with equal misclassification costs. The table provides a comprehensive overview of the selected hyper-parameter values for each model in the study and the same cost matrix.

Predictive Models Performance Metrics

We carried out computational experiments with the three methods. The following tables only report the best performance results. The performance metrics are related to both three-fold cross-validation and five-fold cross-validation. The best values for each performance metric achieved by a ML model are highlighted in bold.
Table 5 summarises the best results related to M e t h o d 1 , i.e., cost-sensitive learning methods. The SVM models with the linear kernel and RBF kernel demonstrated the best overall performances compared to the other models when the misclassification costs were set to c 12 = 1 and c 21 = 30 . The decision tree models exhibited a comparable prediction performance in terms of sensitivity and specificity only when the cost of misclassifying the minority class was set to a high value, that is, c 21 = 200 or c 21 = 300 . Naive Bayes models showed a comparable prediction performance with both c 21 = 30 and c 21 = 40 with five-fold cross-validation. Random forest models showed a comparable prediction performance only when the cost of misclassifying the minority class was set to a high value, that is, c 21 = 500 . Furthermore, it is noteworthy that, in general, hyper-parameters that allow one to find the highest sensitivity have the lowest specificity, and vice versa.
Table 6 reports the performance of the predictive models with M e t h o d 2 and M e t h o d 3 . More specifically, the first row shows the results found by M e t h o d 2 per each ML model reported in the first columns; the rest of the rows show the results found by M e t h o d 3 . Based on the performance results, the MLP with c 12 = 1 , c 21 = 2 has the best performance among the built models with a mean sensitivity and a mean specificity of 65% and 55%, respectively, a mean area under the curve of 0.71, and a G-mean of 0.60. In general, we observed that the performance of the predictive models MLP, naive Bayes, decision tree, and random forest improved when the cost of misclassifying the minority class was higher ( c 21 > c 12 ). Additionally, all constructed models achieved a meaningful prediction performance with a five-fold cross-validation approach. However, the SVM model with a polynomial kernel showed a lower sensitivity performance in certain cases. These findings indicate that the combination of cost-sensitive methods with data over/under-sampling approaches is effective for the meaningful prediction of cardiovascular events. Comparing the three methods, the cost-sensitive learning methods proved to be superior to the sampling approach. They achieved a high performance in terms of G-mean, indicating their efficacy in handling imbalanced data and improving the model’s overall performance.
This model could be particularly useful in CHF patients with NYHA class III or IV, where functional reserves are reduced and each exacerbation leads to a further deterioration of cardiac function that worsens the symptomatology and often requires hospitalisation and sometimes results in death. In this context, the early identification and treatment of a re-exacerbation of HF may improve the patient’s symptoms and prognosis and avoid hospitalisation with reduced healthcare costs [6].

5. Conclusions

Technology is increasingly playing a significant role in clinical practice. In this context, machine learning represents an innovative methodology for managing chronic disorders, empowering clinicians with a key role in predicting cardiovascular events rather than merely being spectators.
The data used in this study were collected during a pilot study in a well-characterised CHF population of patients. The results open up numerous potential perspectives for applying ML approaches to clinical practice. Having a predictive system for CV deterioration events in CHF patients can lead to significant advantages, such as reducing hospitalisations and associated costs. The clinical importance of utilising such a model cannot be understated, as early detection can prevent CV events, improve patient health, and optimise healthcare expenditure.
The findings suggest that cost-sensitive methods can effectively predict CV deterioration events in CHF patients using only a few clinical variables. CHF remains one of the most severe chronic diseases in terms of mortality, hospitalisation rate, and healthcare costs. Successfully addressing even one variable linked to the natural course of this disease would be a significant achievement, benefiting the patients, clinicians, and the National Health System. Our KD process yielded promising results, paving the way for large-scale application. It is interesting that a few clinical variables such as HR, RR, DBP, and SBP had a good performance in the prediction of CV deterioration events. Furthermore, the performance of these models may be enhanced by including other variables of study patients. Further studies, however, are needed to validate this approach in a more large-scaled patient population. As part of future work, we plan to expand the sample size by including more patients, reducing the interval between consecutive visits using remote monitoring, and conduct an in-depth feature selection study to better understand which other features are crucial in diagnosing deterioration events in CHF patients.

Author Contributions

Conceptualisation, M.C.G., R.G., D.C. and A.S.; Formal analysis, M.C.G., R.G., D.C. and A.S.; Investigation, M.C.G., R.G., C.P., G.A., A.F.T., S.M., E.S., M.L.H. and A.S.; Methodology, M.C.G., R.G., D.C., E.S., M.L.H. and A.S.; Software, M.C.G.; Validation, M.C.G.; Writing—original draft, M.C.G., R.G., D.C., E.S., M.L.H. and A.S.; Data acquisition, C.P., G.A., A.F.T., S.M., E.S., M.L.H. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the research project “SI.F.I.PA.CRO.DE.—Sviluppo e industrializzazione farmaci innovativi per terapia molecolare personalizzata PA.CRO.DE. (PON ARS0100568, CUP: B29C20000360005, CONCESSIONE RNA-COR: 4646672), Italian Ministry of University and Research, 2021”.

Informed Consent Statement

The studies involving human participants were reviewed and approved by the institutional review board at the Geriatrics Division of the “Mater Domini” University Hospital, allowing a retrospective review of medical records and granting a waiver of informed consent. Written informed consent was obtained from all subjects involved in the data collection in accordance with the national legislation and the institutional requirements. Patients cannot be identified in this study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet 2018, 392, 1789–1858. [Google Scholar] [CrossRef] [PubMed]
  2. McMurray, J.J. Improving outcomes in heart failure: A personal perspective. Eur. Heart J. 2015, 36, 3467–3470. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, T.; Evans, J.C.; Benjamin, E.; Lévy, D.; Leroy, E.; Vasan, R. Natural History of Asymptomatic Left Ventricular Systolic Dysfunction in the Community. Circ. J. Am. Heart Assoc. 2003, 108, 977–982. [Google Scholar] [CrossRef] [PubMed]
  4. Dunlay, S.M.; Redfield, M.M.; Weston, S.A.; Therneau, T.M.; Long, K.H.; Shah, N.D.; Roger, V.L. Hospitalizations after Heart Failure Diagnosis. J. Am. Coll. Cardiol. 2009, 54, 1695–1702. [Google Scholar] [CrossRef]
  5. Ponikowski, P.; Voors, A.A.; Anker, S.D.; Bueno, H.; Cleland, J.G.F.; Coats, A.J.S.; Falk, V.; González-Juanatey, J.R.; Harjola, V.P.; Jankowska, E.A.; et al. 2016 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure: The Task Force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC)Developed with the special contribution of the Heart Failure Association (HFA) of the ESC. Eur. Heart J. 2016, 37, 2129–2200. [Google Scholar] [CrossRef]
  6. McDonagh, T.A.; Metra, M.; Adamo, M.; Gardner, R.S.; Baumbach, A.; Bohm, M.; Burri, H.; Butler, J.; Celutkien, J.; Chioncel, O.; et al. 2021 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure: Developed by the Task Force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC) With the special contribution of the Heart Failure Association (HFA) of the ESC. Eur. Heart J. 2021, 42, 3599–3726. [Google Scholar]
  7. Stevenson, L.W.; Pagani, F.D.; Young, J.B.; Jessup, M.; Miller, L.; Kormos, R.L.; Naftel, D.C.; Ulisney, K.; Desvigne-Nickens, P.; Kirklin, J.K. INTERMACS profiles of advanced heart failure: The current picture. J. Heart Lung Transplant. Off. Publ. Int. Soc. Heart Transplant. 2009, 28, 535–541. [Google Scholar] [CrossRef]
  8. Ziaeian, B.; Fonarow, G. Epidemiology and aetiology of heart failure. Nat. Rev. Cardiol. 2016, 13, 368–378. [Google Scholar] [CrossRef]
  9. Mehta, P.A.; Dubrey, S.W.; McIntyre, H.F.; Walker, D.M.; Hardman, S.M.C.; Sutton, G.C.; McDonagh, T.A.; Cowie, M.R. Improving survival in the 6 months after diagnosis of heart failure in the past decade: Population-based data from the UK. Heart 2009, 95, 1851–1856. [Google Scholar] [CrossRef]
  10. Tavazzi, L.; Senni, M.; Metra, M.; Gorini, M.; Cacciatore, G.; Chinaglia, A.; Lenarda, A.D.; Mortara, A.; Oliva, F.; Maggioni, A.P. Multicenter Prospective Observational Study on Acute and Chronic Heart Failure. Circ. Heart Fail. 2013, 6, 473–481. [Google Scholar] [CrossRef]
  11. Liao, L.; Allen, L.A.; Whellan, D.J. Economic burden of heart failure in the elderly. PharmacoEconomics 2008, 26, 447–462. [Google Scholar] [CrossRef]
  12. Stewart, S.; Jenkins, A.; Buchan, S.; McGuire, A.; Capewell, S.; McMurray, J. The current cost of heart failure to the National Health Service in the UK. Eur. J. Heart Fail. 2002, 4, 361–371. [Google Scholar] [CrossRef] [PubMed]
  13. Marangoni, E.; Lissoni, F.; Raimondi Cominesi, I.; Tinelli, S. Heart failure: Epidemiology, costs and healthcare programs in Italy. G. Ital. Cardiol. 2012, 13, 139S–144S. [Google Scholar] [CrossRef]
  14. Krumholz, H.M.; Chen, Y.T.; Vaccarino, V.; Wang, Y.; Radford, M.J.; Bradford, W.; Horwitz, R.I. Correlates and impact on outcomes of worsening renal function in patients ≥ years of age with heart failure. Am. J. Cardiol. 2000, 85, 1110–1113. [Google Scholar] [CrossRef] [PubMed]
  15. Brons, M.; Koudstaal, S.; Asselbergs, F.W. Algorithms used in telemonitoring programmes for patients with chronic heart failure: A systematic review. Eur. J. Cardiovasc. Nurs. 2018, 17, 580–588. [Google Scholar] [CrossRef]
  16. Kurtz, B.; Lemercier, M.; Pouchin, S.C.; Benmokhtar, E.; Vallet, C.; Cribier, A.; Bauer, F. Automated home telephone self-monitoring reduces hospitalization in patients with advanced heart failure. J. Telemed. Telecare 2011, 17, 298–302. [Google Scholar] [CrossRef]
  17. Olsen, C.R.; Mentz, R.J.; Anstrom, K.J.; Page, D.; Patel, P.A. Clinical applications of machine learning in the diagnosis, classification, and prediction of heart failure. Am. Heart J. 2020, 229, 1–17. [Google Scholar] [CrossRef] [PubMed]
  18. Averbuch, T.; Sullivan, K.; Sauer, A.; Mamas, M.A.; Voors, A.A.; Gale, C.P.; Metra, M.; Ravindra, N.; Van Spall, H.G.C. Applications of artificial intelligence and machine learning in heart failure. Eur. Heart J.-Digit. Health 2022, 3, 311–322. [Google Scholar] [CrossRef] [PubMed]
  19. Lofaro, D.; Groccia, M.C.; Guido, R.; Conforti, D.; Caroleo, S.; Fragomeni, G. Machine learning approaches for supporting patient-specific cardiac rehabilitation programs. In Proceedings of the 2016 Computing in Cardiology Conference (CinC), Vancouver, BC, Canada, 11–14 September 2016; pp. 149–152. [Google Scholar]
  20. Groccia, M.C.; Lofaro, D.; Guido, R.; Conforti, D.; Sciacqua, A. Predictive Models for Risk Assessment of Worsening Events in Chronic Heart Failure Patients. In Proceedings of the 2018 Computing in Cardiology Conference (CinC), Maastricht, The Netherlands, 23–26 September 2018; Volume 45, pp. 1–4. [Google Scholar] [CrossRef]
  21. Tripoliti, E.E.; Papadopoulos, T.G.; Karanasiou, G.S.; Naka, K.K.; Fotiadis, D.I. Heart failure: Diagnosis, severity estimation and prediction of adverse events through machine learning techniques. Comput. Struct. Biotechnol. J. 2017, 15, 26–47. [Google Scholar] [CrossRef]
  22. Lorenzoni, G.; Sabato, S.S.; Lanera, C.; Bottigliengo, D.; Minto, C.; Ocagli, H.; De Paolis, P.; Gregori, D.; Iliceto, S.; Pisanò, F. Comparison of machine learning techniques for prediction of hospitalization in heart failure patients. J. Clin. Med. 2019, 8, 1298. [Google Scholar] [CrossRef]
  23. Lemaître, G.; Nogueira, F.; Aridas, C.K. Imbalanced-learn: A python toolbox to tackle the curse of imbalanced datasets in machine learning. J. Mach. Learn. Res. 2017, 18, 559–563. [Google Scholar]
  24. Chawla, N.V.; Bowyer, K.W.; Hall, L.O.; Kegelmeyer, W.P. SMOTE: Synthetic minority over-sampling technique. J. Artif. Intell. Res. 2002, 16, 321–357. [Google Scholar] [CrossRef]
  25. Beckmann, M.; Ebecken, N.F.; de Lima, B.S.P. A KNN undersampling approach for data balancing. J. Intell. Learn. Syst. Appl. 2015, 7, 104. [Google Scholar] [CrossRef]
  26. Akbani, R.; Kwek, S.; Japkowicz, N. Applying support vector machines to imbalanced datasets. In Proceedings of the European Conference on Machine Learning, Pisa, Italy, 20–24 September 2004; Springer: Berlin/Heidelberg, Germany, 2004; pp. 39–50. [Google Scholar]
  27. Cheng, F.; Zhang, J.; Wen, C. Cost-sensitive large margin distribution machine for classification of imbalanced data. Pattern Recognit. Lett. 2016, 80, 107–112. [Google Scholar] [CrossRef]
  28. Veropoulos, K.; Campbell, C.; Cristianini, N. Controlling the sensitivity of support vector machines. In Proceedings of the International Joint Conference on AI, Stockholm, Sweden, 31 July–6 August 1999; Volume 55, p. 60. [Google Scholar]
  29. Qi, Z.; Tian, Y.; Shi, Y.; Yu, X. Cost-sensitive support vector machine for semi-supervised learning. Procedia Comput. Sci. 2013, 18, 1684–1689. [Google Scholar] [CrossRef]
  30. Tao, X.; Li, Q.; Guo, W.; Ren, C.; Li, C.; Liu, R.; Zou, J. Self-adaptive cost weights-based support vector machine cost-sensitive ensemble for imbalanced data classification. Inf. Sci. 2019, 487, 31–56. [Google Scholar] [CrossRef]
  31. Thai-Nghe, N.; Gantner, Z.; Schmidt-Thieme, L. Cost-sensitive learning methods for imbalanced data. In Proceedings of the The 2010 International Joint Conference on Neural Networks (IJCNN), Barcelona, Spain, 18–23 July 2010; IEEE: Piscataway, NJ, USA, 2010; pp. 1–8. [Google Scholar]
  32. Kohavi, R.; Provost, F. Glossary of terms. Mach. Learn. 1998, 30, 271–274. [Google Scholar]
  33. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017.
  34. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  35. Burges, C. A Tutorial on Support Vector Machines for Pattern Recognition. Data Min. Knowl. Discov. 1998, 2, 121–167. [Google Scholar] [CrossRef]
  36. Hofmann, T.; Scholkopf, B.; Smola, A.J. Kernel Methods in Machine Learning. Ann. Statist. 2008, 36, 1171–1220. [Google Scholar] [CrossRef]
  37. Krenker, A.; Bešter, J.; Kos, A. Introduction to the artificial neural networks. In Artificial Neural Networks: Methodological Advances and Biomedical Applications; InTech: Houston, TX, USA, 2011; pp. 1–18. [Google Scholar]
  38. Rumelhart, D.E.; Hinton, G.E.; Williams, R.J. Learning Internal Representations by Error Propagation; Technical Report; California Univ San Diego La Jolla Inst for Cognitive Science: La Jolla, CA, USA, 1985. [Google Scholar]
  39. Rish, I. An Empirical Study of the Naive Bayes Classifier. Technical Report. 2001. Available online: https://www.cc.gatech.edu/home/isbell/classes/reading/papers/Rish.pdf (accessed on 29 September 2023).
  40. Shalev-Shwartz, S.; Ben-David, S. Decision Trees. In Understanding Machine Learning; Cambridge University Press: Cambridge, UK, 2014; Chapter 18. [Google Scholar]
  41. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  42. Alam, M.Z.; Rahman, M.S.; Rahman, M.S. A Random Forest based predictor for medical data classification using feature ranking. Inform. Med. Unlocked 2019, 15, 100180. [Google Scholar] [CrossRef]
  43. Yang, F.; Wang, H.; Mi, H.; Lin, C.; Cai, W. Using random forest for reliable classification and cost-sensitive learning for medical diagnosis. BMC Bioinform. 2009, 10, S22. [Google Scholar] [CrossRef] [PubMed]
  44. Elkan, C. The Foundations of Cost-Sensitive Learning. In Proceedings of the 17th International Joint Conference on Artificial Intelligence—Volume 2, Seattle, WA, USA, 4–10 August 2001; IJCAI’01. Morgan Kaufmann Publishers Inc.: San Francisco, CA, USA, 2001; pp. 973–978. [Google Scholar]
  45. Santos, M.S.; Soares, J.P.; Abreu, P.H.; Araujo, H.; Santos, J. Cross-Validation for Imbalanced Datasets: Avoiding Overoptimistic and Overfitting Approaches [Research Frontier]. IEEE Comput. Intell. Mag. 2018, 13, 59–76. [Google Scholar] [CrossRef]
  46. Kong, J.; Kowalczyk, W.; Menzel, S.; Bäck, T. Improving Imbalanced Classification by Anomaly Detection. In Proceedings of the 16th International Conference, PPSN 2020, Leiden, The Netherlands, 5–9 September 2020; pp. 512–523. [Google Scholar] [CrossRef]
  47. Mienye, I.D.; Sun, Y. Performance analysis of cost-sensitive learning methods with application to imbalanced medical data. Inform. Med. Unlocked 2021, 25, 100690. [Google Scholar] [CrossRef]
  48. Guido, R.; Groccia, M.C.; Conforti, D. A hyper-parameter tuning approach for cost-sensitive support vector machine classifiers. Soft Comput. 2023, 27, 12863–12881. [Google Scholar] [CrossRef]
  49. Guido, R.; Groccia, M.C.; Conforti, D. Hyper-Parameter Optimization in Support Vector Machine on Unbalanced Datasets Using Genetic Algorithms. In Optimization in Artificial Intelligence and Data Sciences; Amorosi, L., Dell’Olmo, P., Lari, I., Eds.; Springer International Publishing: Cham, Switzerland, 2022; pp. 37–47. [Google Scholar]
  50. Zhang, F.; Petersen, M.; Johnson, L.; Hall, J.; O’Bryant, S.E. Hyperparameter Tuning with High Performance Computing Machine Learning for Imbalanced Alzheimer Disease Data. Appl. Sci. 2022, 12, 6670. [Google Scholar] [CrossRef]
  51. Bergstra, J.; Bengio, Y. Random search for hyper-parameter optimization. J. Mach. Learn. Res. 2012, 13, 281–305. [Google Scholar]
  52. Rumelhart, D.E.; McClelland, J.L.; PDP Research Group, C. (Eds.) Parallel Distributed Processing: Explorations in the Microstructure of Cognition, Vol. 1: Foundations; MIT Press: Cambridge, MA, USA, 1986. [Google Scholar]
  53. Sun, Y.; Wong, A.K.C.; Kamel, M.S. Classification of Imbalanced Data: A Review. Int. J. Pattern Recognit. Artif. Intell. 2009, 23, 687–719. [Google Scholar] [CrossRef]
  54. Branco, P.; Torgo, L.; Ribeiro, R. A Survey of Predictive Modelling under Imbalanced Distributions. arXiv 2015, arXiv:cs.LG/1505.01658. [Google Scholar]
  55. Eibe, F.; Hall, M.A.; Witten, I.H. The WEKA Workbench. In Online Appendix for Data Mining: Practical Machine Learning Tools and Techniques; Morgan Kaufmann: Burlington, MA, USA, 2016. [Google Scholar]
Figure 1. Our knowledge discovery process.
Figure 1. Our knowledge discovery process.
Information 14 00542 g001
Figure 2. Distribution of the events based on chronic heart failure aetiology.
Figure 2. Distribution of the events based on chronic heart failure aetiology.
Information 14 00542 g002
Figure 3. Graphical representation of the k-fold cross-validation method.
Figure 3. Graphical representation of the k-fold cross-validation method.
Information 14 00542 g003
Table 1. Demographic and clinical characteristics of the patients. NYHA: New York Heart Association; CHF: chronic heart failure; PTCA: percutaneous transluminal coronary angioplasty; ICD: implantable cardioverter defibrillator; TIA: transient ischaemic attack; COPD: chronic obstructive pulmonary disease.
Table 1. Demographic and clinical characteristics of the patients. NYHA: New York Heart Association; CHF: chronic heart failure; PTCA: percutaneous transluminal coronary angioplasty; ICD: implantable cardioverter defibrillator; TIA: transient ischaemic attack; COPD: chronic obstructive pulmonary disease.
Characteristic
Age(years ± SD)72.5 ± 14.2
SexMale36 (72%)
Female14 (28%)
NYHA ClassI3 (6%)
II38 (76%)
III9 (18%)
CHF aetiologyIschemic heart disease23 (46%)
Idiopathic dilatation9 (18%)
Hypertension4 (8%)
Valvular diseases8 (16%)
Valvular diseases + hypertension4 (8%)
Alcoholic habit2 (4%)
Cardiovascular historyInstable angina1 (2%)
PTCA1 (2%)
By-pass7 (14%)
Atrial flutter13 (26%)
Pacemaker3 (6%)
Cardiac resynchronisation1 (2%)
ICD2 (4%)
Mitral insufficiency21 (42%)
Aortic insufficiency4 (8%)
Hypertension28 (56%)
TIA2 (4%)
Other diseasesDiabetes11 (22%)
Hypothyroidism1 (2%)
Renal failure4 (8%)
COPD5 (10%)
Asthma1 (2%)
Sleep apnea4 (8%)
Pulmonary fibrosis1 (2%)
Gastrointestinal diseases4 (8%)
Hepatic diseases3 (6%)
Pharmacological treatmentACE-I/ARB47 (94%)
Diuretic therapy29 (58%)
Beta-blockers44 (88%)
Corticosteroids/NSAIDs0 (0%)
Table 2. Cost matrix for a binary classification problem.
Table 2. Cost matrix for a binary classification problem.
Actual Class
PositiveNegative
Predicted classPositive c 11 c 12
Negative c 21 c 22
Table 3. Hyper-parameter tuning for SVM and MLP and best values.
Table 3. Hyper-parameter tuning for SVM and MLP and best values.
ML ModelParameterSearch SpaceStepSelected Value
SVMd[1, 5]12
(polynomial kernel)C[1, 10]0.55
SVM γ [0.01, 1.00]0.010.82
(RBF kernel)C[1, 10]0.53
MLPLearning rate[0.1, 1.0]0.10.6
Momentum[0.1, 1.0]0.10.1
Number of epochs[400, 600]100500
Table 4. Optimal hyper-parameters values.
Table 4. Optimal hyper-parameters values.
ModelParameterValue
SVM with linear kernelC1
SVM with polynomial kerneld2
C5
SVM with RBF kernel γ 0.82
C3
MLPLearning rate0.6
Momentum0.1
Number of epochs500
Number of neurons in input layer4
Number of neurons in hidden layer3
Number of neurons in output layer2
Naive BayesuseKernelEstimatorFalse
useSupervisedDiscretisationFalse
Decision treeConfidence factor0.25
Random forestBag size percent100
Number of iterations100
Table 5. Predictive cost-sensitive learning model performance metrics by the cost-sensitive method.
Table 5. Predictive cost-sensitive learning model performance metrics by the cost-sensitive method.
3-Fold Cross-Validation5-Fold Cross-Validation
Model c12c21 AUC Sens Spec B-acc G-mean AUC Sens Spec B-acc G-mean
SVM with
linear kernel
1100.5350.1300.9400.5360.2780.6110.2570.9640.6110.498
1200.5970.3940.8010.5980.5420.7030.5860.8200.7030.693
1300.5760.5580.5940.5760.5570.6700.7190.6200.6700.668
1350.4990.6210.3760.4990.4290.6010.7520.4490.6010.581
1400.4820.7880.1760.4820.1820.5350.8810.1890.5350.408
SVM with
polynomial
kernel
1100.5420.130.9540.5420.2780.5830.1950.9720.5840.435
1200.6090.330.8890.6100.5150.6350.3570.9130.6350.571
1300.5440.3970.6920.5450.4880.5860.8430.3290.5860.527
1350.5440.7450.3410.5440.4850.5650.7190.4120.5660.544
1400.4720.6210.3220.4720.4140.5440.7860.3030.5450.488
SVM with
RBF kernel
1100.5570.1610.9520.5570.3170.5830.190.9760.5830.431
1200.5460.2270.8650.5460.4280.6610.4240.8980.6610.617
1300.5780.5580.5990.5790.5590.6770.7190.6350.6770.676
1350.5510.5550.5480.5520.5340.5890.690.4870.5890.580
1400.5050.5580.4520.5050.4700.5480.7240.3710.5480.518
MLP1100.5790.4640.6410.5530.2790.720.3620.830.5960.548
1150.5280.4640.6520.5580.2800.7160.6330.3820.5080.492
1200.6091.0000.0000.5000.0000.5591.0000.0000.5000.000
Naive
Bayes
1100.6190.3240.8860.6060.5280.6670.2570.910.5840.484
1200.6190.3910.8070.5990.5550.6670.3900.8260.6080.568
1300.6190.4550.7240.5900.5650.6670.4900.7460.6180.605
1400.6190.4850.6210.5530.5390.6670.5860.6550.6210.620
1500.6190.5180.5100.5150.5090.6670.6190.5620.5910.590
1600.6190.6150.4480.5320.519-----
Decision
tree
1100.5530.1970.9160.5570.4120.5260.1330.9210.5270.35
1200.560.1970.9240.5610.4150.5270.1330.9190.5260.350
1300.5630.2300.8980.5640.4340.5530.2000.8970.5490.424
1400.5330.1970.8700.5340.4020.5680.2330.8920.5630.456
1500.5430.2300.8630.5470.4250.5770.2670.8740.5710.483
11000.5720.3640.7730.5680.4770.5780.3240.8050.5650.511
12000.5540.4300.6720.5510.4840.6180.5900.6110.6010.600
13000.5380.5240.5150.5200.4960.6420.6520.4730.5630.555
13500.5160.6910.2960.4930.292-----
14000.5160.6910.2960.4930.292-----
Random
forest
1100.5930.1640.9610.5630.3800.6390.1000.9720.5360.312
1200.5790.1640.9610.5630.3800.6170.1000.9640.5320.310
1300.5960.1640.9580.5610.3790.6180.2000.9580.5790.438
1400.5850.1970.9540.5760.4050.6260.1670.9570.5620.400
1500.5720.1970.9520.5750.4050.6240.2000.9490.5750.436
11000.5870.1640.9290.5470.3740.6290.2620.9300.5960.494
12000.5950.2300.8790.5550.4130.6010.3240.8650.5950.529
13000.5680.3640.7870.5760.4960.5970.3900.7590.5750.544
14000.5680.3300.6840.5080.4290.5990.4860.6900.5880.579
15000.5620.3970.6100.5080.4290.5880.5520.6020.5770.576
Table 6. Predictive model performance metrics related to the data resampling method, and the cost-sensitive method combined with the data resampling method.
Table 6. Predictive model performance metrics related to the data resampling method, and the cost-sensitive method combined with the data resampling method.
3-Fold Cross-Validation5-Fold Cross-Validation
Model c12c21 AUC Sens Spec B-acc G-mean AUC Sens Spec B-acc G-mean
SVM with
linear kernel
110.5730.3940.7520.5730.5220.6800.5570.8020.6800.668
11.50.5220.5270.5170.5730.5220.6700.7810.5600.6710.661
120.4590.6850.2320.4590.2790.5450.8760.2140.5450.433
150.5001.0000.0000.5000.0000.5001.0000.0000.5000.000
SVM with
polynomial
kernel
110.5360.3550.7170.5360.4980.5780.3240.8320.5780.519
11.50.5690.8120.3250.5360.4980.5770.8430.3100.5770.511
120.5410.8120.2710.5420.4470.5200.8140.2260.5200.429
150.5150.8420.1870.5150.3740.4710.8430.0990.4710.289
1100.4760.9360.0170.4770.0990.5051.0000.0090.5050.095
1200.5001.0000.0000.5000.0000.5001.0000.0000.5000.000
1300.4730.9060.0410.4740.1450.5031.0000.0060.5030.077
SVM with
RBF kernel
110.5960.4270.7630.5960.5500.6440.4520.8350.6440.614
11.50.5190.4580.5800.5960.5500.6380.6570.6190.6380.638
120.5210.5580.4840.5210.4800.5860.7480.4250.5870.564
150.4660.8120.1200.4660.2340.4650.8760.0540.4650.218
1100.5001.0000.0000.5000.0000.5001.0000.0000.5000.000
MLP110.6340.3970.6800.5380.4750.6590.5570.7640.6610.652
120.6060.6970.2930.4950.0940.7150.6520.5550.6040.602
150.5980.6520.5740.6130.5910.6950.7810.3960.5890.556
1100.5480.7670.3510.5590.3850.6630.7240.3480.5360.502
1200.4980.6520.3850.5190.4460.6290.7900.1360.4630.328
1300.5330.6320.4160.5240.5000.5700.7520.1800.4660.368
Naive
Bayes
110.6190.3910.7890.5900.5480.6800.4240.7970.6110.581
11.50.6190.4880.6810.5900.5480.6800.5520.7160.6340.629
120.6190.5480.5880.5690.5600.6800.6190.6070.6130.613
150.6190.8790.2010.5400.3880.6800.8430.2020.5230.413
1100.6190.8790.0930.4860.1700.6800.9100.04340.4770.199
1200.6190.9390.0570.4980.1250.6801.0000.0120.5060.110
1300.6190.9390.0390.4890.1030.6801.0000.0060.5030.078
Decision
tree
110.5240.1640.8840.5240.2930.6170.3570.8640.6110.555
120.5400.1970.8740.5360.3890.6020.3240.8710.5980.531
150.5150.2970.7540.5260.4450.6150.4570.7410.5990.582
1100.5000.4640.5200.4920.4180.6210.5860.6180.6020.602
1200.5090.5550.4300.4930.4600.6130.6480.4480.5480.539
1300.4750.6210.3390.4800.3040.4730.7480.1550.4520.340
Random
forest
110.5950.1640.9210.5420.3720.6400.1950.9310.5630.426
120.5700.2300.9160.5730.4220.6380.2570.9330.5950.490
150.5990.2300.8790.5550.4310.6380.2570.8950.5760.480
1100.5800.3300.8070.5690.4850.6130.3900.8280.6090.568
1200.5680.4300.6300.5310.4730.6340.5520.6750.6140.610
1300.5720.5580.5050.5310.5000.6010.6190.5030.5610.558
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Groccia, M.C.; Guido, R.; Conforti, D.; Pelaia, C.; Armentaro, G.; Toscani, A.F.; Miceli, S.; Succurro, E.; Hribal, M.L.; Sciacqua, A. Cost-Sensitive Models to Predict Risk of Cardiovascular Events in Patients with Chronic Heart Failure. Information 2023, 14, 542. https://doi.org/10.3390/info14100542

AMA Style

Groccia MC, Guido R, Conforti D, Pelaia C, Armentaro G, Toscani AF, Miceli S, Succurro E, Hribal ML, Sciacqua A. Cost-Sensitive Models to Predict Risk of Cardiovascular Events in Patients with Chronic Heart Failure. Information. 2023; 14(10):542. https://doi.org/10.3390/info14100542

Chicago/Turabian Style

Groccia, Maria Carmela, Rosita Guido, Domenico Conforti, Corrado Pelaia, Giuseppe Armentaro, Alfredo Francesco Toscani, Sofia Miceli, Elena Succurro, Marta Letizia Hribal, and Angela Sciacqua. 2023. "Cost-Sensitive Models to Predict Risk of Cardiovascular Events in Patients with Chronic Heart Failure" Information 14, no. 10: 542. https://doi.org/10.3390/info14100542

APA Style

Groccia, M. C., Guido, R., Conforti, D., Pelaia, C., Armentaro, G., Toscani, A. F., Miceli, S., Succurro, E., Hribal, M. L., & Sciacqua, A. (2023). Cost-Sensitive Models to Predict Risk of Cardiovascular Events in Patients with Chronic Heart Failure. Information, 14(10), 542. https://doi.org/10.3390/info14100542

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop