SVM ensemble training for imbalanced data classification using multi-objective optimization techniques

One of the main problems with classifier training for imbalanced data is defining the correct learning criterion. On the one hand, we want the minority class to be correctly recognized, and on the other hand, we do not want to make too many mistakes in the majority class. Commonly used metrics focus either on the predictive quality of the distinguished class or propose an aggregation of simple metrics. The aggregate metrics, such as Gmean or AUC, are primarily ambiguous, i.e., they do not indicate the specific values of errors made on the minority or majority class. Additionally, improper use of aggregate metrics results in solutions selected with their help that may favor the majority class. The authors realize that a solution to this problem is using overall risk. However, this requires knowledge of the costs associated with errors made between classes, which is often unavailable. Hence, this paper will propose the semoos algorithm - an approach based on multi-objective optimization that optimizes criteria related to the prediction quality of both minority and majority classes. semoos returns a pool of non-dominated solutions from which the user can choose the model that best suits him. Automatic solution selection formulas with a so-called Pareto front have also been proposed to compare state-of-the-art methods. The proposed approach will train a svm classifier ensemble dedicated to the imbalanced data classification task. The experimental evaluations carried out on a large number of benchmark datasets confirm its usefulness.


Introduction
Among imbalanced data classification methods, one of the most promising directions is using models based on classifier ensembles. In the case of ensemble learning, great emphasis is placed, on the one hand, on good prediction quality and, on the other hand, on appropriate diversification of base classifiers. Additionally, for the task in which there are unequal proportions between classes, and thus the cost of misclassification of one class (usually minority class) is higher than an incorrect object from the other class, a significant problem is the selection of a correct classification quality metric. This metric should serve as an optimization criterion when learning the model. The canonical classification model assumes knowledge of the mentioned misclassification costs, which should be provided by the user in the form of a loss function. Then a conceptually simple criterion would be the expectation value of the mentioned function, i.e., the overall risk [20].
Unfortunately, such information is lacking in real-world problems, and determining the cost of an erroneous decision is very difficult. One may suggest that a misclassification cost should be inversely proportional to the frequency of objects in a given class. However, it is easy to show that this is not necessarily true for many practical problems, and especially in the case of highly imbalanced data sets, this would lead to ignoring errors made on the majority class. Hence, approaches are constantly being developed to acquire the misclassification costs mentioned above, such as those given in the work of Branco et al. [9] utility-based learning. Most of the work on classifying imbalanced data analyzes simple metrics such as Recall / Published online: 17 November 2022 Applied Intelligence (2023) 53:15424-15441 and Precision (or Specificity). However, aggregate metrics such as Gmean, AUC, or F-measure are adopted due to the desire to express the quality assessment of a method by a single value. Their undoubted disadvantage is that they assume a fixed relationship between simple criteria, e.g., in the case of Gmean, it is the geometric mean of Precision and Recall. It is also worth noting that these criteria ignore the imbalance ratio or the misclassification costs on objects of different classes. Brzezinski et al. [10] showed that such metrics are biased towards the majority class and suggested using parametric metrics as F β , where β should be proportional to the imbalance ratio. Nevertheless, such a recommendation still does not guarantee that the best-evaluated classifier using this criterion will properly optimize the cost of misclassification because it assumes that the misclassification cost of minority class is proportional to the imbalance ratio, but it is not always true.
Additionally, it should also be noted that the values of these aggregated criteria are ambiguous. Based on their values, we de facto do not know how the given model behaves, i.e., what values the simple criteria take since a given value of the aggregated criterion can be achieved for many pairs of precision and recall values.
Considering the above, it seems interesting to consider the problem of classifier learning for imbalanced data as a multi-criteria optimization task. As a result of such algorithms, we should obtain not a single solution but a set of solutions (Pareto front), from which the selection of a single solution can either be automated or left to the user to decide.
In this paper, we intend to answer whether it is possible to propose classifier learning algorithms using multi-criteria optimization to train a set of Pareto-optimal SVM classifiers to obtain the ensemble as good as the state-of-the-art methods. The classical SVM classifier has difficulties classifying imbalanced data because determining the optimal hyperplane is not as simple as for more balanced data. The higher the imbalance ratio, the weaker the classification ability is. Therefore data sampling mechanisms (undersampling or oversampling) or changes at the classifier architecture level are used [12].
The main contributions of this paper can be summarized as follows: • Proposing the SEMOOS ensemble method using multicriteria optimization in the model learning phase. • Designing a strategy for building the ensemble consisting of Pareto optimal SVM classifiers. • Developing a new bootstrapping variant for sample subspace selection from imbalanced data. • Analysing the impact of hyperparameters setting on the SEMOOS behavior.
• Estimating the computational complexity of the SEMOOS method. • Experimental evaluation of the proposed approach based on diverse benchmark datasets and a detailed comparison with the state-of-the-art methods.
The paper is organized as follows. The following section discusses the related works. Section 3 provides the details of the SEMOOS algorithm and its run-time complexity. Then we discuss the experimental setup, report the results along with our analysis, and finally present the concluding remarks and possible future research directions.

Related works
This section presents the main works related to our research. We shortly introduce the main topics related to multi-objective optimization, imbalanced data analysis, and classifier ensemble learning.

Multi-objective optimization
Multi-criteria (multi-objective) optimization algorithms optimize more than one objective function. Hence unlike singleoptimization methods they return many solutions instead of a single one. The returned solutions may be divided into dominated and non-dominated ones [21,60]. A non-dominated (Pareto-optimal) solution cannot be improved due to any criterion without compromising the value of the others. Most multi-objective optimization algorithms are global methods [17], where approximate solutions are returned. They may be grouped into: (i) Weighted Objectives Methods, (ii) Hierarchical Optimization Methods, (iii) Trade-Off Methods, (iv) Methods of Distance Functions, Min-Max Methods, (v) Goal Programming Methods, (vi) Genetic Methods.
Among the latter methods, NSGA, NSGA-II, NPGA, and FFGA [3] should be mentioned, as well as a multiobjective evolutionary algorithm based on decomposition (MOEA/D) that decomposes an optimization problem into a number of scalar optimization subproblems and optimizes them simultaneously [79]. Nguyen et al. [54] modified the MOEA/D algorithm and proposed two decomposition methods based on multiple reference points, static and dynamic. The developed method gave satisfactory results in the feature selection in the classification. Currently, multiobjective evolutionary algorithms (MOEAs) are eagerly developed and bring significant benefits, especially in largescale multi-objective optimization problems. In this case, the decision space is grouped into several subspaces, is reduced, or novel search strategies are used [69].
Pareto solutions make it possible to reconcile several criteria but simultaneously make it difficult to choose one, the best solution, which is subjective and should depend on the user's preference. The more critical criterion is for a user, the more likely the solution with a better score for that criterion will be selected. The Pareto MTL method was developed to facilitate selecting a Pareto front solution [45]. Thanks to it, solutions in the objective space are welldistributed, which ensures a compromise between various, often opposite, criteria. Another approach that provides well-distributed solutions on the Pareto front is the COSMOS method [61] developed for deep learning.
There are many propositions on how to solve optimization. One of them is the use of the surrogate-assisted Particle Swarm Optimization with the Pareto active learning algorithm [48], which has a relatively low computational cost, fast convergence, and good diversity. Convergence is good when the solutions are close to the Pareto front, and good diversity in the objective space means that the solutions are evenly distributed in space and are not just at one point. The Ma et al. [49] proposed the LSMOEA/D method using reference vectors in the control variable analysis. The experiments confirmed high quality for the huge test problems with 2-10 objectives and 200-1000 variables.
Multi-criteria optimization is a rapidly growing domain, primarily due to the increasing complexity of modeled processes. Increasingly, such optimization is used in real problems such as crude oil price forecasting [35], efficient use of energy in agriculture [36], renewable energy of the building sector [46] and many more. There are also some applications of multi-criteria optimization in classification, e.g., Jin et al. [77] solved the problem of multi-criteria optimization of the structure and parameters of the neural network using evolutionary methods (NSGA-II). The MO-SELM method [75] has been tested for classification and regression. It relies on the optimization of parameters and structure learning of the Extreme Learning Machine network to cope with the overfitting problem. GEMONN [76] employed a gradient-guided evolutionary approach containing the advantages of gradient descent and evolutionary algorithms to train deep neural networks. Optimization is used to determine weights for the network, while multi-criteria optimization is used simultaneously for the network sparsity and the training loss.
Many researchers combine SVM classifiers with optimization algorithms that solve the feature selection problem. The article [78] used NSGA-III for the feature subset selection and CNN-SVM (Convolutional Neural Network -Support Vector Machine) for software defect prediction with an imbalance problem. Mierswa [51] showed the possibility of using multi-criteria optimization techniques in SVM learning, pointing out that thanks to this approach, it is possible to turn away from aggregated optimization criteria as a combination of opposing criteria. Additionally, Paretooptimal solutions allow complexity analysis of the solution so the user can easily see which solutions are overfitted. The combination of SVM classifier to detect malicious traffic and a Genetic Algorithm to optimize hyperparameters is used in the article [62].

Imbalanced data classification
In imbalanced data classification, the disproportion among the different classes is not the sole issue of learning difficulties. One may easily come up with an example where the instance distributions from different classes are well-separated. Napierała and Stefanowski observed that the minority class samples often may form scattered clusters of an unknown structure [53]. An additional complication arises from the possibility that there may be an insufficient number of minority class samples for a classifier learning algorithm to achieve an adequate level of generalization, resulting in overfitting [16].
One may divide imbalanced data classification algorithms into three groups [47].
Data preprocessing methods concentrate on decreasing the number of examples from the majority class (undersampling) or generating new minority class samples (oversampling). These mechanisms aim to balance the number of objects from considered classes. Oversampling randomly replicates existing samples or generates new samples in a guided manner. SMOTE is the most recognized algorithm [14] that generates new minority samples between two randomly selected objects. Unfortunately, methods such as SMOTE may lead to changes in the characteristic of the minority class and, as a result, overfit the classifier. Therefore several modifications of SMOTE have been proposed that can identify the samples to be copied more intelligently, such as BorderlineSMOTE [29] that generates new minority class samples close to the decision border. Safe-Level SMOTE [11] and LN-SMOTE [50] reduce the probability of generating synthetic minority instances in areas where the predominant objects are that of the majority class. Koziarski et al. proposed RBO [39] and CCR that enforce instances from the majority-class to be relocated from the areas where the minority-class instances are present [38]. The alternative preprocessing approach is undersampling. Such methods employ randomly removing the instances from the majority class or removing them from the areas so that the classifier's quality is not disrupted using neighborhood analysis.
Inbuilt mechanisms modify existing classification algorithms for imbalanced tasks ensuring balanced accuracy for both classes. One approach is one-class classification [34], which aims to learn the decision areas associated with one class. Initially, an approach based on building models for the majority class was proposed due to the sufficiently large number of objects representing it, and the minority class was treated as so-called outliers. A different approach was proposed in [41], where a one-class classifier was trained on minority class. In turn cost-sensitive classification considers the asymmetrical loss function that assigns a higher misclassification cost of minority class [30,40,47,80]. Unfortunately, such methods may cause a reverse bias towards the minority class.
Hybrid methods combine the advantages of methods using data preprocessing with the different classification methods. One of the most popular approaches is the hybridization of under-and oversampling with ensemble classifiers [27]. This approach allows the data to be independently processed for each base model. It is also worth noting methods based on ensemble classification [74], such as SMOTEBoost [15] and AdaBoost.NC [71].

Classifier ensembles
The purpose of the classifier ensemble is to make a joint decision by the pool of base classifiers to obtain more significant predictive power [42]. The main factors affecting the predictive quality of the ensemble are the individual predictive quality of base classifiers and their diversity [52]. These factors are the opposite. In the case of an ideal classifier, making a correct prediction in each case, increasing ensemble diversification must lead to the addition of worse quality classifiers. Therefore the ensemble forming process may be treated as a multi-objective optimization problem.
Chandra and Yao proposed DIVACE (DIVerse and Accurate Ensemble Learning Algorithm) [13], which employs multi-objective optimization to the ensemble learning task to find a trade-off between diversity and accuracy. Abbas developed Memetic Pareto Artificial Neural Network (MPANN), which optimizes similar criteria [1]. In [73], the authors optimize the weights of the models in the ensemble, and select the solution from the Pareto front using the PROMETHEE method. Other works focus on ensembles of decision trees [37], recursive networks [63], or fuzzy rules [33] which prove the possibility of using this approach. However, the number of works devoted to employing multi-objective optimization to classifier ensemble design is relatively small, especially compared to the common problem of designing the ensemble using a single criterion. Answers to key questions require further research, i.e., which models are best suited to the multiobjective task, and whether it is possible to develop algorithms that effectively combine different decision models. It is also necessary to develop methods for forming the base classifier and creating combination rules that connect them in the multi-objective optimization task.
Fletcher et al. [26] proposed a non-specific ensemble classification algorithm that uses multi-objective optimization to set the required. Ribero and Reynoso-Meza [59] developed a two-stage ensemble learning framework, where firstly, a set of diverse classifiers is generated, and secondly, a pool of classifiers is pruned. The same authors also [4] analyze different MOOD approaches for ensemble learning and propose a taxonomy of multi-objective ensemble learning. Olivera et al. [56] employed multi-objective optimization to select the base classifier's valuable features and then choose the best ensemble line-op. Onan et al. [57] developed an ensemble method that employs a static classifier selection involving majority voting error and forward search, as well as a multi-objective differential evolution algorithm. Liang et al. [44] described an ensemble learning model based on multimodal multi-objective optimization. Asadi and Roshan [5] formulated an interesting proposition that focuses on the bagging procedure, considering the number of bags and simultaneously the diversity of the trained classifiers. Gu et al. [28] focused on classifier ensemble generation, where they proposed a solution that is a trade-off between accuracy and diversity ensemble. They also showed that the proposed solution could outperform single-objective methods.

Multi-objective optimization for imbalanced data classification
One could find only a few works on imbalanced data classification employing multi-objective optimization methods. It is worth mentioning the work of Bhowan et al. [6], who proposed to build a classifier ensemble based on Pareto-optimal classifiers. In turn, Soda [64] suggested training the classification system named Reliability-based Balancing, using multi-objective optimization methods and maximizing two criteria related to the frequency of correct decisions and Gmean. They used the data preprocessing technique based simultaneously on feature selection and prototype selection obtained from multi-objective optimization. Li et al. [43] proposed a data preprocessing method Adaptive Multiobjective Swarm Crossover Optimization, which used both over-and under-sampling at the same time. This approach selected the best proportion between majority and minority samples by multi-objective optimization.
Also, several ensemble techniques have been employed for imbalanced data classification tasks. Bhowan et al. [7] proposed a two-step approach to evolving ensembles using genetic programming for imbalanced data classification. The Pareto front optimal classifiers form an initial pool of classifiers, then ensemble pruning methods based on genetic programming were employed. Fernandez et al. [25]  employed a decision tree ensemble and multi-objective optimization to find the best combination between feature and instance selections for the multi-class imbalanced task. Felicioni et al. [22] developed an algorithm that took fourth place in the ACM RecSys Challenge 2020, organized by Twitter. The challenge aimed to predict the probability of user engagement based on past interactions on the Twitter platform. Authors employ feature extraction and gradient boosting for decision trees and neural networks and build the ensemble using multi-objective optimization.

Proposition of the method
Let us propose the SVM Ensemble with Multi-Objective Optimization Selection (SEMOOS) method dedicated to training an SVM classifier ensemble for imbalanced data. Its pseudocode is presented in Algorithm 1 and in Fig. 1. 1 The main idea is to find such a pool of SVMs that gives a diversified ensemble. To achieve that, we look for the setting of two parameters C and γ of the SVM Radial Basis Function (RBF) kernel. Additionally, the feature selection for each base classifier is performed to ensure the high diversity of the ensemble. The multi-objective optimization is used to select the best SVM parameter setting, including feature selection [67]. This method depends on several parameters, but we distinguish two main versions of SEMOOS -SEMOOSb and SEMOOSbp. SEMOOSb employs an original bootstrapping method to increase the diversity of the ensemble (Bootstrapping is true). SEMOOSbp employs also pruning to remove similar models from the ensemble (Pruning is true).
Bootstrapping generates data subspaces. Optimization based on this data results in differentiated SVM parameters from a given range and determines which features will be selected. The best non-dominated solutions were used to create a model, and then such a model was added to the ensemble of classifiers. This iteration of center-based Bootstrapping was repeated several times to ensure stability, higher quality, as well as avoid overfitting.

Algorithm
Let us quickly analyze Algorithm 1, which starts with the learning set LS as an input. Suppose the center-based Bootstrapping mechanism is enabled (Bootstrapping = T rue). In that case, it is repeated n-times, and it divides each fold of the dataset into subsets S.  (1) c + x (1) kx (2) c + x (2) . . .
For example, in two-dimensional space, after the second iteration, there are two roots, and the center of them is in the middle of the line that connects these roots. After the third iteration, three roots are selected creating a triangle, so the new root is a centroid of this figure. This process continues until the iterator reaches the value of the nrepetition parameter.
In each iteration, the subset S i is used as the input to the multi-objective optimization fulfilled by the NSGA-II algorithm [18]. Equation 2 presents two fitness functions, F 1 and F 2 that the optimization algorithm uses to maximize Precision and Recall The metrics are calculated during the validation process inside the optimization with the base estimator SVM and different values of its hyperparameters C and γ that have various lower and upper limits in the set of real numbers: C ∈ [1e 6, 1e 9] and γ ∈ [1e-7, 1e-4].x is a binary vector containing selected features. These three parameters (C, γ ,x) form an initial population. The optimization is repeated until the maximum number of evaluations (m) is reached. It returns the Pareto optimal set containing results from the Pareto front, such as C, γ , and selected features.
Then, the estimator is trained and added to the ensemble for each result. The version SEMOOSbp with the Pruning does not add all models to the ensemble. It finds unique values of fitness functions (Precision and Recall) and, based on these indexes of solutions, trains the estimator on (C, γ ,x). Finally, the algorithm returns the ensemble of classifiers. The prediction is based on Average Support Vectors.
The basic version of SEMOOS is much simpler. It starts from the multi-objective optimization described above but the input is LS. Subsequently solutions from the Pareto Set P S such as (C, γ ,x) are used to train the model, and all models are added to the ensemble.

Computational complexity analysis
The computational complexity depends on a few aspects of the proposed method. Firstly, the time complexity of the base classifier SVM is O(N 3 ) [2], where N denotes the size of the dataset, so it is the number of examples in the learning set. Let us assume that M is the number of objectives, and K is the population size. The computational complexity of NSGA-II is O(MK 2 ) [18]. In our case, M = 2, the complexity is O(2K 2 ). The last part of the method is Bootstrapping, and its complexity is given by O(N)

Experimental evaluation
The experiments described in this section are used to test the proposed methods and answer the research questions posed below.
RQ1: What is the impact of the SEMOOS's parameters (especially Bootstrapping and Pruning) on its quality? RQ2: How do variants of the SEMOOS affect classification quality? RQ3: Can SEMOOS methods outperform state-of-the-art algorithms? RQ4: What is diversity of SEMOOS ensemble compared to the reference methods?
As the proposed method aims to train an ensemble of SVMs, the Random Subspace method is also based on SVM, which randomly performs feature selection for base classifiers. By choosing the right features for each base classifier, diversity is assured. The other reference algorithms are not using the ensemble classifier paradigm, i.e., simple SVM classifiers without feature selection and two SVM models based on feature selection. A description of the SEMOOS algorithm's variants and the reference methods used in the experiments are presented below.
• SEMOOS -SVM Ensemble with Multi Objective Optimization Selection Three versions of the proposed method (SEMOOS, SEMOOSb, SEMOOSbp) will be compared with the following benchmark solutions: SVM, an ensemble RS, and two classifiers with feature selection (FS and FSIRSVM). FS is a classical approach to Feature Selection. It is based on the Chi-square statistic [68] and the K-best function, which chooses K-best features. FSIRSVM is almost the same as FS, but it has an additional parameter, IR -Imbalance Ratio of each fold, which is applied to the SVM as the class weight parameter. FS and FSIRSVM select 75% of features from each dataset. RS creates random subspaces, and it has 100 models in the ensemble. The Support Vector Machines (SVM) classifier was used as a base classifier for all methods except all variants of SEMOOS. SVM has default parameters such as the regularization parameter C = 1, and the kernel is RBF, kernel coefficient γ scaled, probability estimates are set to True, and the rest parameters had default values from the scikit-learn.
As mentioned in Section 3, all SEMOOS variants use SVM and optimize its parameters C and γ . The rest parameters of SVM are default. Our methods use the optimization algorithm NSGA-II (Non-dominated Sorting Genetic Algorithm) [18] by 100 populations with the diverse representation, because C and γ are real values, andx is a binary vector of selected features. After pre-experiments genetic operators were used, such as Random Sampling, Polynomial Mutation (for the real representation) and Bitflip Mutation (for the binary representation), Simulated Binary Crossover (for the real representation), and Two Point Crossover (for the binary representation). An eta parameter of Simulated Binary Crossover and Polynomial Mutation was set to 5. A constraint's representationx is 75% meaning that 75% of features are selected and their value in the vector is 1. The process of optimizing parameters is done by 1000 evaluations. The size of the ensemble depends on the variant of our methods. The parameter determines the ensemble size in the SEMOOS methods, which is ten models. In SEMOOSb, this value is multiplied by a number of iterations (five iterations) of the Bootstrapping, and the final ensemble consists of 50 classifiers. Due to the Pruning in SEMOOSbp, we cannot specify its size in advance.
As the experimental protocol, the Repeated Stratified K-Fold (5 repeats x 2 splits) cross-validation was chosen [66]. Such cross-validation was also applied inside the optimization to avoid overfitting. Results were saved as csv files for each dataset fold and metric. Then, we performed Wilcoxon statistical rank-sum tests to see if one method was statistically significantly better than the other in pairwise rankings. Five metrics were used to measure the quality of methods: BAC -Balanced Accuracy, Gmean, Recall, and Precision.
All imbalanced datasets used in the experiments are listed in Table 1, where I D -is the dataset identifier, Dataset -is the name of the dataset, I R -is the Imbalance Ratio, Ex. -is the number of examples, Attr. -is the number of attributes. They were loaded from the KEEL dataset repository [70], where these datasets were divided according to IR. The first 19 datasets are separated with the line and have IR lower than 9. The rest of the datasets has IR higher than 9. We keep up this division in our results to check the effectiveness of our methods on data with low and high imbalance. Most of the datasets are multi-class problems but they have already been prepared as binary classification problems. For example, the glass-0-1-2-3 vs 4-5-6 dataset combines 0, 1, 2, 3 classes of the original dataset as a negative class (majority) and 4, 5, 6 as a positive class (minority). The extended description with class names is characterized by Fernandez et al. [23,24].

Experiments
We carried out four groups of experiments to answer the research questions:

Experiment 1: Setting hyperparameters
Our methods have a few parameters, so we decided to conduct pre-experiments to choose values of these parameters that get the best quality. We tested SEMOOSb on four datasets (ecoli-0-3-4-7 vs 5-6, glass5, vehicle3, yeast-2 vs 4) and averaged results to present figures.  glass-0-6 vs 5 11.00 108 9 40 ecoli-0-  Firstly, we check the eta parameter of crossover and mutation for each metric for values [2,5,10,20]. Deb et al. in [19] point out that small eta values for crossover provide a diverse search among solutions.  Next, the number of iterations of bootstraps and the percent of selected features were tested, and ecoli-0-3-4-7 vs 5-6 dataset is shown in Fig. 3. This situation is similar to the previous one, BAC, Gmean, and Recall indicate the highest metrics for bootstrap = 5 and f eatures = 75%.

Experiment 2: Comparison of three variants of SEMOOS
Based on the parameters selected in the previous experiment, tests are conducted comparing the proposed SEMOOS, SEMOOSb, and SEMOOSbp methods. The results from the folds and all datasets are averaged, and Wilcoxon rank-sum statistical tests are performed. In the figure of the Wilcoxon test, rows are different metrics (BAC, Gmean, Recall, Precision), and columns are labeled with methods. The green color means that the method wins, yellow -ties, and red -loses to the method located at the bar on the chart. The black dashed line indicates the statistical significance of the method as winning. Figure 4 shows a test for datasets with the imbalance ratio of less than 9. The results do not indicate any of the methods as statistically significantly better than the others. However, it can be noticed that SEMOOS and SEMOOSbp win with the SEMOOSb method, especially in the first three metrics, i.e. BAC, Gmean, and Recall. In Fig. 5, there are 59 datasets where I R > 9. Similar conclusions can be drawn from it that SEMOOS is better than SEMOOSb, but the other methods are at a similar level.

Experiment 3: Comparison with reference methods
A vital element of each method is to compare it with state-of-the-art methods to check whether the proposed new method is statistically significantly better than the others. This experiment will compare three proposed variants of the SEMOOS method and the reference methods RS, SVM, FS, FSIRSVM. In Appendix A is a link to tables with the exact results of the metrics averaged over the folds together with the standard deviation for each dataset. The best result for a given dataset is marked in bold. Each table contains a different metric. The results of the proposed methods are compared with the reference ones. Sometimes the quality is slightly higher or lower, but it is difficult to pinpoint a winning method.
Therefore, Wilcoxon statistical tests showed how many datasets each proposed method won. As in the previous section, the figures are grouped by imbalance ratio above or below 9. Each SEMOOS variant is compared to all reference methods for the five metrics. All SEMOOS variants work similarly in Fig. 6, and they are statistically significantly better than FSIRSVM for BAC, Gmean, and Recall metrics. They also achieved many wins compared to the RS method. However, only SEMOOS for the Recall metric shows a statistically significant win.
The results in Fig. 7 are slightly different for data with high imbalance. There are many more datasets in this case. All variants of the SEMOOS method win with static significance with the RS method for BAC, Gmean, and Recall metrics. There are also fewer draws between the methods for these three metrics. An important goal for us was to correct identify the minority class that is represented by Recall.
Analysis of the Pareto front in the objective functions space is presented in Fig. 8. Triangles represent the reference methods. Our method proposals fall into two Fig. 4 Wilcoxon rank-sum test of proposed methods for datasets with I R < 9 (green -win, yellow -tie, red -loss) 1 3

Fig. 5
Wilcoxon rank-sum test of proposed methods for datasets with I R > 9 (green -win, yellow -tie, red -loss) categories, the method name with PF (Pareto Front), and the method name alone is the quality of the final constructed ensemble. These scatter plots show how the Pareto front solutions locate relative to the final built ensembles and the reference methods. Figure 8a shows an exemplary broad Pareto front, which provides more diverse models. The reference methods are behind the Pareto front; when comparing them to ensembles, they obtain a lower Precision value. Figure 8b shows more focused PF solutions than Fig. 8a. SEMOOS variants obtain a similar value for both metrics and win over the reference methods. In these examples, the methods score is both high Precision and Recall, but this is not a rule for all datasets.

Experiment 4: Evaluating SEMOOS ensemble diversity
The last experiment focuses on the evaluation of the diversity of models in the ensemble. We conducted tests for four diversity metrics proposed in [42] for three SEMOOS variants and one reference ensemble RS. Figure 9 shows the Fig. 6 Wilcoxon rank-sum test for datasets with I R < 9 (green -win, yellow -tie, red -loss) SVM ensemble training for imbalanced... 15435

Fig. 7
Wilcoxon rank-sum test for datasets with I R > 9 (green -win, yellow -tie, red -loss) Q-statistic metric, while all other metrics are available on the remote GitHub repository. 5 , 6 Q-statistic is a pairwise diversity measure. It can assess two classifier outputs and return the decision on their similarity. Q is in the range from −1 to 1, where Q = 0 means that classifiers are statistically independent; Q < 0 -classifiers make mistakes on different objects; Q > 0classifiers correctly recognize the same objects.
From Fig. 9a, it may be concluded that for data with a small imbalance, the method diversity hardly differs between each other. However, for a greater imbalance (Fig. 9b), these differences are significant, and the result of the SEMOOS method is closest to the value 0, which means that the models in this method are the most diverse, compared to RS the difference is 0.3.

Lessons learned
Let us try to answer the asked research questions considering the obtained results of the experiments. None of the proposed methods shows a statistically significant difference compared to the others. Therefore, we included all three variants in further experiments. However, based on the presented statistical test results, SEMOOS and SEMOOSbp are better than the SEMOOSb for all datasets. RQ3: Can SEMOOS methods outperform state-of-theart algorithms? Each variant of the SEMOOS method has been compared with four state-of-the-art methods, and the statistical test results for each SEMOOS variant are similar if we consider the number of wins. However, depending on the level of dataset imbalance, the results vary. For datasets with the imbalance ratio below 9, the proposed methods exceed the FSIRSVM method with statistical significance and have a high number of wins with the RS method. However, for the imbalance ratio above 9, statistical significance is obtained for the RS algorithm. RQ4: What is diversity of SEMOOS ensemble compared to the reference methods?
The last element of the research shows that all proposed methods are more diversified than the Random Subspace (RS) method. When analyzing the figure with a greater imbalance ratio with a larger number of datasets, it is noticeable that the SEMOOS method works best.

Conclusion
The main goal of this work was to use multi-objective optimization (NSGA-II) in the form of two independent fitness functions, Precision and Recall, to classify imbalanced data. Three ensemble classifiers using SVM as a base model were proposed: SEMOOS, SEMOOSb, and SEMOOSbp. SEMOOS is a basic version that adds all SVM models obtained from optimization. SEMOOSb has the additional option of Bootstrapping, i.e., resampling a dataset with replacement is performed to get more samples from one dataset. The last version is SEMOOSbp, and the method includes both the Bootstrapping and Pruning needed to remove redundant models from the ensemble.
The experiments were performed according to the experimental protocol on 78 imbalanced datasets. First, SEMOOS hyperparameters were set to select the parameter values that produce the best results for this kind of data. Then, the three versions of the SEMOOS method were compared with each other, but none of the methods showed statistical significance. Therefore all versions were selected for further research. The main experiment was to compare the proposed methods with various state-ofthe-art methods, such as single SVM classifier, ensemble classifier (RS), and two methods using feature selection (FS and FSIRSVM). Statistical tests showed that the SEMOOS variants outperform the RS and FSIRSVM methods with statistical significance. The last element of the work was to compare the diversity of models in ensemble methods. The presented results show that the proposed methods are more diversified than the SOTA solutions.
Our methods could also be tested for multi-class problems because data binarization is only a simplification of this problem. Due to feature selection in the proposed methods and potentially good results for datasets with more features, it would be good to check more such datasets. In case of too little real data, it is worth considering artificially generated data, which will also ensure controlled conditions for testing methods. The comparison could also be made with other reference methods that would have more elements in common with the SEMOOS, SEMOOSb or SEMOOSbp methods. The last of future works is to consider the optimization itself, that is to use other multi-objective optimization methods, e.g. MOEA/D, other metrics in fitness functions, or to increase their number. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.