Ebola is an infectious virus that causes Ebola hemorrhagic fever in primates and humans, which was first found in 1976. The Ebola virus outbreak in West Africa in 2014 was the largest ever. A lot of researchers use mathematical models to analyze the characteristics of infectious diseases. However, many parameters in the model cannot be estimated completely. To ease the difficulty, we proposed an approach to estimate the parameter based on genetic algorithm (GA). GA uses the natural selection method of the fittest to find the optimal solution of the model. The least residual squares sum is used as fitness function to measure the performance of GA in parameter estimation. Moreover, we used a dynamical model and the real data of Ebola in Sierra Leone to verify the validity of GA. The experimental results indicate that the GA has strong competitiveness compared with the classical method, and it is a feasible method for estimating the parameters of infectious disease models.
Program for the Outstanding Innovative Teams (OIT) of Higher Learning Institutions of ShanxiNatural Science Foundation of Shanxi Province201801D221003China Postdoctoral Science Foundation2017M6211101. Introduction
Ebola virus belongs to the family Filoviridae and is considered a prototype pathogen of viral hemorrhagic fever [1]. The virus was first detected in the Ebola river basin in southern Sudan and Congo in 1976 [2–9]. Since the discovery of Ebola virus, only four species of this virus cause human disease, namely, Zaire ebolavirus, Tai Forest ebolavirus, Sudan ebolavirus, and Bundibugyo ebolavirus [10]. The Reston virus causes only animal disease, not human disease. Therefore, the source of the Ebola virus is unknown. Researchers found evidence of asymptomatic infection of Ebola virus in three species of the fruit bats, which suggested that the bats are most likely to be the source of the deadly virus [11]. The bats could carry Ebola virus to other animals and even humans [12–14].
Ebola virus is transmitted through the saliva, the urine, and other body fluids [15, 16]. People can cause infection by direct contact with body fluids which carry the virus, with the virus entering the body through the nose, the mouth, the eyes, and the damaged skin [17]. Humans become infected after contact with the blood, the body fluids, and the infected fruit bats, as well as through the sexual contact [18].
Since there were no good treatments and approved vaccines at the time, the management of Ebola virus was limited to the use of obstacles and palliative care to suppress transmission [19]. A large-scale Ebola outbreak occurred in West Africa in 2014, mainly in Guinea, Liberia, and Sierra Leone. The number of confirmed cases is far greater than that in the past [10]. The lack of effective preventive measures at the time resulted in more people being infected with the Ebola virus [20]. In [21], the authors investigated the effectiveness of small interfering RNAs treatments for Ebola-infected patients. RNA interference can suppress the expression of viral genes; thus it is effective in suppressing Ebola virus replication, and the authors developed monoclonal antibodies against Ebola glycoprotein for the treatment of Ebola-virus-infected people [22]. In addition, some researchers used Sierra Leone’s disease data to study mathematical models of Ebola virus, predict the progress of the epidemic, and propose preventive control measures and recommendations [5, 23–25].
Research on infectious diseases using dynamic models has become one of the important methods [26–32]. The propagation coefficient of the disease in the model affects the prediction results directly, and, consequently, it is important to estimate the propagation coefficient correctly. Classical parameter estimation methods are the Markov Chain Monte Carlo (MCMC) method, the least-squares method, and so on. The basic principle of the MCMC method is to construct a Markov chain by using the joint posterior probability distribution of the model propagation coefficients and assign any initial value to the simulation until it converges to a stable distribution. This determines the propagation coefficient [33–35]. There are many improved MCMC methods, such as using sequential Monte Carlo (SMC) filter techniques to estimate the propagation coefficients in the model [36]. However, firstly, the traditional method is limited by the calculation cost of the high-dimensional nonlinear model, which may take a lot of calculation time; it is usually not easy to obtain high-precision results, and it is not possible to get all the propagation coefficients at once [37, 38]. Secondly, the numerical estimation of the marginal probability distribution is difficult to achieve in the high-dimensional inversion model [39]. The least-squares method is performed by convolving the simulated data with the real data [40]. Although it has low calculation cost and generality, it does not consider the uncertainty of the inverse problem solution, and the initial value of the propagation coefficient will affect the efficiency of the algorithm. The least-squares method has certain flaws in determining the initial value. If it is set close to the optimal propagation coefficient, the result will be obtained quickly. If it is set far from the optimal propagation coefficient, it will increase the time of the algorithm [41].
In this paper, we present a method to solve inverse problems of differential equations based on GA. The GA is a method that is widely used in parameter estimation and other fields, and it has been proven to be a reliable method for estimating parameters based on nonlinear functions [42, 43]. It has a powerful adaptive search technology and uses the natural selection method of the fittest to simulate the evolution process, and thereby it can effectively solve the optimization problem [44]. When searching in high-dimensional models, GAs are superior to the other traditional search techniques due to their simplicity, effectiveness, versatility, and robustness [45, 46]. We have used the GA of adaptive mutation operator to estimate the parameters of differential equations. The advantage of this method is that the parameters in the high-dimensional model can be completely estimated by a small amount of data and all parameters combinations can be quickly obtained in a limited evolution process. In addition, an effective combination of multiple propagation coefficients can be obtained by GA for reference in studying the propagation dynamics.
The remainder of this paper is organized as follows. In Section 2, we introduced the transmission dynamical model of Ebola in Sierra Leone and the theories and processes of GA. In Section 3, we estimated the values of the parameters in a dynamical model based on GA. What is more, we validate the accuracy of the experiment results. Finally, we give the discussion and conclusions in Section 4.(1)dSdt=−β1SISN+β2SIPN+βHSHN+βFSFN+λ1pIS,dEdt=β1SISN+β2SIPN+βHSHN+βFSFN−m1ηE−m21−ηE,dISdt=m21−ηE−λ1pIS−λ2q1−pIS−λ31−q1−pIS,dIPdt=λ31−q1−pIS+m1ηE−rhθIP−δ11−θrdIP−ri1−θ1−δ1IP,dHdt=rhθIP−rdhδ2H−rih1−δ2H,dFdt=λ2q1−pIS+δ11−θrdIP+rdhδ2H−rfF,dRdt=ri1−θ1−δ1IP+rih1−δ2H+rfF.
2. Main Method2.1. A Dynamical Model of Ebola Virus Transmission in Sierra Leone
The time series of Ebola-confirmed case reports were collected from the World Health Organization (WHO) and the Ministry of Health of Sierra Leone. The data include the Ebola outbreak in 14 regions of Sierra Leone, including the suspected cases IS, the probable cases IP, and the hospital-confirmed cases H, which are thought to represent the best available data of the Ebola epidemic. Due to the fact that hospitalization is a result of real infections, while suspected and probable cases may not be completely converted into hospitalized cases, it is more accurate to use hospitalization cases to indicate the actual number of the Ebola infections. We collected the newly infected cases for 34 weeks from May 19th, 2014, to January 11th, 2015. More detailed data can be found in [47].
We used a GA based on adaptive mutation operator. Think of the propagation coefficient of the Ebola virus model as a genetic target. It is binary-encoded, and then the genetic operators of random selection with elitism, multipoint crossover, and gene site mutation are used to simulate evolution to find the optimal solution. Set evaluation index for parameter genetic process as fitness function which is a sum of variances of fitted data and real data. We estimated the parameters based on the dynamical model established in [6]. Record the optimal parameter set of each generation and perform comparison with the optimal parameter set of the next generation, always save the optimal set, and wait until the evolution has completed obtaining all the parameters in the model.
Based on other literature analyses, this article divides Ebola virus transmission into seven categories, namely, susceptible S, exposed E, suspected individuals may be misdiagnosed IS, probable individuals IP, hospital-confirmed cases H, the individuals who may infect others at a funeral F, and removed R [6]. Figure 1 depicted the transmission mechanism of the Ebola virus.
Ebola’s transmission mechanism in Sierra Leone.
Consequently, we have used the following system of equation (1) to simulate the transmission dynamics of the Ebola virus in Sierra Leone and the biological meanings of parameters can be obtained in Table 1. We quantified the uncertainty of parameter estimates, and we give the 95% confidence intervals in Table 1.
Description of parameters in model (1).
Parameter
Description
Value
95% confidence intervals
β1
Transmission coefficient between communities S and Is
0.9162
[0.1173, 1.1247]
β2
Transmission coefficient between communities S and IP
0.8929
[0.1341, 1.1803]
βH
Transmission coefficient at the hospital
0.1674
[0.0623, 1.0713]
βF
Transmission coefficient during the funerals
0.2736
[0.0887, 1.0665]
p
Misdiagnosed proportion in IS
0.657
[0.1638, 0.6891]
q
Proportion of correctly diagnosed suspected cases that move to F
0.1447
[0.0531, 1.1805]
η
Proportion of exposed cases entering the IP compartment
0.619
[0.1008, 1.1189]
θ
Proportion of infected and hospitalized
0.0023
[0.1476, 0.3507]
λ1
Time from suspected to susceptible
1.59 days
[0.0641, 1.002]
λ2
Average time from suspected to death
4.88 days
[0.0418, 1.0139]
λ3
Time from suspected to probable
2.38 days
[0.0661, 1.0172]
m1
Time from exposed to probable
1.0007 days
[0.0611, 1.2091]
m2
Time from exposed to suspected
20.9 days
[0.1053, 1.2612]
rh
Time from probable to hospitalized
3.56 days
[0.2733, 0.5884]
rdh
Average time from hospitalization to death
1.2 days
[0.1109, 1.187]
ri
Time from probable to removed
5.69 days
[0.1221, 1.1318]
rih
Time from hospitalized to removed
2.56 days
[0.0844, 1.1501]
rf
Average time from death to burial
1.21 days
[0.1323, 0.9899]
rd
Time from probable to death
1.11 days
[0.0099, 1.183]
δ1
Mortality rate from probable to death
0.5141
[0.1219, 1.0691]
δ2
Mortality rate from hospitalized to death
0.6214
[0.0455, 1.1318]
N
Size of the Sierra Leone population
6348350
[47]
2.2. Genetic Algorithm
GA is a method for finding solutions based on biological evolution process [42]. The process includes random selection, crossover, and mutation for an individual with the best combination of genes. GA begins with initializing the propagation coefficient in model (1) for binary encoding. The encoding length is determined by parameter range and accuracy. We used par to represent the parameters set: par=β1,β2,βH,βF,λ1,p,m1,η,m2,λ2,q,λ3,rh,θ,δ1,rd,ri,rdh,δ2,rih,rf. Assuming that all propagation coefficients are within [0, 1], the accuracy is 4 digits after decimal point, and thus the coding length can be determined by the following formula:(2)213<rangemax−rangeminaccuracy=1−00.0001=10000<214.
Therefore, a parameter can be represented by 14 bits of binary, and par consisting of all parameters needs to be encoded with 14×21=294 bits of binary. We collected the disease status of Ebola for 34 weeks, and the detailed data are shown above, so we set the initial population of parameters to 34 parameters set. Afterwards, we need to determine the fitness function; we take the minimum residual sum of squares between the solution of the infected case in the model and the actual infected case: (3)fitness=∑t=1I˜t−It2,t=1,2,…,34,where I˜t represents the infected cases in the model at time t and It represents the actual infected cases at time t. We use fitness function to evaluate the initial population and give the initial fitness value of parameters.
GA mainly includes three genetic operators: selection, crossover, and mutation. Selection is to apply the selection operator to the group. The purpose of selection is to inherit the optimized parameters directly to the next generation or to generate a new par to the next generation through pairing and crossover. Selection operations are based on the fitness evaluation of parameters in the population. Here we adopt a random selection combined with elitism, which means that we will copy the parameters with higher fitness once, replace the ones with the least fitness, and retain the best parameters of each generation. This is the elite strategy; then we randomly select parameters for crossover and mutation.
Crossover operator plays a key role in GA. It is mainly divided into single-point crossover, two-point crossover, and multiple-point crossover. The commonly used is the single-point crossover; that is, a cross point is randomly set in the parameter string, when the intersection is performed, in front of this point or partial structures of two parameters sets exchange and thus two new parameters sets are generated. An example of a single-point crossover is shown in Figure 2. Because our parameter binary is too long, single-point crossover cannot meet our needs, so we chose multiple-point crossover to increase the diversity of parameters.
The two parameters sets A and B exchange the part after the cross point and reconstitute two new parameters sets.
The mutation operator is to change certain gene positions of parameter strings in the population; for example, 0 becomes 1, and 1 becomes 0. An example of mutation is shown in Figure 3. We adopted the gene locus mutation with adaptive mutation operator in GA; that is, each gene position was mutated with a certain mutation probability. Moreover, the mutation probability can be adaptively adjusted by the parameters set fitness. When the difference between parameters set fitness and the average fitness of the population is small, it means that parameters are close to each other, which is not conducive to the next crossover. Therefore, it is necessary to increase the mutation probability and reduce the mutation probability when the difference is large, because the mutation probability is usually between 0.001 and 0.1, and the probability of mutation is small, so that it is not easy to destroy the genes of the dominant parameters, and it can jump out when the algorithm falls into the local optimal solution. The mutation characteristic of GA can make the solution process randomly search the entire space where the solution may exist and ensure the diversity of the population, so the global optimal solution can be obtained to a certain extent.
Parameters sets gene mutations generate new parameters sets.
Next, we solve model (1), decode the parameters into decimalism, and substitute them into the model solution to get the estimated value of confirmed-infection cases, that is, hospitalization cases. Fitness function (3) is used to fit the real infected cases. Then we find the parameters with the highest fitness in this generation, that is, the parameters with the least error, and keep them in the next-generation genetic process. Afterwards, we cyclically execute selection, crossover, mutation, and evaluation of new parameters set until the maximum number of iterations is reached.
The general steps of using GA to estimate parameters of Ebola model are explained as follows:
Step 1: consider the parameter to be estimated as a gene chromosome, define the parameter using binary coding, and then initialize the population.
Step 2: assign a fitness value to each parameter using equation (3). Starting from the second generation of parameters, parameters are ranked from small to large according to fitness values, and the first two parameters sets with the greatest fitness are duplicated once to replace the two parameters sets with the smallest fitness.
Step 3 (selection process): add randomly initialized parameters to the population to increase population diversity, and then two parameters sets are randomly selected as paternal parameters sets.
Step 4 (crossover): two new offspring parameters are generated by crossing the two parents at multiple points.
Step 5 (mutation): use the gene locus mutation described above in combination with adaptive mutation operator.
Step 6: solve equation (1) using the ode function.
Step 7: convert the types of parameters sets from binary to decimal, substitute it into the solution of equation (1), obtain the predicted value of the diseased cases, and use equation (3) to evaluate the fitness value of the new parameters set to obtain an optimal parameters set.
Step 8: when the fitness of the new offspring produced by genetic manipulation is higher than that of the parent, the new parameters sets replace the parents and are inserted into the parent population for the next genetic manipulation. If the optimal individual remains unchanged for 30 consecutive times, multiple-point crossover and mutation are carried out.
Step 9: save the best parameters set of this generation.
Step 10: if the number of iterations is not reached, proceed with step 2.
The above steps are executed iteratively until the termination condition is reached. Parameter estimation is finished whenever the genetic operation completed. The parameter values to be used in the model are the optimal solution of the last generation of parameters, and the specific values are shown in Table 1. We can study the propagation dynamics and preventive measures of Ebola.
3. Main Results
There are many methods for parameter estimation, including Markov Chain Monte Carlo (MCMC) method and least-squares method. It is not easy to analyze and estimate the parameters in the infectious disease model because there are many parameters in the model which cannot be estimated fully. Consequently, we propose to estimate the parameters in model (1) using GA, which are described in Section 2. Algorithm 1 shows the scheme of GA used on parameter estimation of the Ebola model. In this algorithm, par denotes a set of parameters, constants a and b represent two parameters sets randomly selected from par set and crossed according to the crossover probability px to get a′ and b′. In mutation, a′ and b′ are changed to a˜ and b˜ by gene locus mutations, according to mutation probability pm. IC denotes the initial value of each variable in equation (1) and t indicates the time of virus transmission. “count” indicates that successive generations of optimal values have not changed.
Algorithm 1: Parameter estimation of GA.
Parameter length CL is determined according to parameter range.
Each parameter is represented by a chromosome of CL length to get a complete parameters set par=β1,β2,…,rf.
Initialize all parameters set, parc=β1,β2,…,rf,c=34.
Evaluate the fitness of each par using f.
for i=1,2,…, Constant do
ifi=2do
Sort par allocation fitness in an ascending order.
par34=β1,β2,…,rf=par1=β1,β2,…,rf.
par33=β1,β2,…,rf=par2=β1,β2,…,rf.
end if
Randomly add an initialization parameter to replace one of the par
if max fpar close to average fpardo
Mutation probability pm=0.1.
end if
Randomly select a,b∈para≠b.
Cross over a and b according to px, a,b⟹a′,b′.
Perform mutation on a′ and b′ according to pm,a′,b′⟹a˜,b˜.
D=ebolaSolpar,IC,t.
Evaluate fitness f=sumD−Data2.
iffa˜<fado
para=para˜.
end if
ifcount>30do
Multiple crossover and variation.
end if
Find the best parameters set in this generation bestpi that satisfies bestpi=minfpar.
end for
Find the best parameters set parbest that satisfies parbest=bestpConstant.
In this study, we used different genetic operators and fitness functions (3) to conduct data fitting for the real hospital-diagnosed cases. We conducted 70 experiments and selected a set of parameters which performed well. The results are listed in Table 1. After experimental verification, we chose the crossover probability px to be 0.8 and the initial mutation probability pm to be 0.01. The genetic algebra Constant is 3000 times. The fitting result for the cumulative number of cases is shown in Figure 4.
The fitting result for actual data of Ebola virus disease in Sierra Leone for 34 weeks with the deterministic model (1).
Figure 5 represents the evolutionary process of the optimal value of each generation in GA. With the increase of genetic algebra, the error between the model solution and the real data is gradually decreasing, which means that the fitness of the model is increasing. Until the maximum genetic algebra is reached, a set of near-optimal parameters are obtained. Since the error persists, we can regard the suboptimal solution as the optimal solution. The subgraph in Figure 5 is an enlarged view of early inheritance. Although this set of parameters performs well and converges quickly, we can see from the figure that it converges around 200 generations, but, because of the instability of GA, sometimes it takes a long time to converge to the optimal value. Therefore, we unified the genetic algebra to 3000 generations during the experiment. It can be seen from the figure that the GA can quickly reduce the error in the early stage, and the convergence speed is very fast, which shows the effectiveness of our algorithm.
Evolutionary process of model optimal solution. Here, the yellow solid line represents the evolutionary process of the fitness values corresponding to the optimal parameters of each generation, while the blue solid line represents the early amplification part of the model solution.
In [5], the parameters βI=0.0498,βH=0.0225, and βF=0.0013 are given, and, in [23], the parameters βI=0.128,βH=0.08,βF=0.111, and δ1=δ2=0.75 are given. Meanwhile, the parameters βF=0.489,δ1=0.8,δ2=0.4 are given in [6]. Although some of the propagation coefficients are very different from those in other papers, because the deterministic models are different, we have that cabins are more, and we cannot directly compare them with those in other papers. We are uncertain about patients in the exposed period, so the propagation coefficient associated with it is uncertain. We can only use GA to solve each parameter value. Thus, we only perform comparison with some important parameters. In this paper, δ1=0.5141,δ2=0.6214 are close to the above literature to some extent. We obtain the basic reproduction number of Ebola virus in Sierra Leone R0≈1.895 which is calculated in the same way as in Xia et al.’s work [6]; it is basically consistent with R0≈1.7 given in [5] and R0≈1.78 given in [23]. The results show that the GA can accurately estimate all the parameters in the model, and the data are fitted well. Another advantage of the GA is that you are free to set the parameters precision, but you need to consider the length of the parameters set, which is very useful for getting high-precision parameters. It can be seen that GA can be used as a feasible method for parameter estimation.
Due to the fact that we are using a certain mathematical model, we only need to use the actual diseased data as test data and apply them in model (1) to get a set of near-optimal parameters. We set each of the parameter ranges to [0, 1], randomly generate the initial parameters as the input of the model, and use the GA to obtain a set of parameters that fit the model better. We can use GA to give a variety of parameter combinations for researchers’ reference. However, some parameters may have overfitting problems, so we conducted 70 experiments to select a set of parameters that are more realistic. We calculated the 95% confidence intervals for the best set of parameters in the 70 experiments, and almost all parameters are within the confidence intervals, which also shows the validity of the GA for parameter estimation. The confidence intervals for all the parameters are shown in Table 1.
4. Discussion and Conclusion
This work proposed a GA parameter estimation method based on adaptive mutation operator, which could be applied to biomathematical models and differential equations in other fields. Through GA’s adaptive search parameters, the parameters in various models can be effectively found, and multiple parameter combination schemes can be given, which reduces the process of manual adjustment of parameters by researchers and provides an effective reference for scientific research. In addition, in the study of infectious diseases, due to the complexity of the model, all of the parameters are often difficult to obtain, and the basic reproduction number R0 for evaluating whether there is an outbreak of an infectious disease needs to be obtained by calculating parameters; and the parameters can indicate the transmission dynamics of infectious diseases and the scale of transmission visually. Thus, parameter estimation is the most important, and the GA can be used to effectively find all the parameters in the model, which make up for the shortcomings of the traditional method, such as long calculation time and slow convergence speed. The method we have proposed has been evaluated in the experimental process, where the performance has reached the desired level. Finally, the GA can be applied not only to infectious disease models but also to other mathematical and physical models, and it proposes a new idea for parameter estimation.
Since the initial population of GA needs to assign values randomly, it will lead to output instability; we cannot quantify the uncertainty due to the discrete distribution of the output parameters, which is a common problem of GA and needs to be improved. On the other hand, fitness function is an important factor to determine the pros and cons of genetic evolution, so it is very important to select a suitable fitness function. Finally, a very important step in GA is to find the solution of the model, so it may not be suitable for the unsolvable equation, but it is applicable for most models. There are many improved genetic algorithms, and we can improve the existing algorithms and expand their scope of application.
The GA can solve specific problems with only a small amount of data, and the corresponding fitness function can be used for searching. Therefore, it is a general-purpose algorithm used in many fields. In function optimization, GA can estimate not only the parameters of biomathematical model but also the kinetic parameters of microorganisms. It can also solve the performance parameters of nonlinear physical problems. In addition, GA performs well in path planning, cloud computing task scheduling, communication network design, image feature extraction, and other fields. Furthermore, there are some studies that combine GA with machine learning methods such as neural network.
In this paper, we introduce the basic process of adaptive mutation genetic algorithm, introduce how to use GA to estimate the parameters of Ebola virus model in Sierra Leone, and give the data curve and genetic iterative process for fitting actual infections. In addition, we offer a new idea for parameter estimation in other research fields, such as dynamical model of disease transmission [48–50] or predator-prey interactions [51–53] with spatial effects [32, 54–56] in the form of reaction-diffusion equations.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by Program for the Outstanding Innovative Teams (OIT) of Higher Learning Institutions of Shanxi, Natural Science Foundation of Shanxi Province (Grant no. 201801D221003), and China Postdoctoral Science Foundation (Grant no. 2017M621110).
BrayM.MahantyS.Ebola hemorrhagic fever and septic shock2003188111613161710.1086/3797272-s2.0-0348147492BorchertM.MutyabaI.KerkhoveM. D. V.Ebola haemorrhagic fever outbreak in Masindi district, Uganda: Outbreak description and lessons learned20111135710.1186/1471-2334-11-3572-s2.0-84455162509ReynardO.Escudero-PerezB.VolchkovV.Dérégulation de l’hémostase dans les infections à filovirus201531214315010.1051/medsci/201531020102-s2.0-84924282962KhanA. S.The reemergence of ebola hemorrhagic fever, democratic Republic of the Congo19951791999S76S86LiL.Transmission dynamics of ebola virus disease with human mobility in Sierra Leone201710457557910.1016/j.chaos.2017.09.0222-s2.0-85029675412XiaZ.-Q.WangS.-F.LiS.-L.Modeling the transmission dynamics of Ebola virus disease in Liberia2015511385710.1038/srep138572-s2.0-84941073373Infectious Diseases Society of America2014Arlington, VA, USAInfectious Diseases Society of AmericaBwakaM. A.BonnetM. J.CalainP.Ebola hemorrhagic fever in Kikwit, democratic Republic of the Congo: Clinical observations in 103 patients1999179s1S1S710.1086/514308RollinP. E.KsiazekT. G.Ebola haemorrhagic fever2011377147148ShearsP.O’DempseyT. J. D.Ebola virus disease in Africa: Epidemiology and nosocomial transmission20159011910.1016/j.jhin.2015.01.0022-s2.0-84926980952LeroyE. M.KumulunguiB.PourrutX.Fruit bats as reservoirs of Ebola virus2005438706857557610.1038/438575a2-s2.0-28444494766Centers for Disease Control and Prevention, Diseases and ConditionsGrardG.BiekR.Muyembe TamfumJ.-J.Emergence of divergent Zaire Ebola virus strains in democratic Republic of the Congo in 2007 and 20082011204suppl 3S776S78410.1093/infdis/jir3642-s2.0-80054754287RouquetP.FromentJ. M.BermejoM.Wild animal mortality monitoring and human ebola outbreaks, Gabon and Republic of Congo, 2001–200320051128329010.3201/eid1102.0405332-s2.0-13744256174LeroyE. M.BaizeS.VolchkovV. E.Human asymptomatic Ebola infection and strong inflammatory response200035592222210221510.1016/s0140-6736(00)02405-32-s2.0-0034709761LuoX.ShaoN.ChengJ.ChenW.Modeling the trend of outbreak of COVID-19 in the diamond princess cruise ship based on a time-delay dynamic system202091522GeorgesA. J.LeroyE. M.RenautA. A.Ebola hemorrhagic fever outbreaks in Gabon, 1994–1997: Epidemiologic and health control issues1999179s1S65S7510.1086/514290MbonyeA.WamalaJ.Winyi-KaboyoW. K. V.TugumizemoVAcengJ.MakumbiI.Repeated outbreaks of viral hemorrhagic fevers in Uganda201212457910.4314/ahs.v12i4.312-s2.0-84873264330QiuX.WongG.AudetJ.Reversion of advanced Ebola virus disease in nonhuman primates with ZMapp20145147520475310.1038/nature137772-s2.0-84907263545MeltzerM. I.AtkinsC. Y.SantibanezS.Estimating the future number of cases in the Ebola epidemic-Liberia and Sierra Leone, 2014-2015201463suppl 3114SunZ.Analysis for the process of preventing and controlling plague20209914WilsonJ. A.Epitopes involved in antibody-mediated protection from Ebola virus200028754581664166610.1126/science.287.5458.16642-s2.0-0034051595EllardL.KatznelsonR.WasowiczM.Modeling the impact of interventions on an epidemic of ebola in Sierra Leone and Liberia20146e62HelleringerS.NoymerA.Assessing the direct effects of the ebola outbreak on life expectancy in Liberia, Sierra Leone and Guinea20157e6210.1371/currents.outbreaks.01a99f8342b42a58d806d7d1749574ea2-s2.0-84958554743KhanA.NaveedM.Dur-E-AhmadM.ImranetM.Estimating the basic reproductive ratio for the Ebola outbreak in Liberia and Sierra Leone2015411310.1186/s40249-015-0043-32-s2.0-85006208450GuoZ.-G.SongL.-P.SunG.-Q.LiC.JinZ.Pattern dynamics of an SIS epidemic model with nonlocal delay2019292195002710.1142/s02181274195002752-s2.0-85062645502LiL.Patch invasion in a spatial epidemic model201525834234910.1016/j.amc.2015.02.0062-s2.0-84923818700SunG.-Q.JusupM.JinZ.WangY.WangZ.Pattern transitions in spatial epidemics: Mechanisms and emergent properties201619437310.1016/j.plrev.2016.08.0022-s2.0-84994099913SunG.-Q.XieJ.-H.HuangS.-H.JinZ.LiM.-T.LiuL.Transmission dynamics of cholera: Mathematical modeling and control strategies20174523524410.1016/j.cnsns.2016.10.0072-s2.0-84992520676LiL.ZhangJ.LiuC.ZhangH.-T.WangY.WangZ.Analysis of transmission dynamics for Zika virus on networks201934756657710.1016/j.amc.2018.11.0422-s2.0-85057786623XingY.SongL.SunG.-Q.JinZ.ZhangJ.Assessing reappearance factors of H7N9 avian influenza in China201730919220410.1016/j.amc.2017.04.0072-s2.0-85017598066SunG.-Q.WangC.-H.WuZ.-Y.Pattern dynamics of a Gierer-Meinhardt model with spatial effects20178821385139610.1007/s11071-016-3317-92-s2.0-85009826116RannalaB.Identifiability of parameters in MCMC bayesian inference of phylogeny200251575476010.1080/106351502901024292-s2.0-0036803348MosselE.VigodaE.Phylogenetic MCMC algorithms are misleading on mixtures of trees200530957442207220910.1126/science.11154932-s2.0-25844516325AndrieuC.de FreitasN.DoucetA.JordanM. I.An introduction to MCMC for machine learning20035054310.1023/a:10202813271162-s2.0-0037262814ShahtoriN. M.ScoglioC.PourhabibA.Sequential Monte Carlo filtering estimation of Ebola progression in West AfricaProceedimgs of the American Control Conference (ACC)July 2016Boston, MA, USA12771282ConradP. R.MarzoukY. M.PillaiN. S.Accelerating asymptotically exact MCMC for computationally intensive models via local approximations201411115911607HadfieldJ. D.MCMC methods for multi-response generalised linear mixed models:The MCMCglmm R package201033212210.18637/jss.v033.i022-s2.0-77749249761GelfandA. E.SmithA. F. M.Sampling-based approaches to calculating marginal densities199085398409SavitzkyA.GolayM. J. E.Smoothing and differentiation of data by simplified least squares procedures19643681627163910.1021/ac60214a0472-s2.0-0042553279SuykensJ. A. K.VandewalleJ.Least squares support vector machine classifiers19999329330010.1023/a:10186286097422-s2.0-0032638628WangL.TangF.WuH.Hybrid genetic algorithm based on quantum computing for numerical optimization and parameter estimation200517121141115610.1016/j.amc.2005.01.1152-s2.0-31144437014ShetaA. F.De JongK.Time-series forecasting using GA-tuned radial basis functions20011333-422122810.1016/s0020-0255(01)00086-x2-s2.0-0035307307GoldbergD. E.HollandJ. H.Genetic algorithms and machine learning198832-3959910.1007/bf001138922-s2.0-34250092221WuR.PainumkalJ. T.VolkJ. M.Parameter estimation of nonlinear nitrate prediction model using genetic algorithmProceedings of the IEEE Congress on Evolutionary ComputationJune 2017Donostia-San Sebastian, SpainLucasiusC. B.KatemanG.Understanding and using genetic algorithms Part 1. Concepts, properties and context199319113310.1016/0169-7439(93)80079-w2-s2.0-0027284874Data for the 2014 Ebola outbeak in West Africa, http://github.com/cmrivers/ebola, 2014van BallegooijenW. M.BoerlijstM. C.Emergent trade-offs and selection for outbreak frequency in spatial epidemics200410152182461825010.1073/pnas.04056821012-s2.0-11144237747GuoZ.-G.SunG.-Q.WangZ.JinZ.LiL.LiC.Spatial dynamics of an epidemic model with nonlocal infection202037712515810.1016/j.amc.2020.125158OstfeldR.GlassG.KeesingF.Spatial epidemiology: An emerging (or re-emerging) discipline200520632833610.1016/j.tree.2005.03.0092-s2.0-20144376492ReichenbachT.MobiliaM.FreyE.Mobility promotes and jeopardizes biodiversity in rock-paper-scissors games200744871571046104910.1038/nature060952-s2.0-34548305141SunG.-Q.WangS.-L.RenQ.JinZ.WuY.-P.Effects of time delay and space on herbivore dynamics: Linking inducible defenses of plants to herbivore outbreak2015511124610.1038/srep112462-s2.0-84935010450ReichenbachT.MobiliaM.FreyE.Noise and correlations in a spatial population model with cyclic competition2007992323810510.1103/physrevlett.99.2381052-s2.0-36849066945SunG.-Q.WangC.-H.ChangL.-L.WuY.-P.LiL.JinZ.Effects of feedback regulation on vegetation patterns in semi-arid environments20186120021510.1016/j.apm.2018.04.0102-s2.0-85046749053LiL.JinZ.LiJ.Periodic solutions in a herbivore-plant system with time delay and spatial diffusion2016407-84765477710.1016/j.apm.2015.12.0032-s2.0-84953248005LiQ.LiuZ.YuanS.Cross-diffusion induced turing instability for a competition model with saturation effect2019347647710.1016/j.amc.2018.10.0712-s2.0-85056640125