Next Article in Journal
Kinematics in the Information Age
Next Article in Special Issue
A Surface of Section for Hydrogen in Crossed Electric and Magnetic Fields
Previous Article in Journal
L-Fuzzy Sets and Isomorphic Lattices: Are All the “New” Results Really New?
Previous Article in Special Issue
Bifurcation Analysis of a Certain Hodgkin-Huxley Model Depending on Multiple Bifurcation Parameters
 
 
Correction published on 15 November 2018, see Mathematics 2018, 6(11), 252.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of an Age-Structured SIR Epidemic Model with a Reduction Method to ODEs

Graduate School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
Mathematics 2018, 6(9), 147; https://doi.org/10.3390/math6090147
Submission received: 21 July 2018 / Revised: 22 August 2018 / Accepted: 22 August 2018 / Published: 23 August 2018
(This article belongs to the Special Issue Applied Analysis of Ordinary Differential Equations 2018)

Abstract

:
In this paper, we are concerned with the asymptotic stability of the nontrivial endemic equilibrium of an age-structured susceptible-infective-recovered (SIR) epidemic model. For a special form of the disease transmission function, we perform the reduction of the model into a four-dimensional system of ordinary differential equations (ODEs). We show that the unique endemic equilibrium of the reduced system exists if the basic reproduction number for the original system is greater than unity. Furthermore, we perform the stability analysis of the endemic equilibrium and obtain a fourth-order characteristic equation. By using the Routh–Hurwitz criterion, we numerically show that the endemic equilibrium is asymptotically stable in some epidemiologically relevant parameter settings.

1. Introduction

The mathematical modeling of epidemics in human populations has been studied for a long time [1]. In 1760, Bernoulli used a mathematical model of differential equations to discuss the benefit of smallpox inoculation [2]. In 1911, Ross claimed that malaria could be eradicated by reducing the number of mosquitoes, and constructed a mathematical model of differential equations to theoretically support his claim [3]. In 1927, Kermack and McKendrick constructed the first susceptible-infective-recovered (SIR) epidemic model, in which the total population is divided into three classes called susceptible, infective, and recovered [4]. Since their work, the theory of various epidemic models such as a susceptible-infective-susceptible (SIS) epidemic model [5], a susceptible-exposed-infective-recovered (SEIR)epidemic model [6], and a susceptible-infective-recovered-susceptible (SIRS) epidemic model [7] with various structures such as the age structure [8,9], the space structure [10], and the network structure [11] has been developed from both mathematical and epidemiological points of view.
Epidemiologically, the basic reproduction number R 0 for an infectious disease is defined by the expected number of secondary cases produced by a typical infective individual in a completely susceptible population ([9], Chapter 5). Mathematically, R 0 is defined by the spectral radius of a linear operator called the next-generation operator [12], and it determines the complete global dynamics of each equilibrium for some basic epidemic models: if R 0 < 1 , then the trivial disease-free equilibrium is globally asymptotically stable, whereas if R 0 > 1 , then the nontrivial endemic equilibrium is globally asymptotically stable [13]. However, for some epidemic models, the endemic equilibrium can be stable even if R 0 < 1 due to the backward bifurcation [14], and it can be unstable even if R 0 > 1 , which leads to a periodic solution due to the Hopf bifurcation [15].
In [16], some conjectures on the threshold property of R 0 for an age-structured SIR epidemic model were proposed, and they were proved in [17]: if R 0 < 1 , then the disease-free equilibrium is globally asymptotically stable and no endemic equilibrium exists, whereas if R 0 > 1 , then the endemic equilibrium uniquely exists and it is locally asymptotically stable under some additional conditions. However, in general, it is known that the endemic equilibrium cannot always be unique and stable for R 0 > 1 . Several authors have studied some special cases where the endemic equilibrium is unstable, and periodic solutions occur through the Hopf bifurcation for R 0 > 1 [18,19,20,21]. From the viewpoint of application, it is important to investigate when the endemic equilibrium of an age-structured SIR epidemic model is stable and when it is not, as the age distribution of infective individuals in such a model should be stable if one tries to estimate the basic reproduction number R 0 for an endemic disease by using statistical data that exhibit an almost unchanged age distribution of infective individuals. The purpose of this study is to obtain a new condition under which the endemic equilibrium of an age-structured SIR epidemic model is (locally) asymptotically stable.
Age-structured SIR epidemic models as stated above are systems of partial differential equations (PDEs), and hence the stability analysis of them often requires a relatively difficult method such as the spectral theory of positive operators ([17], §5). In this paper, we make some assumptions on the parameters of an age-structured SIR epidemic model, under which we can perform the reduction of the model into a four-dimensional system of ordinary differential equations (ODEs). We can then apply the standard method of characteristic equations for the stability analysis of the endemic equilibrium.
This paper is organized as follows. In Section 2, we formulate an age-structured SIR epidemic model, and perform the reduction of it into a four-dimensional system of ODEs. In Section 3, we prove that the reduced system has a unique endemic equilibrium if the basic reproduction number R 0 for the original system is greater than unity. Moreover, we investigate the asymptotic stability of the endemic equilibrium, and obtain a fourth-order characteristic equation. As its coefficients have quite complex forms, we only prove their positivity, and numerically show by using the Routh–Hurwitz criterion that the endemic equilibrium is asymptotically stable in some epidemiologically relevant parameter settings in Section 4. In the parameter settings, the essential supremum of the demographic mortality rate is determined based on a dataset for Japan in 2015, and the recovery rate γ and the basic reproduction number R 0 are varied for the cases of an influenza-like disease ( γ = 52 and 2 R 0 3 ), a chlamydia-like sexually transmitted disease ( γ = 1 and 1 < R 0 1.5 ), and a wider range of realistic values of them ( 1 / 50 γ 365 and 1 < R 0 50 ). Finally, Section 5 is devoted to the discussion.

2. Reduction of an Age-Structured SIR Epidemic Model into ODEs

We first formulate an age-structured SIR epidemic model. Let S ( t , a ) , I ( t , a ) , and R ( t , a ) denote the susceptible, infective, and recovered populations of age a 0 at time t 0 , respectively. Let b > 0 denote the birth rate, let μ ( · ) L + ( 0 , + ) denote the age-specific mortality rate such that 0 + μ ( a ) d a = + , and let γ > 0 denote the recovery rate. As in [18], we focus on the case where the disease transmission function is only dependent on the age of infective individuals: κ = κ ( · ) L + ( 0 , + ) . In this case, the age-structured SIR epidemic model is formulated as follows:
t + a S ( t , a ) = S ( t , a ) 0 + κ ( a ) I ( t , a ) d a μ ( a ) S ( t , a ) , t > 0 , a > 0 , t + a I ( t , a ) = S ( t , a ) 0 + κ ( a ) I ( t , a ) d a μ ( a ) + γ I ( t , a ) , t > 0 , a > 0 , t + a R ( t , a ) = γ I ( t , a ) μ ( a ) R ( t , a ) , t > 0 , a > 0 , S ( t , 0 ) = b , I ( t , 0 ) = R ( t , 0 ) = 0 , t > 0 , S ( 0 , a ) = S 0 ( a ) , I ( 0 , a ) = I 0 ( a ) , R ( 0 , a ) = R 0 ( a ) , a 0 ,
where S 0 ( · ) , I 0 ( · ) , R 0 ( · ) L + 1 ( 0 , + ) × L + 1 ( 0 , + ) × L + 1 ( 0 , + ) denotes the initial age distributions of each population. It is easy to see that (1) has the demographic steady state P ( a ) = b e 0 a μ ( σ ) d σ , a 0 . Let
s ( t , a ) = S ( t , a ) P ( a ) , i ( t , a ) = I ( t , a ) P ( a ) , r ( t , a ) = R ( t , a ) P ( a ) , t 0 , a 0 .
We then can normalize (1) as follows.
t + a s ( t , a ) = s ( t , a ) 0 + κ ( a ) P ( a ) i ( t , a ) d a , t > 0 , a > 0 , t + a i ( t , a ) = s ( t , a ) 0 + κ ( a ) P ( a ) i ( t , a ) d a γ i ( t , a ) , t > 0 , a > 0 , t + a r ( t , a ) = γ i ( t , a ) , t > 0 , a > 0 , s ( t , 0 ) = 1 , i ( t , 0 ) = r ( t , 0 ) = 0 , t > 0 , s ( 0 , a ) = S 0 ( a ) P ( a ) , i ( 0 , a ) = I 0 ( a ) P ( a ) , r ( 0 , a ) = R 0 ( a ) P ( a ) , a 0 .
As shown in [18], for a specific form of κ ( · ) such that κ ( · ) P ( · ) is sufficiently concentrated in one particular age class, the endemic equilibrium of (2) can be destabilized even if it uniquely exists. Thus, our interest in this paper is when it is stable. In this paper, we assume that κ ( · ) has the following form:
κ ( a ) = β b a e k a e 0 a μ ( σ ) d σ , β > 0 , k > μ = ess . sup a 0 μ ( a ) > 0 , a 0 .
Note that κ ( · ) L + ( 0 , + ) is satisfied under this assumption. In this case, the equations of s and i in (2) can be rewritten as follows.
t + a s ( t , a ) = s ( t , a ) β 0 + a e k a i ( t , a ) d a , t > 0 , a > 0 , t + a i ( t , a ) = s ( t , a ) β 0 + a e k a i ( t , a ) d a γ i ( t , a ) , t > 0 , a > 0 , s ( t , 0 ) = 1 , i ( t , 0 ) = 0 , t > 0 , s ( 0 , a ) = S 0 ( a ) P ( a ) , i ( 0 , a ) = I 0 ( a ) P ( a ) , a 0 .
Note that we can omit the equation of r, as it does not affect the dynamics of (4). As in [17], without loss of generality, we can assume that 0 s ( t , a ) 1 and 0 i ( t , a ) 1 for all t 0 and a 0 .
As seen in [17,18], the stability analysis of age-structured PDE systems such as (4) requires complex calculation. In this paper, we perform the reduction of (4) into a four-dimensional system of ODEs. Let
X ( t ) = 0 + e k a s ( t , a ) d a , Y ( t ) = 0 + e k a i ( t , a ) d a , L ( t ) = β 0 + a e k a s ( t , a ) d a , Λ ( t ) = β 0 + a e k a i ( t , a ) d a , t 0 .
Note that these variables have no specific epidemiological implications except Λ ( t ) , which implies the force of infection at time t 0 . By differentiating X ( · ) , we have
d X ( t ) d t = 0 + e k a s ( t , a ) t d a = 0 + e k a s ( t , a ) a s ( t , a ) Λ ( t ) d a = e k a s ( t , a ) 0 + k 0 + e k a s ( t , a ) d a X ( t ) Λ ( t ) = 1 k + Λ ( t ) X ( t ) , t > 0 .
Similarly, we have
d Y ( t ) d t = 0 + e k a i ( t , a ) t d a = 0 + e k a i ( t , a ) a + s ( t , a ) Λ ( t ) γ i ( t , a ) d a = e k a i ( t , a ) 0 + k 0 + e k a i ( t , a ) d a + X ( t ) Λ ( t ) γ Y ( t ) = X ( t ) Λ ( t ) k + γ Y ( t ) , t > 0 ,
d L ( t ) d t = β 0 + a e k a s ( t , a ) t d a = β 0 + a e k a s ( t , a ) a s ( t , a ) Λ ( t ) d a = β a e k a s ( t , a ) 0 + + β 0 + e k a s ( t , a ) d a k 0 + a e k a s ( t , a ) d a L ( t ) Λ ( t ) = β X ( t ) k + Λ ( t ) L ( t ) , t > 0 ,
and
d Λ ( t ) d t = β 0 + a e k a i ( t , a ) t d a = β 0 + a e k a i ( t , a ) a + s ( t , a ) Λ ( t ) γ i ( t , a ) d a = β a e k a i ( t , a ) 0 + + β 0 + e k a i ( t , a ) d a k 0 + a e k a i ( t , a ) d a + L ( t ) Λ ( t ) γ Λ ( t ) = β Y ( t ) + L ( t ) Λ ( t ) k + γ Λ ( t ) , t > 0 .
Hence, combining (5)–(8), we obtain the following new four-dimensional system of ODEs:
d X ( t ) d t = 1 k + Λ ( t ) X ( t ) , t > 0 , d Y ( t ) d t = X ( t ) Λ ( t ) k + γ Y ( t ) , t > 0 , d L ( t ) d t = β X ( t ) k + Λ ( t ) L ( t ) , t > 0 , d Λ ( t ) d t = β Y ( t ) + L ( t ) Λ ( t ) k + γ Λ ( t ) , t > 0 , X ( 0 ) = X 0 , Y ( 0 ) = Y 0 , L ( 0 ) = L 0 , Λ ( 0 ) = Λ 0 .
Note that X 0 , Y 0 , L 0 , Λ 0 R + 4 since S 0 ( · ) , I 0 ( · ) L + 1 ( 0 , + ) × L + 1 ( 0 , + ) . In what follows, we perform the stability analysis of system (9).

3. Existence, Uniqueness, and Stability of the Endemic Equilibrium

Following the theory in [12], the basic reproduction number R 0 for the original (normalized) system (4) can be calculated as follows.
R 0 = β 0 + a e k a 0 a e γ ( a σ ) d σ d a = β γ 0 + a e k a 1 e γ a d a = β γ a e k a k 0 + + 1 k 0 + e k a d a a e ( k + γ ) a k + γ 0 + 1 k + γ 0 + e ( k + γ ) a d a = β γ 1 k 2 1 k + γ 2 = β 2 k + γ k 2 k + γ 2 .
We now prove that system (9) has the unique endemic equilibrium if R 0 > 1 . Let E : X , Y , L , Λ R + { 0 } 4 denote the endemic equilibrium of system (9). The following equations are satisfied:
1 k + Λ X = 0 , X Λ k + γ Y = 0 , β X k + Λ L = 0 , β Y + L Λ k + γ Λ = 0 .
We prove the following theorem.
Theorem 1.
Let R 0 be defined by (10). If R 0 > 1 , then system (9) has the unique endemic equilibrium E .
Proof. 
By rearranging the first three equations in (11), we have
X = 1 k + Λ , Y = X Λ k + γ = Λ k + γ k + Λ , L = β X k + Λ = β k + Λ 2 .
By substituting the equations of Y and L into the last equation in (11), we have
0 = β Λ k + γ k + Λ + β Λ k + Λ 2 k + γ Λ .
Dividing both sides of this equation by Λ and rearranging the equation, we have
1 = β k + γ 1 k + γ k + Λ + 1 k + Λ 2 .
Let F ( Λ ) be a function defined by the right-hand side of this equation. Since F ( Λ ) is monotonically decreasing to 0 as Λ + and
F ( 0 ) = β k + γ 1 k + γ k + 1 k 2 = β 2 k + γ k 2 k + γ 2 = R 0 > 1 ,
there exists the unique positive root Λ > 0 of Equation (13). Substituting it into the three equations in (12), we obtain the unique endemic equilibrium E . This completes the proof. ☐
To investigate the asymptotic stability of the endemic equilibrium E , we consider the following Jacobian matrix J E around E :
J E = ( k + Λ ) 0 0 X Λ ( k + γ ) 0 X β 0 ( k + Λ ) L 0 β Λ L ( k + γ ) .
From (12), we have
k + Λ = 1 X , L = β X 2 .
Using (14), we derive the characteristic polynomial for E as follows:
λ I J E = λ + k + Λ 0 0 X Λ λ + k + γ 0 X β 0 λ + k + Λ L 0 β Λ λ L + k + γ = λ + k + Λ 0 0 X λ + k λ + k + γ 0 0 β 0 λ + k + Λ β X 2 β β λ + k λ + k + γ = λ + 1 X 0 0 X γ λ + k + γ 0 0 β X λ + 2 0 λ + 1 X 0 0 β λ + k λ + k + γ = λ + 1 X 2 λ + k + γ 2 β X γ λ + 1 X + λ + k λ + k + γ X λ + 2 = λ 4 + a 3 λ 3 + a 2 λ 2 + a 1 λ + a 0 ,
where I denotes the identity matrix and
a 3 = 2 X + 2 k + γ β X 2 , a 2 = 1 X 2 + 4 k + γ X + k + γ 2 β X 2 + ( 2 k + γ ) X , a 1 = 2 k + γ X 2 + 2 k + γ 2 X β X γ + k k + γ X + 2 2 k + γ , a 0 = k + γ 2 X 2 β X γ X + 2 k k + γ .
To apply the Routh–Hurwitz criterion, we prove the following proposition.
Proposition 1.
Let R 0 and a i ( i = 0 , 1 , 2 , 3 ) be defined by (10) and (16), respectively. If R 0 > 1 , then a i > 0 for all i = 0 , 1 , 2 , 3 .
Proof. 
By Theorem 1, the unique endemic equilibrium E exists. From (12) and (13), we have
k + γ = β X k + γ + β X 2 , k + γ 2 = β X + k + γ β X 2 , 1 k X = 1 k k + Λ > 0 .
We then have
a 3 = 2 X + 2 β X k + γ + β X 2 β X 2 = 2 X + 2 β X k + γ + β X 2 > 0 , a 2 = 1 X 2 + 4 X β X k + γ + β X 2 + β X + k + γ β X 2 β X 2 + ( 2 k + γ ) X = 1 X 2 + 4 β k + γ + 3 β X k β X 2 = 1 X 2 + 4 β k + γ + 2 β X + 1 k X β X > 0 , a 1 = 2 β k + γ X + 2 β + 2 β + 2 k + γ β X + γ β X k k + γ β X 2 2 k β X 2 k + γ β X = 2 β k + γ X + 2 β 1 k X + 2 β + k + γ β X k β X k k + γ β X 2 = 2 β k + γ X + 3 β 1 k X + β + k + γ β X 1 k X > 0 , a 0 = β X + k + γ β + γ β 2 k k + γ β X = β X + 2 k + γ β k β 2 k k + γ β X = β X 1 k X + 2 k + γ β 1 k X > 0 .
This completes the proof. ☐
By Proposition 1, it follows from the Routh–Hurwitz criterion [15] (Proof of Theorem 3.1) that the endemic equilibrium E is asymptotically stable if and only if
Δ = a 1 a 2 a 3 a 0 a 3 2 a 1 2 > 0 .
However, it seems that quite a long calculation is needed to show (17) analytically. Instead, in the next section, we show (17) numerically in some epidemiologically relevant parameter settings.

4. Numerical Results

Let the unit time be 1 year. By the definition in (3), k should satisfy inequality k > μ = ess . sup a 0 μ ( a ) > 0 . For the sake of simplicity, we regard a = 100 (years old) as the maximum age of individuals. In a dataset available in [22], the mortality rate μ ( a ) is at most 0.39954 (at a = 100 ) for males in Japan, 2015. Hence, we fix k = 0.4 . By (10), we can determine β for chosen R 0 and γ as follows:
β = k 2 k + γ 2 2 k + γ R 0 .
We first consider an influenza-like disease which has an infectious period of about 1 week (see [23]). Therefore, let γ = 52 so that the average infectious period is 1 / γ = 1 / 52 year = 1 week. Following the estimation result in [24], we vary the value of R 0 from 2 to 3. In Figure 1a, we can confirm that criterion Δ defined in (17) is always positive, and hence, the endemic equilibrium E is asymptotically stable. In fact, in Figure 1b, the force of infection Λ ( t ) converges to the positive equilibrium value Λ = 0.2339 > 0 as time evolves for R 0 = 2.5 .
We next consider a chlamydia-like sexually transmitted disease which has an infectious period of about 1 year (see [25]). Therefore, γ = 1 so that the average infectious period is 1 / γ = 1 year. Based on the estimation result in [25], we vary the value of R 0 from 1 to 1.5 . In Figure 2a, we can confirm that criterion Δ is always positive, and hence the endemic equilibrium E is asymptotically stable. In fact, in Figure 2b, the force of infection Λ ( t ) converges to the positive equilibrium value Λ = 0.0539 > 0 as time evolves for R 0 = 1.25 .
Finally, we vary the values of γ from 1 / 50 to 365 (i.e., the average infectious period 1 / γ is varied from 1 / 365 year = 1 day to 50 years) and R 0 from 1 to 50. In Figure 3, we can confirm that criterion Δ is always positive in this parameter region. Hence, we can conclude that the endemic equilibrium E is asymptotically stable for epidemiologically relevant values of γ and R 0 .

5. Discussion

In this paper, we formulated an age-structured SIR epidemic model and performed its reduction into a four-dimensional system of ODEs under an additional assumption on the disease transmission function κ ( · ) . We proved that if the basic reproduction number R 0 is greater than 1, then the system has the unique endemic equilibrium E (Theorem 1). Moreover, we obtained a fourth-order characteristic equation, and proved that all of its coefficients are positive (Proposition 1). By the Routh–Hurwitz criterion, the endemic equilibrium E is asymptotically stable if and only if Δ > 0 . As it seems difficult to show Δ > 0 analytically, we showed it numerically in Section 4 for some epidemiologically relevant parameters. We showed that Δ > 0 holds for parameters for an influenza-like disease ( γ = 52 and 2 R 0 3 ) and a chlamydia-like disease ( γ = 1 and 1 < R 0 1.5 ). Furthermore, we showed that Δ > 0 holds for a wider region of epidemiologically relevant parameters γ 1 / 50 , 365 and R 0 1 , 50 .
The results in this paper contribute to enlarge the stability region of the endemic equilibrium E for R 0 > 1 to the set where we can perform our reduction method of the PDEs system into the ODEs system. Epidemiologically, our results imply that E can be stable for R 0 > 1 with some realistic parameters, and have broadened the possibilities of application of an age-structured SIR epidemic model for the estimation of R 0 based on the real data of endemic diseases.
In this study, we restricted our attention to the case where the disease transmission function κ = κ ( a ) is only dependent on the age a of infective individuals. Of course, the case where it depends on the ages of both susceptible and infective individuals is more general and epidemiologically realistic. Nevertheless, the stability of the endemic equilibrium of age-structured SIR epidemic models has not been clarified well enough even for the former case. In fact, the global asymptotic stability of the endemic equilibrium E of such models does not generally hold even if the basic reproduction number R 0 is greater than unity, as a special case where E becomes unstable for R 0 > 1 was obtained in [18]. As an important future work, we will seek other forms of κ ( · ) for which we can apply a similar reduction method to ODEs as in this study.

Funding

This research was funded by Grant-in-Aid for Young Scientists (B) of Japan Society for the Promotion of Science (grant number 15K17585).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Bacaër, N. A Short History of Mathematical Population Dynamics; Springer-Verlag: London, UK, 2011. [Google Scholar]
  2. Bernoulli, D. Réflexions sur les avantages de l’inoculation. Mercue de France 1760, 173–190. [Google Scholar]
  3. Ross, R. The Prevention of Malaria, 2nd ed.; John Murray: London, UK, 1911. [Google Scholar]
  4. Kermack, W.O.; McKendrick, A.G. Contributions to the mathematical theory of epidemics—I. Proc. R. Soc. 1927, 115, 700–721. [Google Scholar] [CrossRef]
  5. Busenberg, S.N.; Iannelli, M.; Thieme, H.R. Global behavior of an age-structured epidemic model. SIAM J. Math. Anal. 1991, 22, 1065–1080. [Google Scholar] [CrossRef]
  6. Röst, G.; Wu, J. SEIR epidemiological model with varying infectivity and infinite delay. Math. Biosci. Eng. 2008, 5, 389–402. [Google Scholar] [PubMed]
  7. Enatsu, Y.; Nakata, Y.; Muroya, Y. Global stability of SIRS epidemic models with a class of nonlinear incidence rates and distributed delays. Acta Math. Sci. 2012, 32, 851–865. [Google Scholar] [CrossRef]
  8. Iannelli, M. Mathematical Theory of Age-Structured Population Dynamics; Giardini Editori e Stampatori: Pisa, France, 1995. [Google Scholar]
  9. Inaba, H. Age-Structured Population Dynamics in Demography and Epidemiology; Springer: Singapore, 2017. [Google Scholar]
  10. Ruan, S. Spatial-temporal dynamics in nonlocal epidemiological models. In Mathematics for Life Science and Medicine; Takeuchi, Y., Iwasa, Y., Sato, K., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 97–122. [Google Scholar]
  11. Frasca, M.; Buscarino, A.; Rizzo, A.; Fortuna, L.; Boccaletti, S. Dynamical network model of infective mobile agents. Phys. Rev. E 2006, 74, 036110. [Google Scholar] [CrossRef] [PubMed]
  12. Diekmann, O.; Heesterbeek, J.A.J.; Metz, J.A.J. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 1990, 28, 365–382. [Google Scholar] [CrossRef] [PubMed]
  13. Hethcote, H.W. Qualitative analyses of communicable disease models. Math. Biosci. 1976, 28, 335–356. [Google Scholar] [CrossRef]
  14. Gumel, A.B. Causes of backward bifurcations in some epidemiological models. J. Math. Anal. Appl. 2012, 395, 355–365. [Google Scholar] [CrossRef]
  15. Greenhalgh, D. Hopf bifurcation in epidemic models with a latent period and nonpermanent immunity. Math. Comput. Model. 1997, 25, 85–107. [Google Scholar] [CrossRef]
  16. Greenhalgh, D. Threshold and stability results for an epidemic model with an age-structured meeting rate. IMA J. Math. Med. Biol. 1988, 5, 81–100. [Google Scholar] [CrossRef]
  17. Inaba, H. Threshold and stability results for an age-structured epidemic model. J. Math. Biol. 1990, 28, 411–434. [Google Scholar] [CrossRef] [PubMed]
  18. Thieme, H. Stability change of the endemic equilibrium in age-structured models for the spread of S-I-R type infectious diseases. In Differential Equations Models in Biology, Epidemiology and Ecology; Busenberg, S., Martelli, M., Eds.; Springer: Berlin/Heidelberg, Germany, 1991; Volume 92, pp. 139–158. [Google Scholar]
  19. Andreasen, V. Instability in an SIR-model with age-dependent susceptibility. In Mathematical Population Dynamics: Analysis of Heterogeneity; Arino, O., Axelrod, D., Kimmel, M., Langlais, M., Eds.; Wuerz Publishing: Winnipeg, MB, Canada, 1995; pp. 3–14. [Google Scholar]
  20. Cha, Y.; Iannelli, M.; Milner, F.A. Stability change of an epidemic model. Dyn. Syst. Appl. 2000, 9, 361–376. [Google Scholar]
  21. Franceschetti, A.; Pugliese, A.; Breda, D. Multiple endemic states in age-structured SIR epidemic models. Math. Biosci. Eng. 2012, 9, 577–599. [Google Scholar] [CrossRef] [PubMed]
  22. Ministry of Health, Labour and Welfare, the 22nd Life Tables. Available online: https://www.mhlw.go.jp/english/database /db-hw/lifetb22nd/index.html (accessed on 12 July 2018).
  23. Centers for Disease Control and Prevention, Clinical Signs and Symptoms of Influenza. Available online: https://www. cdc.gov/flu/professionals/acip/clinical.htm (accessed on 13 July 2018).
  24. Mills, C.E.; Robins, J.M.; Lipsitch, M. Transmissibility of 1918 pandemic influenza. Nature 2004, 432, 904–906. [Google Scholar] [CrossRef] [PubMed]
  25. Heijne, J.C.M.; Herzog, S.A.; Althaus, C.L.; Low, N.; Kretzschmar, M. Case and partnership reproduction numbers for a curable sexually transmitted infection. J. Theor. Biol. 2013, 331, 38–47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Numerical confirmation of condition (17) and asymptotic stability of the endemic equilibrium E for the case of an influenza-like disease: (a) Variation of Δ defined in (17) for 2 R 0 3 ; (b) Time variation of Λ ( t ) for R 0 = 2 . 5 .
Figure 1. Numerical confirmation of condition (17) and asymptotic stability of the endemic equilibrium E for the case of an influenza-like disease: (a) Variation of Δ defined in (17) for 2 R 0 3 ; (b) Time variation of Λ ( t ) for R 0 = 2 . 5 .
Mathematics 06 00147 g001
Figure 2. Numerical confirmation of condition (17) and asymptotic stability of the endemic equilibrium E for the case of a chlamydia-like sexually transmitted disease: (a) Variation of Δ defined in (17) for 1 < R 0 1 . 5 ; (b) Time variation of Λ ( t ) for R 0 = 1 . 25 .
Figure 2. Numerical confirmation of condition (17) and asymptotic stability of the endemic equilibrium E for the case of a chlamydia-like sexually transmitted disease: (a) Variation of Δ defined in (17) for 1 < R 0 1 . 5 ; (b) Time variation of Λ ( t ) for R 0 = 1 . 25 .
Mathematics 06 00147 g002
Figure 3. Numerical confirmation of condition (17) for realistic values of recovery rate γ 1 / 50 , 365 and the basic reproduction number R 0 1 , 50 .
Figure 3. Numerical confirmation of condition (17) for realistic values of recovery rate γ 1 / 50 , 365 and the basic reproduction number R 0 1 , 50 .
Mathematics 06 00147 g003

Share and Cite

MDPI and ACS Style

Kuniya, T. Stability Analysis of an Age-Structured SIR Epidemic Model with a Reduction Method to ODEs. Mathematics 2018, 6, 147. https://doi.org/10.3390/math6090147

AMA Style

Kuniya T. Stability Analysis of an Age-Structured SIR Epidemic Model with a Reduction Method to ODEs. Mathematics. 2018; 6(9):147. https://doi.org/10.3390/math6090147

Chicago/Turabian Style

Kuniya, Toshikazu. 2018. "Stability Analysis of an Age-Structured SIR Epidemic Model with a Reduction Method to ODEs" Mathematics 6, no. 9: 147. https://doi.org/10.3390/math6090147

APA Style

Kuniya, T. (2018). Stability Analysis of an Age-Structured SIR Epidemic Model with a Reduction Method to ODEs. Mathematics, 6(9), 147. https://doi.org/10.3390/math6090147

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop