Machine Learning Modeling of Aerobic Biodegradation for Azo Dyes and Hexavalent Chromium

: The present study emphasizes the e ﬃ cacy of a biosurfactant-producing bacterial strain Klebsiella sp. KOD36 in biodegradation of azo dyes and hexavalent chromium individually and in a simultaneous system. The bacterial strain has exhibited a considerable potential for biodegradation of chromium and azo dyes in single and combination systems (maximum 97%, 94% in an individual and combined system, respectively). Simultaneous aerobic biodegradation of azo dyes and hexavalent chromium (SBAHC) was modeled using machine learning programming, which includes gene expression programming, random forest, support vector regression, and support vector regression-fruit ﬂy optimization algorithm. The correlation coe ﬃ cient includes the dispersion index, and the Willmott agreement index was employed as statistical metrics to assess the performance of each model separately. In addition, the Taylor diagram was used to further investigate the methods used. The ﬁndings of the present study were that the support vector regression-fruitﬂy optimization algorithm (SVR-FOA) with correlation coe ﬃ cient (CC) of 0.644, (scattered index) SI of 0.374, and (Willmott’s index of agreement) WI of 0.607 performed better than the autonomous support vector regression (SVR), gene expression programming (GEP), and random forest (RF) methods. In addition, the standalone SVR model with CC of 0.146, SI of 0.473, and WI of 0.408 ranked the second best. In summary, the SBAHC can be accurately estimated using the hybrid SVR-FOA method. In other words, FOA has proven to be a powerful optimization algorithm for increasing the accuracy of the SVR method. and A.K.; methodology, Z.A., A.M. and H.S.; project administration, H.S. and S.M.; resources, Z.A., H.Z. and M.S.; software, A.M.; supervision, A.M.; validation, N.N.; visualization, M.S. and N.N.; writing—original draft, Z.A., and A.M.; writing—review and editing, Z.A., A.M. and N.N. All authors


Introduction
Many salts, in particular azo dyes and chromium sulphate, are the most commonly used chemicals in various processes in the leather tanning industry [1]. Azo dyes are regarded as xenobiotics because they are very resistant. Their chemical structure is designed to resist discoloration when exposed to and enhances solubility of compounds in aqueous system. Moreover, application of SVR-FOA for estimation of simultaneous biodegradation of azo dyes and chromium (VI) in the current study is a novel strategy as the technique is superior compared to previously employed algorithm programs, as SVR-FOA has fewer parameters and is easy to program. Hence, it optimizes the complex linear regression problem inspired by the fruit fly food searching phenomena by a specialized way of smell and vision. The main contribution of the study was to develop a novel modeling approach based on machine learning programing for estimating and predicting optimal parameters for SBAHC by the newly isolated bacterial strain Klebsiella sp. KOD36.

Culture Medium, Chemicals, and Microorganism
Reactive black dyes and di-potassium chromate were used in this study to determine possible degradation by bacterial stain. All solutions were prepared in distilled water at 121 • C for 20 min. A working solution of 15, 100, and 150 mg L −1 azo dyes and 2, 5, and 10 mg L −1 chromium was prepared from the stock solution. The diphenyl carbide reagent was used to determine the chromium concentration in the spectrophotometer. Degradation studies using mineral salts (MSM), NaCl (1.0 g/L), CaCl 2 (0.1 g/L), KH 2 PO 4 (1.0 g/L), MgSO 4 . 7H 2 O (0.5 g/L), Na 2 HPO 4 (1.0 g/L), and yeast extract (4.0 g/l) amended with various concentration of Cr(VI) and azo dyes were used for the purification and streaking of bacterial strains already reported to produce biosurfactants. The strain used was isolated previously [33] and identified as Klebsiella sp. through 16 s RNA. The sequence of identified strain was submitted with accession number KT364873 to gene bank [33].

Evaluation of Klebsiella sp. KOD36 Potential for Biodegradation of Azo Dyes and Chromium (VI) in Single and Combined System
The strain Klebsiella sp. KOD36, which produces biosurfactants, was investigated for capability to utilize as bioremediation agent in biodegradation of azo dyes and hexavalent chromium in single and combined aerobic systems. For this purpose, the MSM broth was supplemented with various concentrations of chromium (0, 5, 10 mg L −1 ) and azo dyes (0, 100, 150 mg L −1 ) in single or combination systems (100 mg L −1 azo dyes and 5 mg L −1 Cr, 100 mg L −1 azo dyes and 10 mg L −1 Cr, 150 mg L −1 azo dyes and 5 mg L −1 Cr, 150 mg L −1 azo dyes and 10 mg L −1 Cr) inoculated with Klebsiella sp. KOD36 uniform bacterial community 2.3 × 10 7 colony forming unit (CFU) mL −1 . The glass vials were placed at 150 rpm in a stirred, temperature-controlled incubator for 24 h. Aliquots were regularly taken alternately from each vial to determine their concentration of reduced species. The degradation potential was measured using the technique as described by Desai et al. [34].

Azo Dye and Chromium Concentration Measurement
Chromium and azo dye concentrations were measured using the technique described by Desai et al. [34]. In short, 1 mL of desired sample was drawn from a tube and subjected to centrifugation at 10,000 rev/min for 20 min. The chromium concentration was determined by complexing the chromium with a diphenyl carbide reagent. The pink color was developed, which was determined at wavelengths of 540 nm. The concentration of the azo dye at 570 nm was also calculated for each dye with a spectrophotometer.

Machine Learning Methodolgies
In the present study, our major aim was to design a novel methodology to estimate and predict the simultaneous biodegradation of azo dyes and chromium (VI) in a single and combine system by means of machine learning programing. The acquired data were subjected to use in training GEP, RF, SVR, and SVR-FOA procedures. The gene expression program was up-to-date algorithm methodology that described the relationship between given input and resulting output variables through software development [35].
In contrast to genetic algorithm (GA) and genetic programming (GP), the GEP is actually the combination of these two [35,36]. GEP is a specialized program for solving regression problems. Support vector regression (SVR) is specifically important for its capability of management and required performance in solving problems, particularly nonlinear regression problems [37]. The selection of basic parameters is an important aspect in SVR technique. The main parameters of SVR technique include c, ε, and kernel functions [38]. The drawback of this technique is, in certain cases, the incorrect values of SVR basic parameters may result in either under-or over-tuning. Therefore, an optimal value of each basic parameter should be chosen while in the SVR training phase. Various methodologies have been adopted to select the optimal value of these basic parameter of SVR. Fruit fly optimization algorithm (FOA) is such a technique, introduced by Pan [30], which employs the foraging strategy of a fruit fly to search the optimal position of each basic value in support vector regression. These methods are described as following.

Gene Expression Programming (GEP)
Genetic expression programming is a specialized genetic programing technique that is capable of solving optimal problems by developing an expression tree (ET). The baseline of this technique is a specialized tree-like structure which is first trained as living organisms by changing shape, size, composition, and other affecting factors. Genetic expression programs, similar to living organisms, are coded as fixed traits on chromosomes. Therefore, the GEP is just like a genotype-phenotype expression system that uses a simple genetic information code to exhibit the biologically inspired phenotype. The particular fixed-length chromosome has genetic information similar to actual information stored in a chromosome part. Each chromosome contains several genes that are called sub-entity types. In genetic expression programing, all the other sub-entity types are connected to a root, similar to a tree, and make a connection to each other. These sub-connections in genetic expression technique include division, multiplication, subtraction, and addition [39]. These particular genes, irrespective of their fixed length, have varying size and shape. So, the varying length of different genes allows these genes in GEP to progress adoption and evolution. Each specific area in a gene is known as open reader frame (ORF), which provides solution while exhibiting as code in expression tree [40]. Similar to other evolutionary programs, the genetic expression program is based on scattered information on chromosomes. In a general set of data (population), fitness function is used to evaluate each chromosome and designate a specific value. Different fitness functions in a genetic expression program have been used previously [35]. Suitable chromosomes are picked up in the next generation. These chromosomes are further controlled by particular gene operators after being selected. The process of selection continues until a proper optimal value is attained [35,41].

Support Vector Regression
SVR is normally used for solving regression problems and has been used as an estimation and prediction system in biological systems [42][43][44]. This is a supervised program which utilizes structural risk minimization (SRM), in contrast to empirical risk minimization (ERM) used in conventional neural network. ERM basically reduces the error of a training dataset while SRM is helpful in reducing errors at higher extent. Therefore, SVM has the potential to reduce errors of commonly practiced neural networks [37,38]. The main background of SVR is to design the given input data into more precise dimensions and improve the efficiency of linear regression problems by adopting kernel function; yet these kernel functions are not adequate specifically in more complex linear regression problems. Each function (kernel) must have two features, symmetricity and compliance, with Cauchy-Schwarz criteria to address these issues. These features ensure that new space is capable of being defined by these functions (kernel). In general, the support vector regression performance depends on its parameters, which include ε (intensive zone, which is normally used to fire the training dataset), C (trade-off), and γ (determine of relative error and smoothness). To optimize the various SVR parameters, different algorithms have been developed.

FOA (Fruit Fly Optimization Algorithm)
The FOA is a biological program inspired by the behavior of Drosophila insect for food search (Pan, 2012). The Drosophila insect has a unique quality of superior smell and vision, which differentiate and make it superior from other insects. The insect employs a sense of vision when approaching a food source. This information is transmitted to the insect body, which determines the route leading to the food source identified [45]. The schematic diagram of SVR-FOA methodology, adapted by Nabipour et al. [46], is shown in Figure 2.

Parameters for Evaluation of Models' Performance
The performance (predictive) of the recommended model was evaluated as CC, SI, and WI. These statistics are presented as follows [47,48]: I: CC expressed as:

Modeling Methodologies
The simultaneous biodegradation of azo dyes and chromium was estimated and modeled using DT techniques (e.g., GEP, SVR, RF, and SVR-FOA). DT techniques present simple interpretation of results by handling nonlinear and nonparametric variables. Evaluating the performance of the mentioned techniques was performed by comparing results with the empirical relationships presented by other researchers.
Simultaneous degradation of azo dyes and chromium (%) = −3.87091 Nitrogen + Incubation − Nitrogen 28408Nitrogen In the above formulation, set the mentioned following values for nitrogen as yeast, ammonium sulfate, and urea, while for carbon source glucose set sucrose and starch.

Parameters for Evaluation of Models' Performance
The performance (predictive) of the recommended model was evaluated as CC, SI, and WI. These statistics are presented as follows [47,48]: I: CC expressed as: II: SI follows as: O III: WI expressed as: The targeted and predicted values are denoted as O i and P i , respectively.

Results and Discussion
In the present study, the potential of biosurfactant-producing strain Klebsiella sp. KOD36 was tested for its simultaneous reduction of chromium and reactive black-5 azo dyes (RB-5). Additionally, optimization of environmental and nutritional parameters during simultaneous biodegradation of chromium and azo dyes was assessed using GEP, SVR, RF, and SVR-FOA. Table 3 presents the statistical analysis parameters of the used data. Results presented in Figure 3 represent that for chromium at 5 mg L −1 , 98% reduction was observed as compared to control after 24 h of inoculation of bacterial strains Klebsiella sp. KOD36, while for 10 mg L −1 of chromium up to 80% reduction occurred, and for chromium (5 mg L −1 ) 80% reduction of chromium was observed compared to control ( Figure 4).

Results and Discussion
In the present study, the potential of biosurfactant-producing strain Klebsiella sp. KOD36 was tested for its simultaneous reduction of chromium and reactive black-5 azo dyes (RB-5). Additionally, optimization of environmental and nutritional parameters during simultaneous biodegradation of chromium and azo dyes was assessed using GEP, SVR, RF, and SVR-FOA. Table 3 presents the statistical analysis parameters of the used data. Results presented in Figure 3 represent that for chromium at 5 mg L −1 , 98% reduction was observed as compared to control after 24 h of inoculation of bacterial strains Klebsiella sp. KOD36, while for 10 mg L −1 of chromium up to 80% reduction occurred, and for chromium (5 mg L −1 ) 80% reduction of chromium was observed compared to control (Figure 4). Regarding azo dyes at 100 mg L −1 , a maximum 90% degradation of azo dyes was observed, while for 150 mg L −1 azo dyes' concentration, a maximum 87% degradation was observed after 24 h. The inoculation of biosurfactant-producing strain Klebsiella sp. KOD36 significantly enhanced the biodecolorization of reactive black dyes compared with the sample lacking the bacterial strain. Based on results of the previous study, B. circulans BWL1061 decolorizes azo dyes [37], which lead to achieving the biodecolorization by improving the enzymes responsible for degradation and dyes. The significant improvement in biodegradation of azo dyes and chromium (VI) by biosurfactant at critical micelle concentration (CMC) and biosurfactant-producing bacteria indicates that electrostatic attraction forces and hydrophobic part of the biosurfactant play a vital role in biodegradation of dyes. A similar finding and mechanism was described by a previous study conducted by Liu et al. [49], who described the isolated strain BWL1061 exhibted degradation potential for azo dyes, which may likely have been due to interaction of biosurfactant hydrophobic moiety of biosurfactant and dyes. Similarly, in another study conducted, Thacker and Madamwar [50] described the capability of the biosurfactant-producing bacterial strain Ochrobactrum sp. and Bacillus sp. for the reduction of hexavalent chromium in a batch study experiment.  Regarding azo dyes at 100 mg L −1 , a maximum 90% degradation of azo dyes was observed, while for 150 mg L −1 azo dyes' concentration, a maximum 87% degradation was observed after 24 h. The inoculation of biosurfactant-producing strain Klebsiella sp. KOD36 significantly enhanced the biodecolorization of reactive black dyes compared with the sample lacking the bacterial strain. Based on results of the previous study, B. circulans BWL1061 decolorizes azo dyes [37], which lead to achieving the biodecolorization by improving the enzymes responsible for degradation and dyes. The significant improvement in biodegradation of azo dyes and chromium (VI) by biosurfactant at critical micelle concentration (CMC) and biosurfactant-producing bacteria indicates that electrostatic attraction forces and hydrophobic part of the biosurfactant play a vital role in biodegradation of dyes. A similar finding and mechanism was described by a previous study conducted by Liu et al. [49], who described the isolated strain BWL1061 exhibted degradation potential for azo dyes, which may likely have been due to interaction of biosurfactant hydrophobic moiety of biosurfactant and dyes. Similarly, in another study conducted, Thacker and Madamwar [50] described the capability of the biosurfactant-producing bacterial strain Ochrobactrum sp. and Bacillus sp. for the reduction of hexavalent chromium in a batch study experiment.

In Combined System
Microorganisms are important in the way that they are excellent bioremediation agents for heavy metal contamination (soil and water). Microorganisms showing a significant resistance to heavy metals have potential as remediation agent in detoxification of these heavy metals. However, in certain cases, more specifically under co-contamination, their efficiency is hindered, as heavy metals cause toxicity to microorganisms and lower their efficiency for biodegradation of azo dyes. In this scenario, the following investigation was designed to investigate the efficacy of biosurfactant and biosurfactant-producing bacteria for decolorization of azo dyes with various concentrations of chromium (VI).
Regarding simultaneous degradation results (azo dyes at 100 mg L −1 and Cr at 5 mg L −1 concentration), 89% degradation was observed after 24 h of inoculation of bacterial strains Klebsiella sp. KOD36 after 24 h, as compared to control ( Figure 5). While for 100 mg L −1 (azo dyes and 10 mg L −1 Cr), up to 94% simultaneous degradation was observed after 24 h of inoculation of bacterial strains Klebsiella sp. KOD36, as compared to control. The percent degradation for 150 mg L −1 azo dyes and 5 mg L −1 Cr, and 150 mg L −1 azo dyes and 10 mg L −1 Cr were 91% and 82%, respectively, after 24 h, as compared to control. Similar effects could also be observed in the study conducted by Halmi et al. [51], who isolated a novel potential stain for decolorization of four different dyes, namely amaranth dye, Biebrich scarlet, direct blue, and metanil yellow, under aerobic environment. The isolate bacterial strain exhibited decolorization a maximum of 52% of dyes (initial concentration 150 ppm potassium dichromate) in nutrient broth medium after an incubation of 24 h under shaking at 150 rpm. The enhanced biodecolorization effect could be likely have been due to the fact that biosurfactants reduce the toxicity of hexavalent chromium by entrapping it in micelles and reduce their bioavailability to microorganisms and, meanwhile, bacterial decolorization for azo dyes was enhanced. chromium (VI).
Regarding simultaneous degradation results (azo dyes at 100 mg L −1 and Cr at 5 mg L −1 concentration), 89% degradation was observed after 24 h of inoculation of bacterial strains Klebsiella sp. KOD36 after 24 h, as compared to control ( Figure 5). While for 100 mg L −1 (azo dyes and 10 mg L −1 Cr), up to 94% simultaneous degradation was observed after 24 h of inoculation of bacterial strains Klebsiella sp. KOD36, as compared to control. The percent degradation for 150 mg L −1 azo dyes and 5 mg L −1 Cr, and 150 mg L −1 azo dyes and 10 mg L −1 Cr were 91% and 82%, respectively, after 24 h, as compared to control. Similar effects could also be observed in the study conducted by Halmi et al. [51], who isolated a novel potential stain for decolorization of four different dyes, namely amaranth dye, Biebrich scarlet, direct blue, and metanil yellow, under aerobic environment. The isolate bacterial strain exhibited decolorization a maximum of 52% of dyes (initial concentration 150 ppm potassium dichromate) in nutrient broth medium after an incubation of 24 h under shaking at 150 rpm. The enhanced biodecolorization effect could be likely have been due to the fact that biosurfactants reduce the toxicity of hexavalent chromium by entrapping it in micelles and reduce their bioavailability to microorganisms and, meanwhile, bacterial decolorization for azo dyes was enhanced. The presence of organic and non-organic compounds emit mixed pollution in industrial zones [38]. Conventional wastewater contains various types of organic and inorganic contaminants, which require immediate attention. Chang et al. [52] conducted a study using a high salinity-tolerant bacterial strain, A12 and L, for its biodegradability for sulfamethoxazole (SMX). They also found that under aerobic and anaerobic conditions the bacterial strain denoted as A12 and L showed a significant degradation of SMX in milkfish culture pond sediment batch experiments. Biosurfactant plays a vital role in this regard. Micelle formation in biosurfactant entrapped the heavy metal (chromium) in its core, thus reducing the bioavailability to bacterial cell by preventing the cells of chromium [30]. Subsequently, bacterial cell efficiently decolorized azo dyes in the presence of Cr(VI). The above findings showed that Klebsiella sp. KOD36 is a proper choice for reclaiming azo dyes and metal (Cr(VI))-contaminated sites.

Modeling Outcomes
There is not any significant instruction for splitting training and testing data. In this study, data were divided into training (67%) and testing (33%) to develop GEP, RF, SVR, and SVR-FOA models for SBAHC estimation. Moreover, SVR-FOA optimized the default values of SVR for increasing the accuracy of predictions. So, the default and optimized values of SVR and SVR-FOA are presented in Table 4. Also, the GEP modeling functional parameters are shown in Table 5. Therefore, with default and optimized parameters, the defined scenarios for SVR, SVR-FOA, GEP, and RF models' parameters are shown in Table 6. This is clear in Table 6 that SVR-FOA showed the maximum estimation performance. In other words, SVR-FOA with CC value of 0.644, SI value of 0.374, and WI value of 0.607 estimated SBAHC more accurately than other considered models and, hence, chosen as the best model among others studies, followed by SVR with CC value of 0.146, SI value of 0.473, and WI value of 0.408. Although the CC value of SVR was low, due to lower SI error, it may be more appropriate than GEP and RF models. Additionally, among GEP and RF models, GEP showed weak performance with CC value of 0.387, SI value of 0.647, and WI value of 0.456. Furthermore, it can be concluded from Table 6 that SVR-FOA increased CC values of SVR, GEP, and RF by 341.1%, 66.4%, and 57.1%, respectively. Also, it reduced SI values of the mentioned models by 20.9%, 42.2%, and 27.9%, respectively. Finally, SVR-FOA increased the WI values of the mentioned models by 48.8%, 33.1%, and 19.7%, respectively. Although the GEP had lower accuracy in predicting simultaneous biodegradation, the model could be used for SBAHC. The mentioned GEP formulation is presented below. SBAHC = −3.87091 Nitrogen + Incubation − Nitrogen +
These performance parameters of various models used in the present study are also shown as bar chart (Figure 6). It is obvious from the chart that SVR-FOA had highest potential for predicting SBAHC in prediction and estimation.
( ) In the above formulation, the following values should be considered for nitrogen, yeast 1, ammonium 2, urea 3; and for carbon, glucose 1 and sucrose 2. These performance parameters of various models used in the present study are also shown as bar chart (Figure 6). It is obvious from the chart that SVR-FOA had highest potential for predicting SBAHC in prediction and estimation. The SBAHC predictive results of various models is also shown in Figure 7. It can be observed from Figure 5 that SVR-FOA had a higher performance than other considered models. Furthermore, Figure 8 indicates scatter plots of prediction of the SBAHC values with SVR-FOA, SVR, GEP, and RF models. The less-scattered points exhibited by SVR-FOA is a clear indication that the values of SVR-FOA were more accurate than standalone SVR, GEP, and RF models. The SBAHC predictive results of various models is also shown in Figure 7. It can be observed from Figure 5 that SVR-FOA had a higher performance than other considered models. Furthermore, Figure 8 indicates scatter plots of prediction of the SBAHC values with SVR-FOA, SVR, GEP, and RF models. The less-scattered points exhibited by SVR-FOA is a clear indication that the values of SVR-FOA were more accurate than standalone SVR, GEP, and RF models.    Furthermore, Taylor diagrams (TD) were employed to examine standard deviation (SD) and CC values for the SVR-FOA, SVR, GEP, and RF models. Figure 9 presents TD for all models. It can be understood from Figure 7 that SVR-FOA (a point with grey color), due to a shorter distance from the Furthermore, Taylor diagrams (TD) were employed to examine standard deviation (SD) and CC values for the SVR-FOA, SVR, GEP, and RF models. Figure 9 presents TD for all models. It can be understood from Figure 7 that SVR-FOA (a point with grey color), due to a shorter distance from the observed green point, provided relatively precise predictions of SBAHC values. As a conclusive remark, it can be stated that SVR-FOA with optimized values (C = 1.7242, γ= 0.0517m andε = 0.0468) and using input parameters of temperature, pH, incubation period (IP), and shaking is more capable for accurate estimation of SBAHC in comparison to standalone SVR and GEP models, and may be recommended for further implementations [36,53,54]. and using input parameters of temperature, pH, incubation period (IP), and shaking is more capable for accurate estimation of SBAHC in comparison to standalone SVR and GEP models, and may be recommended for further implementations [36,53,54]. The need for a robust model for estimation of large number of input variables is obvious in the recent world. Another study, investigated by Amato et al. [55], illustrated that in a well-defined geographic area the application of social nets (on-line) may increase the detection efficiency (real time) and alert diffusion. They proposed a multicomplex big data system that uses clustering event detection techniques along with multimedia content and biologically inspired programming to develop alerts.

Conclusions
Klebsiella sp. KOD36 significantly boosted the biodegradation of azo dyes. Additionally, the capabilities of SVR, SVR-FOA, GEP, and RF models in estimation of SBAHC values were inspected. Accordingly, the enactments of studied methods were comprehensively examined using CC, SI, and WI parameters. Also, Taylor diagrams were utilized for further assessment. The obtained results indicated that SVR-FOA with CC of 0.644, SI of 0.374, and WI of 0.607 had better performance The need for a robust model for estimation of large number of input variables is obvious in the recent world. Another study, investigated by Amato et al. [55], illustrated that in a well-defined geographic area the application of social nets (on-line) may increase the detection efficiency (real time) and alert diffusion. They proposed a multicomplex big data system that uses clustering event detection techniques along with multimedia content and biologically inspired programming to develop alerts.

Conclusions
Klebsiella sp. KOD36 significantly boosted the biodegradation of azo dyes. Additionally, the capabilities of SVR, SVR-FOA, GEP, and RF models in estimation of SBAHC values were inspected. Accordingly, the enactments of studied methods were comprehensively examined using CC, SI, and WI parameters. Also, Taylor diagrams were utilized for further assessment. The obtained results indicated that SVR-FOA with CC of 0.644, SI of 0.374, and WI of 0.607 had better performance comparing to standalone SVR, GEP, and RF methods. Moreover, standalone SVR model ranked the second best with CC of 0.146, SI of 0.473, and WI of 0.408. This can be verified by the presented Taylor diagram. Because of less-scattered points exhibited by SVR-FOA, it can be concluded that the estimates of SVR-FOA were much more accurate than other studied models. Conclusively, a fruit fly optimization algorithm had remarkable impact in reducing the prediction errors of a standalone SVR method and it can be recommended for SBAHC estimation. However, the main drawback in using the SVR-based model is that, when there is a big dataset or large number of population load, the time consumed as learning/training is very high, thus determination of parameters involves mainly the researcher experience.