
The formation of a biological organism, or an organ within it, can often be regarded as the unfolding of successive equilibria of a mechanical system. In a mathematical model, these changes of equilibria may be considered to be responses of mechanically constrained systems to a change of a reference configuration and of a reference metric, which are in turn driven by genes and their expression. This paper brings together three major threads of research. These are: Lie-type symmetries of equations; models as well as data on growth and pattern formation; and the relation between Lie algebras (and groups) and special functions associated with them. We show that symmetry methods can be generalized to map between solutions to models with different reference metrics. In the case in which we attempt to obtain such equations, they seem too complicated to be of any immediate service to the community of researchers on cortical growth. However, models and data on growth may be used to obtain generators of these Lie algebras empirically and numerically. These generators result in new classes of special functions. The paper is an invitation to develop what we may call empirical Lie algebras and associated functions. The hypothesis that remains to be tested is whether the confluence of ideas described in the paper, namely the Lie algebraic-related consequences of pattern formation and growth, prove useful for deepened understanding of biological growth patterns.
Citation: Raghu Raghavan. Growth and form, Lie algebras and special functions[J]. Mathematical Biosciences and Engineering, 2021, 18(4): 3598-3645. doi: 10.3934/mbe.2021181
[1] | Eunha Shim . Optimal dengue vaccination strategies of seropositive individuals. Mathematical Biosciences and Engineering, 2019, 16(3): 1171-1189. doi: 10.3934/mbe.2019056 |
[2] | Hyun Mo Yang, André Ricardo Ribas Freitas . Biological view of vaccination described by mathematical modellings: from rubella to dengue vaccines. Mathematical Biosciences and Engineering, 2019, 16(4): 3195-3214. doi: 10.3934/mbe.2019159 |
[3] | Vanessa Steindorf, Sergio Oliva, Jianhong Wu . Cross immunity protection and antibody-dependent enhancement in a distributed delay dynamic model. Mathematical Biosciences and Engineering, 2022, 19(3): 2950-2984. doi: 10.3934/mbe.2022136 |
[4] | Haitao Song, Dan Tian, Chunhua Shan . Modeling the effect of temperature on dengue virus transmission with periodic delay differential equations. Mathematical Biosciences and Engineering, 2020, 17(4): 4147-4164. doi: 10.3934/mbe.2020230 |
[5] | Bo Zheng, Lihong Chen, Qiwen Sun . Analyzing the control of dengue by releasing Wolbachia-infected male mosquitoes through a delay differential equation model. Mathematical Biosciences and Engineering, 2019, 16(5): 5531-5550. doi: 10.3934/mbe.2019275 |
[6] | Guangming Qiu, Zhizhong Yang, Bo Deng . Backward bifurcation of a plant virus dynamics model with nonlinear continuous and impulsive control. Mathematical Biosciences and Engineering, 2024, 21(3): 4056-4084. doi: 10.3934/mbe.2024179 |
[7] | Wei You, Jie Ren, Qimin Zhang . Finite-time contraction stability of a stochastic reaction-diffusion dengue model with impulse and Markov switching. Mathematical Biosciences and Engineering, 2023, 20(9): 16978-17002. doi: 10.3934/mbe.2023757 |
[8] | Biao Tang, Weike Zhou, Yanni Xiao, Jianhong Wu . Implication of sexual transmission of Zika on dengue and Zika outbreaks. Mathematical Biosciences and Engineering, 2019, 16(5): 5092-5113. doi: 10.3934/mbe.2019256 |
[9] | Saima Akter, Zhen Jin . Simulations and fractional modeling of dengue transmission in Bangladesh. Mathematical Biosciences and Engineering, 2023, 20(6): 9891-9922. doi: 10.3934/mbe.2023434 |
[10] | Zuohuan Zheng, Kaihua Wang, Daiyu Yang, Feifei Yin, Dingwei Sun, Weiyan Yu, Jialun Lin, Ying Liu, Changkuan Chen, Zehui Yang, Biao Wu . Factors affecting the transmission of dengue fever in Haikou city in 2019. Mathematical Biosciences and Engineering, 2023, 20(9): 16045-16059. doi: 10.3934/mbe.2023716 |
The formation of a biological organism, or an organ within it, can often be regarded as the unfolding of successive equilibria of a mechanical system. In a mathematical model, these changes of equilibria may be considered to be responses of mechanically constrained systems to a change of a reference configuration and of a reference metric, which are in turn driven by genes and their expression. This paper brings together three major threads of research. These are: Lie-type symmetries of equations; models as well as data on growth and pattern formation; and the relation between Lie algebras (and groups) and special functions associated with them. We show that symmetry methods can be generalized to map between solutions to models with different reference metrics. In the case in which we attempt to obtain such equations, they seem too complicated to be of any immediate service to the community of researchers on cortical growth. However, models and data on growth may be used to obtain generators of these Lie algebras empirically and numerically. These generators result in new classes of special functions. The paper is an invitation to develop what we may call empirical Lie algebras and associated functions. The hypothesis that remains to be tested is whether the confluence of ideas described in the paper, namely the Lie algebraic-related consequences of pattern formation and growth, prove useful for deepened understanding of biological growth patterns.
Dengue virus (DENV), a flavivirus primarily transmitted to human via a bites of an infected mosquitoes of the genus Aedes (i.e., Aedes aegypti or Aedes albopictus) [1,2,3,4,5,6,7,8,9,10]. This species is capable of transmitting other mosquito-borne diseases such as Chikungunya, Zika, Yellow fever [3,11,12,13,14]. The DENV is endemic in at least 100 countries in the tropical and subtropical regions of the world, with nearly 4 billion people at risk [7,15]. It is second to malaria amongst the deadly mosquito-borne diseases with fast-spreading ability [7], with an estimate of 390 million new infections of which about 96 million are clinical, up to 1 million are severe cases (i.e., dengue hemorrhagic fever/dengue shock syndrome also known as DHF/DSS), and more than 20 thousand deaths every year globally [6,7,16]. The DHF/DSS have a high mortality risk with death occurring within 24hrs after the onset of shock [17,18,19]).
The DENV are of four serotypes (i.e., DENV-1, DENV-2, DENV-3 and DENV-4) with low cross-immunity among each other which can result in secondary infection after an infection with one serotype has happened [7]. An individual infected with one serotype can be infected with different serotypes about a half year later [3], yet there is no proof of reinfection with the same serotype. DEN infection of all serotypes can cause three distinct syndromes: classic dengue fever, DHF and DSS [3].
The epidemiological cycle occurs as follows. Female Aedes aegypti mosquitoes are infected by the blood meal from biting an infectious human during the viremic stage of the infection. After a period of time called the extrinsic incubation period (EIP), the mosquitoes become infectious [20]. The range of the extrinsic incubation period lies between 4 to 10 days [20]. The infected mosquito does not have immunity and can transmit the virus throughout its lifespan [3,7,20]. An infectious mosquito can transmit the virus to a susceptible individual via bites [21], after a period of time called the intrinsic incubation period (IIP), which varies from 4 to 7 days [20], the virus can evolve from classical dengue fever to the secondary cases of dengue infection which may lead to DHF or DSS [17]. The manifestation of the symptoms varies between 3 to 7 days [3], which corresponds roughly to the infectious period. Thereafter, the individual develops the life-long immunity against the same serotype of the virus [3].
The prevention against DENV is limited to controlling mechanisms applied on the mosquito, because the vaccine is not yet available (even though there is a vaccine called dengvaxia, which is still under clinical trial, released by Sanofi Pasteur in 2015 which was approved by morethan 10 countries [7,22]. The vaccine efficacy varies by the serotypes of the DENV (i.e., 54.7 % for DENV-1, 43.0 % for DENV-2, 71.6 % for DENV-3, and 76.9 % for DENV-4 [7,23]). To describe the overall transmission dynamics of the DENV, we developed a mathematical model that takes into account the female adult mosquito coupled with the human population to analyses and explain the transmission dynamics of the DENV epidemics in Taiwan in 2014–2015.
A number of mathematical modeling studies have been carried out to gain insight into the DENV transmission dynamics in human and/or mosquito populations [3,4,5,6,11,24,25]. Garba et al. [6] developed a deterministic model for the DENV transmission within human and mosquito populations, and extended the model to incorporate an imperfect vaccine against the virus, and find that both models exhibited backward bifurcation phenomenon, and the courses have been explained in detail by Gumel [26]. Yang and Ferreiri [3] proposed a model to describe the transmission dynamics of the DENV in human and mosquito populations, they incorporated aquatic stages (i.e., eggs, larvae, pupae) in the mosquito population, and calculated the basic reproduction number to be less than one (which indicate that the disease can dies out in time). Yang et al. [24] designed a simple model in order to assess the effect of temperature on the population of the mosquito by using the temperature controlled experiments. Yang et al. [25] extended Yang et al. [24] by including human compartments and updated most of the parameters to assess the effect of temperature on the risk of the DENV outbreaks. Yang [5] developed and analyzed a model encompassing four quiescence stages of the mosquito, and used the model to assess the influence of the quiescent eggs on the transmission dynamics of the mosquito.
In the current study, we focus on the design and analysis of a model, which extends some of the aforementioned studies, for assessing the transmission dynamics of the DENV in human and mosquito populations. Our model includes both the human and mosquito populations, and incorporates asymptomatic (Ah) and severe cases (those with DHF/DSS) of the DENV infection (Ih2) in human population. We performed model fitting by using the plug-and-play inference framework for the human population, and our results contribute to better understanding of the transmission dynamics of the DENV epidemics in Taiwan in 2014-15 and provide useful guidelines for the design of control strategies in the future.
The paper is organized as follows. The model is formulated in section 2 and analyzed in section 3. Sensitivity analysis results of the full model is also given in section 4. Finally numerical analysis of the human only sub-model is performed with a plug-and-play inference framework in section 5.
The total human population at time t, denoted by Nh(t), is divided into six compartments of susceptible, Sh(t), exposed (infected but not infectious individuals that can progress to Ih1 class), Eh(t), asymptomatically infected, Ah(t) (i.e., expose individuals that can progress and becomes infectious without showing any clinical symptoms, which are lumped together for computational conveniences), symptomatically infected, Ih1(t), symptomatically infected with severe DENV (DHF/DSS), Ih2(t), and recovered, Rh(t), humans. Hence,
Nh(t)=Sh(t)+Eh(t)+Ah(t)+Ih1(t)+Ih2(t)+Rh(t). |
The total mosquito population at time t, denoted by Nm(t), is split into susceptible Sm(t), exposed (infected but not infectious mosquitoes) Em(t) and symptomatic (infectious mosquitoes) Im(t) classes. Hence,
Nm(t)=Sm(t)+Em(t)+Im(t). |
The model for the DENV transmission in the human and mosquito populations is given by the following deterministic ordinary differential equations (ODE) systems (2.1). The flow diagram of the model (2.1) is depicted in Figure 1. The state variables and parameters of the model are described in Tables 1 and 2, respectively.
dShdt=Πh−λhSh−μhSh,dEhdt=θλhSh−(vh+μh)Eh,dAhdt=(1−θ)λhSh−(γh+μh)Ah,dIh1dt=vhEh−(γ1+μh)Ih1,dIh2dt=(1−η)γ1Ih1−(δh+γ2+μh)Ih2,dRhdt=ηγ1Ih1+γ2Ih2+γhAh−μhRh,dSmdt=Πm−λmSm−μmSm,dEmdt=λmSm−(σm+μm)Em,dImdt=σmEm−μmIm. | (2.1) |
Variable | Interpretation |
Nh(t) | Total population of humans |
Sh(t) | Population of susceptible humans |
Eh(t) | Population of humans exposed to DENV |
Ah(t) | Population of asymptomatic humans |
Ih1(t) | Population of humans with clinical symptoms of the DENV |
Ih2(t) | Population of humans with severe clinical symptoms of the DENV |
Rh(t) | Population of humans recovered from the DENV |
Nm(t) | Total population of adult female mosquitoes |
Sm(t) | population of susceptible adult mosquitoes |
Em(t) | population of exposed adult mosquitoes |
Im(t) | population of infectious adult mosquitoes |
Parameter | Interpretation/Description |
Πh/Πm | Recruitment rate of humans/mosquitoes |
μh | Natural death rate of humans |
μm | Death rate of female adult mosquitoes |
βh | Transmission probability from infectious mosquitoes to susceptible humans |
βm | Transmission probability from infectious humans to susceptible mosquitoes |
a | Mosquito biting rate |
b | Maximum number of bites a human can receive per unit time |
θ | Fraction of infected humans that are exposed |
vh | Progression rate of exposed humans to infectious humans with clinical symptoms |
η | Fraction of infectious humans that will not get a severe symptom |
γh, γ1, γ2 | Recovery rate of infectious humans from Ah, Ih1, Ih2, respectively |
δh | Disease-induced death rate of humans |
σm | Progression rate of exposed mosquitoes to the infectious mosquitoes |
ρ | Reporting rate to actual case ratio |
m | Average mosquito to human ratio |
Here, the infection rate for humans (λh) and vectors (λm) are given by
λh=abβhImNh,andλm=abβm(Ah+Ih1+Ih2)Nh. | (2.2) |
In model (2.1), Πh is the recruitment rate for human by birth, λh is the infection rate of susceptible human from infectious mosquitoes, where a is the mosquito biting rate, b is the maximum number of bites a human can receive per unit time, βh is the transmission probability from infected mosquitoes to the susceptible humans while βm is transmission probability from an infectious humans to the susceptible mosquitoes per bite, and μh is the natural death rate of human. The term θ is the symptomatic ratio, and (1−θ) remains the fraction of asymptomatic ratio. The parameter vh measures the rate at which an individuals in the Eh class develop a primary clinical symptoms of the virus, while η is a fraction of humans who will not get a severe DENV, δh accounts for the disease-induced death rate of humans from the Ih2 class, and γi(i=0,1,2) measures the rate of recovery of individual from Ah, Ih1, and Ih2, respectively.
The susceptible adult mosquitoes acquire the DENV at a rate λm, this population is decreased by natural death at a rate μm. Finally, σm is a progression rate of the DENV from the Em(t) class to the Im(t). It is also worth stating that the model (2.1) accounts for the conservation law of mosquito bites. Thus, the infection rates, λh and λm, are normalized by the total host population, that is Nh(t) [27].
The basic properties of the model (2.1) will now be explored. Consider the following equations for the rate of change of the total human N′h(t) and mosquito N′m(t) populations
N′h(t)=dNhdt=Πh−μhNh−δhIh2⩽Πh−μhNh, | (2.3) |
and
N′m(t)=dNmdt=Πm−μmNm, | (2.4) |
here, the prime, ′, represent a differentiation with respect to time, t.
Furthermore, consider the region,
Ω={(Sh,Eh,Ah,Ih1,Ih2,Rh,Sm,Em,Im)∈R9+:Nh⩽Πhμh,Nm⩽Πmμm}. |
It can be shown by solving for Nh and Nm in Eqs (2.3)-(2.4) so that all solutions of the system starting in the region Ω will remain in Ω for all time t with t⩾0. Thus, the region Ω is positively-invariant, and it is sufficient to consider solutions restricted in Ω. In this region, the usual existence, uniqueness and continuation results hold for the system (2.1) [28,29].
The disease-free equilibrium (DFE) of the model (2.1) obtained at steady state is given by
E0=(S0h,E0h,A0h,I0h1,I0h2,R0h,S0m,E0m,I0m)=(Πhμh,0,0,0,0,0,Πmμm,0,0). |
Using the next generation matrix method [30], we obtained the associated reproduction number of the model (2.1), denoted by R0=ρ(FV−1). The term ρ represents the function to find the spectral radius of the next generation matrix, G=FV−1. The matrices F (for the new infection terms) and V (for the remaining transition terms), associated with the model (2.1), are given by
F=[00000θabβh00000(1−θ)abβh0000000000000abβmS0mNhabβmS0mNhabβmS0mNh00000000], | (3.1) |
V=[n1000000n20000−vh0n300000−(1−η)γ1n4000000n500000−σmμm], | (3.2) |
where n1=vh+μh, n2=γh+μh, n3=γ1+μh, n4=δh+γ2+μh, and n5=σm+μm. Therefore, the basic reproduction number, R0, is given by
R0=Rh⋅Rm, | (3.3) |
where
Rh=√abβhμhΠh⋅[θvhn1n3+θ(1−η)γ1vhn1n3n4+(1−θ)n2], | (3.4) |
and
Rm=√abβmσmS0mn5μm. | (3.5) |
R0 interpretation: The threshold quantity R0 is ecologically and epidemiologically interpreted as follows.
ⅰ. abβhNh accounts for the number of new infected human hosts caused by an infected mosquito over its expected infectious period,
ⅱ. vhn1 is the probability that an exposed human survives the exposed stage and move to the infectious stage, Ih1,
ⅲ. 1n3 is the average duration in the infectious stage,
ⅳ. (1−η)γ1n4 is the probability that an infectious human survives the infectious stage (including severe stage of the DENV) and move to the recovered class,
ⅴ. (1−θ)n2 is the probability that an infected human move to the asymptomatic class,
ⅵ. abβmσmn5S0mμm accounts for the possible number of new infected mosquito caused by an infected human hosts over its exposed infectious period.
The theorem 3.1 below follows the Theorem 2 of Ref [30].
Theorem 3.1. The DFE, E0, of the model (2.1), is locally-asymptotically stable (LAS) in Ω if R0<1, but unstable if R0>1.
The endemic equilibrium (EE) for the system (2.1), which is
E∗=(S∗h,E∗h,A∗h,I∗h1,I∗h2,R∗h,S∗m,E∗m,I∗m), |
in terms of the forces of infection, λ∗h and λ∗m, is given by
S∗h=Πhλh+μh,E∗h=θλhΠhn1(λh+μh),A∗h=(1−θ)λhΠhn2(λh+μh),I∗h1=vhθλhΠhn1n3(λh+μh),I∗h2=(1−η)γ1vhθλhΠhn1n3n4(λh+μh),R∗h=Πhλhμh(λh+μh)((1−θ)γhn2+ηγ1vhθn1n3+γ1γ2(1−η)vhθn1n3n4),S∗m=Πmλm+μm,E∗m=Πmλmn5(λm+μm)andI∗m=σmΠmλmn5μm(λm+μm). | (3.6) |
Substituting the above equilibrium points (equation (3.6)) into equation (2.2), we have
Aλ2h+Bλh+C=0, | (3.7) |
where
A=Πh(σm+μm)μm(−abβmμhθvhn2n4−abβmμhγ1θvhn2+abβmμhγ1θvhn2η−abβmμh(1−θ)n1n3n4−μmθμhn2n3n4−μm(1−θ)μhn1n3n4−μmθvhμhn2n4−μmγ1θvhμhn2+μmγ1θvhμhn2η−μmγh(1−θ)n1n3n4−μmηγ1vhθn2n4−μmγ1γ2vhθn2+μmγ1γ2vhθn2η)(−θμhn2n3n4−(1−θ)μhn1n3n4−θvhμhn2n4−γ1θvhμhn2+γ1θvhμhn2η−γh(1−θ)n1n3n4−ηγ1vhθn2n4−γ1γ2vhθn2+γ1γ2vhθn2η), |
B=n4n2n1μh(−μm(σm+μm)[(((((−2n3−2vh)n4+2γ1vh(−1+η))n2+2n1n3n4)μh−2vh(n4η−γ2(−1+η))γ1n2+2γhn1n3n4)θ−2n1n3n4(μh+γh))μm+a(((−n4+γ1(−1+η))vhn2+n1n3n4)θ−n1n3n4)μhbβm]Πh+a2σm(((−n4+γ1(−1+η))vhn2+n1n3n4)θ−n1n3n4)μhΠmb2βmβh)n3,andC=(1−R20). |
Therefore, one can obtain the positive EE of the model (2.1) by simplifying the quadratic equation (3.7), and substituting only the positive values of the λh into the EE points. The occurrence of the quadratic equation (3.7) shows that the backward bifurcation (BB) for the model (2.1) exists.
Hence, the following Theorem 3.2 is established.
Theorem 3.2. The model (2.1) has
i) a EE if C<0⇔R0>1;
ii) a unique EE if B<0 and C=0 or B2−4AB=0;
iii) two EEs if B<0, C>0 and B2−4AC>0; or
iv) no EE otherwise.
The existence of the phenomenon of the BB in the model (2.1) highlights a co-existence between a stable disease free-equilibrium (DFE) and a stable endemic equilibrium (EE) even if the basic reproduction number (R0) is less than unity, this makes the disease control more difficult. Furthermore, the disease eradication would no longer depends on the basic concept of R0, i.e., when R0<1, the disease dies out in time; and while R0>1, the disease persists in the community. Although the BB has been firstly shown to exist in the DENV from the work of Garba et al. [6] and its courses has been explained extensively by Gumel [26]. In this study, our analysis suggests that the BB can be removed if the R0 can be equal to one (i.e., R0=1), so that equation (3.7) will be Aλh+B=0 which is linear equation.
The following Theorem 3.3 is claimed.
Theorem 3.3. The endemic equilibrium (EE), E∗, is globally-asymptotically stable (GAS) in Ω when R0>1 provided that
(1−λhλ∗h)(1−Ih1λ∗hI∗h1λh)⩾0,(1−λhλ∗h)(1−Ih2λ∗hI∗h2λh)⩾0,(1−λmλ∗m)(1−Imλ∗mI∗mλm)⩾0,and(Ih1I∗h1−lnIh1Ih1+lnAhAh+AhA∗h)⩽0. | (3.8) |
The proof of the above Theorem 3.3 is given in Appendix 7, and the Figure 2 is used to demonstrate the analytical result obtained.
Parameter | Baseline [Range] | Units | Source |
Πh | 2.5 [1,5] | person~day−1 | [6] |
Πm | 5000 [1000,6000] | Day−1 | estimated from [6] |
μh | 3.9×10−5 [3.6×10−5, 4.0×10−5] | Day−1 | [12] |
μm | 0.05714 [0.01,0.1] | Day−1 | [11] |
βh | 0.75 [0.1,0.95] | Day−1 | [6] |
βm | 0.75 [0.1,0.95] | Day−1 | [6] |
a | 0.4997 [0.1,1] | Day−1 | [13,14] |
b | 0.4 [0.1,1] | Day−1 | [13] |
θ | 0.18 [0.1,0.6] | Dimensionless | [31] |
vh | 0.0666 [0,1] | Day−1 | [6,32] |
η | 0.1 [0.05,0.15] | Dimensionless | [19,33] |
γh | 0.1428 [0.1,0.2] | Day−1 | [13] |
γ1 | 0.2 [0.1428,0.3] | Day−1 | [13] |
γ2 | 0.05 [0.0333,0.07143] | Day−1 | [13] |
δh | 10−3 [10−4,2×10−3] | Day−1 | [6] |
σm | 0.5 [0,1] | Day−1 | [6] |
m | 10 [1,20] | Dimensionless | assumed |
Following the previous studies [13,14,34,35], we adopted the partial ranked correlation coefficient (PRCC) for sensitivity analysis. The PRCC of the R0 and IAR of the model (2.1) is estimated. The sensitivity analysis results in Figure 3 suggest the top three ranked parameters (i.e., a, b and Πm) are prioritized in controlling the DENV epidemics in Taiwan.
The analyses in this section will be carried out for the special scenario of the model (2.1) in the absence of mosquito compartments (i.e., a human only sub-model). We assumed that there is no birth and death rate (Πh=μh=0), and no disease-induced mortality rate (δh=0) in a short period. We also assumed that the effective contact rate of the simplified version of the model to be β=β(t)=ab⋅m(t), where m(t) is the time-dependent mosquito to human ratio. This setting allows us to fit the "SIR" based model to a vector-borne disease under the "absence of vector" model structure. The term of the product of (abm) is equivalent to the standard Ross-Macdonald malaria model [36].
The equivalent (simplified) model without mosquito compartment together with the above assumptions is given by
dShdt=−βIh1ShNh,dEhdt=θβIh1ShNh−vhEh,dIh1dt=vhEh−γ1Ih1,dIh2dt=(1−η)γ1Ih1−γ2Ih2,dAhdt=(1−θ)βIh1ShNh−γhAh,dRhdt=ηγ1Ih1+γ2Ih2+γhAh, | (5.1) |
where the basic reproduction ratio of the simplified model (5.1) is given by R0(t)=βθγ1.
The weekly number of cases of the i-th week is
Zi=∫weekiργ1Ih1dt, | (5.2) |
where the term ρ denotes a constant reporting rate of the DENV cases. The Zi denotes the theoretical weekly DENV cases yield from the model (5.1).
We model the observed (reported) DENV cases, Ci for the i-th week during the study period, as a partially observed Markov process (POMP, also know as hidden Markov model, HMM). Instead of implementing Poisson-distributed priors [34], all Cis are assumed to follow over-dispersed Poisson distributions according to the theoretical modelling (Eq (5.2)) outputs, Zis. Since the rate of Poisson distribution is a Gamma random variable, the observed weekly number of people who are confirmed DENV infections (Ci of the i-th week) is a random sample from a Negative-binomial (NB) distribution. Therefore,
Ci∼NB(mean=Zi,variance=Zi(1+ψZi)), | (5.3) |
where ψ is the over-dispersion parameter for NB distribution under estimation. In this model, we set Li(⋅) be the likelihood function for the week, which is the "probability" of Ci, given the real cases from simulations Zi under the NB distribution [37,38]. In this section, the model simulations are conducted by using the R (version 3.4.1) package "POMP" [39].
The overall log-likelihood, l, for the whole-time series is
l(Θ)=T∑i=1ln[Li(Ci∣C0,...,Ci−1;Θ)], | (5.4) |
where Θ is the parameter vector under estimation. The term T denotes the total number of weeks during the study period. We apply the iterated filtering algorithm with the plug-and-play likelihood-based inference framework to estimate the maximum likelihood estimates (MLE) of Θ, (see [37,40,41,42,43,44,45,46]. We use the fixed-time-step Euler-multinomial algorithm to simulate the ODE system (5.1) [14,34,37,38,40], We compared different models using the small-sample-size corrected Akaike's Information Criterion (AICc) [44] as a measurement of the trade-off between model complexity and the goodness-of-fit. The AICc is given by
AICc=−2l(ˆΘ)+2k+2k(k+1)N−k−1, | (5.5) |
where N is the number of data points and k is the number of free parameters.
The (simplified) epidemic model (5.1) is fitted to the dengue cases time series in Kaohsiung and Tainan cities in 2014–2016. We proposed the time-dependent basic reproduction number (i.e., R0(t)) driven by variable mosquito abundance and dengue transmission rate (e.g., mosquito biting rate, etc). Here, we compared two numerical approaches to fit the epidemic data and reconstruct the R0(t) series, namely,
● the cubic spline reconstruction approach, and
● the square wave reconstruction approach.
The cubic spline approach is defined from the previous studies [14,47]. We used pq to denote the number of nodes in the cubic spline of R0(t) to be reconstructed. We tested different possible pq's with the aim to find the value of pq that leads to the lowest AICc. We further restrict that Kaohsiung and Tainan shared exactly the same cubic spline function of R0(t). Hence, the two cities were fitted simultaneously with the same set of parameter estimates. For the cubic spline approach, the fitting is conducted with completely the same time series of R0(t) (i.e., same timing, scale, and all settings are the same) for both cities. And we find the goodness-of-fit (in term of the MLL) can still be improved for both cities. Even though the spline can yield a relatively flexible (in term of the degree of freedom of the parameters to be estimated) R0(t) reconstruction outcome, the fittings are still not satisfied. Thus, we speculate a non-synchronized R0(t) would be more reasonable. The fitting results of the R0(t) reconstructed by the cubic spline function are shown in Figure 4, and the parameter estimates are summarized in Table 5. The square wave is defined as a baseline transmission rate plus a constant additional value for each year. The additional values are different for both Kaohsiung and Tainan in the same year, and also vary in different years. To explore the possible non-synchronization of R0(t)'s in Kaohsiung and Tainan cities, we further allow the changing time of the square wave to have different estimates in the two cities. Thus, the reconstructed R0(t)s in Kaohsiung and Tainan are expected to have different shape which indicates that the R0(t)'s of the two cities are not perfectly synchronized. In Tainan, we fix the time of the occurrence of the R0(t) high values (i.e., with the additional value) to be May 1 and October 31 in both 2014 and 2015. In Kaohsiung, our estimates are that the high values of R0(t) occurred between March 16 and October 31 in 2014, and between May 1 and November 25 in 2015. These time shifts are crucial to achieve good fitting outcomes in Kaohsiung city. Thus, this finding indicates that the DENV outbreaks in Kaohsiung turned on early in 2014, and turned off late in 2015. For the square-wave approach, the fitting is conducted with similar series of R0(t)'s for the two cities. The timing (of change) and the scale of R0(t)'s were slightly different. Although this approach is less flexible as the spline approach (i.e., with 4 degrees of freedom decreased), it leads to better fitting results with MLL improved by 10 units. And together, the AIC improves (i.e., decreases) by 24 units. This indicates that the non-synchronized square-wave approach is a strong and significant improvement from the synchronized spline approach. The fitting results of R0(t) reconstructed by the square wave function are shown in Figure 5, and the parameter estimates are summarized in Table 6. Therefore, for the spline versus square-wave approach, the non-synchronized square-wave approach has better goodness-of-fit (i.e., higher MLL) with less flexibility (i.e., less degree of freedom). In other words, the non-synchronized square-wave approach has better fitting performance (in term of the AIC). Hence, we conclude that non-synchronized R0(t)'s were likely to occur in the two cities.
Our (simplified) DENV epidemic model well-explained the temporal patterns of the DENV outbreaks in the two cities in Taiwan under biologically reasonable conditions (i.e., set of model parameters). We discovered that the epidemics in the two cites can be reconstructed by similar transmission series (i.e., R0(t)), and also estimated with the same model structure. The square wave approach not only reduced the number of free parameters by 4, but also increased the log likelihood by 10 units. Hence, its AICc decreases roughly by 24 units from the cubic spline reconstruction, which is a significant improvement.
Although we find that the estimated transmission rates exhibited similar seasonal patterns in Kaohsiung and Tainan, they are neither completely synchronized between the two cities despite their proximity, nor periodically in-phase perfectly. If they were completely synchronized or perfectly periodic, the cubic spline approach would identify such a pattern reflected by an equivalent level of goodness-of-fit as in the square wave approach. We conclude that the time lags between the seasonal waves of the transmission rates (R0(t)) in the two cities are likely to occurred. Other than the pathogenic factors that can causes the phase difference in the transmissibility [48], these time lags could be due to the heterogeneity in the local conditions in Kaohsiung and Tainan cities.
We formulated and analyzed a mathematical model to study the transmission dynamics of the dengue outbreaks in Kaohsiung and Tainan cities in Taiwan in 2014–2015. We showed that the model (2.1) has two equilibria, of which the disease-free equilibrium (DFE) is locally asymptotically stable (LAS) whenever the basic reproduction number (R0<1), and unstable if the R0>1. And the endemic equilibrium (EE) is globally asymptotically stable (GAS) in the region of attraction, Ω, whenever the R0>1 (numerical examples, see Figure 2, is used to demonstrate the analytical results obtained). Our analysis shows the existence of backward bifurcation (BB) phenomenon of the model (2.1), a situation where a stable disease free-equilibrium co-exists with a stable endemic equilibrium even when the basic reproduction number (R0) is less than one, which makes the disease control even more difficult and is no longer dependent on the basic concept of the R0 (i.e., R0<1 the disease dies out in time, while R0>1 the disease persist to the community). In this paper, the analysis suggest that the backward bifurcation can be removed when the basic reproduction number is equal to one (i.e., R0=1).
We employed the plug-and-play statistical inference framework on the special case of the model (2.1) in the absence of mosquito. We fitted the simplified model (5.1) to the weekly number of cases in the Kaohsiung and Tainan Cities in Taiwan (obtained from the National Infectious Disease Statistics System [49], see Figures. 4 and 5). Our (simplified) dengue epidemic model is able to explained the temporal patterns of the dengue outbreaks in the two cities in Taiwan under biologically reasonable conditions (i.e., set of model parameters). We obtained that the epidemics in the two cites can be reconstructed by the same transmission series (i.e., R0(t)) and also estimated with the same model structure. This findings suggests the similarity of the dengue outbreak in a two different places during the same period of time, which could be due to the similar demographic and/or meteorological conditions. The proposed two reconstruction approaches presents (roughly) equivalent goodness-of-fit in terms of the values of the log-likelihood. However, the AICc improves (decreases) roughly by 24 units from the cubic spline to the square wave reconstruction. We also found that despite the proximity in Kaohsiung and Tainan, the estimated transmission rates were neither completely synchronized, nor periodically in-phase perfectly in the two cities. The time lags between the seasonal waves of the transmission rates in the two cities were likely to occur. These time lags could be due to the differences in their local conditions. Furthermore, we show that the non-synchronized square-wave approach has better goodness-of-fit (i.e., higher MLL) with less flexibility (i.e., less degree of freedom). In other words, the non-synchronized square-wave approach has a better fitting performance (in term of the AIC) than the cubic spline approach. Hence, we concluded that non-synchronized R0(t)s were likely to occur in the two cities.
The sensitivity analysis results shows the top three ranked parameters (i.e., a, b and Πm), that is, the mosquito biting rate, the maximum number of bites a human can receive per unit time, and the mosquito recruitment rate are the key parameters to be prioritized for controlling the Taiwan dengue outbreaks. This further suggests that proper sanitation of mosquito breeding sites and avoiding the mosquito bites are the key control measures to future dengue outbreaks in Taiwan.
List of Abbreviations
We listed the abbreviations used in Table 4.
abbreviation | explanation / full name |
AICc | corrected Akaike's Information Criterion |
BB | Backward bifurcation |
CI | Confidence interval |
DENV | Degue virus |
DENV-1 | Dengue virus serotype 1 |
DENV-2 | Dengue virus serotype 2 |
DENV-3 | Dengue virus serotype 3 |
DENV-4 | Dengue virus serotype 4 |
DFE | Disease-free equilibrium |
DHF | Dengue hemorrhagic fever |
DSS | Dengue shock syndrome |
EE | Endemic equilibrium |
EIP | Extrinsic incubation period |
GAS | Globally-asymptotically stable |
HMM | Hidden Markov Model |
IIP | Intrinsic incubation period |
LAS | Locally-asymptotically stable |
MLE | Maximum likelihood estimates |
MLL | Maximum log-likelihood |
NB | Negative binomial (distribution) |
ODE | Ordinary differential equation |
POMP | partially observed Markov process |
PRCC | Partial rank correlation coefficient |
Zhen Jin is partly supported by Shanxi Scientific and Technology Innovation Team (201705D15111172 & 201805D131012-1). We thank reviewers for helpful comments.
The authors declare that they have no conflict of interest in this paper.
Proof. Define a Lyapunov function as follows
V(t)=J1(Sh−S∗h−S∗hlnShS∗h)+J2(Eh−E∗h−E∗hlnEhE∗h)+J3(Ah−A∗h−A∗hlnAhA∗h)+J4(Ih1−I∗h1−I∗h1lnIh1I∗h1)+J5(Ih2−I∗h2−I∗h2lnIh2I∗h2)+K1(Sm−S∗m−S∗mlnSmS∗m)+K2(Em−E∗m−E∗mlnEmE∗m)+K3(Im−I∗m−I∗mlnImI∗m), | (A1-1) |
where J1=2, J2=11−θ, J3=1(θ), J4=λ∗hS∗hvhE∗h, J5=λ∗hS∗h(1−η)γ1I∗h1, K1=K2=1, and K3=λ∗mS∗mσmE∗m.
Thus, the Lyapunov derivative computed along solutions of the system (2.1) is given by
˙V(t)=J1(1−S∗hSh)˙Sh+J2(1−A∗hAh)˙Ah+J3(1−E∗hEh)˙Eh+J4(1−I∗h1Ih1)˙Ih1+J5(1−I∗h2Ih2)˙Ih2+K1(1−S∗mSm)˙Sm+K2(1−E∗mEm)˙Em+K3(1−I∗mIm)˙Im. | (A1-2) |
Hence,
˙V(t)=−J1μh(Sh−S∗h)2Sh+J1λ∗hS∗h[1−λhShλ∗hS∗h−S∗hSh+λhλ∗h]+J2(1−θ)λ∗hS∗h[λhShλ∗hS∗h−AhA∗h−A∗hλhShAhλ∗hS∗h+1]+J3θλ∗hS∗h[λhShλ∗hS∗h−EhE∗h−E∗hλhShEhλ∗hS∗h+1]+J4vhE∗m[EhE∗h−Ih1I∗h1−I∗h1EhIh1E∗h+1]+J5(1−η)γ1Ihh1∗[Ih1I∗h1−Ih2I∗h2−I∗h2Ih1Ih2I∗h1+1]+−K1μm(Sm−S∗m)2Sm+K1λ∗mS∗m[1−λmSmλ∗mS∗m−S∗mSm+λmλ∗m]+K2λ∗mS∗m[λmSmλ∗mS∗m−EmE∗m−E∗mλmSmEmλ∗mS∗m+1]+K3σmE∗m[EmE∗m−ImI∗m−I∗mEmImE∗m+1]. | (A1-3) |
Following [50,51], the function v(x)=1−x+lnx, if x>0, it leads to v(x)⩽0. And if x=1, we have v(x)=0. Thus, x−1⩾ln(x) for any x>0. Then, we also have that
λ∗hS∗h[2−S∗hSh−EhE∗h−λhShE∗hλ∗hS∗hEh+λhλ∗h]=λ∗hS∗h[−(1−λhλ∗h)(1−Ih1λ∗hI∗λh)+3−S∗hSh−λhShE∗hλ∗hS∗hEh−Ih1λ∗I∗h1λh−EhE∗h+Ih1I∗h1]⩽λ∗hS∗h[−(S∗hSh−1)−(λhShE∗hλ∗hS∗hEh−1)−(Ih1λ∗hI∗h1λh1−1)−EhE∗h+Ih1I∗h1]⩽λ∗hS∗h[−ln(S∗hShλhShE∗hλ∗S∗hEhIh1λ∗hI∗h1λh)−EhE∗h+Ih1I∗h1]=λ∗hS∗h[Ih1I∗h1−ln(Ih1I∗h1)+ln(EhE∗h)−EhE∗h]. | (A1-4) |
Similarly, we have
λ∗hS∗h[2−S∗hSh−AhA∗h−λhShA∗hλ∗hS∗hAh+λhλ∗h]⩽λ∗hS∗h[Ih2I∗h2−ln(Ih2I∗h2)+ln(AhA∗h)−AhA∗h], | (A1-5) |
and
λ∗mS∗m[2−S∗mSm−EmE∗m−λmSmE∗mλ∗mS∗mEm+λmλ∗m]⩽λ∗mS∗m[ImI∗m−ln(ImI∗m)+ln(EmE∗m)−EmE∗m]. | (A1-6) |
Furthermore, we have that
λ∗hS∗h(EhE∗h−Ih1I∗h1−I∗h1EhIh1E∗h+1)⩽λ∗hS∗h(EhE∗h−ln(EhE∗h)+ln(Ih1I∗h1)−Ih1I∗h1), | (A1-7) |
and also
λ∗hS∗h(Ih1I∗h1−Ih2I∗h2−I∗h2Ih1Ih2I∗h1+1)⩽λ∗hS∗h(Ih1I∗h1−ln(Ih1I∗h1)+ln(Ih2I∗h2)−Ih2I∗h2), | (A1-8) |
as well as
λ∗mS∗m(EmE∗m−ImI∗m−I∗mEmImE∗m+1)⩽λ∗mS∗m(EmE∗m−ln(EmE∗m)+ln(ImI∗m)−ImI∗m). | (A1-9) |
Hence, the equations (A1-1)-(A1-9) together with condition (3.8) ensure that dVdt⩽0. Furthermore, the strict inequality dVdt=0 holds only for Sh=S∗h, Eh=E∗h, Ah=A∗h, Ih1=I∗h1, Ih2=I∗h2, Rh=R∗h, Sm=S∗m Em=E∗m and Im=I∗m. Thus, the endemic equilibrium state E∗ is the only positive invariant set to the system (2.1) contained entirely in {(Sh,Eh,Ah,Ih1,Ih2,Rh,Sm,Em,Im)∈Ω:Sh=S∗h,Eh=E∗h,Ah=A∗h,Ih1=I∗h1,Ih2=I∗h2,Rh=R∗h,Sm=S∗m,Em=E∗m,Im=I∗m\bigamma}. Therefore, it follows from the LaSalle's invariance principle [52] that every solutions to the equations in (A1-2) with initial conditions in Ω converge to stable endemic equilibrium point, E∗, as t→∞. Hence, the positive endemic equilibrium is globally asymptotically stable.
Parameter | Value (Kaohsiung) | Value (Tainan) | Source |
Sh(0) | 0.93212 | 0.8774 | estimated |
Eh(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ah(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ih1(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ih1(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Rh(0) | 0.0688 | 0.1226 | estimated |
β | 1×10−20 | 1×10−20 | estimated |
ρ | 0.0631 | 0.0631 | [53,54,55] |
nm | 9 | 9 | estimated |
Parameter | Value (Kaohsiung) | Value (Tainan) | Source |
Sh(0) | 0.9483 | 0.9161 | estimated |
Eh(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ah(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ih1(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ih1(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Rh(0) | 0.0517 | 0.0839 | estimated |
β | 1.9986 | 1.9986 | estimated |
ρ | 0.0837 | 0.0837 | estimated |
[1] | W. T. D'Arcy, On Growth and Form, Dover Publications, 1992. |
[2] | J. Altman, Neural and mental evolution: origins of the human body, brain, behavior, consciousness, and culture, Laboratory of Developmental Neurobiology, Inc., 2013. Available from: https://brainmindevolution.org/. |
[3] |
K. Amunts, C. Lepage, L. Borgeat, H. Mohlberg, T. Dickscheid, M. É. Rousseau, et al., BigBrain: an ultrahigh-resolution 3D human brain model, Science, 340 (2013), 1472–1475. doi: 10.1126/science.1235381
![]() |
[4] |
P. V. Bayly, R. J. Okamoto, G. Xu, Y. Shi, L. A. Taber, A cortical folding model incorporating stress-dependent growth explains gyral wavelengths and stress patterns in the developing brain, Phys. Biol., 10 (2013), 016005. doi: 10.1088/1478-3975/10/1/016005
![]() |
[5] |
P. V. Bayly, L. A. Taber, C. D. Kroenke, Mechanical forces in cerebral cortical folding: a review of measurements and models, J. Mech. Behav. Biomed. Mater., 29 (2014), 568–581. doi: 10.1016/j.jmbbm.2013.02.018
![]() |
[6] | M. A. Biot, Mechanics of Incremental Deformations, John Wiley & Sons, 1965. |
[7] | S. Budday, C. Raybaud, E. Kuhl, A mechanical model predicts morphological abnormalities in the developing human brain, Sci. Rep., 4 (2014), 5644. |
[8] |
R. Chelakkot, L. Mahadevan, On the growth and form of shoots, J. R. Soc., Interface, 14 (2017), 20170001. doi: 10.1098/rsif.2017.0001
![]() |
[9] |
B. Fischl, M. I. Sereno, R. B. H. Tootell, A. M. Dale, High-resolution intersubject averaging and a coordinate system for the cortical surface, Hum. Brain Mapp., 8 (1999), 272–284. doi: 10.1002/(SICI)1097-0193(1999)8:4<272::AID-HBM10>3.0.CO;2-4
![]() |
[10] |
L. M. J. Florack, B. M. T. H. Romeny, J. J. Koenderink, M. A. Viergever, Cartesian differential invariants in scale-space, J. Math. Imaging Vision, 3 (1993), 327–348. doi: 10.1007/BF01664793
![]() |
[11] | L. M. J. Florack, B. M. ter Haar Romeny, J. J. Koenderink, M. A. Viergever, Linear scale-space, J. Math. Imaging Vision, 4 (1994), 325–351. |
[12] |
G. Franceschini, D. Bigoni, P. Regitnig, G. A. Holzapfel, Brain tissue deforms similarly to filled elastomers and follows consolidation theory, J. Mech. Phys. Solids, 54 (2006), 2592–2620. doi: 10.1016/j.jmps.2006.05.004
![]() |
[13] |
K. E. Garcia, E. C. Robinson, D. Alexopoulos, D. L. Dierker, M. F. Glasser, T. S. Coalson, et al., Dynamic patterns of cortical expansion during folding of the preterm human brain, Proc. Natl. Acad. Sci., 115 (2018), 3156–3161. doi: 10.1073/pnas.1715451115
![]() |
[14] |
G. Geng, L. A. Johnston, E. Yan, J. M. Britto, D. W. Smith, D. W. Walker, et al., Biomechanisms for modelling cerebral cortical folding, Med. Image Anal., 13 (2009), 920–930. doi: 10.1016/j.media.2008.12.005
![]() |
[15] |
J. N. Giedd, J. Blumenthal, N. O. Jeffries, F. X. Castellanos, J. L. Rapoport, Brain development during childhood and adolescence: a longitudinal MRI study, Nat. Neurosci., 2 (1999), 861–863. doi: 10.1038/13158
![]() |
[16] |
M. F. Glasser, T. S. Coalson, E. C. Robinson, C. D. Hacker, J. Harwell, E. Yacoub, et al., A multi-modal parcellation of human cerebral cortex, Nature, 536 (2016), 171–178. doi: 10.1038/nature18933
![]() |
[17] | E. Hohlfeld, L. Mahadevan, Unfolding the sulcus, Phys. Rev. Lett., 106 (2011), 105702. |
[18] |
E. Hohlfeld, L. Mahadevan, Scale and nature of sulcification patterns, Phys. Rev. Lett., 109 (2012), 025701. doi: 10.1103/PhysRevLett.109.025701
![]() |
[19] |
M. A. Holland, K. E. Miller, E. Kuhl, Emerging brain morphologies from axonal elongation, Ann. Biomed. Eng., 43 (2015), 1640–1653. doi: 10.1007/s10439-015-1312-9
![]() |
[20] |
R. C. Knickmeyer, S. Gouttard, C. Kang, D. Evans, K. Wilber, J. Keith Smith, et al., A structural MRI study of human brain development from birth to 2 years, J. Neurosci., 28 (2008), 12176–12182. doi: 10.1523/JNEUROSCI.3479-08.2008
![]() |
[21] | T. Lindeberg, Scale-Space Theory in Computer Vision, Springer Science & Business Media, 2013. |
[22] |
K. L. Mills, C. K. Tamnes, Methods and considerations for longitudinal structural brain imaging analysis across development", author="Mills, Kathryn L and Tamnes, Christian K, Dev. Cognit. Neurosci., 9 (2014), 172–190. doi: 10.1016/j.dcn.2014.04.004
![]() |
[23] |
J. Nie, L. Guo, G. Li, C. Faraco, L. S. Miller, T. Liu, A computational model of cerebral cortex folding, J. Theor. Biol., 264 (2010), 467–478. doi: 10.1016/j.jtbi.2010.02.002
![]() |
[24] |
M. Nonaka-Kinoshita, I. Reillo, B. Artegiani, M. A. Martínez-Martínez, M. Nelson, V. Borrell, et al., Regulation of cerebral cortex size and folding by expansion of basal progenitor, EMBO J., 32 (2013), 1817–1828. doi: 10.1038/emboj.2013.96
![]() |
[25] |
I. Reillo, C. de Juan Romero, M. A. García-Cabezas, V. Borrell, A role for intermediate radial glia in the tangential expansion of the mammalian cerebral cortex, Cereb. Cortex, 21 (2011), 1674–1694. doi: 10.1093/cercor/bhq238
![]() |
[26] |
L. Ronan, N. Voets, C. Rua, A. Alexander-Bloc, M. Hough, C. Mackay, et al., Differential tangential expansion as a mechanism for cortical gyrification, Cereb. Cortex, 24 (2014), 2219–2228. doi: 10.1093/cercor/bht082
![]() |
[27] | A. Serag, P. Aljabar, G. Ball, S. J. Counsell, J. P. Boardman, M. A. Rutherford, et al., Construction of a consistent high-definition spatio-temporal atlas of the developing brain using adaptive kernel regression, Cereb. Cortex, 59 (2012), 2255–2265. |
[28] |
K. L. Spalding, O. Bergmann, K. Alkass, S. Bernard, M. Salehpour, H. B. Huttner, et al., Dynamics of hippocampal neurogenesis in adult humans, Cell, 153 (2013), 1219–1227. doi: 10.1016/j.cell.2013.05.002
![]() |
[29] | G. F. Striedter, S. Srinivasan, E. S. Monuki, Cortical folding: when, where, how, and why?, Ann. Rev. Neurosci., 38 (2015), 291–307. |
[30] | T. Sun, R. F. Hevner, Growth and folding of the mammalian cerebral cortex: from molecules to malformations, Nat. Rev. Neurosci., 15 (2014), 217–232. |
[31] | T. Tallinen, J. S. Biggins, L. Mahadevan, Surface sulci in squeezed soft solids, Phys. Rev. Lett., 15 (2013), 024302. |
[32] |
T. Tallinen, J. Y. Chung, J. S. Biggins, L. Mahadevan, Gyrification from constrained cortical expansion, Proc. Natl. Acad. Sci., 111 (2014), 12667–12672. doi: 10.1073/pnas.1406015111
![]() |
[33] |
T. Tallinen, J. Y. Chung, F. Rousseau, N. Girard, J. Lefèvre, L. Mahadevan, On the growth and form of cortical convolutions, Nat. Phys., 12 (2016), 588–593. doi: 10.1038/nphys3632
![]() |
[34] |
R. Toro, Y. Burnod, A morphogenetic model for the development of cortical convolutions, Cereb. Cortex, 15 (2005), 1900–1913. doi: 10.1093/cercor/bhi068
![]() |
[35] |
R. Toro, M. Perron, B. Pike, L. Richer, S Veillette, Z. Pausova, et al., Brain size and folding of the human cerebral cortex, Cereb. Cortex, 18 (2008), 2352–2357. doi: 10.1093/cercor/bhm261
![]() |
[36] |
G. Xu, A. K. Knutsen, K. Dikranian, C. D. Kroenke, P. V. Bayly, L. A. Taber, Axons pull on the brain, but tension does not drive cortical folding, J. Biomech. Eng., 132 (2010), 071013-1–071013-8. doi: 10.1115/1.4001683
![]() |
[37] | E. Armstrong, K.Ziles, M. Man, A. Schleicher, The ontogeny of cortical folding in the human brain, Soc. Neurosci. Abstr., 14 (1988), 1257. |
[38] |
E. Armstong, A. Schleicher, H. Omran, M. Curtis, K. Ziles, The ontogeny of human gyrification, Cereb. Cortex, 5 (1995), 56–63. doi: 10.1093/cercor/5.1.56
![]() |
[39] | J. G. Chi, E. C. Dooling, F. H. Gilles, Gyral Development of the Human Brain, Ann. Neurol., 1 (1977), 83–93. |
[40] |
M. A. Hofman, Size and shape of the cerebral cortex in mammals, Brain, Behav. Evol., 27 (1985), 28–40. doi: 10.1159/000118718
![]() |
[41] | M. S. van der Knaap, G. van Wezel-Meijle, P. G. Barth, Normal gyration and sulcation in preterm and term neonates: Appearance on MR images, Brain, Behav. Evol., 200 (1985), 389–396. |
[42] |
K. Zilles, E. Armstrong, A. Schlcicher, H. Kreuchmann, The human pattern of gyrification in the cerebral cortex, Anat. Embryol., 179 (1988), 173–179. doi: 10.1007/BF00304699
![]() |
[43] |
K. Zilles, A. Schleicher, C. Langemann, K. Amunts, P. Morosan, N. Palomero-Gallagher, et al., Quantitative analysis of sulci in the human cerebral cortex: Development, regional heterogencity, gender difference, asymmetry, intersubject variabilty and cortical architecture, Hum. Brain Mapp., 5 (1997), 218–221. doi: 10.1002/(SICI)1097-0193(1997)5:4<218::AID-HBM2>3.0.CO;2-6
![]() |
[44] | P. Rakic, Corticogenesis in human and nonhuman primate, Cognit. Neurosci., 1994 (1994), 127–145. |
[45] | P. Rakic, Critical cellular events in cortical development: Opportunities for biophysics, Periodical Biol., 100 (1998), 161–169. |
[46] | A. C. Flint, A. R. Kriegstein, Mechanisms underlying neuronal migration disorders and epilepsy, J. Phys. Chem. Ref. Data, 10 (1997), 92–97. |
[47] | B. Nadarajah, J. G. Parnavelas, Modes of neuronal migration in the developing cerebral cortex, Nat. Rev. Neurosci., 10 (2002), 423–432. |
[48] |
V. S. Caviness, D. N. Kennedy, C. Richelme, The human brain age 7–11 years: A volumetric analysis based on magnetic resonance images, Cereb. Cortex, 6 (1996), 726–736. doi: 10.1093/cercor/6.5.726
![]() |
[49] |
S. F. Witelson, D. L. Kigar, T. Harvey, The exceptional brain of Albert Einstein, Lancet, 353 (1999), 2149–2153. doi: 10.1016/S0140-6736(98)10327-6
![]() |
[50] | S. F. Witelson, D. L. Kigar, T. Harvey, The exceptional brain of Einstein, Recherche, 326 (2000), 31–35. |
[51] |
E. Armstrong, K. Zilles, M. Curtis, Cortical folding, the lunate sulcus and the evolution of the human brain, J. Hum. Evol., 20 (1991), 341–348. doi: 10.1016/0047-2484(91)90014-M
![]() |
[52] |
E. Armstrong, K. Zilles, A. Schleicher, Cortical folding and the evolution of the human brain, J. Hum. Evol., 25 (1993), 387–392. doi: 10.1006/jhev.1993.1057
![]() |
[53] |
R. A. Barton, P. H. Harvey, Mosaic evolution of brain structure in mammals, Nature, 405 (2000), 1055–1058. doi: 10.1038/35016580
![]() |
[54] |
T. W. Deacon, Rethinking mammalian brain evolution, Am. Zool., 30 (1990), 629–705. doi: 10.1093/icb/30.3.629
![]() |
[55] | D. Falk, Hominid Paleoneurology, Ann. Rev. Anthropol., 16 (1985), 13–30. |
[56] |
D. Falk, Apples, oranges and the lunate sulcus, Am. J. Phys. Anthropol., 67 (1985), 313–315. doi: 10.1002/ajpa.1330670403
![]() |
[57] |
D. Falk, Ape-like endocast of ape-man taung, Am. J. Phys. Anthropol., 80 (1989), 335–339. doi: 10.1002/ajpa.1330800307
![]() |
[58] |
D. Falk, C. Hildebolt, M. W. Vannier, Reassesment of the taung early hominid from a neurological perspective, J. Hum. Evol., 18 (1989), 485–492. doi: 10.1016/0047-2484(89)90077-8
![]() |
[59] |
D. Falk, Shifting positions on the lunate sulcus-reply, Am. J. Phys. Anthropol., 84 (1991), 89–91. doi: 10.1002/ajpa.1330840109
![]() |
[60] |
R. L. Holloway, The Taung endocast and the lunate sulcus: a rejection of the hypothesis of its anterior position, Am. J. Phys. Anthropol., 64 (1984), 285–287. doi: 10.1002/ajpa.1330640310
![]() |
[61] |
R. L. Holloway, On Falk 1989 accusations regarding Holloway study of the taung endocast a reply, Am. J. Phys. Anthropol., 84 (1991), 87–88. doi: 10.1002/ajpa.1330840108
![]() |
[62] |
R. L. Holloway, The failure of the gyrification index (GI) to account for volumetric reorganization in the evolution of the human brain, J. Hum. Evol., 22 (1992), 163–170. doi: 10.1016/S0047-2484(05)80005-3
![]() |
[63] | T. M. Mayhew, G. L. M. Mwamengele, V. Dantzer, The gyrification of mammalian cerebral cortex: Quantitive evidence of anisomorphic surface expansion during phylogenetic and otogenetic development, J. Anat., 188 (1996), 53–58. |
[64] | J. Prothero, Small brain, large brain-a quest for nature's scale up rules, J. fur Hirnforsch., 39 (1999), 335–347. |
[65] |
E. Bullmore, M. Brammer, I. Harvey, Cerebral hemispheric asymmetry revisited-effects of handedness, gender and schizophrenia measured by radius of gyration in magnetic resonance images, Psychol. Med., 25 (1995), 349–362. doi: 10.1017/S0033291700036254
![]() |
[66] |
A. Ide, C. Dolezal, M. Fernandez, Hemispheric differences in variabilty of fissural patterns in parasylvian and cingulate regions of human brains, J. Comp. Neurol., 410 (1999), 235–242. doi: 10.1002/(SICI)1096-9861(19990726)410:2<235::AID-CNE5>3.0.CO;2-G
![]() |
[67] |
S. F. Walker, Lateralization of functions in the vertebrate brain-a review, Brit. J. Psychol., 71 (1980), 329–367. doi: 10.1111/j.2044-8295.1980.tb01750.x
![]() |
[68] |
M. C. DelaCoste, D. S. Horvath, D. J. Woodward, Possible sex-differences in the developing human fetal brain, J. Clin. Exper. Neuropsychiatry, 13 (1991), 831–846. doi: 10.1080/01688639108405101
![]() |
[69] | E. Armstrong, M. Curtis, D. P. Buxhocveden, C. Fregoe, K. Zilles, M. F. Casanova, et al., Cortical Gyrification in the rhesus monkey: A test of mechanical folding hypothesis, Cereb. Cortex, 2 (1992), 462–432. |
[70] | T. McBride, S. E. Arnold, R. C. Gur, A comparative volumetric analysis of the prefrontal cortex in human and baboon MRI, Brain Behav. Evol., 2 (1999), 159–166. |
[71] |
R. E. Passingham, Rates of brain development in mammals including man, Brain Behav. Evol., 26 (1985), 167–175. doi: 10.1159/000118773
![]() |
[72] |
M. Mimura, T.Tsujikawa, Aggregating pattern dynamics in a chemotaxis model including growth, Phys. A, 230 (1996), 499–543. doi: 10.1016/0378-4371(96)00051-9
![]() |
[73] |
H. G. Othmer, A. Stevens, Aggregation, blowup and collapse: the ABC's of taxis in reinforced random walks, SIAM J. Appl. Math., 57 (1997), 1044–1081. doi: 10.1137/S0036139995288976
![]() |
[74] | R. Albert, H. G. Othmer, The topology of regulatory interactions predicts the expression pattern of the segment polarity genes in rosophilia melonagester, J. Biophys. Soci., 44 (2003), 79–90. |
[75] | V. V. Gursky, J. Jaeger, K. N. Kozlov, J. Reinitz, A. M. Samsonov, Pattern nuclear divisions are uncoupled in Drosophilia segmentation: comparison of spatially discrete and continuous models, Phys. D, 197 (2003), 286–302. |
[76] |
J. Reinitz, D. H. Sharp, Mechanism of eve stripe formation, Mech. Dev., 49 (1995), 133–158. doi: 10.1016/0925-4773(94)00310-J
![]() |
[77] |
D. H. Sharp, J. Reinitz, Prediction of mutant expression patterns using gene circuits, Biosystems, 47 (1998), 79–90. doi: 10.1016/S0303-2647(98)00014-8
![]() |
[78] |
D. M. Holloway, J. Reinitz, A. Spirov, C. E. Vanario-Alonso, Sharp borders from fuzzy gradients, Trends Genet., 18 (2002), 385–386. doi: 10.1016/S0168-9525(02)02724-5
![]() |
[79] |
E. Mjolsness, D. H. Sharp, J. Reinitz, A Connectionist Model of Development, J. Theor. Biol., 152 (1991), 429–453. doi: 10.1016/S0022-5193(05)80391-1
![]() |
[80] | G. B. Muller, S. A. Newman, Origination of Organismal Form: beyond the gene in developmental and evolutionary biology, MIT Press, 2003. |
[81] |
J. L. Rubenstein, P. Rakic, Genetic control of cortical development, Cereb. Cortex, 9 (1999), 521–523. doi: 10.1093/cercor/9.6.521
![]() |
[82] |
J. L. Rubenstein, S. Anderson, L. M. Shi, Genetic control of cortical regionalization and connectivity, Cereb. Cortex, 9 (1999), 524–532. doi: 10.1093/cercor/9.6.524
![]() |
[83] |
S. L. Rutherford, From genotype to phenotype: buffering mechanisms and the storage of genetic information, BioEssays, 22 (2000), 1095–1105. doi: 10.1002/1521-1878(200012)22:12<1095::AID-BIES7>3.0.CO;2-A
![]() |
[84] | J. D. Murray, Mathematical biology II: spatial models and biomedical applications, Springer, 2003. |
[85] |
G. M. Odell, G. Oster, P. Alberch, B. Burnside, The Mechanical Basis of Morphogenesis, Dev. Biol., 85 (1981), 446–462. doi: 10.1016/0012-1606(81)90276-1
![]() |
[86] | G. F. Oster, J. D. Murray, G. M. Odell, The Formation of Microvilli, Mol. Determ. Anim. Form, 1985 (1985), 365–384. |
[87] | L. Y. Cheng, J. D. Murray, G. M. Odell, G. F. Oster, The Cortical Tractor, Lect. Notes Biomath., 71 (1985), 208–216. |
[88] |
T. Bollenbach, K. Kruse, P. Pantazis, M. Gonzalez-Gaitan, F. Julicher, Robust formation of morphogen gradients, Phys. Rev. Lett., 94 (2005), 018103-1–018103-4. doi: 10.1103/PhysRevLett.94.018103
![]() |
[89] |
J. L. England, J. Cardy, Morphogen gradient from a noisy source, Phys. Rev. Lett., 94 (2005), 078101-1–078101-4. doi: 10.1103/PhysRevLett.94.078101
![]() |
[90] | L. W. Swanson, Brain architecture: understanding the basic plan, Oxford University Press, 2003. |
[91] | S. Kim, In A Mathematical Model of Cerebral Cortical Folding Development Based on a Biomechanical Hypothesis, Conference: Society for Mathematical Biology Annual Meeting, Atlanta, 2015. |
[92] |
R. K. Schofield, G. W. S. Blair, The relationship between viscosity, elasticity and plastic strength of soft materials as illustrated by some mechanical properties of flour doughs I, Proc. R. Soc. London, Ser. A, 138 (1932), 707–718. doi: 10.1098/rspa.1932.0211
![]() |
[93] |
R. K. Schofield, G. W. S. Blair, The relationship between viscosity, elasticity and plastic strength of soft materials as illustrated by some mechanical properties of flour doughs II, Proc. R. Soc. London, Ser. A, 139 (1933), 557–566. doi: 10.1098/rspa.1933.0038
![]() |
[94] |
R. K. Schofield, G. W. S. Blair, The relationship between viscosity, elasticity and plastic strength of soft materials as illustrated by some mechanical properties of flour doughs III, Proc. R. Soc. London, Ser. A, 141 (1933), 72–85. doi: 10.1098/rspa.1933.0104
![]() |
[95] |
G. Hunt, H. Muhlhaus, B. Hobbs, Localized folding of viscoelastic layers, Geol. Rundsch., 85 (1996), 58–64. doi: 10.1007/s005310050052
![]() |
[96] |
H. B. Muhlhaus, H. Sakagushi, B. E. Hobbs, Evolution of three-dimensional fods for a non-Newtonian plate in a viscous medium, Proc. R. Soc. London, Ser. A, 454 (1998), 3121–3143. doi: 10.1098/rspa.1998.0294
![]() |
[97] |
S. M. Schmalholz, Y. Y. Podladchikov, Finite amplitude folding: transition from expotential to layer length controlled growth, Earth Planet. Sci. Lett., 179 (2000), 363–377. doi: 10.1016/S0012-821X(00)00116-3
![]() |
[98] |
C. B. Muratov, V. V. Osipov, Scenarios of domain pattern formation in a reaction-diffusion system, Phys. Rev. B, 54 (1996), 4860–4879. doi: 10.1103/PhysRevE.54.4860
![]() |
[99] |
J. H. E. Cartwright, Labyrinthine Turing pattern formation in the cerebral cortex, J. Theor. Biol., 217 (2002), 97–103. doi: 10.1006/jtbi.2002.3012
![]() |
[100] | M. Ono, S. Kunik, C. D. Abernathy, Atlas of the Cerebral Sulci, Georg Thieme Verlag, 1990. |
[101] | L. B. Leopold, A View of the River, Harvard University Press, 1994. |
[102] | J. A. Adam, Mathematics in Nature: Modeling Patterns in the Natural World, Princeton University Press, 2003. |
[103] |
P. S. Dodds, D. H. Rothman, Unified View of Scaling Laws for River Networks, Phys. Rev. E, 59 (1999), 4865–4877. doi: 10.1103/PhysRevE.59.4865
![]() |
[104] | H. C. Patel, Growth Analysis by Non-linear Continuum Theory, Ph.D thesis, Columbia University, 1983. |
[105] | L. D. Landau, E. M. Lifshitz, Theory of Elasticity, Pergamon, 1970. |
[106] | R. W. Lardner, Mathematical Theory of Dislocations and Fracture, University of Toronto Press, 1974. |
[107] |
E. K. Rodriquez, A. Hoger, A. D. McCulloch, Stress-Dependent Finite Growth in Soft Elastic Tissue, J. Biomech., 27 (1994), 455–467. doi: 10.1016/0021-9290(94)90021-3
![]() |
[108] |
S. M. Klisch, T. J. vanDyke, A. Hoger, A theory of volumetric growth for compressible elastic biological materials, Math. Mech. Solids, 6 (2001), 551–575. doi: 10.1177/108128650100600601
![]() |
[109] | E. Kuhl, A. Menzel, P. Steinmann, Computational Modeling of Growth, Comput. Mech., 88 (2003), 71–88. |
[110] | V. A. Lubarda, Constitutive theories based on th multiplicative decomposition of deformation gradient: Themoelasticity, elastoplasticity, and biomechanics, Am. Soc. Mech. Eng., 57 (2004), 95–108. |
[111] | P. Olver, Applications of Lie Groups to Differential Equations, Springer-Verlag, 1993. |
[112] | P. J. Olver, Equivalence, Invariants and Symmetry, Cambridge University Press, 1995. |
[113] |
W. Nowinski, A. Fang, B. Nguyen, L. Jagannathan, J. Raphel, R. Raghavan, et al., Multiple brain atlas database and atlas-based neuroimaging system, Comput. Aided Surg., 2 (1997), 42–66. doi: 10.3109/10929089709149082
![]() |
[114] |
R. Raghavan, S. R. Ranjan, R. Viswanathan, W. Lawton, A continuum mechanical model for cortical growth, J. Theor. Biol., 187 (1997), 285–296. doi: 10.1006/jtbi.1997.0450
![]() |
[115] | L. V. Ovsiannikov, Group Analysis of Differential Equations, Academic Press, 1982. |
[116] |
N. H. Ibragimov, M. H. Torrisi, A. Valenti, Preliminary group classification of the equations vtt=f(x,vx)vxx+g(x,vx), J. Math. Phys., 32 (1991), 2988–2995. doi: 10.1063/1.529042
![]() |
[117] | I. Lisle, Equivalence transformations for classes of differential equations, Ph.D thesis, The University of British Columbia, 1992. |
[118] | S. Bayer, J. Altman, Atlas of Human Central Nervous System Development, in six volumes, CRC Press, 2008. |
[119] | H. J. Kretschmann, G. Kammradt, E. C. Cowart, The Yakovlev collection- a unique resource for brain research and the basis for a multinational data bank, J. Hirnforsch., 23 (1982), 647–656. |
[120] |
P. Feinsilver, R. Schott, Differential Relations and Recurrence Formulas for Lie Groups, Stud. Appl. Math., 96 (1996), 387–406. doi: 10.1002/sapm1996964387
![]() |
[121] | P. Feinsilver, Lie algebras, Representations, and Analytic Semigroups through Dual Vector Fields, Adv. Math., 2006 (2006). |
[122] |
P. Feinsilver, R. Schott, Vector Fields and their Duals, Adv. Math., 149 (2000), 182–192. doi: 10.1006/aima.1999.1850
![]() |
[123] | P. Feinsilver, R. Schott, Volume I: Representations and Probability Theory, in Algebraic Structures and Operator Calculus: , Springer Science & Business Media, 1993. |
[124] | P. Feinsilver, R. Schott, Volume II: Special Functions and Computer Science, in Algebraic Structures and Operator Calculus, Kluwer Academic Publishers, 1994. |
[125] | P. Feinsilver, R. Schott, Volume III, Representations of Lie Groups, in Algebraic Structures and Operator Calculus, Kluwer Academic Publishers, 1996. |
[126] | F. Alshammari, P. S. Isaac, I. Marquette, A differential operator realisation approach for constructing Casimir operators of non-semisimple Lie algebras, J. Phyics A Math. Theor., 51 (2018), 18. |
1. | Salihu S. Musa, Shi Zhao, Daozhou Gao, Qianying Lin, Gerardo Chowell, Daihai He, Mechanistic modelling of the large-scale Lassa fever epidemics in Nigeria from 2016 to 2019, 2020, 493, 00225193, 110209, 10.1016/j.jtbi.2020.110209 | |
2. | Qianying Lin, Salihu S. Musa, Shi Zhao, Daihai He, Modeling the 2014–2015 Ebola Virus Disease Outbreaks in Sierra Leone, Guinea, and Liberia with Effect of High- and Low-risk Susceptible Individuals, 2020, 82, 0092-8240, 10.1007/s11538-020-00779-y | |
3. | Shi Zhao, Salihu S. Musa, Jay T. Hebert, Peihua Cao, Jinjun Ran, Jiayi Meng, Daihai He, Jing Qin, Modelling the effective reproduction number of vector-borne diseases: the yellow fever outbreak in Luanda, Angola 2015–2016 as an example, 2020, 8, 2167-8359, e8601, 10.7717/peerj.8601 | |
4. | Akhil Kumar Srivastav, Pankaj Kumar Tiwari, Mini Ghosh, Modeling the impact of early case detection on dengue transmission: deterministic vs. stochastic, 2020, 0736-2994, 1, 10.1080/07362994.2020.1804403 | |
5. | Abdullahi Yusuf, Bahar Acay, Umar Tasiu Mustapha, Mustafa Inc, Dumitru Baleanu, Mathematical modeling of pine wilt disease with Caputo fractional operator, 2021, 143, 09600779, 110569, 10.1016/j.chaos.2020.110569 | |
6. | Salihu Sabiu Musa, Shi Zhao, Daihai He, Chao Liu, The Long-Term Periodic Patterns of Global Rabies Epidemics Among Animals: A Modeling Analysis, 2020, 30, 0218-1274, 2050047, 10.1142/S0218127420500479 | |
7. | Salihu Sabiu Musa, Sania Qureshi, Shi Zhao, Abdullahi Yusuf, Umar Tasiu Mustapha, Daihai He, Mathematical modeling of COVID-19 epidemic with effect of awareness programs, 2021, 6, 24680427, 448, 10.1016/j.idm.2021.01.012 | |
8. | Rahul Subramanian, Victoria Romeo-Aznar, Edward Ionides, Claudia T. Codeço, Mercedes Pascual, Predicting re-emergence times of dengue epidemics at low reproductive numbers: DENV1 in Rio de Janeiro, 1986–1990, 2020, 17, 1742-5689, 20200273, 10.1098/rsif.2020.0273 | |
9. | Yongli Cai, Shi Zhao, Yun Niu, Zhihang Peng, Kai Wang, Daihai He, Weiming Wang, Modelling the effects of the contaminated environments on tuberculosis in Jiangsu, China, 2021, 508, 00225193, 110453, 10.1016/j.jtbi.2020.110453 | |
10. | Daihai He, Shi Zhao, Qianying Lin, Salihu S. Musa, Lewi Stone, Christopher M. Barker, New estimates of the Zika virus epidemic attack rate in Northeastern Brazil from 2015 to 2016: A modelling analysis based on Guillain-Barré Syndrome (GBS) surveillance data, 2020, 14, 1935-2735, e0007502, 10.1371/journal.pntd.0007502 | |
11. | G. Udhaya Sankar, C. Ganesa Moorthy, 2020, chapter 4, 9781799821977, 64, 10.4018/978-1-7998-2197-7.ch004 | |
12. | Salihu Sabiu Musa, Shi Zhao, Nafiu Hussaini, Abdulrazaq Garba Habib, Daihai He, Mathematical modeling and analysis of meningococcal meningitis transmission dynamics, 2020, 13, 1793-5245, 2050006, 10.1142/S1793524520500060 | |
13. | Salihu Sabiu Musa, Shi Zhao, Nafiu Hussaini, Salisu Usaini, Daihai He, Dynamics analysis of typhoid fever with public health education programs and final epidemic size relation, 2021, 10, 25900374, 100153, 10.1016/j.rinam.2021.100153 | |
14. | Clarisse Lins de Lima, Ana Clara Gomes da Silva, Cecilia Cordeiro da Silva, Giselle Machado Magalhães Moreno, Abel Guilhermino da Silva Filho, Anwar Musah, Aisha Aldosery, Livia Dutra, Tercio Ambrizzi, Iuri Valério Graciano Borges, Merve Tunali, Selma Basibuyuk, Orhan Yenigün, Tiago Lima Massoni, Kate Jones, Luiza Campos, Patty Kostkova, Wellington P. dos Santos, 2022, Chapter 17, 978-3-030-79752-2, 299, 10.1007/978-3-030-79753-9_17 | |
15. | Umar T. Mustapha, Salihu S. Musa, Muhammad A. Lawan, Aliyu Abba, Evren Hincal, Musa D. Mohammed, Bashir D. Garba, Rabiu B. Yunus, Shehu A. Adamu, Daihai He, Mathematical modeling and analysis of schistosomiasis transmission dynamics, 2021, 12, 1793-9623, 2150021, 10.1142/S1793962321500215 | |
16. | Salihu S. Musa, Isa A. Baba, Abdullahi Yusuf, Tukur A. Sulaiman, Aliyu I. Aliyu, Shi Zhao, Daihai He, Transmission dynamics of SARS-CoV-2: A modeling analysis with high-and-moderate risk populations, 2021, 26, 22113797, 104290, 10.1016/j.rinp.2021.104290 | |
17. | Bilgen Kaymakamzade, Evren Hincal, Nezihal Gokbulut, Tamer Sanlidag, 2022, Chapter 8, 978-3-030-96681-2, 65, 10.1007/978-3-030-96682-9_8 | |
18. | Javier Armando Gutierrez, Karina Laneri, Juan Pablo Aparicio, Gustavo Javier Sibona, Meteorological indicators of dengue epidemics in non-endemic Northwest Argentina, 2022, 7, 24680427, 823, 10.1016/j.idm.2022.10.004 | |
19. | Afeez Abidemi, Hammed Olawale Fatoyinbo, Joshua Kiddy K. Asamoah, Sishu Shankar Muni, 2022, Evaluation of the Efficacy of Wolbachia Intervention on Dengue Burden in a Population: A Mathematical Insight, 978-1-6654-9501-1, 1618, 10.1109/DASA54658.2022.9765106 | |
20. | Emmanuelle Sylvestre, Clarisse Joachim, Elsa Cécilia-Joseph, Guillaume Bouzillé, Boris Campillo-Gimenez, Marc Cuggia, André Cabié, Victor S. Santos, Data-driven methods for dengue prediction and surveillance using real-world and Big Data: A systematic review, 2022, 16, 1935-2735, e0010056, 10.1371/journal.pntd.0010056 | |
21. | RUBAYYI T. ALQAHTANI, ABDULLAHI YUSUF, DEVELOPMENT AND ANALYSIS OF A SEIR MODEL FOR COVID-19 EPIDEMIC WITH VACCINATION AND NONSINGULAR KERNEL, 2022, 30, 0218-348X, 10.1142/S0218348X22400400 | |
22. | Salihu S. Musa, Abdullahi Yusuf, Emmanuel A. Bakare, Zainab U. Abdullahi, Lukman Adamu, Umar T. Mustapha, Daihai He, Unravelling the dynamics of Lassa fever transmission with differential infectivity: Modeling analysis and control strategies, 2022, 19, 1551-0018, 13114, 10.3934/mbe.2022613 | |
23. | Abdurrahman Abdulhamid, Nafiu Hussaini, Salihu S. Musa, Daihai He, Mathematical analysis of Lassa fever epidemic with effects of environmental transmission, 2022, 35, 22113797, 105335, 10.1016/j.rinp.2022.105335 | |
24. | Ziyue Yuan, Salihu Sabiu Musa, Shu-Chien Hsu, Clara Man Cheung, Daihai He, Post pandemic fatigue: what are effective strategies?, 2022, 12, 2045-2322, 10.1038/s41598-022-13597-0 | |
25. | Ziyue Yuan, Shu-Chien Hsu, Clara Man Cheung, Vahid Asghari, Effectiveness of Interventions for Controlling COVID-19 Transmission between Construction Workers and Their Close Contacts, 2022, 38, 0742-597X, 10.1061/(ASCE)ME.1943-5479.0001033 | |
26. | Mohammed M. Al-Shomrani, Salihu S. Musa, Abdullahi Yusuf, Unfolding the Transmission Dynamics of Monkeypox Virus: An Epidemiological Modelling Analysis, 2023, 11, 2227-7390, 1121, 10.3390/math11051121 | |
27. | Xin-You Meng, Chong-Yang Yin, DYNAMICS OF A DENGUE FEVER MODEL WITH UNREPORTED CASES AND ASYMPTOMATIC INFECTED CLASSES IN SINGAPORE, 2020, 2023, 13, 2156-907X, 782, 10.11948/20220111 | |
28. | Alhassan Ibrahim, Usa Wannasingha Humphries, Parinya Sa Ngiamsunthorn, Isa Abdullahi Baba, Sania Qureshi, Amir Khan, Modeling the dynamics of COVID-19 with real data from Thailand, 2023, 13, 2045-2322, 10.1038/s41598-023-39798-9 | |
29. | Abdoulaye Sow, Cherif Diallo, Hocine Cherifi, Benny Yong, Interplay between vaccines and treatment for dengue control: An epidemic model, 2024, 19, 1932-6203, e0295025, 10.1371/journal.pone.0295025 | |
30. | Nawaf L. Alsowait, Mohammed M. Al-Shomrani, Ismail Abdulrashid, Salihu S. Musa, Synergistic Impact of Active Case Detection and Early Hospitalization for Controlling the Spread of Yellow Fever Outbreak in Nigeria: An Epidemiological Modeling and Optimal Control Analysis, 2024, 12, 2227-7390, 3817, 10.3390/math12233817 |
Variable | Interpretation |
Nh(t) | Total population of humans |
Sh(t) | Population of susceptible humans |
Eh(t) | Population of humans exposed to DENV |
Ah(t) | Population of asymptomatic humans |
Ih1(t) | Population of humans with clinical symptoms of the DENV |
Ih2(t) | Population of humans with severe clinical symptoms of the DENV |
Rh(t) | Population of humans recovered from the DENV |
Nm(t) | Total population of adult female mosquitoes |
Sm(t) | population of susceptible adult mosquitoes |
Em(t) | population of exposed adult mosquitoes |
Im(t) | population of infectious adult mosquitoes |
Parameter | Interpretation/Description |
Πh/Πm | Recruitment rate of humans/mosquitoes |
μh | Natural death rate of humans |
μm | Death rate of female adult mosquitoes |
βh | Transmission probability from infectious mosquitoes to susceptible humans |
βm | Transmission probability from infectious humans to susceptible mosquitoes |
a | Mosquito biting rate |
b | Maximum number of bites a human can receive per unit time |
θ | Fraction of infected humans that are exposed |
vh | Progression rate of exposed humans to infectious humans with clinical symptoms |
η | Fraction of infectious humans that will not get a severe symptom |
γh, γ1, γ2 | Recovery rate of infectious humans from Ah, Ih1, Ih2, respectively |
δh | Disease-induced death rate of humans |
σm | Progression rate of exposed mosquitoes to the infectious mosquitoes |
ρ | Reporting rate to actual case ratio |
m | Average mosquito to human ratio |
Parameter | Baseline [Range] | Units | Source |
Πh | 2.5 [1,5] | person~day−1 | [6] |
Πm | 5000 [1000,6000] | Day−1 | estimated from [6] |
μh | 3.9×10−5 [3.6×10−5, 4.0×10−5] | Day−1 | [12] |
μm | 0.05714 [0.01,0.1] | Day−1 | [11] |
βh | 0.75 [0.1,0.95] | Day−1 | [6] |
βm | 0.75 [0.1,0.95] | Day−1 | [6] |
a | 0.4997 [0.1,1] | Day−1 | [13,14] |
b | 0.4 [0.1,1] | Day−1 | [13] |
θ | 0.18 [0.1,0.6] | Dimensionless | [31] |
vh | 0.0666 [0,1] | Day−1 | [6,32] |
η | 0.1 [0.05,0.15] | Dimensionless | [19,33] |
γh | 0.1428 [0.1,0.2] | Day−1 | [13] |
γ1 | 0.2 [0.1428,0.3] | Day−1 | [13] |
γ2 | 0.05 [0.0333,0.07143] | Day−1 | [13] |
δh | 10−3 [10−4,2×10−3] | Day−1 | [6] |
σm | 0.5 [0,1] | Day−1 | [6] |
m | 10 [1,20] | Dimensionless | assumed |
abbreviation | explanation / full name |
AICc | corrected Akaike's Information Criterion |
BB | Backward bifurcation |
CI | Confidence interval |
DENV | Degue virus |
DENV-1 | Dengue virus serotype 1 |
DENV-2 | Dengue virus serotype 2 |
DENV-3 | Dengue virus serotype 3 |
DENV-4 | Dengue virus serotype 4 |
DFE | Disease-free equilibrium |
DHF | Dengue hemorrhagic fever |
DSS | Dengue shock syndrome |
EE | Endemic equilibrium |
EIP | Extrinsic incubation period |
GAS | Globally-asymptotically stable |
HMM | Hidden Markov Model |
IIP | Intrinsic incubation period |
LAS | Locally-asymptotically stable |
MLE | Maximum likelihood estimates |
MLL | Maximum log-likelihood |
NB | Negative binomial (distribution) |
ODE | Ordinary differential equation |
POMP | partially observed Markov process |
PRCC | Partial rank correlation coefficient |
Parameter | Value (Kaohsiung) | Value (Tainan) | Source |
Sh(0) | 0.93212 | 0.8774 | estimated |
Eh(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ah(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ih1(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ih1(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Rh(0) | 0.0688 | 0.1226 | estimated |
β | 1×10−20 | 1×10−20 | estimated |
ρ | 0.0631 | 0.0631 | [53,54,55] |
nm | 9 | 9 | estimated |
Parameter | Value (Kaohsiung) | Value (Tainan) | Source |
Sh(0) | 0.9483 | 0.9161 | estimated |
Eh(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ah(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ih1(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ih1(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Rh(0) | 0.0517 | 0.0839 | estimated |
β | 1.9986 | 1.9986 | estimated |
ρ | 0.0837 | 0.0837 | estimated |
Variable | Interpretation |
Nh(t) | Total population of humans |
Sh(t) | Population of susceptible humans |
Eh(t) | Population of humans exposed to DENV |
Ah(t) | Population of asymptomatic humans |
Ih1(t) | Population of humans with clinical symptoms of the DENV |
Ih2(t) | Population of humans with severe clinical symptoms of the DENV |
Rh(t) | Population of humans recovered from the DENV |
Nm(t) | Total population of adult female mosquitoes |
Sm(t) | population of susceptible adult mosquitoes |
Em(t) | population of exposed adult mosquitoes |
Im(t) | population of infectious adult mosquitoes |
Parameter | Interpretation/Description |
Πh/Πm | Recruitment rate of humans/mosquitoes |
μh | Natural death rate of humans |
μm | Death rate of female adult mosquitoes |
βh | Transmission probability from infectious mosquitoes to susceptible humans |
βm | Transmission probability from infectious humans to susceptible mosquitoes |
a | Mosquito biting rate |
b | Maximum number of bites a human can receive per unit time |
θ | Fraction of infected humans that are exposed |
vh | Progression rate of exposed humans to infectious humans with clinical symptoms |
η | Fraction of infectious humans that will not get a severe symptom |
γh, γ1, γ2 | Recovery rate of infectious humans from Ah, Ih1, Ih2, respectively |
δh | Disease-induced death rate of humans |
σm | Progression rate of exposed mosquitoes to the infectious mosquitoes |
ρ | Reporting rate to actual case ratio |
m | Average mosquito to human ratio |
Parameter | Baseline [Range] | Units | Source |
Πh | 2.5 [1,5] | person~day−1 | [6] |
Πm | 5000 [1000,6000] | Day−1 | estimated from [6] |
μh | 3.9×10−5 [3.6×10−5, 4.0×10−5] | Day−1 | [12] |
μm | 0.05714 [0.01,0.1] | Day−1 | [11] |
βh | 0.75 [0.1,0.95] | Day−1 | [6] |
βm | 0.75 [0.1,0.95] | Day−1 | [6] |
a | 0.4997 [0.1,1] | Day−1 | [13,14] |
b | 0.4 [0.1,1] | Day−1 | [13] |
θ | 0.18 [0.1,0.6] | Dimensionless | [31] |
vh | 0.0666 [0,1] | Day−1 | [6,32] |
η | 0.1 [0.05,0.15] | Dimensionless | [19,33] |
γh | 0.1428 [0.1,0.2] | Day−1 | [13] |
γ1 | 0.2 [0.1428,0.3] | Day−1 | [13] |
γ2 | 0.05 [0.0333,0.07143] | Day−1 | [13] |
δh | 10−3 [10−4,2×10−3] | Day−1 | [6] |
σm | 0.5 [0,1] | Day−1 | [6] |
m | 10 [1,20] | Dimensionless | assumed |
abbreviation | explanation / full name |
AICc | corrected Akaike's Information Criterion |
BB | Backward bifurcation |
CI | Confidence interval |
DENV | Degue virus |
DENV-1 | Dengue virus serotype 1 |
DENV-2 | Dengue virus serotype 2 |
DENV-3 | Dengue virus serotype 3 |
DENV-4 | Dengue virus serotype 4 |
DFE | Disease-free equilibrium |
DHF | Dengue hemorrhagic fever |
DSS | Dengue shock syndrome |
EE | Endemic equilibrium |
EIP | Extrinsic incubation period |
GAS | Globally-asymptotically stable |
HMM | Hidden Markov Model |
IIP | Intrinsic incubation period |
LAS | Locally-asymptotically stable |
MLE | Maximum likelihood estimates |
MLL | Maximum log-likelihood |
NB | Negative binomial (distribution) |
ODE | Ordinary differential equation |
POMP | partially observed Markov process |
PRCC | Partial rank correlation coefficient |
Parameter | Value (Kaohsiung) | Value (Tainan) | Source |
Sh(0) | 0.93212 | 0.8774 | estimated |
Eh(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ah(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ih1(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Ih1(0) | 3.229552×10−6 | 3.241949×10−6 | estimated |
Rh(0) | 0.0688 | 0.1226 | estimated |
β | 1×10−20 | 1×10−20 | estimated |
ρ | 0.0631 | 0.0631 | [53,54,55] |
nm | 9 | 9 | estimated |
Parameter | Value (Kaohsiung) | Value (Tainan) | Source |
Sh(0) | 0.9483 | 0.9161 | estimated |
Eh(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ah(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ih1(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Ih1(0) | 2.773339×10−6 | 3.24338×10−6 | estimated |
Rh(0) | 0.0517 | 0.0839 | estimated |
β | 1.9986 | 1.9986 | estimated |
ρ | 0.0837 | 0.0837 | estimated |