COVID-19 Outbreak Prediction with Machine Learning
==================================================

* Sina F. Ardabili
* Amir Mosavi
* Pedram Ghamisi
* Filip Ferdinand
* Annamaria R. Varkonyi-Koczy
* Uwe Reuter
* Timon Rabczuk
* Peter M. Atkinson

## Abstract

Several outbreak prediction models for COVID-19 are being used by officials around the world to make informed-decisions and enforce relevant control measures. Among the standard models for COVID-19 global pandemic prediction, simple epidemiological and statistical models have received more attention by authorities, and they are popular in the media. Due to a high level of uncertainty and lack of essential data, standard models have shown low accuracy for long-term prediction. Although the literature includes several attempts to address this issue, the essential generalization and robustness abilities of existing models needs to be improved. This paper presents a comparative analysis of machine learning and soft computing models to predict the COVID-19 outbreak. Among a wide range of machine learning models investigated, two models showed promising results (i.e., multi-layered perceptron, MLP, and adaptive network-based fuzzy inference system, ANFIS). Based on the results reported here, and due to the highly complex nature of the COVID-19 outbreak and variation in its behavior from nation-to-nation, this study suggests machine learning as an effective tool to model the outbreak.

Keywords
*   COVID-19
*   Coronavirus disease
*   Coronavirus
*   SARS-CoV-2
*   model
*   prediction
*   machine learning

## 1. Introduction

Access to accurate outbreak prediction models is essential to obtain insights into the likely spread and consequences of infectious diseases. Governments and other legislative bodies rely on insights from prediction models to suggest new policies and to assess the effectiveness of the enforced policies [1]. The novel Coronavirus disease (COVID-19) has been reported to infect more than 2 million people, with more than 132,000 confirmed deaths worldwide. The recent global COVID-19 pandemic has exhibited a nonlinear and complex nature [2]. In addition, the outbreak has differences with other recent outbreaks, which brings into question the ability of standard models to deliver accurate results [3]. Besides the numerous known and unknown variables involved in the spread, the complexity of population-wide behavior in various geopolitical areas and differences in containment strategies had dramatically increased model uncertainty [4]. Consequently, standard epidemiological models face new challenges to deliver more reliable results. To overcome this challenge, many novel models have emerged which introduce several assumptions to modeling (e.g., adding social distancing in the form of curfews, quarantines, etc.) [5-7].

To elaborate on the effectiveness of enforcing such assumptions understanding standard dynamic epidemiological (e.g., susceptible-infected-recovered, SIR) models is essential [8]. The modeling strategy is formed around the assumption of transmitting the infectious disease through contacts, considering three different classes of well-mixed populations; susceptible to infection (class *S*), infected (class *I*), and the removed population (class *R* is devoted to those who have recovered, developed immunity, been isolated or passed away). It is further assumed that the class *I* transmits the infection to class *S* where the number of probable transmissions is proportional to the total number of contacts [9-11]. The number of individuals in the class *S* progresses as a time-series, often computed using a basic differential equation as follows: ![Formula][1]</img>  where *I* is the infected population, and *S*is the susceptible population both as fractions. *α*represents the daily reproduction rate of the differential equation, regulating the number of susceptible infectious contacts. The value of *S*in the time-series produced by the differential equation gradually declines. Initially, it is assumed that at the early stage of the outbreak *S*≈1while the number of individuals in class *I* is negligible. Thus, the increment ![Graphic][2]</img> becomes linear and the class *I* eventually can be computed as follows: ![Formula][3]</img>  where *β*regulates the daily rate of new infections by quantifying the number of infected individuals competent in the transmission. Furthermore, the class *R*, representing individuals excluded from the spread of infection, is computed as follows: ![Formula][4]</img> 

Under the unconstrained conditions of the excluded group, Eq. 3, the outbreak exponential growth can be computed as follows: ![Formula][5]</img> 

The outbreaks of a wide range of infectious diseases have been modeled using Eq. 4. However, for the COVID-19 outbreak prediction, due to the strict measures enforced by authorities, the susceptibility to infection has been manipulated dramatically. For example, in China, Italy, France, Hungary and Spain the SIR model cannot present promising results, as individuals committed voluntarily to quarantine and limited their social interaction. However, for countries where containment measures were delayed (e.g., United States) the model has shown relative accuracy [12]. Figure. 1 shows the inaccuracy of conventional models applied to the outbreak in Italy by comparing the actual number of confirmed infections and epidemiological model predictions1. The SEIR models through considering the significant incubation period during which individuals have been infected showed progress in improving the model accuracy for Varicella and Zika outbreak [13,14]. SEIR models assume that the incubation period is a random variable and similarly to the SIR model, there would be a disease-free-equilibrium [15,16]. It is worth mentioning that SEIR model will not work well where the parameters are non-stationary through time [17]. A key cause of non-stationarity is where the social mixing (which determines the contact network) changes through time. Social mixing determines the reproductive number *R* which is the number of susceptible individuals that an infected person will infect. Where *R* is less than 1 the epidemic will die out. Where it is greater than 1 it will spread. *R* for COVID-19 prior to lockdown was estimated as a massive 4 presenting a pandemic. It is expected that lockdown measures should bring *R* down to less than 1. the KEY reason why SEIR models are difficult to fit for COVID-19 is non-stationarity of mixing, caused by nudging (step-by-step) intervention measures.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F1)

Figure 1. Italy’s COVID-19 outbreak: the actual number of confirmed infections vs. epidemiological model.

One can calculate that standard epidemiological models can be effective and reliable only if (a) the social interactions are stationary through time (i.e., no changes in interventions or control measures), or (b) there exists a great deal of knowledge of class *R* with which to compute Eq. 3. Often to acquire information on class *R*, several novel models included data from social media or call data records (CDR), which showed promising results [18-25]. However, observation of the behavior of COVID-19 in several countries demonstrates a high degree of uncertainty and complexity [26]. Thus, for epidemiological models to be able to deliver reliable results, they must be adapted to the local situation with an insight into susceptibility to infection [27]. This imposes a huge limit on the generalization ability and robustness of conventional models. Advancing accurate models with a great generalization ability to be scalable to model both the regional and global pandemic is, thus, essential [28].

A further drawback of conventional epidemiological models is the short lead-time. To evaluate the performance of the models, the median success of the outbreak prediction presents useful information. The median prediction factor can be calculated as follows: ![Formula][6]</img> 

As the lead-time increases, the accuracy of the model declines. For instance, for the COVID-19 outbreak in Italy, the accuracy of the model for more than 5-days-in-the-future reduces from *f*=1for the first five days to *f*=0.86 for day 6 [12].

Due to the complexity and the large-scale nature of the problem in developing epidemiological models, machine learning (ML) has recently gained attention for building outbreak prediction models. ML approaches aim at developing models with higher generalization ability and greater prediction reliability for longer lead-times [29-33].

Although ML methods were used in modeling former pandemics (e.g., Ebola, Cholera, swine fever, H1N1 influenza, dengue fever, Zika, oyster norovirus [8,34-43]), there is a gap in the literature for peer-reviewed papers dedicated to COVID-19. Table 1 represents notable ML methods used for outbreak prediction. These ML methods are limited to the basic methods of random forest, neural networks, Bayesian networks, Naïve Bayes, genetic programming and classification and regression tree (CART). Although ML has long been established as a standard tool for modeling natural disasters and weather forecasting [44,45], its application in modeling outbreak is still in the early stages. More sophisticated ML methods (e.g., hybrids, ensembles) are yet to be explored. Consequently, the contribution of this paper is to explore the application of ML for modeling the COVID-19 pandemic. This paper aims to investigate the generalization ability of the proposed ML models and the accuracy of the proposed models for different lead-times.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T1)

Table 1. 
Notable ML methods for outbreak prediction

The rest of this paper is organized as follows. Section two describes the methods and materials. The results are given in section three. Sections four and five present the discussion and the conclusions, respectively.

## 2. Materials and Methods

Data were collected from [https://www.worldometers.info/coronavirus/country](https://www.worldometers.info/coronavirus/country) for five countries, including Italy, Germany, Iran, USA, and China on total cases over 30 days. Figure 2 presents the total case number (cumulative statistic) for the considered countries. Currently, to contain the outbreak, the governments have implemented various measures to reduce transmission through inhibiting people’s movements and social activities. Although for advancing the epidemiological models information on changes in social distancing is essential, for modeling with machine learning no assumption is required. As can be seen in Figure 2, the growth rate in China is greater than that for Italy, Iran, Germany and the USA in the early weeks of the disease.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F2)

Figure 2. Cumulative number of cases for five countries during thirty days ([https://www.worldometers.info/coronavirus/country](https://www.worldometers.info/coronavirus/country))

The next step is to find the best model for the estimation of the time-series data. Logistic, Linear, Logarithmic, Quadratic, Cubic, Compound, Power and exponential equations (Table 2) are employed to develop the desired model.

View this table:
[Table 2.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T2)

Table 2. 
Models for mathematical forecasting

A, B, C, µ, and L are parameters (constants) that characterize the above-mentioned functions. These constants need to be estimated to develop an accurate estimation model. One of the goals of this study was to model time-series data based on the logistic microbial growth model. For this purpose, the modified equation of logistic regression was used to estimate and predict the prevalence (i.e., *I*/Population at a given time point) of disease as a function of time. Estimation of the parameters was performed using evolutionary algorithms like GA, particle swarm optimizer, and the grey wolf optimizer. These algorithms are discussed in the following.

### Evolutionary algorithms

Evolutionary algorithms (EA) are powerful tools for solving optimization problems through intelligent methods. These algorithms are often inspired by natural processes to search for all possible answers as an optimization problem [46-48]. In the present study, the frequently used algorithms, (i.e., genetic algorithm (GA), particle swarm optimizer (PSO) and grey wolf optimizer (GWO)) are employed to estimate the parameters by solving a cost function.

#### Genetic Algorithm (GA)

GAs are considered a subset of “computational models” inspired by the concept of evolution [49]. These algorithms use “Potential Solutions” or “Candidate Solutions” or “Possible Hypotheses” for a specific problem in a “chromosome-like” data structure. GA maintains vital information stored in these chromosome data structures by applying “Recombination Operators” to chromosome-like data structures [50-53]. In many cases, GAs are employed as “Function Optimizer” algorithms, which are algorithms used to optimize “Objective Functions.” Of course, the range of applications that use the GA to solve problems is very wide [52,54]. The implementation of the GA usually begins with the production of a population of chromosomes generated randomly and bound up and down by the variables of the problem. In the next step, the generated data structures (chromosomes) are evaluated, and chromosomes that can better display the optimal solution of the problem are more likely to be used to produce new chromosomes. The degree of “goodness” of an answer is usually measured by the population of the current candidate’s answers [55-59]. The main algorithm of a GA process is demonstrated in Figure 3.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F3)

Figure 3. GA algorithm

In the present study, GA [59] was employed for estimation of the parameters of Eq. 6 to 13. The population number was selected to be 300 and the maximum generation (as iteration number) was determined to be 500 according to different trial and error processes to reduce the cost function value. The cost function was defined as the mean square error between the target and estimated values according to Eq. 14: ![Formula][7]</img>  where, *Es* refer to estimated values, *T* refers to the target values and *N* refers to the number of data.

#### Particle Swarm Optimization (PSO)

In 1995, Kennedy and Eberhart [60] introduced the PSO as an uncertain search method for optimization purposes. The algorithm was inspired by the mass movement of birds looking for food. A group of birds accidentally looked for food in a space. There is only one piece of food in the search space. Each solution in PSO is called a particle, which is equivalent to a bird in the bird’s mass movement algorithm. Each particle has a value that is calculated by a competency function which increases as the particle in the search space approaches the target (food in the bird’s movement model). Each particle also has a velocity that guides the motion of the particle. Each particle continues to move in the problem space by tracking the optimal particles in the current state [60-62]. The PSO method is rooted in Reynolds’ work, which is an early simulation of the social behavior of birds. The mass of particles in nature represents collective intelligence. Consider the collective movement of fish in water or birds during migration. All members move in perfect harmony with each other, hunt together if they are to be hunted, and escape from the clutches of a predator by moving another prey if they are to be preyed upon [63-65]. Particle properties in this algorithm include [65-67]:

*   Each particle independently looks for the optimal point.

*   Each particle moves at the same speed at each step.

*   Each particle remembers its best position in the space.

*   The particles work together to inform each other of the places they are looking for.

*   Each particle is in contact with its neighboring particles.

*   Every particle is aware of the particles that are in the neighborhood.

*   Every particle is known as one of the best particles in its neighborhood.

The PSO implementation steps can be summarized as: the first step establishes and evaluates the primary population. The second step determines the best personal memories and the best collective memories. The third step updates the speed and position. If the conditions for stopping are not met, the cycle will go to the second step.

The PSO algorithm is a population-based algorithm [68,69]. This property makes it less likely to be trapped in a local minimum. This algorithm operates according to possible rules, not definite rules. Therefore, PSO is a random optimization algorithm that can search for unspecified and complex areas. This makes PSO more flexible and durable than conventional methods. PSO deals with non-differential target functions because the PSO uses the information result (performance index or target function to guide the search in the problem area). The quality of the proposed route response does not depend on the initial population. Starting from anywhere in the search space, the algorithm ultimately converges on the optimal answer. PSO has great flexibility to control the balance between the local and overall search space. This unique PSO property overcomes the problem of improper convergence and increases the search capacity. All of these features make PSO different from the GA and other innovative algorithms [61,65,67].

In the present study, PSO was employed for estimation of the parameters of Eq. 6 to 13. The population number was selected to be 1000 and the iteration number was determined to be 500 according to different trial and error processes to reduce the cost function value. The cost function was defined as the mean square error between the target and estimated values according to Eq. 14.

#### Grey Wolf Optimizer (GWO)

One recently developed smart optimization algorithm that has attracted the attention of many researchers is the grey wolf algorithm. Like most other intelligent algorithms, GWO is inspired by nature. The main idea of the grey wolf algorithm is based on the leadership hierarchy in wolf groups and how they hunt [70]. In general, there are four categories of wolves among the herd of grey wolves, alpha, beta, delta and omega. Alpha wolves are at the top of the herd’s leadership pyramid, and the rest of the wolves take orders from the alpha group and follow them (usually there is only one wolf as an alpha wolf in each herd). Beta wolves are in the lower tier, but their superiority over Delta and omega wolves allows them to provide advice and help to alpha wolves. Beta wolves are responsible for regulating and orienting the herd based on alpha movement. Delta wolves, which are next in line for the power pyramid in the wolf herd, are usually made up of guards, elderly population, caregivers of damaged wolves, and so on. Omega wolves are also the weakest in the power hierarchy [70]. Eq. 15 to 18 are used to model the hunting tool: ![Formula][8]</img>  ![Formula][9]</img>  ![Formula][10]</img>  ![Formula][11]</img>  where *t* is represents repetition of the algorithm. ![Graphic][12]</img> and ![Graphic][13]</img> are vectors of the prey site and the ![Graphic][14]</img> vectors represent the locations of the grey wolves. ![Graphic][15]</img> is linearly reduced from 2 to 0 during the repetition. ![Graphic][16]</img> and ![Graphic][17]</img> are random vectors where each element can take on realizations in the range [0.1]. The GWO algorithm flowchart is shown in Figure 4.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F4)

Figure 4. GWO algorithm

In the present study, GWO [70] was employed for estimation of the parameters of Eq.1 to 8. The population number was selected to be 500 and the iteration number was determined to be 1000 according to different trial and error processes to reduce the cost function value. The cost function was defined as the mean square error between the target and estimated values according to Eq. 14.

### Machine learning (ML)

ML is regarded as a subset of AI. Using ML techniques, the computer learns to use patterns or “training samples” in data (processed information) to predict or make intelligent decisions without overt planning [71,72]. In other words, ML is the scientific study of algorithms and statistical models used by computer systems that use patterns and inference to perform tasks instead of using explicit instructions [73,74].

Time-series are data sequences collected over a period of time [75], which can be used as inputs to ML algorithms. This type of data reflects the changes that a phenomenon has undergone over time. Let Xt be a time-series vector, in which *x**t* is the outbreak at time point *t* and *T* is the set of all equidistant time points. To train ML methods effectively, we defined two scenarios, listed in Table 3.

View this table:
[Table 3.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T3)

Table 3. 
Input and output variables for training ML methods by time-series data

As can be seen in Table 3, scenario 1 employs data for three weeks to predict the outbreak on day *t* and scenario 2 employs outbreak data for five days to predict the outbreak for day *t*. Both of these scenarios were employed for fitting the ML methods. In the present research, two frequently used ML methods, the multi-layered perceptron (MLP) and adaptive network-based fuzzy inference system (ANFIS) are employed for the prediction of the outbreak in the five countries.

#### Multi-layered perceptron (MLP)

ANN is an idea inspired by the biological nervous system, which processes information like the brain. The key element of this idea is the new structure of the information processing system [76-78]. The system is made up of several highly interconnected processing elements called neurons that work together to solve a problem [78,79]. ANNs, like humans, learn by example. The neural network is set up during a learning process to perform specific tasks, such as identifying patterns and categorizing information. In biological systems, learning is regulated by the synaptic connections between nerves. This method is also used in neural networks [80]. By processing experimental data, ANNs transfer knowledge or a law behind the data to the network structure, which is called learning. Basically, learning ability is the most important feature of such a smart system. A learning system is more flexible and easier to plan, so it can better respond to new issues and changes in processes [81].

In ANNs, with the help of programming knowledge, a data structure is designed that can act like a neuron. This data structure is called a node [82,83]. In this structure, the network between these nodes is trained by applying an educational algorithm to it. In this memory or neural network, the nodes have two active states (on or off) and one inactive state (off or 0), and each edge (synapse or connection between nodes) has a weight. Positive weights stimulate or activate the next inactive node, and negative weights inactivate or inhibit the next connected node (if active) [78,84]. In the ANN architecture, for the neural cell *c*, the input *b**p* enters the cell from the previous cell *p. w**pc* is the weight of the input *b**p* with respect to cell *c* and *a**c* is the sum of the multiplications of the inputs and their weights [85]: ![Formula][18]</img> 

A non-linear function *Ө**c* is applied to *a**c*. Accordingly, *b**c* can be calculated as Eq. 20 [85]: ![Formula][19]</img> 

Similarly, *w**cn* is the weight of the *b**cn* which is the output of *c* to *n. W* is the collection of all the weights of the neural network in a set. For input *x* and output *y, h**w**(x)* is the output of the neural network. The main goal is to learn these weights for reducing the error values between *y* and *h**w**(x)*. That is, the goal is to minimize the cost function *Q(W)*, Eq. 21 [85]: ![Formula][20]</img> 

In the present research, one of the frequently used types of ANN called the MLP [76] was employed to predict the outbreak. MLP was trained using a dataset related to both scenarios (according to Table 2). For the training of the network, 8, 12, and 16 inner neurons were tried to achieve the best response. Results were evaluated by RMSE and correlation coefficient to reduce the cost function value. Figure 5 presents the architecture of the MLP.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F5)

Figure 5. Architecture of MLP

#### Adaptive neuro fuzzy inference system (ANFIS)

An adaptive neuro fuzzy inference system is a type of ANN based on the Takagi-Sugeno fuzzy system [86]. This approach was developed in the early 1990s. Since this system integrates the concepts of neural networks and fuzzy logic, it can take advantage of both capabilities in a unified framework. This technique is one of the most frequently used and robust hybrid ML techniques. It is consistent with a set of fuzzy if-then rules that can be learned to approximate nonlinear functions [87,88]. Hence, ANFIS was proposed as a universal estimator. An important element of fuzzy systems is the fuzzy partition of the input space [89,90]. For input *k*, the fuzzy rules in the input space make a *k* faces fuzzy cube. Achieving a flexible partition for nonlinear inversion is non-trivial. The idea of this model is to build a neural network whose outputs are a degree of the input that belongs to each class [91-93]. The membership functions (MFs) of this model can be nonlinear, multidimensional and, thus, different to conventional fuzzy systems [94-96]. In ANFIS, neural networks are used to increase the efficiency of fuzzy systems. The method used to design neural networks is to employ fuzzy systems or fuzzy-based structures. This model is a kind of division and conquest method. Instead of using one neural network for all the input and output data, several networks are created in this model:

*   A fuzzy separator to cluster input-output data within multiple classes.

*   A neural network for each class.

*   Training neural networks with output input data in the corresponding classes. Figure 6 presents a simple architecture for ANFIS.

In the present study, ANFIS is developed to tackle two scenarios described in table 3. Each input included by two MFs with the Tri. shape; Trap. shape and Gauss. shape MFs. The output MF type was selected to be linear with a hybrid optimizer type.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F6)

Figure 6. ANFIS architecture

#### Evaluation criteria

Evaluation was conducted using the root mean square error (RMSE) and correlation coefficient. These statistics compare the target and output values and calculate a score as an index for the performance and accuracy of the developed methods [87,97]. Table 4 presents the evaluation criteria equations.

View this table:
[Table 4.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T4)

Table 4. 
Model Evaluation metrics

Where, *N* is the number of data, *P* and *A* are, respectively, the predicted (output) and desired (target) values.

## 3. Results

Tables 5 to 12 present the results of the accuracy statistics for the logistic, linear, logarithmic, quadratic, cubic, compound, power and exponential equations, respectively. The coefficients of each equation were calculated by the three ML optimizers; GA, PSO and GWO. The table contains country name, model name, population size, number of iterations, processing time, RMSE and correlation coefficient.

View this table:
[Table 5.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T5)

Table 5. 
Accuracy statistics for the logistic model

View this table:
[Table 6.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T6)

Table 6. 
Accuracy statistics for the linear model

View this table:
[Table 7.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T7)

Table 7. 
Accuracy statistics for the logarithmic model

View this table:
[Table 8.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T8)

Table 8. 
Accuracy statistics for the quadratic model

View this table:
[Table 9.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T9)

Table 9. 
Accuracy statistics for the cubic model

View this table:
[Table 10.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T10)

Table 10. 
Accuracy statistics for the compound model

View this table:
[Table 11.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T11)

Table 11. 
Accuracy statistics for the power model

View this table:
[Table 12.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T12)

Table 12. 
Accuracy statistics for the exponential model

According to Tables 5 to 12, GWO provided the highest accuracy (smallest RMSE and largest correlation coefficient) and smallest processing time compared to PSO and GA for fitting the logistic, linear, logarithmic, quadratic, cubic, power, compound, and exponential-based equations for all five countries. It can be suggested that GWO is a sustainable optimizer due to its acceptable processing time compared with PSO and GA. Therefore, GWO was selected as the best optimizer by providing the highest accuracy values compared with PSO and GA. In general, it can be claimed that GWO, by suggesting the best parameter values for the functions presented in Table 2, increases outbreak prediction accuracy for COVID-19 in comparison with PSO and GA. Therefore, the functions derived by GWO were selected as the best predictors for this research.

Tables 13 to 17 present the description and coefficients of the linear, logarithmic, quadratic, cubic, compound, power, exponential and logistic equations estimated by GWO. Tables 13 to 17 also present the RMSE and *r*-square values for each equation fitted to data for China, Italy, Iran, Germany and USA, respectively.

View this table:
[Table 13.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T13)

Table 13. 
Model description for China fitted by GWO **Model**

View this table:
[Table 15.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T14)

Table 15. 
Model description for Iran fitted by GWO

View this table:
[Table 16.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T15)

Table 16. 
Model description for Germany fitted by GWO

View this table:
[Table 17.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T16)

Table 17. 
Model description for USA fitted by GWO

As is clear from Tables 13 to 17, in general, the logistic equation followed by the quadratic and cubic equations provided the smallest RMSE and the largest *r*-square values for the prediction of COVID-19 outbreak. The claim can also be considered from Figure 7 to 11, which presents the capability and trend of each model derived by GWO in the prediction of COVID-19 cases for China, Italy, Iran, Germany, and the USA, respectively.

![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F7.medium.gif)

[Figure 7.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F7)

Figure 7. Fitness graph for China fitted by GWO

Figures 7 to 11 illustrate the fit of the models investigated in this paper. The best fit for the prediction of COVID-19 cases was achieved for the logistic model followed by cubic and quadratic models for China (Figure 7), logistic followed by cubic models for Italy (Figure 8), cubic followed by logistic and quadratic models for Iran (Figure 9), the logistic model for Germany (Figure 10), and logistic model for the USA (Figure 11).

![Figure 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F8.medium.gif)

[Figure 8.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F8)

Figure 8. Set of models for Italy fitted by GWO

![Figure 9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F9.medium.gif)

[Figure 9.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F9)

Figure 9. Set of models for Iran fitted by GWO

![Figure 10.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F10.medium.gif)

[Figure 10.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F10)

Figure 10. Set of models for Germany fitted by GWO

![Figure 11.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F11.medium.gif)

[Figure 11.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F11)

Figure 11. Set of models for USA fitted by GWO

### Machine learning results

This section presents the results for the training stage of ML methods. MLP and ANFIS were employed as single and hybrid ML methods, respectively. ML methods were trained using two datasets related to scenario 1 and scenario 2. Table 14 presents the results of the training phase.

View this table:
[Table 14.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T17)

Table 14. 
Model description for Italy fitted by GWO

According to Table 18, the dataset related to scenarios 1 and 2 have different performance values. Accordingly, for Italy, the MLP with 16 neurons provided the highest accuracy for scenario 1 and ANFIS with tri. MF provided the highest accuracy for scenario 2. By considering the average values of the RMSE and correlation coefficient, it can be concluded that scenario 1 is more suitable for modeling outbreak cases in Italy, as it provides a higher accuracy (the smallest RMSE and the largest correlation coefficient) than scenario 2.

View this table:
[Table 18.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T18)

Table 18. 
Results for the training phase of the ML methods

For the dataset related to China, for both scenarios, MLP with 12 and 16 neurons, respectively for scenarios 1 and 2, provided the highest accuracy compared with the ANFIS model. By considering the average values of RMSE and correlation coefficient, it can be concluded that scenario 2 with a larger average correlation coefficient and smaller average RMSE than scenario 1 is more suitable for modeling the outbreak in China.

For the dataset of Iran, MLP with 12 neurons in the hidden layer for scenario 1 and ANFIS with Gaussian MF type for scenario 2 provided the best performance for the prediction of the outbreak. By considering the average values of the RMSE and correlation coefficient, it can be concluded that scenario 1 provided better performance than scenario 2. Also, in general, the MLP has higher prediction accuracy compared with the ANFIS method.

In Germany, MLP with 12 neurons in its hidden layer provided the highest accuracy (smallest RMSE and largest correlation coefficient). By considering the average values of the RMSE and correlation coefficient, it can be concluded that scenario 1 is more suitable for the prediction of the outbreak in Germany than scenario 2.

In the USA, the MLP with 8 and 12 neurons, respectively, for scenarios 1 and 2, provided higher accuracy (the smallest RMSE and the largest correlation coefficient values) than the ANFIS model. By considering the average values of the RMSE and correlation coefficient values, it can be concluded that scenario 1 is more suitable than scenario 2, and MLP is more suitable than ANFIS for outbreak prediction.

Figures 12 to 16 present the model fits for Italy, China, Iran, Germany, and the USA, respectively. By comparing Figure 12 to 16 with Figures 7 to 11, it can be concluded that the MLP and the logistic model fitted by GWO provided a better fit than the other models. In addition, the ML methods provided better performance compared with other models.

![Figure 12.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F12.medium.gif)

[Figure 12.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F12)

Figure 12. Set of models for Italy fitted by ML methods

![Figure 13.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F13.medium.gif)

[Figure 13.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F13)

Figure 13. Set of models for China fitted by ML methods

![Figure 14.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F14.medium.gif)

[Figure 14.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F14)

Figure 14. Set of models for Iran fitted by ML methods

![Figure 15.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F15.medium.gif)

[Figure 15.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F15)

Figure 15. Set of models for Germany fitted by ML methods

![Figure 16.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F16.medium.gif)

[Figure 16.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F16)

Figure 16. Set of models for USA fitted by ML methods

#### Comparing the fitted models

This section presents a comparison of the accuracy and performance of the selected models for the prediction of 30 days’ outbreak. Figure 17 to 21 shows the deviation from the target values for the selected models.

![Figure 17.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F17.medium.gif)

[Figure 17.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F17)

Figure 17. Deviation from target value for models related to Italy

![Figure 18.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F18.medium.gif)

[Figure 18.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F18)

Figure 18. Deviation from target value for models related to China

![Figure 19.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F19.medium.gif)

[Figure 19.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F19)

Figure 19. Deviation from target value for models related to Iran

![Figure 20.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F20.medium.gif)

[Figure 20.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F20)

Figure 20. Deviation from target value for models related to Germany

![Figure 21.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F21.medium.gif)

[Figure 21.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F21)

Figure 21. Deviation from target value for models related to USA

As is clear from Figure 17 to 21, the smallest deviation from the target values is related to the MLP for scenario 1 followed by MLP for scenario 2. This indicates the highest performance of the MLP method for the prediction of the outbreak.

Figures 22 to 26 present the outbreak prediction for 75 days and Tables 19 to 23 present the outbreak prediction for 150 days.

View this table:
[Table 19.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T19)

Table 19. 
The outbreak prediction for Italy through 150 days

View this table:
[Table 20.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T20)

Table 20. 
The outbreak prediction for China through 150 days

View this table:
[Table 21.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T21)

Table 21. 
The outbreak prediction for Iran through 150 days

View this table:
[Table 22.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T22)

Table 22. 
The outbreak prediction for Germany through 150 days

View this table:
[Table 23.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T23)

Table 23. 
The outbreak prediction for the USA for 150 days

![Figure 22.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F22.medium.gif)

[Figure 22.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F22)

Figure 22. The outbreak prediction for Italy through 75 days

![Figure 23.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F23.medium.gif)

[Figure 23.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F23)

Figure 23. The outbreak prediction for China through 75 days

![Figure 24.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F24.medium.gif)

[Figure 24.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F24)

Figure 24. The outbreak prediction for Iran through 75 days

![Figure 25.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F25.medium.gif)

[Figure 25.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F25)

Figure 25. The outbreak prediction for Germany through 75 days

![Figure 26.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.17.20070094/F26.medium.gif)

[Figure 26.](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/F26)

Figure 26. The outbreak prediction for the USA through 75 days

## 4. Discussion

The parameters of several simple mathematical models (i.e., logistic, linear, logarithmic, quadratic, cubic, compound, power and exponential) were fitted using GA, PSO, and GWO. The logistic model outperformed other methods and showed promising results on training for 30 days. Extrapolation of the prediction beyond the original observation range of 30-days should not be expected to be realistic considering the new statistics. The fitted models generally showed low accuracy and also weak generalization ability for the five countries. Although the prediction for China was promising, the model was insufficient for extrapolation, as expected. In turn, the logistic GWO outperformed the PSO and GA and the computational cost for GWO was reported as satisfactory. Consequently, for further assessment of the ML models, the logistic model fitted with GWO was used for comparative analysis.

In the next step, for introducing the machine learning methods for time-series prediction, two scenarios were proposed. Scenario 1 considered four data samples from the progress of the infection from previous days, as reported in table 3. The sampling for data processing was done weekly for scenario 1. However, scenario 2 was devoted to daily sampling for all previous consecutive days. Providing these two scenarios expanded the scope of this study. Training and test results for the two machine learning models (MLP and ANFIS) were considered for the two scenarios. A detailed investigation was also carried out to explore the most suitable number of neurons. For the MLP, the performances of using 8, 12 and 16 neurons were analyzed throughout the study. For the ANFIS, the membership function (MF) types of Tri, Trap, and Gauss were analyzed throughout the study. The five counties of Italy, China, Iran, Germany, and USA were considered. The performance of both ML models for these countries varied amongst the two different scenarios. Given the observed results, it is not possible to select the most suitable scenario. Therefore, both daily and weekly sampling can be used in machine learning modeling. Comparison between analytical and machine learning models using the deviation from the target value (figures 17 to 21) indicated that the MLP in both scenarios delivered the most accurate results. Extrapolation for long-term prediction of up to 150 days using the ML models was tested. The actual prediction of MLP and ANFIS for the five countries was reported which showed the progression of the outbreak.

## 5. Conclusions

The global pandemic of the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) has become the primary national security issue of many nations. Advancement of accurate prediction models for the outbreak is essential to provide insights into the spread and consequences of this infectious disease. Due to the high level of uncertainty and lack of crucial data, standard epidemiological models have shown low accuracy for long-term prediction. This paper presents a comparative analysis of ML and soft computing models to predict the COVID-19 outbreak. The results of two ML models (MLP and ANFIS) reported a high generalization ability for long-term prediction. With respect to the results reported in this paper and due to the highly complex nature of the COVID-19 outbreak and differences from nation-to-nation, this study suggests ML as an effective tool to model the outbreak.

For the advancement of higher performance models for long-term prediction, future research should be devoted to comparative studies on various ML models for individual countries. Due to the fundamental differences between the outbreak in various countries, advancement of global models with generalization ability would not be feasible. As observed and reported in many studies, it is unlikely that an individual outbreak will be replicated elsewhere [1].

Although the most difficult prediction is to estimate the maximum number of infected patients, estimation of the *n*(deaths) / *n*(infecteds) is also essential. The mortality rate is particularly important to accurately estimate the number of patients and the required beds in intensive care units. For future research, modeling the mortality rate would be of the utmost importance for nations to plan for new facilities.

## Nomenclatures

View this table:
[Table24](http://medrxiv.org/content/early/2020/04/22/2020.04.17.20070094/T24)

## Data Availability

Data is available on request.

[https://www.worldometers.info/coronavirus/country](https://www.worldometers.info/coronavirus/country) 

## Conflicts of Interest

The authors declare no conflict of interest.

## Footnotes

*   s.ardabili{at}ieee.org

*   filipf{at}ujs.sk

*   p.ghamisi{at}hzdr.de

*   pma{at}lancaster.ac.uk

*   Received April 17, 2020.
*   Revision received April 17, 2020.
*   Accepted April 22, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Remuzzi, A.; Remuzzi, G. COVID-19 and Italy: what next? Lancet 2020.
    
    

2.  2.Ivanov, D. Predicting the impacts of epidemic outbreaks on global supply chains: A simulation-based analysis on the coronavirus outbreak (COVID-19/SARS-CoV-2) case. Transp. Res. Part E Logist. Transp. Rev. 2020, 136, doi:10.1016/j.tre.2020.101922.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tre.2020.101922&link_type=DOI) 

3.  3.Koolhof, I.S.; Gibney, K.B.; Bettiol, S.; Charleston, M.; Wiethoelter, A.; Arnold, A.L.; Campbell, P.T.; Neville, P.J.; Aung, P.; Shiga, T.; et al. The forecasting of dynamical Ross River virus outbreaks: Victoria, Australia. Epidemics 2020, 30, doi:10.1016/j.epidem.2019.100377.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2019.100377&link_type=DOI) 

4.  4.Darwish, A.; Rahhal, Y.; Jafar, A. A comparative study on predicting influenza outbreaks using different feature spaces: application of influenza-like illness data from Early Warning Alert and Response System in Syria. BMC Res. Notes 2020, 13, 33, doi:10.1186/s13104-020-4889-5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13104-020-4889-5&link_type=DOI) 

5.  5.Rypdal, M.; Sugihara, G. Inter-outbreak stability reflects the size of the susceptible pool and forecasts magnitudes of seasonal epidemics. Nat. Commun. 2019, 10, doi:10.1038/s41467-019-10099-y.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-10099-y&link_type=DOI) 

6.  6.Scarpino, S.V.; Petri, G. On the predictability of infectious disease outbreaks. Nat. Commun. 2019, 10, doi:10.1038/s41467-019-08616-0.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-08616-0&link_type=DOI) 

7.  7.Zhan, Z.; Dong, W.; Lu, Y.; Yang, P.; Wang, Q.; Jia, P. Real-Time Forecasting of Hand-Foot-and-Mouth Disease Outbreaks using the Integrating Compartment Model and Assimilation Filtering. Sci. Rep. 2019, 9, doi:10.1038/s41598-019-38930-y.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-019-38930-y&link_type=DOI) 

8.  8.Koike, F.; Morimoto, N. Supervised forecasting of the range expansion of novel non-indigenous organisms: Alien pest organisms and the 2009 H1N1 flu pandemic. Global Ecol. Biogeogr. 2018, 27, 991–1000, doi:10.1111/geb.12754.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/geb.12754&link_type=DOI) 

9.  9.Dallas, T.A.; Carlson, C.J.; Poisot, T. Testing predictability of disease outbreaks with a simple model of pathogen biogeography. R. Soc. Open Sci. 2019, 6, doi:10.1098/rsos.190883.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rsos.190883&link_type=DOI) 

10. 10.de Groot, M.; Ogris, N. Short-term forecasting of bark beetle outbreaks on two economically important conifer tree species. For. Ecol. Manage. 2019, 450, doi:10.1016/j.foreco.2019.117495.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.foreco.2019.117495&link_type=DOI) 

11. 11.Kelly, J.D.; Park, J.; Harrigan, R.J.; Hoff, N.A.; Lee, S.D.; Wannier, R.; Selo, B.; Mossoko, M.; Njoloko, B.; Okitolonda-Wemakoy, E.; et al. Real-time predictions of the 2018–2019 Ebola virus disease outbreak in the Democratic Republic of the Congo using Hawkes point process models. Epidemics 2019, 28, doi:10.1016/j.epidem.2019.100354.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2019.100354&link_type=DOI) 

12. 12.Maier, B.F.; Brockmann, D. Effective containment explains sub-exponential growth in confirmed cases of recent COVID-19 outbreak in Mainland China. medRxiv 2020, 10.1101/2020.02.18.20024414, 2020.2002.2018.20024414, doi:10.1101/2020.02.18.20024414.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMi4xOC4yMDAyNDQxNHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMjIvMjAyMC4wNC4xNy4yMDA3MDA5NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

13. 13.Pan, J.R.; Huang, Z.Q.; Chen, K. [Evaluation of the effect of varicella outbreak control measures through a discrete time delay SEIR model]. Zhonghua Yu Fang Yi Xue Za Zhi 2012, 46, 343–347.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22800634&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.17.20070094.atom) 

14. 14.Zha, W.T.; Pang, F.R.; Zhou, N.; Wu, B.; Liu, Y.; Du, Y.B.; Hong, X.Q.; Lv, Y. Research about the optimal strategies for prevention and control of varicella outbreak in a school in a central city of China: Based on an SEIR dynamic model. Epidemiol. Infect. 2020, 10.1017/S0950268819002188, doi:10.1017/S0950268819002188.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268819002188&link_type=DOI) 

15. 15.Dantas, E.; Tosin, M.; Cunha, A.; Jr.. Calibration of a SEIR–SEI epidemic model to describe the Zika virus outbreak in Brazil. Appl. Math. Comput. 2018, 338, 249-259, doi:10.1016/j.amc.2018.06.024.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.amc.2018.06.024&link_type=DOI) 

16. 16.Leonenko, V.N.; Ivanov, S.V. Fitting the SEIR model of seasonal influenza outbreak to the incidence data for Russian cities. Russ J Numer Anal Math Modell 2016, 31, 267-279, doi:10.1515/rnam-2016-0026.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1515/rnam-2016-0026&link_type=DOI) 

17. 17.Imran, M.; Usman, M.; Dur-e-Ahmad, M.; Khan, A. Transmission Dynamics of Zika Fever: A SEIR Based Model. Differ. Equ. Dyn. Syst. 2020, 10.1007/s12591-017-0374-6, doi:10.1007/s12591-017-0374-6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12591-017-0374-6&link_type=DOI) 

18. 18.Miranda, G.H.B.; Baetens, J.M.; Bossuyt, N.; Bruno, O.M.; De Baets, B. Real-time prediction of influenza outbreaks in Belgium. Epidemics 2019, 28, doi:10.1016/j.epidem.2019.04.001.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2019.04.001&link_type=DOI) 

19. 19.Sinclair, D.R.; Grefenstette, J.J.; Krauland, M.G.; Galloway, D.D.; Frankeny, R.J.; Travis, C.; Burke, D.S.; Roberts, M.S. Forecasted Size of Measles Outbreaks Associated With Vaccination Exemptions for Schoolchildren. JAMA Netw Open 2019, 2, e199768, doi:10.1001/jamanetworkopen.2019.9768.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2019.9768&link_type=DOI) 

20. 20.Zhao, S.; Musa, S.S.; Fu, H.; He, D.; Qin, J. Simple framework for real-time forecast in a data-limited situation: The Zika virus (ZIKV) outbreaks in Brazil from 2015 to 2016 as an example. Parasites Vectors 2019, 12, doi:10.1186/s13071-019-3602-9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13071-019-3602-9&link_type=DOI) 

21. 21.Fast, S.M.; Kim, L.; Cohn, E.L.; Mekaru, S.R.; Brownstein, J.S.; Markuzon, N. Predicting social response to infectious disease outbreaks from internet-based news streams. Ann. Oper. Res. 2018, 263, 551–564, doi:10.1007/s10479-017-2480-9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10479-017-2480-9&link_type=DOI) 

22. 22.McCabe, C.M.; Nunn, C.L. Effective network size predicted from simulations of pathogen outbreaks through social networks provides a novel measure of structure-standardized group size. Front. Vet. Sci. 2018, 5, doi:10.3389/fvets.2018.00071.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fvets.2018.00071&link_type=DOI) 

23. 23.Bragazzi, N.L.; Mahroum, N. Google trends predicts present and future plague cases during the plague outbreak in Madagascar: Infodemiological study. J. Med. Internet Res. 2019, 21, doi:10.2196/13142.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2196/13142&link_type=DOI) 

24. 24.Jain, R.; Sontisirikit, S.; Iamsirithaworn, S.; Prendinger, H. Prediction of dengue outbreaks based on disease surveillance, meteorological and socio-economic data. BMC Infect. Dis. 2019, 19, doi:10.1186/s12879-019-3874-x.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-019-3874-x&link_type=DOI) 

25. 25.Kim, T.H.; Hong, K.J.; Shin, S.D.; Park, G.J.; Kim, S.; Hong, N. Forecasting respiratory infectious outbreaks using ED-based syndromic surveillance for febrile ED visits in a Metropolitan City. Am. J. Emerg. Med. 2019, 37, 183–188, doi:10.1016/j.ajem.2018.05.007.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajem.2018.05.007&link_type=DOI) 

26. 26.Zhong, L.; Mu, L.; Li, J.; Wang, J.; Yin, Z.; Liu, D. Early Prediction of the 2019 Novel Coronavirus Outbreak in the Mainland China Based on Simple Mathematical Model. IEEE Access 2020, 8, 51761–51769, doi:10.1109/ACCESS.2020.2979599.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/ACCESS.2020.2979599&link_type=DOI) 

27. 27.Wu, J.T.; Leung, K.; Leung, G.M. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. Lancet 2020, 395, 689–697, doi:10.1016/S0140-6736(20)30260-9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30260-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32014114&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.17.20070094.atom) 

28. 28.Reis, J.; Yamana, T.; Kandula, S.; Shaman, J. Superensemble forecast of respiratory syncytial virus outbreaks at national, regional, and state levels in the United States. Epidemics 2019, 26, 1–8, doi:10.1016/j.epidem.2018.07.001.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2018.07.001&link_type=DOI) 

29. 29.Burke, R.M.; Shah, M.P.; Wikswo, M.E.; Barclay, L.; Kambhampati, A.; Marsh, Z.; Cannon, J.L.; Parashar, U.D.; Vinjé, J.; Hall, A.J. The Norovirus Epidemiologic Triad: Predictors of Severe Outcomes in US Norovirus Outbreaks, 2009-2016. J. Infect. Dis. 2019, 219, 1364–1372, doi:10.1093/infdis/jiy569.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiy569&link_type=DOI) 

30. 30.Carlson, C.J.; Dougherty, E.; Boots, M.; Getz, W.; Ryan, S.J. Consensus and conflict among ecological forecasts of Zika virus outbreaks in the United States. Sci. Rep. 2018, 8, doi:10.1038/s41598-018-22989-0.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-018-22989-0&link_type=DOI) 

31. 31.Kleiven, E.F.; Henden, J.A.; Ims, R.A.; Yoccoz, N.G. Seasonal difference in temporal transferability of an ecological model: near-term predictions of lemming outbreak abundances. Sci. Rep. 2018, 8, doi:10.1038/s41598-018-33443-6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-018-33443-6&link_type=DOI) 

32. 32.Rivers-Moore, N.A.; Hill, T.R. A predictive management tool for blackfly outbreaks on the Orange River, South Africa. River Res. Appl. 2018, 34, 1197–1207, doi:10.1002/rra.3357.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/rra.3357&link_type=DOI) 

33. 33.Yin, R.; Tran, V.H.; Zhou, X.; Zheng, J.; Kwoh, C.K. Predicting antigenic variants of H1N1 influenza virus based on epidemics and pandemics using a stacking model. PLoS ONE 2018, 13, doi:10.1371/journal.pone.0207777.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0207777&link_type=DOI) 

34. 34.Agarwal, N.; Koti, S.R.; Saran, S.; Senthil Kumar, A. Data mining techniques for predicting dengue outbreak in geospatial domain using weather parameters for New Delhi, India. Curr. Sci. 2018, 114, 2281–2291, doi:10.18520/cs/v114/i11/2281-2291.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18520/cs/v114/i11/2281-2291&link_type=DOI) 

35. 35.Anno, S.; Hara, T.; Kai, H.; Lee, M.A.; Chang, Y.; Oyoshi, K.; Mizukami, Y.; Tadono, T. Spatiotemporal dengue fever hotspots associated with climatic factors in taiwan including outbreak predictions based on machine-learning. Geospatial Health 2019, 14, 183–194, doi:10.4081/gh.2019.771.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4081/gh.2019.771&link_type=DOI) 

36. 36.Chenar, S.S.; Deng, Z. Development of artificial intelligence approach to forecasting oyster norovirus outbreaks along Gulf of Mexico coast. Environ. Int. 2018, 111, 212–223, doi:10.1016/j.envint.2017.11.032.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2017.11.032&link_type=DOI) 

37. 37.Chenar, S.S.; Deng, Z. Development of genetic programming-based model for predicting oyster norovirus outbreak risks. Water Res. 2018, 128, 20–37, doi:10.1016/j.watres.2017.10.032.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.watres.2017.10.032&link_type=DOI) 

38. 38.Iqbal, N.; Islam, M. Machine learning for dengue outbreak prediction: A performance evaluation of different prominent classifiers. Informatica 2019, 43, 363–371, doi:10.31449/inf.v43i1.1548.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.31449/inf.v43i1.1548&link_type=DOI) 

39. 39.Liang, R.; Lu, Y.; Qu, X.; Su, Q.; Li, C.; Xia, S.; Liu, Y.; Zhang, Q.; Cao, X.; Chen, Q.; et al. Prediction for global African swine fever outbreaks based on a combination of random forest algorithms and meteorological data. Transboundary Emer. Dis. 2020, 67, 935–946, doi:10.1111/tbed.13424.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/tbed.13424&link_type=DOI) 

40. 40.Mezzatesta, S.; Torino, C.; De Meo, P.; Fiumara, G.; Vilasi, A. A machine learning-based approach for predicting the outbreak of cardiovascular diseases in patients on dialysis. Comput. Methods Programs Biomed. 2019, 177, 9–15, doi:10.1016/j.cmpb.2019.05.005.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmpb.2019.05.005&link_type=DOI) 

41. 41.Raja, D.B.; Mallol, R.; Ting, C.Y.; Kamaludin, F.; Ahmad, R.; Ismail, S.; Jayaraj, V.J.; Sundram, B.M. Artificial intelligence model as predictor for dengue outbreaks. Malays. J. Public Health Med. 2019, 19, 103–108.
    
    

42. 42.Tapak, L.; Hamidi, O.; Fathian, M.; Karami, M. Comparative evaluation of time series models for predicting influenza outbreaks: Application of influenza-like illness data from sentinel sites of healthcare centers in Iran. BMC Res. Notes 2019, 12, doi:10.1186/s13104-019-4393-y.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13104-019-4393-y&link_type=DOI) 

43. 43.Titus Muurlink, O.; Stephenson, P.; Islam, M.Z.; Taylor-Robinson, A.W. Long-term predictors of dengue outbreaks in Bangladesh: A data mining approach. Infect. Dis. Modelling 2018, 3, 322–330, doi:10.1016/j.idm.2018.11.004.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.idm.2018.11.004&link_type=DOI) 

44. 44.Nádai, L.; Imre, F.; Ardabili, S.; Gundoshmian, T.M.; Gergo, P. and Mosavi, A.; 2020. Performance Analysis of Combine Harvester using Hybrid Model of Artificial Neural Networks Particle Swarm Optimization. arXiv preprint arxiv:2002.11041.
    
    

45. 45.Ardabili, S.; Beszédes, B.; Nádai, L.; Széll, K.; Mosavi, A. and Imre, F.; 2020. Comparative Analysis of Single and Hybrid Neuro-Fuzzy-Based Models for an Industrial Heating Ventilation and Air Conditioning Control System. arXiv preprint arxiv:2002.11042.
    
    

46. 46.Back, T. Evolutionary algorithms in theory and practice: evolution strategies, evolutionary programming, genetic algorithms; Oxford university press: 1996.
    
    

47. 47.Deb, K. Multi-objective optimization using evolutionary algorithms; John Wiley & Sons: 2001; Vol. 16.
    
    

48. 48.Zitzler, E.; Deb, K.; Thiele, L. Comparison of multiobjective evolutionary algorithms: Empirical results. Evolutionary computation 2000, 8, 173–195.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1162/106365600568202&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10843520&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.17.20070094.atom) 

49. 49.Ghamisi, P.; Benediktsson, J.A. Feature selection based on hybridization of genetic algorithm and particle swarm optimization. IEEE Geoscience and remote sensing letters 2014, 12, 309–313.
    
    

50. 50.Ardabili, S.; Beszédes, B.; Nádai, L.; Széll, K.; Mosavi, A. and Imre, F.; 2020. Comparative Analysis of Single and Hybrid Neuro-Fuzzy-Based Models for an Industrial Heating Ventilation and Air Conditioning Control System. arXiv preprint arxiv:2002.11042.
    
    

51. 51.Ardabili, S.; Mosavi, A. and Várkonyi-Kóczy, A.R.; 2019, September. Advances in machine learning modeling reviewing hybrid and ensemble methods. In International Conference on Global Research and Education (pp. 215–227). Springer, Cham.
    
    

52. 52.Mühlenbein, H.; Schomisch, M.; Born, J. The parallel genetic algorithm as function optimizer. Parallel computing 1991, 17, 619–632.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0167-8191(05)80052-3&link_type=DOI) 

53. 53.Whitley, D. A genetic algorithm tutorial. Statistics and computing 1994, 4, 65–85.
    
    

54. 54.Horn, J.; Nafpliotis, N.; Goldberg, D.E. A niched Pareto genetic algorithm for multiobjective optimization. In Proceedings of Proceedings of the first IEEE conference on evolutionary computation. IEEE world congress on computational intelligence; pp. 82–87.
    
    

55. 55.Reeves, C.R. A genetic algorithm for flowshop sequencing. Computers & operations research 1995, 22, 5–13.
    
    

56. 56.Jones, G.; Willett, P.; Glen, R.C.; Leach, A.R.; Taylor, R. Development and validation of a genetic algorithm for flexible docking. Journal of molecular biology 1997, 267, 727–748.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/jmbi.1996.0897&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9126849&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.17.20070094.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WT33400020&link_type=ISI) 

57. 57.Ardabili, S.; Mosavi, A.; Varkonyi-Koczy, A.R. Advances in machine learning modeling reviewing hybrid and ensemble methods. 2019.
    
    

58. 58.Houck, C.R.; Joines, J.; Kay, M.G. A genetic algorithm for function optimization: a Matlab implementation. Ncsu-ie tr 1995, 95, 1–10.
    
    

59. 59.Whitley, D.; Starkweather, T.; Bogart, C. Genetic algorithms and neural networks: Optimizing connections and connectivity. Parallel computing 1990, 14, 347–361.
    
    

60. 60.Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of Proceedings of ICNN’95- International Conference on Neural Networks; pp. 1942–1948.
    
    

61. 61.Poli, R.; Kennedy, J.; Blackwell, T. Particle swarm optimization. Swarm intelligence 2007, 1, 33–57.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11721-007-0002-0&link_type=DOI) 

62. 62.Clerc, M. Particle swarm optimization; John Wiley & Sons: 2010; Vol. 93.
    
    

63. 63.Angeline, P.J. Using selection to improve particle swarm optimization. In Proceedings of 1998 IEEE International Conference on Evolutionary Computation Proceedings. IEEE World Congress on Computational Intelligence (Cat. No. 98TH8360); pp. 84–89.
    
    

64. 64.Sun, J.; Feng, B.; Xu, W. Particle swarm optimization with particles having quantum behavior. In Proceedings of Proceedings of the 2004 congress on evolutionary computation (IEEE Cat. No. 04TH8753); pp. 325–331.
    
    

65. 65.Parsopoulos, K.E.; Vrahatis, M.N. Recent approaches to global optimization problems through particle swarm optimization. Natural computing 2002, 1, 235–306.
    
    

66. 66.Bai, Q. Analysis of particle swarm optimization algorithm. Computer and information science 2010, 3, 180.
    
    

67. 67.Parsopoulos, K.E.; Vrahatis, M.N. Particle swarm optimization method in multiobjective problems. In Proceedings of Proceedings of the 2002 ACM symposium on Applied computing; pp. 603–607.
    
    

68. 68.Ghamisi, P.; Couceiro, M.S.; Benediktsson, J.A.; Ferreira, N.M. An efficient method for segmentation of images based on fractional calculus and natural selection. Expert Systems with Applications 2012, 39, 12407–12417.
    
    

69. 69.Ghamisi, P.; Couceiro, M.S.; Martins, F.M.; Benediktsson, J.A. Multilevel image segmentation based on fractional-order Darwinian particle swarm optimization. IEEE Transactions on Geoscience and Remote sensing 2013, 52, 2382–2394.
    
    

70. 70.Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Advances in engineering software 2014, 69, 46–61.
    
    

71. 71.Ardabili, S.; Mosavi, A.; Varkonyi-Koczy, A. Building Energy information: demand and consumption prediction with Machine Learning models for sustainable and smart cities. Engineering for Sustainable Future, Lecture Notes in Networks and Systems, Springer Nature Switzerland.: 2019.
    
    

72. 72.Ardabili;, S.; Mosavi;, A.; Varkonyi-Koczy;, A. Systematic Review of Deep Learning and Machine Learning Models in Biofuels Research,. Engineering for Sustainable Future, Lecture Notes in Networks and Systems, Springer Nature Switzerland 2019.
    
    

73. 73.Gundoshmian, T.M.; Ardabili, S.; Mosavi, A.; Varkonyi-Koczy, A.R. Prediction of combine harvester performance using hybrid machine learning modeling and response surface methodology. 2019.
    
    

74. 74.Gundoshmian, T.M.; Ardabili, S.; Mosavi, A. and Várkonyi-Kóczy, A.R.; 2019, September. Prediction of combine harvester performance using hybrid machine learning modeling and response surface methodology. In International Conference on Global Research and Education (pp. 345–360). Springer, Cham.
    
    

75. 75.Hamilton, J.D. Time series analysis; Princeton university press Princeton, NJ: 1994; Vol. 2.
    
    

76. 76.Ardabili, S.; Mosavi, A.; Mahmoudi, A.; Gundoshmian, T.M.; Nosratabadi, S.; Varkonyi-Koczy, A.R. Modelling Temperature Variation of Mushroom Growing Hall Using Artificial Neural Networks. Engineering for Sustainable Future, Lecture Notes in Networks and Systems, Springer Nature Switzerland. 2019.
    
    

77. 77.Khosravi, A.; Koury, R.; Machado, L.; Pabon, J. Prediction of wind speed and wind direction using artificial neural network, support vector regression and adaptive neuro-fuzzy inference system. Sustainable Energy Technologies and Assessments 2018, 25, 146–160.
    
    

78. 78.Amid, S.; Mesri Gundoshmian, T. Prediction of output energies for broiler production using linear regression, ANN (MLP, RBF), and ANFIS models. Environmental Progress & Sustainable Energy 2017, 36, 577–585.
    
    

79. 79.Hossain, M.A.; Ayodele, B.V.; Cheng, C.K.; Khan, M.R. Artificial neural network modeling of hydrogen-rich syngas production from methane dry reforming over novel Ni/CaFe 2 O 4 catalysts. International Journal of Hydrogen Energy 2016, 41, 11119–11130.
    
    

80. 80.Taghavifar, H.; Mardani, A. Wavelet neural network applied for prognostication of contact pressure between soil and driving wheel. Information Processing in Agriculture 2014, 1, 51–56.
    
    

81. 81.Sharabiani, V.; Kassar, F.; Gilandeh, Y.; Ardabili, S. APPLICATION OF SOFT COMPUTING METHODS AND SPECTRAL REFLECTANCE DATA FOR WHEAT GROWTH MONITORING. The Iraqi Journal of Agricultural Science 2019, 50, 1064–1076.
    
    

82. 82.Mirzazadeh, A.; Abdollahpour, S.; Mahmoudi, A.; Bukat, A. Intelligent modeling of material separation in combine harvester’s thresher by ANN. International Journal of Agriculture and Crop Sciences 2012, 4, 1767–1777.
    
    

83. 83.Khalesi, S.; Mahmoudi, A.; Hosainpour, A.; Alipour, A. Detection of walnut varieties using impact acoustics and artificial neural networks (ANNs). Modern Applied Science 2012, 6, 43.
    
    

84. 84.Reshadsedghi, A.; Mahmoudi, A. Detection of almond varieties using impact acoustics and artificial neural networks. International Journal of Agriculture and Crop Sciences 2013, 6, 1008.
    
    

85. 85.Hassoun, M.H. Fundamentals of artificial neural networks; MIT press: 1995.
    
    

86. 86.Jang, J.-S. ANFIS: adaptive-network-based fuzzy inference system. IEEE transactions on systems, man, and cybernetics 1993, 23, 665–685.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/21.256541&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1993MF08900003&link_type=ISI) 

87. 87.Faizollahzadeh_Ardabili; Mahmoudi A., Mesri Gundoshmian, T. Modeling and simulation controlling system of HVAC using fuzzy and predictive (radial basis function, RBF) controllers. Journal of Building Engineering 2016, 6, 301–308.
    
    

88. 88.Amid, S.; Mesri Gundoshmian, T. Prediction of output energy based on different energy inputs on broiler production using application of adaptive neural-fuzzy inference system. Agri. Sci. Dev 2016, 5, 14–21.
    
    

89. 89.Ardabili, S.; Mosavi, A.; Dehghani, M.; Várkonyi-Kóczy, A.R. Deep learning and machine learning in hydrological processes climate change and earth systems a systematic review. In Proceedings of International Conference on Global Research and Education; pp. 52–62.
    
    

90. 90.Nosratabadi, S.; Mosavi, A.; Keivani, R.; Ardabili, S.; Aram, F. State of the art survey of deep learning and machine learning models for smart cities and urban sustainability. In Proceedings of International Conference on Global Research and Education; pp. 228–238.
    
    

91. 91.Güler, I.; Übeyli, E.D. Adaptive neuro-fuzzy inference system for classification of EEG signals using wavelet coefficients. Journal of neuroscience methods 2005, 148, 113–121.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneumeth.2005.04.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16054702&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.17.20070094.atom) 

92. 92.Chang, J.-R.; Wei, L.-Y.; Cheng, C.-H. A hybrid ANFIS model based on AR and volatility for TAIEX forecasting. Applied Soft Computing 2011, 11, 1388–1395.
    
    

93. 93.Chang, F.-J.; Chang, Y.-T. Adaptive neuro-fuzzy inference system for prediction of water level in reservoir. Advances in water resources 2006, 29, 1–10.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.advwatres.2005.04.015&link_type=DOI) 

94. 94.Moghaddamnia, A.; Gousheh, M.G.; Piri, J.; Amin, S.; Han, D. Evaporation estimation using artificial neural networks and adaptive neuro-fuzzy inference system techniques. Advances in Water Resources 2009, 32, 88–97.
    
    

95. 95.Çaydaş, U.; Hasçalik, A.; Ekici, S. An adaptive neuro-fuzzy inference system (ANFIS) model for wire- EDM. Expert Systems with Applications 2009, 36, 6135–6139.
    
    

96. 96.Yadav, D.; Chhabra, D.; Gupta, R.K.; Phogat, A.; Ahlawat, A. Modeling and analysis of significant process parameters of FDM 3D printer using ANFIS. Materials Today: Proceedings 2020, 21, 1592–1604.
    
    

97. 97.Faizollahzadeh_Ardabili; Mahmoudi A., Gundoshmian, T.M.; Roshanianfard, A. Modeling and comparison of fuzzy and on/off controller in a mushroom growing hall. Measurement 2016, 90, 127–134.

 [1]: /embed/graphic-1.gif
 [2]: /embed/inline-graphic-1.gif
 [3]: /embed/graphic-2.gif
 [4]: /embed/graphic-3.gif
 [5]: /embed/graphic-4.gif
 [6]: /embed/graphic-6.gif
 [7]: /embed/graphic-11.gif
 [8]: /embed/graphic-12.gif
 [9]: /embed/graphic-13.gif
 [10]: /embed/graphic-14.gif
 [11]: /embed/graphic-15.gif
 [12]: /embed/inline-graphic-2.gif
 [13]: /embed/inline-graphic-3.gif
 [14]: /embed/inline-graphic-4.gif
 [15]: /embed/inline-graphic-5.gif
 [16]: /embed/inline-graphic-6.gif
 [17]: /embed/inline-graphic-7.gif
 [18]: /embed/graphic-18.gif
 [19]: /embed/graphic-19.gif
 [20]: /embed/graphic-20.gif