<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>The Fifth International Workshop on Computer Modeling and Intelligent Systems, Zaporizhzhia, Ukraine, May</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Feature  Engineering  and  Missing  Data  Imputation  Method  of  Medical Data Analysis </article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Nataliya Shakhovska</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Nataliia Melnykova</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Lviv Polytechnic National University</institution>
          ,
          <addr-line>S.Bandera str,12, Lviv, 79013</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2022</year>
      </pub-date>
      <volume>12</volume>
      <issue>2022</issue>
      <fpage>0000</fpage>
      <lpage>0002</lpage>
      <abstract>
        <p>   This work provides an alternative way to preprocessing procedure for consolidated data. Two methods are proposed. The first one is used for feature selection based on ensemble of machine learning algorithms. And the second one organizes missing data imputation based on combination of functional dependencies and associative rules. Ensemble methods for processing multimodal data based on a hierarchical classifier, a set of weak classifiers and a number of methods for selecting important characteristics with a much higher value of accuracy on unbalanced data sets compared to existing machine learning methods are developed. The methods are validated on medical dataset. The percentage of recovery data is on 1.2% comparing with associative rules. The proposed missing data imputation method creates additional data values operating a based domain and functional dependencies and includes these values to available training data. The correctness of the filled-in values is proved on the predictor built on the original dataset. The proposed PPD method conducts 12% better than RF and EM models for 30% missing data.</p>
      </abstract>
      <kwd-group>
        <kwd> 1  feature selection</kwd>
        <kwd>missing data</kwd>
        <kwd>machine learning</kwd>
        <kwd>ensemble</kwd>
        <kwd>data preprocessing</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction </title>
      <p>are not used inefficiently. Thus, there is a problem of developing an effective strategy of data analysis
that can be applied to dispersed databases of various domains.</p>
      <p>Data analysis methods are effectively used on clear and previously processed data. That is why the
paper aimed to analyses and develop methods for data preprocessing. Particularly, we have taken into
account feature selection (feature engineering) and missing data imputation.</p>
      <p>The main contribution of the paper is as follows:
1. A new hybrid ensemble feature selection model for a machine learning-based
post-COVID19 prediction system is proposed as an automatic feature cut-off rank identifier.</p>
      <p>2. A method based on probabilistic dependence is developed and tested. The percentage of
recovery data is on 1.2% comparing with associative rules.</p>
      <p>3. Development of a method for finding dependencies in large data sets with gaps and
uncertainties based on an ensemble of clustering methods and auto-associative dependencies. Boruta,
Decision Tree and Random Forest are used to select an object. The importance of variables is
different for different methods (logistic regression, Support vector machine, Naive Bayes, XGBoost,
Random Forest, neural network, decision tree). This means that the relationship between the
parameters is supported only for part of the data set. That is why we propose to find the dependency
for separate clusters and use this dependence for classification.</p>
    </sec>
    <sec id="sec-2">
      <title>2. Literature Review </title>
      <p>
        Selecting the appropriated features can be a more significant task than lessening computation time
or enhancing classification or predictive accuracy. For example, in medicine [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ], finding the optimal
set of optimal features for the classification or predictive problem can help develop a diagnostic test.
      </p>
      <p>Selecting important features (for example, determining genes appropriate for a specific type of
cancer) can help decipher the mechanisms underlying the welfare problem for research. A complete
enumeration of features can implement the selection method, that is, having checked all possible sets,
selecting those signs for which the error is minimal. This method is simple to implement, but it is
entirely ineffective on big data. Therefore, in this case, other algorithms are most often used.</p>
      <p>
        The three primary classes of feature extract algorithms – filters, wrappers, and built-in algorithms
are used [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ].
      </p>
      <p>
        Filters are based on some metrics that are independent of the classification method. For example,
the correlation of features with the target vector and information content criteria. They are applied
before classification. The most significant benefit of filtering is that it can be used as preprocessing to
reduce space dimensionality and overcome overfitting. Filtering methods are generally fast. Filters are
used to choose features in clustering or to build an initial approximation [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. Unfortunately, such
methods are not designed to detect complex relationships between elements and, as a rule, are not
sensitive enough to specify all dependences in the data.
      </p>
      <p>
        Embedded algorithms organize feature extraction during the classifier training approach, and it is
they explicitly optimize the set of features used to achieve better accuracy [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]. The benefits of built-in
algorithms are that, as a rule, they discover solutions quickly, decrease the possibility of retraining
data, eliminate the need to separate data into training and test subsamples. Nevertheless, these
algorithms are not ubiquitous.
      </p>
      <p>
        Wrappers rely on feature significance information from several classifications or regression
models and can thus find deeper patterns in dataset than filters. Wrappers can be built on any
classifier that defines the degree of significance of the features [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ].
      </p>
      <p>Imputation is a procedure for assessing unknown or missing values based on available data, which
allows you to form a complete set of data with some plausible estimates.</p>
      <p>
        Methods of imputation are divided into non-model and model-based approaches. There are
approaches based on single and multiple imputation methods in terms of the quantity of values
received from imputation methods [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. One-time filling algorithms provide a single complete set of
data, where a deal replaces each space. The benefits of this method are the use of methods of
comprehensive data analysis in the subsequent stages of processing. Replenishment algorithms form
several complete data sets analyzed separately and later combined according to specific rules. This
minimizes standard errors in the following steps by processing comprehensive datasets. Nevertheless,
mentioned methods require multiple resources to create more data sets, spend more time performing
the analysis, and have more memory to store the results [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ].
      </p>
      <p>
        The method based on mean substitution solves the problem of incomplete data by replacing each
missing variable with an average value. There are the following types of substitutions: median value,
mean value for subgroup [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ], value with the highest frequency, and replacement with minimum /
maximum value. This method can lead to undesirable results [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ], such as variance change, negative
correlation shift, and misrepresentation of the population.
      </p>
      <p>
        Hot-deck imputation approach replaces each gap with a random value taken from an existing
dataset [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ]. Its significant disadvantage is the warping of correlations and covariates.
      </p>
      <p>
        Cold-deck imputation (CD) approach implements the replacement of each gap by some constant
value from an external source [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ]. The specific case of CD is zero replacement. It has the same
disadvantages as hot-deck imputation.
      </p>
      <p>The regression model applies to replace data gaps with expected values emanated from a
regression equation constructed from a complete dataset. Disadvantages of regression completion
include the need to accurately define regression models, exaggerate correlation and covariance, the
likelihood of moving to predicted values outside the logical sequence, and the need for extensive
quantities of data to obtain consistent estimates.</p>
      <p>
        The association rules (AR) mining method uses the constructed associative rules for the data
imputation [
        <xref ref-type="bibr" rid="ref10 ref11">10, 11</xref>
        ]. However, the temporal complexity of this method needs to be improved [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ].
      </p>
      <p>Diabetes increases the likelihood of severe COVID19. New clinical data and investigations show
that this might work in the opposite direction: scientists are recording new cases where COVID-19
has sharply provoked type 1 diabetes in humans. The World Health Organization views diabetes as
one of the existing diseases, on a par with respectable age, making someone more vulnerable to severe
COVID 19 infection. Cellular immunity is essential for protecting against viral diseases, and its
effectiveness decreases with age. In particular, a decrease in T-cell receptors explains the significant
increase in mortality of COVID19 with age.</p>
      <p>Therefore, another important factor is the search for possible relationships between biomarkers of
aging and COVID resistance.</p>
      <p>
        Existing research and the data sets collected for this use standard machine learning methods, while
demonstrating not very high accuracy of prediction. Thus, in the paper [
        <xref ref-type="bibr" rid="ref13">13</xref>
        ], there are used feature
selection, XGBoost and decision tree to determine COVID biomarkers, F1-score does not rise above
0.7. In the paper [
        <xref ref-type="bibr" rid="ref14">14</xref>
        ], markers of CH4 and CH8 immunodeficiency and their association with
coronavirus infections were analyzed using statistical models, including the Cox model. Accordingly,
it is impossible to prevent negative situations. In paper [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ], there is used the empirical mode
decomposition (EEMD) and the artificial neural network (ANN) to predict the COVID-19 epidemic.
Thus, in order to prolong the active period of life, it is necessary to track the dynamics of changes in
molecular and biochemical markers, anthropometric indicators, behavioral factors, environmental
parameters and habitat, and so on. As a result, it is necessary to use a big data-based approach to
collect information from disparate datasets, process them, and further analyze them. At the same time,
it is necessary to analyze small data samples, which will include time multimodal series of changes in
human parameters.
      </p>
      <p>The analysis of literature sources showed the lack of a comprehensive approach to solving the
problem of prolonging the active period of life and preventing exacerbation of chronic diseases. At
the same time, as we see in the case of COVID, chronic diseases (diabetes and obesity) can not only
reduce the ability to work, but also increase the likelihood of severe course of other diseases.</p>
    </sec>
    <sec id="sec-3">
      <title>3. Materials and Methods </title>
      <p>Biomarkers of aging are used to predict possible changes in the body that lead to disability due to
functional age-related changes. Biomarkers of aging are markers that can predict the functional
capacity of an organism at a certain age better than chronological age. Since the immune system is a
change the language here for English version</p>
      <p>leading factor in aging, the main impact of which is realized through increased inflammation and
reduced effectiveness of cellular immunity, it becomes clear the need to involve relevant markers to
develop interventions to increase the duration of healthy longevity.</p>
      <p>Thus, control of chronic age-related diseases (diabetes and obesity) and biological markers can be
used to predict functional changes in the body, and analysis of other personal indicators will
determine how to reduce the negative effects of such changes by extending the period of active
longevity.</p>
      <p>Sociological research also shows that people in certain regions remain active for a long time and
there are far fewer people who are obese. Therefore, it is also advisable to analyze the parameters of
the environment and habitat and its impact on the parameters of the organism.</p>
      <p>The baseline of the hybrid ensemble feature selection model looks like the following:
 Several selectors using,
 Aggregation of the results.</p>
      <p>Several wrapper algorithms will be used in the preprocessing stage for the first stage.</p>
      <p>The correlation matrix shows the numerical value of the correlation coefficient for all possible
combinations of variables. It is mainly used to find out the relationship between more than two
variables. The decision tree returns the feature weight as the criterion for evaluating features. It allows
building a ranked list of selected features using different measures. Classification and Regression
Trees (CART) was used for feature selection with Gini-index as a measure in our case.</p>
      <p>Random Forest is an ensemble of numerous training-sensitive algorithms (decision trees).
Mentioned approach has a slight compensation. The bias of the training method is the deviation of the
average response of the trained algorithm from the reaction of the ideal algorithm. Each of these
classifiers is built on a random subset of entities and a random subset of characters.</p>
      <p>
        Boruta is a heuristic algorithm for choosing important features based on Random Forest approach
[
        <xref ref-type="bibr" rid="ref16">16</xref>
        ]. Features with the Z-measure less than the maximum Z-measure among the added features are
removed at each iteration. To calculate the Z-measure, we need to calculate its importance, obtained
using the built-in algorithm in Random Forest, and divide it by the standard deviation of the feature
importance. Added features are obtained as follows: the characteristics available in the selection are
copied, and then every new attribute is filled by mixing its values. This procedure is repeated several
times to get statistically significant results, and variables are generated independently at each
iteration.
      </p>
      <p>
        The Jaccard index [
        <xref ref-type="bibr" rid="ref17">17</xref>
        ] measures the similarity of the feature subsets chosen by separated feature
selectors (each selector is organized as a separated iteration):
 , … , 
| ∩…∩ |,
| ∪…∪ |
(1)
where Si is the subset of features at the i-th iteration, for i=1,..,n. The value of the Jaccard index varies
from 0 to 1, where 1 implies the absolute similarity of subsets.
      </p>
      <p>The schema of the hybrid ensemble feature selection model is given in Fig. 1.</p>
      <p>Next, the method of missing data imputation is developed. This method is based on classical
functional dependencies in relational databases and association rules from non-relational databases. It
consists of two parts:
 Probabilistic Production Dependencies mining;
 the Probabilistic Production Dependencies usage for missing data imputation.</p>
      <p>
        Investigating extensive data requires specifying attribute clusters that form functional
dependencies (FD) [
        <xref ref-type="bibr" rid="ref18">18</xref>
        ]. However, datasets are extensively standard in the real world, with essential
dependencies defined only on a subset of crucial attribute group values. Moreover, the reliance can
appear not only for tuples in relational data sources but also between subsets of values in different
tuples. We will name them Probabilistic Production Dependency (PPD).
      </p>
      <p>Figure 1: The hybrid ensemble feature selection model. 
 </p>
      <p>
        Probabilistic Production Dependency is a dependence similar to associative rule in the primary
ratio selection that is proper for many entities. The significance threshold should be specified expertly
or based on estimation of the probability of erroneous selection of this dependence. The main
distinction between associative rules and PPD is that PPD will be generated from existing functional
dependencies (FD) in the dataset [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ].
      </p>
      <p>FI : K  ai, ai  A, D  aj,
aj  A,:P k  K  d  D  p
where k and d are the tuples of groups of attributes K and D , respectively.</p>
      <p>The gaps and missing data presented among the values of the attribute Y of the relation r are
classified using PPD. The following algorithm for PPD mining is proposed.</p>
      <p>Algorithm 1. PPD mining algorithm</p>
      <sec id="sec-3-1">
        <title>1. Entities with the same X-values will be grouped;</title>
        <sec id="sec-3-1-1">
          <title>2. To choose attributes from FD with same X and add them to Y;</title>
          <p>3. To calculate the Support and Confidence, Imputation of the
tuple selected in step 2);</p>
        </sec>
        <sec id="sec-3-1-2">
          <title>4. To identify the tuples with the highest value of Confidence;</title>
        </sec>
        <sec id="sec-3-1-3">
          <title>5. To add X→Y to PPDset.</title>
          <p>To fill in missing data, the PPD should be built.</p>
          <p>Next, the novel algorithm for missing data imputation is developed.</p>
        </sec>
      </sec>
      <sec id="sec-3-2">
        <title>Algorithm 2. Data imputation algorithm</title>
        <p>Completeness=0</p>
        <sec id="sec-3-2-1">
          <title>While Completeness/100&lt; Imputation</title>
        </sec>
        <sec id="sec-3-2-2">
          <title>Arrange all attributes from PPDset by Confidence level</title>
        </sec>
      </sec>
      <sec id="sec-3-3">
        <title>For each group:</title>
        <sec id="sec-3-3-1">
          <title>If percentage of non-empty Y-value is higher or equal to Support</title>
          <p>fill in empty values using PPDset</p>
        </sec>
      </sec>
      <sec id="sec-3-4">
        <title>Completeness++</title>
      </sec>
      <sec id="sec-3-5">
        <title>Else</title>
      </sec>
      <sec id="sec-3-6">
        <title>Merge PPD using Armstrong rules</title>
        <p>
          Next, a hierarchical classifier as a two-stage data prediction algorithm is developed. The first stage
is clustering; the next step is to build a classification model for each separate cluster. K-means [
          <xref ref-type="bibr" rid="ref19">19</xref>
          ]
together with the random forest do not dominate other cluster models. The hierarchical classifier is
constructed as follows:
 the appropriate number of clusters was found using gap statistics;
 the density of distribution is calculated;
 XGboost and Random Forest are used for each cluster separately;
 hard voting for the results. Based on it, the class with the highest number of votes will be
selected. If the voices are the same, the result of the classifier with the minimum value of
depth will be selected.
        </p>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Experiments and results </title>
      <p>To validate the proposed methods, dataset with medical data is used.</p>
      <p>Dataset consists of 35 features and 122 instances collected from Lviv regional rehabilitation center
for post-COVID patients with short- and long-term (more than 20 days) treatment and rehabilitation.</p>
      <p>The personal data were removed from the dataset and replaced with unique random identifiers.
The next feature, sex, is processed using one-hot encoding technics and in the final dataset is given in
two components – female and male. Features like age, weight, height, BMI, CAT, pulse, the function
of external respiration are taken as physiological parameters measured before inpatient treatment. The
rest of the features were immune-based biomarkers as described below.</p>
      <p>Zero cells (0-lymphocytes) do not carry T- and B-cells markers. Zero cells make 10–20 % of the
total lymphocytes in human peripheral blood. Some researchers consider them immature or overripe
T- or B-lymphocytes because they have a small number of antigens common to B- and T-cells. Zero
cells include K-cells and NK-cells.</p>
      <p>CD3+ is a surface marker specific to all T-lymphocyte subpopulation cells. By function, it belongs
to the family of proteins that form a complex of membrane signaling associated with the T-cell
receptor. Mature T-lymphocytes are "responsible" for cellular immune reactions and conduct
immunological monitoring of antigenic homeostasis in the body.</p>
      <p>
        CD4+ is a characteristic of helper T-cells; also represented on monocytes, macrophages, dendritic
cells. It binds to class II MHC molecules expressed on antigen-presenting cells, facilitating the
recognition of peptide antigens. Helper T-lymphocytes (CD4+) are helpers (inducers) of the immune
response, cells that regulate the strength of the body's immune response to a foreign antigen, as well
as control the stability of the body's internal environment (antigenic homeostasis) and cause increased
antibody synthesis [
        <xref ref-type="bibr" rid="ref20">20</xref>
        ].
      </p>
      <p>CD8+ is a characteristic of suppressor and cytotoxic T-cells, NK-cells, mostly thymocytes. It is a
T-cell activation receptor that facilitates the recognition of cell-bound class I MHC antigens.</p>
      <p>CD16+ natural killers are part of innate immunity; they are involved in early response against viral
infections and intracellular bacteria. Compared with cells of specific immunity (T- and
Blymphocytes), they have the advantage that they do not require long-term activation. Besides,
NKcells complement the action of T-cytotoxic cells can also regulate the immune response by producing
various cytokines, including interferon-γ. They are the primary cells of antitumor protection. Their
role is vital in manifesting cellular immunity in viral, protozoan, fungal and bacterial diseases caused
by intracellular parasites. Their action is enhanced by interferon. The functions performed by natural
killers can be divided into two main types: the production of cytokines that regulate the work of other
cells of the immune system and the direct destruction of damaged cells.</p>
      <p>Mature B lymphocytes express CD22+ markers. B-lymphocytes are responsible for the humoral
adaptive immune response, primarily at removing extracellular infectious agents. After binding to a
specific antigen, B-lymphocytes, in cooperation with T-lymphocytes and T-helpers proliferate,
differentiate into plasma cells that secrete antibodies/immunoglobulins and memory cells. Defects of
humoral immunity associated with the B-cells are sporadic, so a common hypoimmunoglobulinemia
is mainly caused by other reasons.</p>
      <p>
        CD4/CD8 immunoregulatory index reflects the ratio of CD4+ cells (T-helpers) to CD8+ cells
(Tcytotoxic cells). It is a relative indicator that has an indicative value. Its small increase or decrease has
no independent diagnostic value. Changes in the index significance the clinician to focus on the
reasons for the deviation of this index. The immunoregulatory index is assessed relative to the phase
of the immune response. In the period of exacerbation and remission of clinical manifestations, the
immunoregulatory index reaches high values due to the high percentage of T-helpers (CD4+ T-cells).
During the recovery period, the indicator's value decreases due to the increase in the level of CD8+
Tcells (killers). Violation of this pattern indicates the inadequacy of the immune response and the
possibility of chronic infection due to incomplete removal of the pathogen [
        <xref ref-type="bibr" rid="ref21">21</xref>
        ].
      </p>
      <p>We implemented our approach in Rstudio. The essential packages we used were caret, rpart,
Metrics, Boruta, Random forest, rules and ggplot2 for visualization. To generate PPD, the minimal
support threshold equal to 0.0001 in the a priori algorithm is chosen. In addition, all rules with
confidence below 0.001 are filtered out.</p>
      <p>First, the predictive accuracy for whole dataset and selected features was analyzed (Table 1).</p>
      <p>The much higher accuracy is obtained with SVM (Support vector machine, polynomial kernel) and
artificial neural network (one hidden layer with 12 neurons in it, sigmoid activation function) for
selected features [22]. The used measures are the following [23]:
 mean squared error MSE,
 mean absolute error MAE,
 R2.</p>
      <p>Next, missing data imputation method is used. The developed method was compared with the
existing ones: associative rules (AR), random forest (RF), support vector machine (SVM),
multilayered perceptron (MLP), expectation-maximization (EM) and k-nearest neighbor (KNN) (Fig.
2). The recovery error is presented using normalized root-mean-square error (NRMSE).</p>
      <p>The proposed missing data imputation method creates additional data values operating a based
domain and functional dependencies and includes these values to available training data. The
correctness of the filled-in values is proved on the predictor built on the original dataset. The
proposed PPD method conducts 12% better than RF and EM models for 30% missing data.</p>
    </sec>
    <sec id="sec-5">
      <title>5. Conclusion </title>
      <p>Current trends in the development of information technology and databases are analyzed. As a
result, unsolved problems in the field of dependency search in large databases of several subject areas,
in particular, in medicine, were revealed. The analysis of existing methods and means of detection of
dependences in data is carried out. This made it possible to identify a new subclass of dependencies
Probabilistic Production dependencies. A method for deriving PPD in relational databases has been
developed.</p>
      <p>A new hybrid ensemble feature selection model for a machine learning-based post-COVID
prediction system is proposed as an automatic feature cut-off rank identifier. Ensemble methods for
processing of multimodal data based on a hierarchical classifier, a set of weak classifiers and a
number of methods for selecting important characteristics with a much higher value of accuracy on
unbalanced data sets compared to existing machine learning methods are developed.</p>
      <p>The associative rules are found together with weak predictors usage to improve the classification
quality.</p>
      <p>The proposed missing data imputation method creates additional data values operating a based
domain and functional dependencies and includes these values to available training data.</p>
      <p>The correctness of the filled-in values is proved on the predictor built on the original dataset. The
proposed PPD method conducts 12% better than RF and EM models for 30% missing data. The EM
method looks the best for more additional missing data (the range of about 40%-50% missing data),
and the PPD has equivalent results with the SVM (support vector machine).</p>
      <p>Comparison of approaches for investigating and modeling the statistical processes by qualitative
criteria verified the proposed method has the subsequent benefits:
• retaining the characteristics of resistance to errors in the data;
• allowing the parallel implementation in distributed databases;
• automating and performing the analysis of various data types.</p>
    </sec>
    <sec id="sec-6">
      <title>6. References </title>
      <p>[22] D.A. Otchere, et al., Application of supervised machine learning paradigms in the prediction of
petroleum reservoir properties: Comparative analysis of ANN and SVM models, Journal of
Petroleum Science and Engineering 200 (2021): 108182.
[23] D. Chicco, M. J. Warrens, G. Jurman, The coefficient of determination R-squared is more
informative than SMAPE, MAE, MAPE, MSE and RMSE in regression analysis evaluation.
PeerJ Computer Science, 7, e623 (2021).</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>M.B.</given-names>
            <surname>Kursa</surname>
          </string-name>
          ,
          <string-name>
            <given-names>W.R.</given-names>
            <surname>Rudnicki</surname>
          </string-name>
          ,
          <article-title>The all relevant feature selection using random forest</article-title>
          ,
          <source>arXiv preprint arXiv:1106.5112</source>
          ,
          <year>2011</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>G.</given-names>
            <surname>Chandrashekar</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Sahin</surname>
          </string-name>
          ,
          <article-title>A survey on feature selection methods</article-title>
          .
          <source>Computers Electrical Engineering</source>
          ,
          <volume>40</volume>
          (
          <issue>1</issue>
          ),
          <fpage>16</fpage>
          -
          <lpage>28</lpage>
          (
          <year>2014</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>A.</given-names>
            <surname>Bommert</surname>
          </string-name>
          ,
          <string-name>
            <given-names>X.</given-names>
            <surname>Sun</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Bischl</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Rahnenführer</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Lang</surname>
          </string-name>
          ,
          <article-title>Benchmark for filter methods for feature selection in high-dimensional classification data</article-title>
          .
          <source>Computational Statistics Data Analysis</source>
          ,
          <volume>143</volume>
          ,
          <issue>106839</issue>
          (
          <year>2020</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>B.</given-names>
            <surname>Venkatesh</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Anuradha</surname>
          </string-name>
          ,
          <article-title>A review of feature selection and its methods</article-title>
          .
          <source>Cybernetics and Information Technologies</source>
          ,
          <volume>19</volume>
          (
          <issue>1</issue>
          ),
          <fpage>3</fpage>
          -
          <lpage>26</lpage>
          (
          <year>2019</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>L. N.</given-names>
            <surname>Sanchez-Pinto</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L. R.</given-names>
            <surname>Venable</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Fahrenbach</surname>
          </string-name>
          ,
          <string-name>
            <surname>M. M. Churpek</surname>
          </string-name>
          ,
          <article-title>Comparison of variable selection methods for clinical predictive modeling</article-title>
          .
          <source>International journal of medical informatics</source>
          ,
          <volume>116</volume>
          ,
          <fpage>10</fpage>
          -
          <lpage>17</lpage>
          (
          <year>2018</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>P.</given-names>
            <surname>Hayati Rezvan</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K. J.</given-names>
            <surname>Lee</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. A.</given-names>
            <surname>Simpson</surname>
          </string-name>
          ,
          <article-title>The Rise of Multiple Imputation: A Review of the Reporting and Implementation of the Method in Medical Research</article-title>
          .
          <source>BMC Med Res Methodol</source>
          , ,
          <volume>15</volume>
          ,
          <issue>30</issue>
          (
          <year>2015</year>
          ). https://doi.org/10.1186/s12874-015-0022-1
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>N. Ahmat</given-names>
            <surname>Zainuri</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A. A.</given-names>
            <surname>Jemain</surname>
          </string-name>
          ,
          <string-name>
            <surname>N.</surname>
          </string-name>
          <article-title>A Muda, Comparison of Various Imputation Methods for Missing Values in Air Quality Data</article-title>
          . JSM,
          <volume>44</volume>
          ,
          <fpage>449</fpage>
          -
          <lpage>456</lpage>
          (
          <year>2015</year>
          ). https://doi.org/10.17576/jsm2015-4403-17
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>C. A.</given-names>
            <surname>Leke</surname>
          </string-name>
          , T. Marwala,
          <article-title>Introduction to Missing Data Estimation</article-title>
          .
          <source>Deep Learning and Missing Data in Engineering Systems</source>
          ,
          <volume>1</volume>
          -
          <fpage>20</fpage>
          (
          <year>2019</year>
          ). https://doi.org/10.1117/12.2053057
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>C.</given-names>
            <surname>Wang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.</given-names>
            <surname>Shakhovska</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Sachenko</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Komar</surname>
          </string-name>
          ,
          <article-title>A New Approach for Missing Data Imputation in Big Data Interface</article-title>
          .
          <source>Information Technology and Control</source>
          ,
          <volume>49</volume>
          (
          <issue>4</issue>
          ),
          <fpage>541</fpage>
          -
          <lpage>555</lpage>
          (
          <year>2020</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>M.</given-names>
            <surname>Azmi</surname>
          </string-name>
          ,
          <string-name>
            <given-names>G. C.</given-names>
            <surname>Runger</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Berrado</surname>
          </string-name>
          ,
          <article-title>Interpretable regularized class association rules algorithm for classification in a categorical data space</article-title>
          .
          <source>Information Sciences</source>
          ,
          <volume>483</volume>
          ,
          <fpage>313</fpage>
          -
          <lpage>331</lpage>
          (
          <year>2019</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>F.</given-names>
            <surname>Thabtah</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Cowling</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Peng</surname>
          </string-name>
          ,
          <article-title>MCAR: multi-class classification based on association rule</article-title>
          .
          <source>In The 3rd ACS/IEEE International Conference on Computer Systems and Applications</source>
          , (
          <year>2005</year>
          )/
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <given-names>K.</given-names>
            <surname>Mittal</surname>
          </string-name>
          , G. Aggarwal,
          <string-name>
            <given-names>P.</given-names>
            <surname>Mahajan</surname>
          </string-name>
          ,
          <article-title>A comparative study of association rule mining techniques and predictive mining approaches for association classification</article-title>
          .
          <source>International Journal of Advanced Research in Computer Science</source>
          ,
          <volume>8</volume>
          (
          <issue>9</issue>
          ) (
          <year>2017</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <given-names>L.</given-names>
            <surname>Yan</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H. T.</given-names>
            <surname>Zhang</surname>
          </string-name>
          , J. Goncalves,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Xiao</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Wang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Guo</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Yuan</surname>
          </string-name>
          ,
          <article-title>An interpretable mortality prediction model for COVID-19 patients</article-title>
          ,
          <source>Nature machine intelligence</source>
          ,
          <volume>2</volume>
          (
          <issue>5</issue>
          ),
          <fpage>283</fpage>
          -
          <lpage>288</lpage>
          (
          <year>2020</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <given-names>A.</given-names>
            <surname>Trickey</surname>
          </string-name>
          , M. T. May,
          <string-name>
            <given-names>P.</given-names>
            <surname>Schommers</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Tate</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S. M.</given-names>
            <surname>Ingle</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. L.</given-names>
            <surname>Guest</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J. A.</given-names>
            <surname>Sterne</surname>
          </string-name>
          ,
          <article-title>CD4: CD8 ratio and CD8 count as prognostic markers for mortality in human immunodeficiency virus - infected patients on antiretroviral therapy: the Antiretroviral Therapy Cohort Collaboration (ART-CC)</article-title>
          ,
          <source>Clinical Infectious Diseases</source>
          ,
          <volume>65</volume>
          (
          <issue>6</issue>
          ),
          <fpage>959</fpage>
          -
          <lpage>966</lpage>
          (
          <year>2017</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [15]
          <string-name>
            <given-names>N.</given-names>
            <surname>Hasan</surname>
          </string-name>
          ,
          <article-title>A Methodological Approach for Predicting COVID-19 Epidemic Using EEMD-ANN Hybrid Model</article-title>
          .
          <source>Internet of Things</source>
          <volume>11</volume>
          ,
          <issue>100228</issue>
          (
          <year>2020</year>
          ). https://doi.org/10.1016/j.iot.
          <year>2020</year>
          .100228
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          [16]
          <string-name>
            <given-names>L. N.</given-names>
            <surname>Sanchez-Pinto</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L. R.</given-names>
            <surname>Venable</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Fahrenbach</surname>
          </string-name>
          ,
          <string-name>
            <surname>M. M. Churpek</surname>
          </string-name>
          ,
          <article-title>Comparison of variable selection methods for clinical predictive modeling</article-title>
          .
          <source>International journal of medical informatics</source>
          ,
          <volume>116</volume>
          ,
          <fpage>10</fpage>
          -
          <lpage>17</lpage>
          (
          <year>2018</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          [17]
          <string-name>
            <given-names>A.</given-names>
            ,
            <surname>Dariush</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Bastanfard</surname>
          </string-name>
          ,
          <article-title>An objective method to evaluate exemplar‐based inpainted images quality using Jaccard index</article-title>
          ,
          <source>Multimedia Tools and Applications</source>
          <volume>80</volume>
          .17,
          <fpage>26199</fpage>
          -
          <lpage>26212</lpage>
          (
          <year>2021</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          [18]
          <string-name>
            <given-names>A.</given-names>
            ,
            <surname>Federico</surname>
          </string-name>
          , et al.,
          <article-title>Association rules extraction for the identification of functional dependencies in complex technical infrastructures</article-title>
          ,
          <source>Reliability Engineering System Safety</source>
          <volume>209</volume>
          (
          <year>2021</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref19">
        <mixed-citation>
          [19]
          <string-name>
            <surname>I.-D. Borlea</surname>
          </string-name>
          , et al.,
          <article-title>A unified form of fuzzy C-means and K-means algorithms and its partitional implementation</article-title>
          ,
          <source>Knowledge-Based Systems</source>
          <volume>214</volume>
          (
          <year>2021</year>
          ):
          <fpage>106731</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref20">
        <mixed-citation>
          [20]
          <string-name>
            <surname>A.M. Miggelbrink</surname>
          </string-name>
          , et al.,
          <article-title>CD4 T-cell exhaustion: Does it exist and what are its roles in cancer?</article-title>
          ,
          <source>Clinical Cancer Research</source>
          <volume>27</volume>
          .
          <volume>21</volume>
          ,
          <fpage>5742</fpage>
          -
          <lpage>5752</lpage>
          (
          <year>2021</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref21">
        <mixed-citation>
          [21]
          <string-name>
            <given-names>L. I. U.</given-names>
            <surname>Huan-xia</surname>
          </string-name>
          , et al.,
          <article-title>Analysis of CD4/CD8 ratio in HIV-infected patients who accepted initial antiretroviral therapy for 48 weeks</article-title>
          ,
          <source>China Tropical Medicine 21.3</source>
          (
          <year>2021</year>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>