A sharp threshold for disease persistence in host metapopulations
|
|
- Charla Strickland
- 5 years ago
- Views:
Transcription
1 A sharp threshold for disease persistence in host metapopulations Thanate Dhirasakdanon, Horst R. Thieme, and P. van den Driessche Department of Mathematics and Statistics, Arizona State University, Tempe, AZ , USA Department of Mathematics and Statistics, University of Victoria, Victoria, BC, Canada V8W 3P4 May 7, 2007 A sharp threshold is established that separates disease persistence from the extinction of small disease outbreaks in an S E I R S type metapopulation model. The travel rates between patches depend on disease prevalence. The threshold is formulated in terms of a basic replacement ratio (disease reproduction number), R 0, and, equivalently, in terms of the spectral bound of a transmission and travel matrix. Since frequency-dependent (standard) incidence is assumed, the threshold results do not require knowledge of a disease-free equilibrium. As a trade-off, for R 0 > 1, only uniform weak disease persistence is shown in general, while uniform strong persistence is proved for the special case of constant recruitment of susceptibles into the patch populations. For R 0 < 1, Lyapunov s direct stability method shows that small disease outbreaks do not spread much and eventually die out. Keywords: Persistence, extinction, basic reproduction number 2000MSC: 92D30, 37B25 1
2 2 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche 1 Introduction Many compartmental epidemic models focus on how infectious diseases develop with time and ignore the other category that shapes reality, space. Spatial disease spread can be included with distance as a continuous variable leading to partial differential equation models as in [2, 7, 17, 18], to integrodifferential equations as in [25, 26], or to equations of mixed type [21]. (For the spatial spread of epidemic outbreaks see [20, 21, 31] and the references therein.) Alternatively, space can be included as a discrete variable leading to (spatially explicit) metapopulation models as in [24] that consist of a (possibly very large) system of ordinary differential equations (ODEs). The underlying idea is to subdivide the spatial region under consideration into a number of discrete patches, which may represent districts, cities, countries, ponds (for the study of amphibian decline) etc. The disease is carried from one patch to another by individuals traveling between patches. As a simplifying assumption, individuals only change disease status when they are on a patch and not when they travel. (See [5] for a model with disease transmission during travel.) Work on S I S type models [1,23,33 36] and on models for influenza [13, 14] indicates that a patchy environment and travel between patches can influence disease spread in a complicated way and, depending on parameters values that include characteristics of the specific disease, can enhance or stifle disease spread. Our aim is to add to the threshold analysis of the S E I R S type metapopulation model in [23] by discussing dynamic disease persistence and the existence of an endemic equilibrium. We also want to extend the threshold analysis to travel rates that depend on the disease prevalence on the patches and to recruitment rates that depend in addition on the patch population densities (see Section 2). The variable travel rates take into account behavioral changes that may occur as individuals adjust to such factors as the severity of the disease or travel restrictions. We prove the existence of a sharp threshold that separates disease persistence (Section 4 and 5) from extinction of small disease outbreaks (Section 6). This threshold can be formulated in terms of a basic replacement ratio (basic parasite reproduction number), R 0, or, equivalently, in terms of the spectral bound of a transmission and travel matrix (Section 3). By contrast with [23], our results (except for Theorem 6.3) do not depend on the existence of a disease-free equilibrium because our methods fully exploit the crucial choice that disease transmission is modeled by frequency-dependent (standard) incidence. The threshold
3 Disease persistence in host metapopulations 3 condition for density-dependent (mass action) incidence would involve the patch population densities at a unique disease-free equilibrium (see the Discussion). As a trade-off, our approach struggles with the fact that, under frequency-dependent incidence, the disease can drive both host and parasite to extinction if it induces fatalities or reduces fertility [9,28,38]. This is why, under general recruitment, we can only establish uniform weak persistence of the disease if R 0 > 1 (Section 4). For the same reason, we formulate disease persistence in terms of the frequencies rather than the densities of the infectives on the patches. So disease persistence does not necessarily imply persistence of the causative agent (parasite, pathogen). For uniform strong disease persistence, we assume constant recruitment into the patch populations (Section 5) as in [23]. Under this strong assumption, not only the disease but also the host and the parasite persist, and there exists an endemic equilibrium (a question left unanswered in [23]). Alternatively, the assumption could be made that the disease neither induces fatalities nor reduces fertility as has been done in [34, 36] for an S I S type model. The threshold condition R 0 > 1 is sharp as we show that small disease outbreaks do not spread much and eventually die out if R 0 < 1 (Section 6). Since the disease-free dynamics are unclear, we do not use a linearized stability analysis, but Lyapunov s direct stability method. The epidemic metapopulation model introduced in [24] not only keeps track of the patch where an individual is currently located (as in the model analyzed in this paper) but also of the patch on which an individual is born and usually resides. For m epidemiological classes on n patches such models lead to a system of mn 2 rather than mn ODEs. Under the assumption that travel is independent of disease prevalence, they are analyzed for the S I S case in [3], for the S E I R S case in [4], and for a model of SARS that includes quarantine in [22]. An explicit expression for the basic reproduction number is given, and numerical simulations for the S I S model indicate that this number acts as a threshold between extinction and persistence of the disease. It remains to be explored whether the methods we develop in this paper will also work for these even more complex models with frequency-dependent incidence.
4 4 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche 2 The model We consider a host metapopulation that is geographically distributed over n patches (districts, countries). The disease divides each patch population into four classes (or compartments): susceptible individuals (represented by letter S), exposed individuals (infected but not yet infective, represented by letter E), infective individuals (represented by letter I), and removed (or recovered) individuals (represented by letter R). We include the possibility that recovered individuals can become susceptible again, so our model will be of S E I R S type. We assume that disease transmission and transition between disease classes only occurs on patches and not during travel. See [5] for a model with disease transmission during travel. Let S i (t), E i (t), I i (t), R i (t) denote the respective numbers of susceptible, exposed, infective, and removed individuals and N i (t) denote the total population size in patch i at time t. We introduce the following notation for the various epidemiological classes, S =(S 1,..., S n ), E = (E 1,..., E n ), I = (I 1,..., I n ), R =(R 1,..., R n ) and J = (E, I, R). (1) The letter C {S, E, I, R} =: C is used for an epidemiological class. The 3n-tuple J represents the state of the non-susceptible part of the host population. The dynamics of the host population and the disease are described by the following system, for
5 Disease persistence in host metapopulations 5 i = 1,..., n, ds i dt = de i dt = di i dt = dr i dt Λ S i I i i(s, J) κ i µ S i S i + γi R R i N i ( ) + m S ik (J)S k m S ki(j)s i, κ S i I i i ( ) γi E + µ E i Ei + N i γe i E i ( ) γi I + µ I i Ii + = γ I ii i ( γ R i + µ R i N i = S i + E i + I i + R i. ) Ri + ( ) m E ik (J)E k m E ki(j)e i, ( ) m I ik (J)I k m I ki(j)i i, ( ) m R ik (J)R k m R ki(j)r i, (2) Individuals are recruited into the local population of patch i at a rate Λ i (S, J). This may happen by immigration or birth, or, in the case of a sexually transmitted disease, by entering the sexually active part of the patch population. Individuals on patch i and in the class C C die at a per capita rate µ C i and move to the next epidemiological class at a per capita rate γi C. Disease transmission is modeled by frequency-dependent (standard) incidence with κ i being the per capita infection rate on patch i (consult the Discussion for the consequences of this choice). Individuals in epidemiological class C travel from patch k to patch i at a per capita rate m C ik (J). Without loss of generality we can assume that m C kk(j) = 0, k = 1,..., n, C C, J R 3n +, (3) where R + = [0, ). More generally than in [23], the recruitment rate may depend on the state of the population for some of our results. Furthermore individuals may adjust their travel rates to the disease status of the population. The following assumptions are made throughout the paper. Assumption 2.1. The functions Λ i (S, J) are defined for S R n +, J R 3n The functions m C ik (J), C C, are defined for J R3n +. These functions are non-negative and are locally Lipschitz continuous on their respective domains. +.
6 6 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche t i S i (t) E i (t) I i (t) R i (t) N i (t) time patch number number of susceptibles in patch i at time t number of exposed individuals in patch i at time t number of infectives in patch i at time t number of removed individuals in patch i at time t population size of patch i at time t Table 1: Model variables All parameters are positive with the possible exception of γi R of which may be 0. Furthermore some or all (a) (b) Λ i (S, J) lim sup (S,J) (S, J) 0, Λ i (S, J) lim inf > µ S i, (S,J) 0,S 0 S i (c) the matrix (m S ik (0)) is irreducible, (d) the matrix (m E ik (0) + mi ik (0)) is irreducible. All norms are sum-norms, x = x j. So J is the total number of non-susceptibles and (S, J) is the size of the total population. S 0 means that S (0, ) n, i.e., the vector S has all its coordinates positive. Assumption 2.1 (a) is satisfied, e.g., if all recruitment rates are bounded. Assumption 2.1 (b) guarantees that each local population survives in the absence of the disease. By Assumption 2.1 (c) and (d), every patch can be reached from every other patch by susceptible and infected individuals. For N i = 0, we define S ii i N i = 0. The expression S ii i N i is then a locally Lipschitz continuous function of (S i, E i, I i, R i ) R 4 +. We define the total population size as N(t) = N i. (4) Notice that N(t) = (S(t), E(t), I(t), R(t)) = (S(t), J(t)) if all vectors are positive.
7 Disease persistence in host metapopulations 7 µ C i per capita mortality rate of class C in patch i γi C p. cap. transition rate in patch i from class C to the next κ i per capita infection rate on patch i Λ i (S, J) recruitment rate of patch i m C ik (J) p. cap. travel rate of class C from patch k to patch i Table 2: Parameters and parameter functions Theorem 2.2. For all S, E, I, R R n +, there exists a unique solution (S(t), E(t), I(t), R(t)) of (2) with initial data (S, E, I, R ) that is defined for all t 0 and takes values in R 4n +. Further there exists some c > 0 such that, for all non-negative solutions of (2), lim sup t N(t) c and N(t) max{c, N(0)} for all t 0. Proof. By Theorem [30; Thm.A.4], there exists a unique solution with values in R 4n + that is defined on some interval [0, b) with b (0, ]. If b <, then lim sup t b N(t) =. We add all equations in system (2), d N dt = (Λ i (S, J) µ Si S i µ Ei E i µ IiI ) i µ Ri R i. (5) Set µ = min n {µ C i ; C C}, then µ > 0 and d N dt By Assumption 2.1 (a), we can choose c > 0 such that Λ i (S, J) µ N. (6) Λ i (S, J) µ 2n N, whenever (S, J) = N c. d N(t) By (6), µ N(t) whenever t [0, b) and N(t) c. This implies dt 2 N(t) max{c, N(0)} for all t [0, b). So b =, and lim sup N(t) c. t
8 8 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche E= (E 1,..., E n ) I = (I 1,..., I n ) R = (R 1,..., R n ) J = (E, I, R) C {S, E, I, R} =: C C= C i = C N = N i = (S, J) generic epidemiological class total size of epidemiological class C total population size Table 3: Vector and other notation 3 The basic replacement ratio, R 0 To establish a threshold parameter that separates disease persistence from the extinction of small disease outbreaks, the following linear operator B on R n R n = R 2n plays a crucial role: for i = 1,..., n, [B(E, I)] i = κ i I i ( ) γi E + µ E i Ei + [B(E, I)] n+i = γ E i E i (γ I i + µ I i)i i + ( ) m E ik E k m E kie i, ( m I iki k m I kii i ), where m C ik = mc ik (0). The operator B can be represented by a quasi-positive matrix which is irreducible by Assumption 2.1 (d) and can be rewritten as (7) B(x, y) = (D I y P x, D E x Qy) (8) with appropriate linear operators D I, D E, P, Q on R n. The positive operators D I and D E are represented by diagonal matrices with all diagonal entries being positive. The diagonal entries of D I are the per capita infection rates κ i, while the diagonal entries of D E are the per capita transition rates γi E from the exposed to the infectious disease stage. The operator P is represented by a matrix with off-diagonal entries m E ik and diagonal entries n j=1 me ji + µ E i + γi E. For Q, the entries are the same as those for P except that the superscript E is replaced by the superscript I. Let s(b) denote the spectral bound of a linear (bounded) operator B, i.e. the largest real part of its eigenvalues, and let r(b) denote its spectral
9 Disease persistence in host metapopulations 9 radius, i.e. the largest absolute value of its eigenvalues. The operators P and Q have the form D M where D is represented by a diagonal matrix with all diagonal entries being positive and M is represented by ( m C ik ) with C = E for P and C = I for Q. Thus P is represented by a quasi-positive matrix. Further the column sums of P are negative, which implies that s(p ) > 0 and r(d 1 M) < 1 ([6; Thm.6.13], [32], [30; Thm.A.44]). With I denoting the identity operator, P 1 = (D(I D 1 M)) 1 = (D 1 M) j D 1, j=0 and P 1 (and also Q 1 ) is a positive operator represented by a matrix with strictly positive diagonal. Proposition 3.1. With B given by (8), s(b) and r(d I Q 1 D E P 1 ) 1 have the same sign. Proof. Let s(b) > 0. By the Perron-Frobenius theory (see, for example, [27; A] or [30; Sec.A.8]), there exists an eigenvector (v, w) (0, 0) of B with B(v, w) (0, 0). By (8), D I w P v, D E v Qw. Since P 1 and Q 1 are positive operators represented by matrices with strictly positive diagonals, 0 v P 1 D I w, 0 w Q 1 D E v. We substitute one inequality into the other, 0 w Q 1 D E P 1 D I w. We set z = D I w. Then 0 z D I Q 1 D E P 1 z. This implies that the spectral radius r of D I Q 1 D E P 1 strictly exceeds 1. Now let s(b) = 0. By the Perron-Frobenius theory, there exists an eigenvector (v, w) (0, 0) of B with B(v, w) = (0, 0). By (8), D I w = P v, D E v = Qw. Proceeding similarly as before, 0 w = Q 1 D E P 1 D I w. We set z = D I w. Then 0 z = D I Q 1 D E P 1 z. Since z 0, the spectral radius r of D I Q 1 D E P 1 is one. Similarly we show that s(b) < 0 implies r < 1.
10 10 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche As in [23], R 0 = r(d I Q 1 D E P 1 ) (9) can be interpreted as the basic reproduction number (basic replacement ratio) of the disease [6, 32]. R 0 will turn out to be the threshold parameter that separates disease persistence (if R 0 > 1) from the extinction of small outbreaks of the disease (if R 0 < 1). An alternative threshold parameter is s(b), the spectral bound of the travel and transmission operator B given by (7). 4 Uniform weak disease persistence if R 0 > 1 Since our choice of frequency-dependent incidence makes it possible that the disease drives the host (and the causative disease agent) to extinction [9, 28, 38], we formulate disease persistence in terms of the frequencies rather than the densities of infective individuals on the patches. Definition 4.1. The disease is said to be uniformly weakly persistent if there exists some ɛ > 0 such that lim sup t I i (t) N i (t) ɛ (10) for all non-negative solutions of (2) with n (I i(0)+e i (0)) > 0 and S(0) 0. We will show that the disease is uniformly weakly persistent if R 0 > 1. This condition is sharp, as we will show in Section 6 that, if R 0 < 1, small disease outbreaks die out. Uniform weak disease persistence is a stronger concept than instability of the disease-free state which is proved in [23] for R 0 > 1. In our model, which is more general than the one in [23], there is even no clear candidate for the disease-free state. We will use the following Lemma several times. Lemma 4.2. There exist δ 1, δ 2 > 0 such that the following holds for all non-negative solutions of (2) with S(0) > 0: S i (t) If lim sup J(t) δ 1, then lim inf t t S(t 1) δ 2, i = 1,..., n.
11 Disease persistence in host metapopulations 11 Proof. We obtain the following differential inequality from the first subsystem in (2), ds i dt ( ) (µs i + κ i )S i + m S ik (J)S k m S ki(j)s i. For δ 1 > 0 small enough, Assumption 2.1 (c) provides a positive operator A represented by an irreducible n n matrix (α ik ) such that α ik m S ik(j) whenever J 2δ 1. By Theorem 2.2, we can choose c > 0 and t 2 1 large enough such that S = (S 1,..., S n ) satisfies the differential vector inequality ds dt cs + AS, t t 2 1. Here c > 0 can be chosen independently of the solution while t 2 may depend on the solution. We integrate, S(t) e c e A S(t 1), t t 2. Since A is an irreducible non-negative matrix, all entries of e A are positive. Let ξ be their minimum. Then ξ > 0 and S i (t) e c ξ S(t 1) for all t t 2. Set δ 2 = e c ξ. Theorem 4.3. If R 0 > 1, the disease is uniformly weakly persistent. Proof. Suppose the statement is false. Choose an arbitrarily small ɛ > 0. By Definition 4.1, there exists a solution of (2) such that n (E i(0)+i i (0)) > 0 and S(0) 0, but I i (t) lim sup t N i (t) < ɛ. Then S(t) 0 for all t 0 and n E i(t) > 0 and n I i(t) > 0 for all t > 0. Shifting forward in time we can assume that I i (t) N i (t) < ɛ t 0. (11)
12 12 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche By Theorem 2.2, there exists some c > 0 such that lim sup t N(t) c. Shifting forward in time again and increasing c, we can assume that N(t) c for all t 0. For C C, we define Then C = C i = C. We complete the proof in three steps. Ī(t) < ɛc t 0. Step 1 There exists some c > 0 such that lim sup t C(t) < ɛ c, C = E, I, R. To prove step 1, we obtain the following differential inequality from the second subsystem in (2), de i dt κ ii i µ E i E i + ( ) m E ik (J)E k m E ki(j)e i. We add over i and set κ = max n κ i, µ = min n µ E i, This implies For R, the proof is similar. dē dt κ Ī µ Ē. lim sup Ē < κ ɛc. t µ Step 2 There exists some δ 2 > 0 such that lim inf t S i (t) δ 2, i = 1,..., n. To prove step 2, we derive the following differential inequality from the first subsystem in (2) and (11), with S = n S i = S, d S dt ) (Λ i (S, J) (µ Si + κ i ɛ)s i. (12)
13 Disease persistence in host metapopulations 13 By Assumption 2.1 (b), there exist some δ 1, ɛ 1 > 0 such that Λ i (S, J) (µ S i + 2ɛ 1 )S i whenever (S, J) 2δ 1, i = 1,..., n. Let S(t) δ 1. We choose ɛ > 0 so small that, by Step 1, there exists some t 1 > 0 such that J(t) δ 1 for all t t 1. So (S(t), J(t)) 2δ 1 for all t t 1. Choose ɛ > 0 so small that κ i ɛ < ɛ 1. By (12), This implies that d S(t) dt ɛ 1 S(t) whenever t t1, S(t) δ1. S(t) > δ 1, t t 1. (13) Step 2 now follows from Lemma 4.2. Step 3 The contradiction Combining step 1 and step 2 and shifting forward in time, we can assume that S i (t) δ, C i (t) cɛ, and 3 C i(t) N i ĉɛ for i = 1,..., n and C = E, I, R. (t) The positive constants δ, c and ĉ = 3 c/δ do not depend on ɛ. Then S i(t) N i (t) 1 ĉɛ and, from the second and third subsystem in (2), de i dt κ i(1 ĉɛ)i i ( γ E i di i dt = γe i E i ( ) γi I + µ I i Ii + ) + µ E i Ei + ( ) m E ik (J)E k m E ki(j)e i, (14) ( ) m I ik (J)I k m I ki(j)i i. By step 2 and by Theorem 2.2, again shifting forward in time if necessary, we can assume that δ S i (t) N(t) c for all t 0, i = 1,..., n with δ > 0 and c > 0 not depending on ɛ. Since the functions m C ik are continuous, m C ik (J) mc ik (0) = mc ik as J 0. So, for any η (0, 1) we can choose ɛ > 0 small enough such that, by step 1, (1 η) m C ik m C ik(j(t)) (1 + η) m C ik, t 0. By (14), we have a differential vector inequality d dt (E, I) B ɛ(e, I)
14 14 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche with a linear operator B ɛ that is associated with an irreducible quasi-positive matrix by Assumption 2.1 (d). Moreover B ɛ B as ɛ 0, where B is the operator in (8). Since R 0 > 1, s(b) > 0 by Proposition 3.1. The eigenvalues depend continuously on the operator (or representing matrix) and the spectral bound is an eigenvalue, therefore s(b ɛ ) > 0 for sufficiently small ɛ > 0. By the Perron-Frobenius theory (see, for example, [27; A] or [30; Sec.A.8]), we can choose a vector v 0 such that B ɛ v = s(b ɛ )v where denotes the dual operator. Then d dt (E(t), I(t)), v B ɛ(e(t), I(t)), v = (E(t), I(t)), B ɛ v =s(b ɛ ) (E(t), I(t)), v. Since v 0, (E(t), I(t)), v > 0 first for t = 0 and then for all t 0. Since s(b ɛ ) > 0, (E(t), I(t)), v as t. Without loss of generality we can choose v such that v j 1 for j = 1,..., 2n, which implies N(t) (E(t), I(t)), v as t. This contradiction to Theorem 2.2 completes the proof. 5 Uniform strong disease persistence and existence of endemic equilibria Typically, uniform weak implies uniform strong persistence (and the existence of an endemic equilibrium) if the dynamical system has a compact attractor. Unfortunately, Theorem 4.3 has the proviso that S(0) 0 (see Definition 4.1), which must be built into the state space. So the existence of a compact attractor would imply strong host persistence. However, since our model assumes standard (frequency-dependent) incidence, it is possible that the disease drives the host to extinction [9, 28, 38]. A suitable assumption that rules this out can be found in [23], namely that all recruitment rates Λ i are positive constants. We adopt it for this section. Under the assumption of constant recruitment, it is easy to see that the host is uniformly strongly persistent. Lemma 5.1. If all Λ i are positive constants, then, for i = 1,..., n, S i (t) > 0 for all t > 0, and there exist constants δ i > 0 such that lim inf t S i (t) δ i for all non-negative solutions of (2).
15 Disease persistence in host metapopulations 15 Proof. By the first subsystem in (2), ds ( i dt Λ i µ S i + κ i + By the fluctuation method ([12], [30; Prop.A.22]), ) m S ki(j) S i. lim inf t S i(t) µ S i + κ i + Λ i lim sup m S ki(j(t)) t. By Theorem 2.2, there exist ζ i (0, ), independent of the solution, such that lim inf S Λ i i(t) t µ S i + κ =: δ i. i + ζ i If R 0 > 1, the parasite is uniformly strongly persistent as well, at least in total abundance. Since S(t) 0 for t > 0 by Lemma 5.1, the subsequent persistence results do not need the solutions of (2) to satisfy S(0) 0. Theorem 5.2. If R 0 > 1 and all Λ i are positive constants, there exists some ɛ > 0 such that lim inf t C i (t) ɛ, C = E, I, R, for all non-negative solutions of (2) with (E i (0) + I i (0)) > 0. Proof. We apply [30; Thm.A.32]. Let X = {(S, E, I, R); S (0, ) n, E, I, R R n +}. By Lemma 5.1, the solution takes its values in X for t > 0. Define ρ : X R + by I i ρ(s, E, I, R) =. N i
16 16 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche In the language of [30; Sec.A.5], the semiflow Φ induced by the solutions of (2) is uniformly weakly ρ-persistent by Theorem 4.3. The compactness condition in [30; Sec.A.5] follows from Theorem 2.2 and Lemma 5.1. By [30; Thm.A.32], Φ is uniformly strongly ρ-persistent. The claim of this theorem for C = I now follows from Lemma 5.1 and the fact that n (E i(0) + I i (0)) > 0 implies ρ(s(t), E(t), I(t), R(t)) > 0 for all t > 0. Once we have the statement for C = I, it follows for C = E by Lemma 5.1 and also for C = R. Corollary 5.3. If R 0 > 1 and all Λ i are positive constants, there exists an equilibrium state of (2) with n C i > 0 for C = E, I, R. Proof. We apply [37; Thm.1.3.7]. Define X as in the proof of Theorem 5.2 and X 0 = {(S, E, I, R) X; n I i > 0}. Then X 0 is convex and relatively open in X. By Theorems 2.2 and 4.3 and Lemma 5.1, the semiflow induced by the solutions of (2) is uniformly persistent with respect to (X 0, X 0 ) in the language of [37]. By [37; Thm.1.3.7], (2) has an equilibrium in X 0. It easily follows that n C i > 0 not only for C = I, but also for C = E, R. Under a further irreducibility assumption concerning the migration of exposed or infective individuals, the parasite persists uniformly strongly on each patch. Theorem 5.4. Let R 0 > 1 and all Λ i be positive constants. Further assume that for each J R 3n + the matrix (m E ik (J)+mI ik (J)) is irreducible. Then there exists some ɛ > 0 such that lim inf t C i(t) ɛ, C C, i = 1,..., n, for all solutions with n (E i(0) + I i (0)) > 0. Further there exists an (endemic) equilibrium of (2) in (0, ) 4n. Proof. Let the state space X be as in Theorem 5.2. Fix i {1,..., n} and define ρ : X R + by ρ(s, E, I, R) = I i. By Theorem 2.2 and Lemma 5.1, the compactness condition of [30; Thm.A.34] is satisfied. Notice that every total orbit w : R X of Φ is associated with a solution of (2) that is defined for all times and takes value in X. By our irreducibility assumption, ρ(w(0)) > 0 whenever ρ(w(t)) > 0 for all t R. The claim for C = I now follows from [30; Thm.A.34]. For C {E, R}, modify ρ(s, E, I, R) = C i. For C = S, the statement has already been shown in Lemma 5.1. The existence of an equilibrium in (0, ) 4n follows from [37; Thm.1.3.7].
17 Disease persistence in host metapopulations 17 6 Extinction of small outbreaks if R 0 < 1 In order to illustrate that the condition R 0 > 1 is sharp for disease persistence, we show that small disease outbreaks do not spread much and eventually die out if R 0 < 1. Recall that J = (E, I, R) and all norms are sum norms. So J is the total number of individuals that are not susceptible. Theorem 6.1. If R 0 < 1, the following local stability results hold for nonnegative solutions of (2). For every ɛ > 0 there exists some δ > 0 such that J(t) ɛ for all t 0 provided J(0) δ. Further there exists some δ 0 > 0 such that J(t) 0 as t whenever J(0) < δ 0. Proof. We use Lyapunov s direct stability method. The second and third subsystems of (2) imply the inequalities de i dt κ ii i ( γ E i ) + µ E i Ei + di i dt =γe i E i ( ) γi I + µ I i Ii + which can be written as ( ) m E ik (J)E k m E ki(j)e i, (15) ( ) m I ik (J)I k m I ki(j)i i, d dt (E, I) B J(E, I) (16) with the linear operator B J being given by the right hand sides of the system (15). By Assumption 2.1 (d), B J is represented by a quasi-positive irreducible matrix for J = 0. Since R 0 < 1, s(b J ) < 0 for J = 0 by Proposition 3.1. Since the eigenvalues of a matrix depend continuously on the entries of the matrix and s(b J ) is an eigenvalue, we find some δ 1 > 0 and a quasi-positive irreducible matrix B such that s( B) < 0 and B J B whenever J R 3n + and J δ 1. By the Perron-Frobenius theory, there exists an eigenvector v (0, ) 2n such that B v = s( B)v. The eigenvector v serves to construct the local Lyapunov function V = (E, I), v. Let ɛ > 0. Choose δ (0, δ 1 /2). Let J(0) < δ and τ = sup { t > 0; J(s) δ 1 s [0, t] }. (17)
18 18 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche We claim that τ = if δ is chosen small enough. Suppose τ <. Then J(t) δ 1 for all t [0, τ] and d dt (E, I) B(E, I). So d d dt (E, I), v = dt (E, I), v B(E, I), v = (E, I), B v and d dt (E, I), v s( B) (E, I), v 0. (18) So (E, I), v is a decreasing function of t [0, τ] and in particular (E(t), I(t)), v (E(0), I(0)), v. Since v 0, there exists some c 1, only dependent on v, such that (E(t), I(t)) c (E(0), I(0)) c J(0). Let γ I = max n γ I i and ξ = min n (γ R i + µ R i ). Then, for t [0, τ], d dt R γi I ξ R γ I c (E(0), I(0)) ξ R. (19) We integrate this inequality, R(t) R(0) + γi c (E(0), I(0)) ĉ J(0) ξ with ĉ = max { } 1, γi c ξ. So, for t [0, τ], J(t) = (E(t), I(t)) + R(t) (c + ĉ) J(0) (c + ĉ)δ. We choose δ so small that (c + ĉ)δ < δ 1. But then we find some τ > τ such that J(t) δ 1 for all t [0, τ], a contradiction to (17). So τ = and J(t) (c + ĉ)δ < δ 1 for all t 0. To obtain the first part of our stability statement, we choose δ > 0 small enough that (c + ĉ)δ < ɛ. For the second statement we choose δ 0 > 0 such that (c + ĉ)δ 0 = δ 1. Let J(0) < δ 0. Then J(t) < δ 1 for all t 0. By (18), which now holds for all t 0, (E(t), I(t)) 0 as t. By (19) and the fluctuation method, lim sup R(t) γi t ξ lim sup I(t) = 0. t
19 Disease persistence in host metapopulations 19 The preceding result does not preclude that both parasite and host go extinct, which is a possibility as our model uses standard (frequency-dependent) incidence [9, 28, 38]. The next result shows that, if R 0 < 1, the disease dies out and the host persists if the initial numbers of exposed, infective or removed individuals are small in comparison to the initial number of susceptible individuals. Theorem 6.2. If R 0 < 1, there exist ɛ 1 ɛ 0 > 0 such that the following holds for all non-negative solutions of (2): for any ɛ (0, ɛ 1 ] there exists some δ > 0 such that S(t) ɛ t 0, S(0) ɛ J(0) δ } = E i (t) + I i (t) + R i (t) N i (t) t 0, lim inf t S(t) ɛ 1, lim inf t S i(t) ɛ 0, i = 1,..., n. Proof. Let η > 0 to be chosen later. By Theorem 6.1, there exists some δ > 0 such that J(t) η for all t 0 whenever J(0) < δ. Let J(0) δ. We add the first subsystem in (2), d dt S ( Λi (S, J) µ S i S i κ i S i I i N i ). (20) Then d dt S ( Λi (S, J) µ S S i I i ) i S i κ i S i + I i ( Λi (S, J) µ S S I ) i S i κ i S + I ( Λi (S, J) µ S i S i ) κ S η S + η with κ = n κ i. By Assumption 2.1 (b), there exist ɛ 1, ɛ 2 > 0 such that Λ i (S, J) µ S i S i ɛ 2 S i whenever (S, J) < 2ɛ 1. (21)
20 20 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche We choose η < ɛ 1. Then ( Λi (S, J) µ S i S i ) ɛ 2 S whenever S ɛ 1, J η. Let ɛ (0, ɛ 1 ] and S(0) > ɛ. Let t be a time at which S(t) = ɛ if such a time exists. Then d ( dt S(t) S(t) ɛ 2 κ η ). ɛ + η We now choose η small enough (in dependence of ɛ) so that the expression in the parentheses above is positive. Then d S(t) > 0 whenever S(t) = ɛ. dt So S(t) ɛ for all t 0 and, by Lemma 4.2, lim inf t S i (t) > 0, i = 1,..., n. By Theorem 6.1, By (20) and (21), E i (t) + I i (t) + R i (t) N i (t) lim inf t S(t) ɛ 1. t 0. If η > 0 is small enough (recall J(t) η for all t 0), we obtain from Lemma 4.2 that lim inf t S i (t) ɛ 0, i = 1,..., n, for some ɛ 0 > 0 that does not depend on the solution. We now assume as in Section 5 that the recruitment rates Λ i are positive constants for i = 1,..., n. Then the disease-free system, E i = I i = R i = 0 and d dt S i =Λ i µ S i S + ( m S S ik k m S S ki i ), i = 1,..., n, m S ik = m S ik(0), (22) has a globally stable equilibrium S (0, ) n. Indeed, (22) can be rewritten as d S = Λ + AS with the operator A having a spectral bound s(a) < 0. By dt the variation of constants formula, S(t) = e ta S(0) + t 0 e sa Λds A 1 Λ 0, t.
21 Disease persistence in host metapopulations 21 Theorem 6.3. If R 0 < 1 and all Λ i are positive constants, there exists some ɛ 1 > 0 with the following property: for any ɛ (0, ɛ 1 ] there exists some δ > 0 such that } { S(0) ɛ J(t) δ, S(t) ɛ t 0, J(0) δ = J(t) t 0, S(t) t S. Proof. By Theorems 6.1 and 6.2, the first subsystem (the one for S) of system (2) is asymptotically autonomous with limit system (22). The claim follows from the global stability of S for (22) and results on asymptotically autonomous systems in [16] (or [29]). If the travel matrices (m C ik ) are independent of J, the disease-free equilibrium (S, 0) is globally asymptotically stable if R 0 < 1 [23; Thm.2.2]. 7 Discussion The sharp threshold results in this paper that separate disease persistence (if R 0 > 1) from the extinction of small disease outbreaks (if R 0 < 1) crucially depend on modeling disease transmission by frequency-dependent (standard) incidence. If frequency-dependent incidence were replaced by density-dependent (mass action) incidence, R 0 would involve the population densities on all patches at a disease-free equilibrium, which would need to be unique and globally attracting for the disease-free population dynamics. As a trade-off, for frequency-dependent incidence, it is possible that the disease drives the host and the causative disease agent to extinction [9, 28, 38] while host persistence is typically automatic for density-dependent incidence. For a complex model like this, it seems very hard to establish a threshold for host persistence unless one makes assumptions such as the one made in Section 5 and [23], constant recruitment into each patch population, or the one made in [34, 36], no disease fatalities and no fertility reduction. Strangely enough, this problem has also repercussions for disease persistence: we can only prove uniform weak disease persistence in general, and host persistence is required in order to go from uniform weak to uniform strong disease persistence under the assumptions that we have imposed. We refer to [10, 15] and the references therein for the discussion as to whether and when frequency-dependent or density-dependent incidence is more adequate. The choice of one over the other has not only consequences
22 22 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche for whether or not the disease can eradicate the host, but also for the occurrences of undamped oscillations: they seem to occur more often for densitydependent incidence than for frequency-dependent incidence [8, 11, 19]. Acknowledgements The authors thank two anonymous referees for helpful comments. This study was partially supported by NSF grant DMS (T.D. and H.R.T.) and by NSERC and MITACS (P.v.d.D). References [1] Allen, L.J.S., B.M. Bolker, Y. Lou, A.L. Nevai, Asymptotic profiles of the steady states for an SIS epidemic patch model, SIAM J. Appl. Math., in press. [2] Allen, L.J.S., B.M. Bolker, Y. Lou, A.L. Nevai, Asymptotic profiles of the steady states for an SIS epidemic reaction-diffusion model, Discrete Contin. Dyn. Syst., in press. [3] Arino, J., P. van den Driessche, A multi-city epidemic model, Mathematical Population Studies 10 (2003), [4] Arino, J., P. van den Driessche, The basic reproduction number in a multi-city compartmental epidemic model, LN Control and Information Sciences 294 (2003), [5] Cui, J., Y. Takeuchi, Y. Saito, Spreading disease with transport-related infection. J. Theor. Biol. 239 (2006), [6] Diekmann, O., J.A.P. Heesterbeek, Mathematical Epidemiology of Infectious Diseases. Model Building, Analysis, and Interpretation, Wiley, Chichester [7] Fitzgibbon, W. E., M. Langlais, A diffusive S.I.S. model describing the propagation of F.I.V., Commun. Appl. Anal. 7 (2003), [8] Gao, L.Q., J. Mena-Lorca, H.W. Hethcote, Variations on a theme of SEI endemic models, Differential Equations and Applications to Biology
23 Disease persistence in host metapopulations 23 and Industry (M. Martelli, C.L. Cooke, E. Cumberbatch, B. Tang, H.R. Thieme, eds.), , World Scientific, Singapore [9] Greenhalgh, D., R. Das, Modelling epidemics with variable contact rates, Theor. Pop. Biol. 47 (1995), [10] Hethcote, H.W., The mathematics of infectious diseases, SIAM Review 42 (2000), [11] Hethcote, H.W., W. Wang, Y. Li, Species coexistence and periodicity in host-host-pathogen models, J. Math. Biol. 51 (2005), [12] Hirsch, W.M., H. Hanisch, J.-P. Gabriel, Differential equation models for some parasitic infections: methods for the study of asymptotic behavior, Comm. Pure Appl. Math. 38 (1985), [13] Hsieh, Y.-H., P. van den Driessche, L. Wang, Impact of travel between patches for spatial spread of disease, Bull. Math. Biol. (to appear), Online First (DOI) /s [14] Hyman, J.M, T. LaForce, Modeling the spread of influenza among cities, Bioterrorism: Mathematical Modeling Applications in Homeland Security, (H.T. Banks, C. Castillo-Chavez, eds.), , SIAM, Philadelphia [15] McCallum, H., N. Barlow, J. Hone, How should pathogen transmission be modelled?, Trends Ecol. Evol. 16 (2001), [16] Markus, L., Asymptotically autonomous differential systems, Contributions to the Theory of Nonlinear Oscillations III (S. Lefschetz, ed.), Ann. Math. Stud. 36, 17-29, Princeton University Press, Princeton [17] Murray, J.D., Mathematical Biology, Springer, Berlin Heidelberg [18] Murray, J.D., Mathematical Biology II: Spatial Models and Biomedical Application, 3rd edition, Springer, Berlin Heidelberg [19] Pugliese, A., An S E I epidemic model with varying population size, Differential Equations Models in Biology, Epidemiology and Ecology, (S. Busenberg, M. Martelli, eds.), , Lecture Notes in Biomathematics 92, Springer, Berlin Heidelberg 1991.
24 24 T. Dhirasakdanon, H.R. Thieme, P. van den Driessche [20] Rass, L., J. Radcliffe, Spatial Deterministic Epidemics, AMS, [21] Ruan, S., Spatial-temporal dynamics in non-local epidemiological models, Mathematics for Life Science and Medicine (Y. Takeuchi, Y. Iwasa, K. Sato, eds.), , Springer, Berlin Heidelberg [22] Ruan, S., W. Wang, S.A. Levin, The effect of global travel on the spread of SARS, Math. Biosci. Eng. 3 (2006), [23] Salmani, M., P. van den Driessche, A model for disease transmission in a patchy environment, Discrete Contin. Dyn. Syst. Ser. B 6 (2006), [24] Sattenspiel, L., K. Dietz, A structured epidemic model incorporating geographic mobility among regions, Math. Biosci. 128 (1995), [25] Schumacher, K., Travelling-front solutions for integro-differential equations. I, J. Reine Angew. Math. 316 (1980), [26] Schumacher, K., Travelling-front solutions for integro-differential equations. II, Biological Growth and Spread: Mathematical Theories and Applications (Proc. Conf., Heidelberg 1979, W. Jäger, H. Rost, P. Tautu, eds.), , Lecture Notes in Biomath. 38, Springer, Berlin New York [27] Smith, H.L., P. Waltman, The Theory of the Chemostat. Dynamics of Microbial Competition, Cambridge University Press, Cambridge [28] Thieme, H.R., Epidemic and demographic interaction in the spread of potentially fatal diseases in growing populations, Math. Biosci. 111 (1992), [29] Thieme, H.R., Convergence results and a Poincaré-Bendixson trichotomy for asymptotically autonomous differential equations, J. Math. Biol. 30 (1992), [30] Thieme, H.R., Mathematics in Population Biology, Princeton University Press, Princeton [31] Thieme, H.R. (book report), Spatial Deterministic Epidemics, by Linda Rass and John Radcliffe, Amer. Math. Soc. 2003, Math. Biosci. 202 (2006),
25 Disease persistence in host metapopulations 25 [32] van den Driessche, P., J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci. 180 (2002), [33] Wang, W., Population dispersal and disease spread, Discrete Contin. Dyn. Syst. Ser. B 4 (2004), [34] Wang, W., Epidemic models with population dispersal, Mathematics for Life Science and Medicine (Y. Takeuchi, Y. Iwasa, K. Sato, eds.), 67-95, Springer, Berlin Heidelberg [35] Wang, W., G. Mulone, Threshold of disease transmission in a patch environment, J. Math. Anal. Appl. 285 (2003), [36] Wang, W., X.-Q. Zhao, An epidemic model in a patchy environment, Math. Biosci. 190 (2004), [37] Zhao, X.-Q., Dynamical Systems in Population Biology, Springer, New York [38] Zhou, J., H.W. Hethcote, Population size dependent incidence in models for diseases without immunity, J. Math. Biol. 32 (1994),
Apparent paradoxes in disease models with horizontal and vertical transmission
Apparent paradoxes in disease models with horizontal and vertical transmission Thanate Dhirasakdanon, Stanley H.Faeth, Karl P. Hadeler*, Horst R. Thieme School of Life Sciences School of Mathematical and
More informationGlobal Stability of SEIRS Models in Epidemiology
Global Stability of SRS Models in pidemiology M. Y. Li, J. S. Muldowney, and P. van den Driessche Department of Mathematics and Statistics Mississippi State University, Mississippi State, MS 39762 Department
More informationIntroduction to SEIR Models
Department of Epidemiology and Public Health Health Systems Research and Dynamical Modelling Unit Introduction to SEIR Models Nakul Chitnis Workshop on Mathematical Models of Climate Variability, Environmental
More informationCompetitive Exclusion in a Discrete-time, Size-structured Chemostat Model
Competitive Exclusion in a Discrete-time, Size-structured Chemostat Model Hal L. Smith Department of Mathematics Arizona State University Tempe, AZ 85287 1804, USA E-mail: halsmith@asu.edu Xiao-Qiang Zhao
More informationA REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 A REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES KING-YEUNG LAM
More informationThe effect of population dispersal on the spread of a disease
J. Math. Anal. Appl. 308 (2005) 343 364 www.elsevier.com/locate/jmaa The effect of population dispersal on the spread of a disease Yu Jin, Wendi Wang Department of Mathematics, Southwest China Normal University,
More informationGLOBAL STABILITY OF SIR MODELS WITH NONLINEAR INCIDENCE AND DISCONTINUOUS TREATMENT
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 304, pp. 1 8. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY
More informationThursday. Threshold and Sensitivity Analysis
Thursday Threshold and Sensitivity Analysis SIR Model without Demography ds dt di dt dr dt = βsi (2.1) = βsi γi (2.2) = γi (2.3) With initial conditions S(0) > 0, I(0) > 0, and R(0) = 0. This model can
More informationStability of SEIR Model of Infectious Diseases with Human Immunity
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 6 (2017), pp. 1811 1819 Research India Publications http://www.ripublication.com/gjpam.htm Stability of SEIR Model of Infectious
More informationGLOBAL STABILITY OF A VACCINATION MODEL WITH IMMIGRATION
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 92, pp. 1 10. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY
More informationUNIFORM WEAK IMPLIES UNIFORM STRONG PERSISTENCE FOR NON-AUTONOMOUS SEMIFLOWS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 127, Number 8, Pages 2395 2403 S 0002-9939(99)05034-0 Article electronically published on April 15, 1999 UNIFORM WEAK IMPLIES UNIFORM STRONG PERSISTENCE
More informationMathematical models on Malaria with multiple strains of pathogens
Mathematical models on Malaria with multiple strains of pathogens Yanyu Xiao Department of Mathematics University of Miami CTW: From Within Host Dynamics to the Epidemiology of Infectious Disease MBI,
More informationGlobal compact attractors and their tripartition under persistence
Global compact attractors and their tripartition under persistence Horst R. Thieme (joint work with Hal L. Smith) School of Mathematical and Statistical Science Arizona State University GCOE, September
More informationImpact of Travel Between Patches for Spatial Spread of Disease
Impact of Travel Between Patches for Spatial Spread of Disease Ying-Hen Hsieh Department of Applied Mathematics National Chung Hsing University Taichung, Taiwan P. van den Driessche Department of Mathematics
More informationPermanence Implies the Existence of Interior Periodic Solutions for FDEs
International Journal of Qualitative Theory of Differential Equations and Applications Vol. 2, No. 1 (2008), pp. 125 137 Permanence Implies the Existence of Interior Periodic Solutions for FDEs Xiao-Qiang
More informationGLOBAL DYNAMICS OF A MATHEMATICAL MODEL OF TUBERCULOSIS
CANADIAN APPIED MATHEMATICS QUARTERY Volume 13, Number 4, Winter 2005 GOBA DYNAMICS OF A MATHEMATICA MODE OF TUBERCUOSIS HONGBIN GUO ABSTRACT. Mathematical analysis is carried out for a mathematical model
More informationDisease Spread in Metapopulations
Fields Institute Communications Volume 00, 0000 Disease Spread in Metapopulations Julien Arino Department of Mathematics, University of Manitoba, Winnipeg, MB, Canada R3T 22, arinoj@cc.umanitoba.ca P.
More informationGlobal Analysis of an Epidemic Model with Nonmonotone Incidence Rate
Global Analysis of an Epidemic Model with Nonmonotone Incidence Rate Dongmei Xiao Department of Mathematics, Shanghai Jiaotong University, Shanghai 00030, China E-mail: xiaodm@sjtu.edu.cn and Shigui Ruan
More informationThe dynamics of disease transmission in a Prey Predator System with harvesting of prey
ISSN: 78 Volume, Issue, April The dynamics of disease transmission in a Prey Predator System with harvesting of prey, Kul Bhushan Agnihotri* Department of Applied Sciences and Humanties Shaheed Bhagat
More informationMATHEMATICAL MODELS Vol. III - Mathematical Models in Epidemiology - M. G. Roberts, J. A. P. Heesterbeek
MATHEMATICAL MODELS I EPIDEMIOLOGY M. G. Roberts Institute of Information and Mathematical Sciences, Massey University, Auckland, ew Zealand J. A. P. Heesterbeek Faculty of Veterinary Medicine, Utrecht
More informationResearch Article A Delayed Epidemic Model with Pulse Vaccination
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume 2008, Article ID 746951, 12 pages doi:10.1155/2008/746951 Research Article A Delayed Epidemic Model with Pulse Vaccination
More informationHETEROGENEOUS MIXING IN EPIDEMIC MODELS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 2, Number 1, Spring 212 HETEROGENEOUS MIXING IN EPIDEMIC MODELS FRED BRAUER ABSTRACT. We extend the relation between the basic reproduction number and the
More informationA GRAPH-THEORETIC APPROACH TO THE METHOD OF GLOBAL LYAPUNOV FUNCTIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 136, Number 8, August 2008, Pages 2793 2802 S 0002-993908)09341-6 Article electronically published on March 27, 2008 A GRAPH-THEORETIC APPROACH TO
More informationGlobal stability of an SEIS epidemic model with recruitment and a varying total population size
Mathematical Biosciences 170 2001) 199±208 www.elsevier.com/locate/mbs Global stability of an SEIS epidemic model with recruitment and a varying total population size Meng Fan a, Michael Y. Li b, *, Ke
More informationReproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission
Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission P. van den Driessche a,1 and James Watmough b,2, a Department of Mathematics and Statistics, University
More informationCan multiple species of Malaria co-persist in a region? Dynamics of multiple malaria species
Can multiple species of Malaria co-persist in a region? Dynamics of multiple malaria species Xingfu Zou Department of Applied Mathematics University of Western Ontario London, Ontario, Canada (Joint work
More informationThe Dynamic Properties of a Deterministic SIR Epidemic Model in Discrete-Time
Applied Mathematics, 05, 6, 665-675 Published Online September 05 in SciRes http://wwwscirporg/journal/am http://dxdoiorg/046/am056048 The Dynamic Properties of a Deterministic SIR Epidemic Model in Discrete-Time
More informationGLOBAL DYNAMICS OF A TIME-DELAYED DENGUE TRANSMISSION MODEL
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 2, Number 1, Spring 212 GLOBAL DYNAMICS OF A TIME-DELAYED DENGUE TRANSMISSION MODEL Dedicated to Herb Freedman on the occasion of his 7th birthday ZHEN WANG
More informationStability Analysis of a SIS Epidemic Model with Standard Incidence
tability Analysis of a I Epidemic Model with tandard Incidence Cruz Vargas-De-León Received 19 April 2011; Accepted 19 Octuber 2011 leoncruz82@yahoo.com.mx Abstract In this paper, we study the global properties
More informationA multi-species epidemic model with spatial dynamics
A multi-species epidemic model with spatial dynamics Julien Arino Jonathan R. Davis David Hartley Richard Jordan Joy M. Miller P. van den Driessche Version of January 6, 25 Abstract A model is formulated
More informationMathematical Modeling and Analysis of Infectious Disease Dynamics
Mathematical Modeling and Analysis of Infectious Disease Dynamics V. A. Bokil Department of Mathematics Oregon State University Corvallis, OR MTH 323: Mathematical Modeling May 22, 2017 V. A. Bokil (OSU-Math)
More informationGLOBAL STABILITY OF THE ENDEMIC EQUILIBRIUM OF A TUBERCULOSIS MODEL WITH IMMIGRATION AND TREATMENT
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 19, Number 1, Spring 2011 GLOBAL STABILITY OF THE ENDEMIC EQUILIBRIUM OF A TUBERCULOSIS MODEL WITH IMMIGRATION AND TREATMENT HONGBIN GUO AND MICHAEL Y. LI
More informationLAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC
LAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC R. G. DOLGOARSHINNYKH Abstract. We establish law of large numbers for SIRS stochastic epidemic processes: as the population size increases the paths of SIRS epidemic
More informationBehavior Stability in two SIR-Style. Models for HIV
Int. Journal of Math. Analysis, Vol. 4, 2010, no. 9, 427-434 Behavior Stability in two SIR-Style Models for HIV S. Seddighi Chaharborj 2,1, M. R. Abu Bakar 2, I. Fudziah 2 I. Noor Akma 2, A. H. Malik 2,
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Centre for Research on Inner City Health St Michael s Hospital Toronto On leave from Department of Mathematics University of Manitoba Julien Arino@umanitoba.ca
More information6. Age structure. for a, t IR +, subject to the boundary condition. (6.3) p(0; t) = and to the initial condition
6. Age structure In this section we introduce a dependence of the force of infection upon the chronological age of individuals participating in the epidemic. Age has been recognized as an important factor
More informationA comparison of delayed SIR and SEIR epidemic models
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 2, 181 190 181 A comparison of delayed SIR and SEIR epidemic models Abdelilah Kaddar a, Abdelhadi Abta b, Hamad Talibi Alaoui b a Université
More informationA Mathematical Model for the Spatial Spread of HIV in a Heterogeneous Population
A Mathematical Model for the Spatial Spread of HIV in a Heterogeneous Population Titus G. Kassem * Department of Mathematics, University of Jos, Nigeria. Abstract Emmanuel J.D. Garba Department of Mathematics,
More informationMODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS
Journal of Biological Systems, Vol. 11, No. 3 2003 325 335 c World Scientific Publishing Company MODELING AND ANALYSIS OF THE SPREAD OF CARRIER DEPENDENT INFECTIOUS DISEASES WITH ENVIRONMENTAL EFFECTS
More informationDynamics of Disease Spread. in a Predator-Prey System
Advanced Studies in Biology, vol. 6, 2014, no. 4, 169-179 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/asb.2014.4845 Dynamics of Disease Spread in a Predator-Prey System Asrul Sani 1, Edi Cahyono
More informationThe death of an epidemic
LECTURE 2 Equilibrium Stability Analysis & Next Generation Method The death of an epidemic In SIR equations, let s divide equation for dx/dt by dz/ dt:!! dx/dz = - (β X Y/N)/(γY)!!! = - R 0 X/N Integrate
More informationSTABILITY ANALYSIS OF A GENERAL SIR EPIDEMIC MODEL
VFAST Transactions on Mathematics http://vfast.org/index.php/vtm@ 2013 ISSN: 2309-0022 Volume 1, Number 1, May-June, 2013 pp. 16 20 STABILITY ANALYSIS OF A GENERAL SIR EPIDEMIC MODEL Roman Ullah 1, Gul
More informationQualitative Analysis of a Discrete SIR Epidemic Model
ISSN (e): 2250 3005 Volume, 05 Issue, 03 March 2015 International Journal of Computational Engineering Research (IJCER) Qualitative Analysis of a Discrete SIR Epidemic Model A. George Maria Selvam 1, D.
More informationTowards a theory of ecotone resilience: coastal vegetation on a salinity gradient
Towards a theory of ecotone resilience: coastal vegetation on a salinity gradient Jiang Jiang, Daozhou Gao and Donald L. DeAngelis Appendix A 1.1 The model without diffusion dn 1 = N 1 ρ 1 fs 11 N 1 1
More informationPersistence theory applied to Keen s model a link between mathematical biology and mathematical economics
Persistence theory applied to Keen s model a link between mathematical biology and mathematical economics Jianhong Wu and Xiang-Sheng Wang Mprime Centre for Disease Modelling York University, Toronto Persistence
More informationDelay SIR Model with Nonlinear Incident Rate and Varying Total Population
Delay SIR Model with Nonlinear Incident Rate Varying Total Population Rujira Ouncharoen, Salinthip Daengkongkho, Thongchai Dumrongpokaphan, Yongwimon Lenbury Abstract Recently, models describing the behavior
More informationSTABILITY AND PERSISTENCE IN A MODEL FOR BLUETONGUE DYNAMICS
STABILITY AND PERSISTENCE IN A MODEL FOR BLUETONGUE DYNAMICS STEPHEN A. GOURLEY, HORST R. THIEME, AND P. VAN DEN DRIESSCHE Abstract. A model for the time evolution of bluetongue, a viral disease in sheep
More informationGlobal Analysis of an SEIRS Model with Saturating Contact Rate 1
Applied Mathematical Sciences, Vol. 6, 2012, no. 80, 3991-4003 Global Analysis of an SEIRS Model with Saturating Contact Rate 1 Shulin Sun a, Cuihua Guo b, and Chengmin Li a a School of Mathematics and
More informationGlobal analysis of multi-strains SIS, SIR and MSIR epidemic models
Global analysis of multi-strains I, IR and MIR epidemic models Derdei Bichara, Abderrahman Iggidr, Gauthier allet To cite this version: Derdei Bichara, Abderrahman Iggidr, Gauthier allet. Global analysis
More informationSIS and SIR Epidemic Models Under Virtual Dispersal
Bull Math Biol DOI 10.1007/s11538-015-0113-5 ORIGINAL ARTICLE SIS and SIR Epidemic Models Under Virtual Dispersal Derdei Bichara 1 Yun Kang 2 Carlos Castillo-Chavez 1 Richard Horan 3 Charles Perrings 4
More informationGlobal attractors, stability, and population persistence
Global attractors, stability, and population persistence Horst R. Thieme (partly joint work with Hal L. Smith*) School of Mathematical and Statistical Sciences Arizona State University partially supported
More informationSpectral Properties of Matrix Polynomials in the Max Algebra
Spectral Properties of Matrix Polynomials in the Max Algebra Buket Benek Gursoy 1,1, Oliver Mason a,1, a Hamilton Institute, National University of Ireland, Maynooth Maynooth, Co Kildare, Ireland Abstract
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a ournal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationGlobal Dynamics of an SEIRS Epidemic Model with Constant Immigration and Immunity
Global Dynamics of an SIRS pidemic Model with Constant Immigration and Immunity Li juan Zhang Institute of disaster prevention Basic Course Department Sanhe, Hebei 065201 P. R. CHIA Lijuan262658@126.com
More informationThe effect of emigration and immigration on the dynamics of a discrete-generation population
J. Biosci., Vol. 20. Number 3, June 1995, pp 397 407. Printed in India. The effect of emigration and immigration on the dynamics of a discrete-generation population G D RUXTON Biomathematics and Statistics
More informationStability of a Numerical Discretisation Scheme for the SIS Epidemic Model with a Delay
Stability of a Numerical Discretisation Scheme for the SIS Epidemic Model with a Delay Ekkachai Kunnawuttipreechachan Abstract This paper deals with stability properties of the discrete numerical scheme
More informationPERSISTENCE AND PERMANENCE OF DELAY DIFFERENTIAL EQUATIONS IN BIOMATHEMATICS
PERSISTENCE AND PERMANENCE OF DELAY DIFFERENTIAL EQUATIONS IN BIOMATHEMATICS PhD Thesis by Nahed Abdelfattah Mohamady Abdelfattah Supervisors: Prof. István Győri Prof. Ferenc Hartung UNIVERSITY OF PANNONIA
More informationSIR Epidemic Model with total Population size
Advances in Applied Mathematical Biosciences. ISSN 2248-9983 Volume 7, Number 1 (2016), pp. 33-39 International Research Publication House http://www.irphouse.com SIR Epidemic Model with total Population
More informationGLOBAL DYNAMICS OF A TICK IXODES SCAPULARIS MODEL
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 19, Number 1, Spring 2011 GLOBAL DYNAMICS OF A TICK IXODES SCAPULARIS MODEL YIJUN LOU AND JIANHONG WU ABSTRACT. Lyme disease remains the world s most frequently
More informationAustralian Journal of Basic and Applied Sciences
AENSI Journals Australian Journal of Basic and Applied Sciences ISSN:1991-8178 Journal home page: www.ajbasweb.com A SIR Transmission Model of Political Figure Fever 1 Benny Yong and 2 Nor Azah Samat 1
More informationA Note on the Spread of Infectious Diseases. in a Large Susceptible Population
International Mathematical Forum, Vol. 7, 2012, no. 50, 2481-2492 A Note on the Spread of Infectious Diseases in a Large Susceptible Population B. Barnes Department of Mathematics Kwame Nkrumah University
More informationLOCAL STABILITY IMPLIES GLOBAL STABILITY IN SOME ONE-DIMENSIONAL DISCRETE SINGLE-SPECIES MODELS. Eduardo Liz. (Communicated by Linda Allen)
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 7, Number 1, January 2007 pp. 191 199 LOCAL STABILITY IMPLIES GLOBAL STABILITY IN SOME ONE-DIMENSIONAL DISCRETE
More informationPERMANENCE IN LOGISTIC AND LOTKA-VOLTERRA SYSTEMS WITH DISPERSAL AND TIME DELAYS
Electronic Journal of Differential Equations, Vol. 2005(2005), No. 60, pp. 1 11. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) PERMANENCE
More informationBifurcation Analysis in Simple SIS Epidemic Model Involving Immigrations with Treatment
Appl. Math. Inf. Sci. Lett. 3, No. 3, 97-10 015) 97 Applied Mathematics & Information Sciences Letters An International Journal http://dx.doi.org/10.1785/amisl/03030 Bifurcation Analysis in Simple SIS
More informationGLOBAL STABILITY FOR A CLASS OF DISCRETE SIR EPIDEMIC MODELS. Yoichi Enatsu and Yukihiko Nakata. Yoshiaki Muroya. (Communicated by Yasuhiro Takeuchi)
MATHEMATICAL BIOSCIENCES doi:10.3934/me.2010.7.347 AND ENGINEERING Volume 7, Numer 2, April 2010 pp. 347 361 GLOBAL STABILITY FOR A CLASS OF DISCRETE SIR EPIDEMIC MODELS Yoichi Enatsu and Yukihiko Nakata
More informationAARMS Homework Exercises
1 For the gamma distribution, AARMS Homework Exercises (a) Show that the mgf is M(t) = (1 βt) α for t < 1/β (b) Use the mgf to find the mean and variance of the gamma distribution 2 A well-known inequality
More informationA simple two-patch epidemiological model with Allee effects and disease-modified fitness
Contemporary Mathematics A simple two-patch epidemiological model with Allee effects and disease-modified fitness Yun Kang and Carlos Castillo-Chavez This paper is dedicated to Professor Ronald Mickens
More informationA NONLOCAL REACTION-DIFFUSION POPULATION MODEL WITH STAGE STRUCTURE
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 11, Number 3, Fall 23 A NONLOCAL REACTION-DIFFUSION POPULATION MODEL WITH STAGE STRUCTURE Dedicated to Professor Paul Waltman on the occasion of his retirement
More informationGlobal Qualitative Analysis for a Ratio-Dependent Predator Prey Model with Delay 1
Journal of Mathematical Analysis and Applications 266, 401 419 (2002 doi:10.1006/jmaa.2001.7751, available online at http://www.idealibrary.com on Global Qualitative Analysis for a Ratio-Dependent Predator
More informationLOTKA-VOLTERRA SYSTEMS WITH DELAY
870 1994 133-140 133 LOTKA-VOLTERRA SYSTEMS WITH DELAY Zhengyi LU and Yasuhiro TAKEUCHI Department of Applied Mathematics, Faculty of Engineering, Shizuoka University, Hamamatsu 432, JAPAN ABSTRACT Sufftcient
More informationResearch Article Modeling Computer Virus and Its Dynamics
Mathematical Problems in Engineering Volume 213, Article ID 842614, 5 pages http://dx.doi.org/1.1155/213/842614 Research Article Modeling Computer Virus and Its Dynamics Mei Peng, 1 Xing He, 2 Junjian
More informationDENSITY DEPENDENCE IN DISEASE INCIDENCE AND ITS IMPACTS ON TRANSMISSION DYNAMICS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 19, Number 3, Fall 2011 DENSITY DEPENDENCE IN DISEASE INCIDENCE AND ITS IMPACTS ON TRANSMISSION DYNAMICS REBECCA DE BOER AND MICHAEL Y. LI ABSTRACT. Incidence
More informationAnalysis of SIR Mathematical Model for Malaria disease with the inclusion of Infected Immigrants
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 14, Issue 5 Ver. I (Sep - Oct 218), PP 1-21 www.iosrjournals.org Analysis of SIR Mathematical Model for Malaria disease
More informationStochastic modelling of epidemic spread
Stochastic modelling of epidemic spread Julien Arino Department of Mathematics University of Manitoba Winnipeg Julien Arino@umanitoba.ca 19 May 2012 1 Introduction 2 Stochastic processes 3 The SIS model
More informationModelling of the Hand-Foot-Mouth-Disease with the Carrier Population
Modelling of the Hand-Foot-Mouth-Disease with the Carrier Population Ruzhang Zhao, Lijun Yang Department of Mathematical Science, Tsinghua University, China. Corresponding author. Email: lyang@math.tsinghua.edu.cn,
More informationAGE OF INFECTION IN EPIDEMIOLOGY MODELS
24 Conference on Diff. Eqns. and Appl. in Math. Biology, Nanaimo, BC, Canada. Electronic Journal of Differential Equations, Conference 12, 25, pp. 29 37. ISSN: 172-6691. URL: http://ejde.math.txstate.edu
More informationarxiv: v2 [q-bio.pe] 3 Oct 2018
Journal of Mathematical Biology manuscript No. (will be inserted by the editor Global stability properties of renewal epidemic models Michael T. Meehan Daniel G. Cocks Johannes Müller Emma S. McBryde arxiv:177.3489v2
More information(mathematical epidemiology)
1. 30 (mathematical epidemiology) 2. 1927 10) * Anderson and May 1), Diekmann and Heesterbeek 3) 7) 14) NO. 538, APRIL 2008 1 S(t), I(t), R(t) (susceptibles ) (infectives ) (recovered/removed = βs(t)i(t)
More informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationT.8. Perron-Frobenius theory of positive matrices From: H.R. Thieme, Mathematics in Population Biology, Princeton University Press, Princeton 2003
T.8. Perron-Frobenius theory of positive matrices From: H.R. Thieme, Mathematics in Population Biology, Princeton University Press, Princeton 2003 A vector x R n is called positive, symbolically x > 0,
More informationSI j RS E-Epidemic Model With Multiple Groups of Infection In Computer Network. 1 Introduction. Bimal Kumar Mishra 1, Aditya Kumar Singh 2
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.13(2012) No.3,pp.357-362 SI j RS E-Epidemic Model With Multiple Groups of Infection In Computer Network Bimal Kumar
More informationFISHER WAVES IN AN EPIDEMIC MODEL. Xiao-Qiang Zhao. Wendi Wang. (Communicated by Hal Smith)
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 4, Number 4, November 2004 pp. 1117 1128 FISHER WAVES IN AN EPIDEMIC MODEL Xiao-Qiang Zhao Department of Mathematics
More information1. Introduction: time-continuous linear population dynamics. Differential equation and Integral equation. Semigroup approach.
Intensive Programme Course: Lecturer: Dates and place: Mathematical Models in Life and Social Sciences Structured Population Dynamics in ecology and epidemiology Jordi Ripoll (Universitat de Girona, Spain)
More informationA Time Since Recovery Model with Varying Rates of Loss of Immunity
Bull Math Biol (212) 74:281 2819 DOI 1.17/s11538-12-978-7 ORIGINAL ARTICLE A Time Since Recovery Model with Varying Rates of Loss of Immunity Subhra Bhattacharya Frederick R. Adler Received: 7 May 212
More informationEpidemics in Two Competing Species
Epidemics in Two Competing Species Litao Han 1 School of Information, Renmin University of China, Beijing, 100872 P. R. China Andrea Pugliese 2 Department of Mathematics, University of Trento, Trento,
More informationGlobal Properties for Virus Dynamics Model with Beddington-DeAngelis Functional Response
Global Properties for Virus Dynamics Model with Beddington-DeAngelis Functional Response Gang Huang 1,2, Wanbiao Ma 2, Yasuhiro Takeuchi 1 1,Graduate School of Science and Technology, Shizuoka University,
More informationMathematical and Numerical Comparisons of Five Single-Population Growth Models
Atlanta University Center DigitalCommons@Robert W. Woodruff Library, Atlanta University Center Clark Atlanta University Faculty Publications Clark Atlanta University 2016 Mathematical and Numerical Comparisons
More informationIn particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with
Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient
More informationDynamical models of HIV-AIDS e ect on population growth
Dynamical models of HV-ADS e ect on population growth David Gurarie May 11, 2005 Abstract We review some known dynamical models of epidemics, given by coupled systems of di erential equations, and propose
More informationA BINOMIAL MOMENT APPROXIMATION SCHEME FOR EPIDEMIC SPREADING IN NETWORKS
U.P.B. Sci. Bull., Series A, Vol. 76, Iss. 2, 2014 ISSN 1223-7027 A BINOMIAL MOMENT APPROXIMATION SCHEME FOR EPIDEMIC SPREADING IN NETWORKS Yilun SHANG 1 Epidemiological network models study the spread
More informationMathematical Epidemiology Lecture 1. Matylda Jabłońska-Sabuka
Lecture 1 Lappeenranta University of Technology Wrocław, Fall 2013 What is? Basic terminology Epidemiology is the subject that studies the spread of diseases in populations, and primarily the human populations.
More informationStability Analysis of an SVIR Epidemic Model with Non-linear Saturated Incidence Rate
Applied Mathematical Sciences, Vol. 9, 215, no. 23, 1145-1158 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.41164 Stability Analysis of an SVIR Epidemic Model with Non-linear Saturated
More informationA NEW SOLUTION OF SIR MODEL BY USING THE DIFFERENTIAL FRACTIONAL TRANSFORMATION METHOD
April, 4. Vol. 4, No. - 4 EAAS & ARF. All rights reserved ISSN35-869 A NEW SOLUTION OF SIR MODEL BY USING THE DIFFERENTIAL FRACTIONAL TRANSFORMATION METHOD Ahmed A. M. Hassan, S. H. Hoda Ibrahim, Amr M.
More informationThe Fractional-order SIR and SIRS Epidemic Models with Variable Population Size
Math. Sci. Lett. 2, No. 3, 195-200 (2013) 195 Mathematical Sciences Letters An International Journal http://dx.doi.org/10.12785/msl/020308 The Fractional-order SIR and SIRS Epidemic Models with Variable
More informationStability and bifurcation analysis of epidemic models with saturated incidence rates: an application to a nonmonotone incidence rate
Published in Mathematical Biosciences and Engineering 4 785-85 DOI:.3934/mbe.4..785 Stability and bifurcation analysis of epidemic models with saturated incidence rates: an application to a nonmonotone
More informationA periodic epidemic model in a patchy environment
J Math Anal Appl 325 (2007) 496 516 wwwelseviercom/locate/jmaa A periodic epidemic model in a patchy environment Fang Zhang, Xiao-Qiang Zhao Department of Mathematics and Statistics, Memorial University
More informationc 2006 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 66, No. 5, pp. 1563 1587 c 2006 Society for Industrial and Applied Mathematics DISCRETE-TIME SIS EPIDEMIC MODEL IN A SEASONAL ENVIRONMENT JOHN E. FRANKE AND ABDUL-AZIZ YAKUBU Abstract.
More informationASYMPTOTIC SPEED OF SPREAD AND TRAVELING WAVES FOR A NONLOCAL EPIDEMIC MODEL. Dashun Xu and Xiao-Qiang Zhao. (Communicated by Hal Smith)
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 5, Number 4, November 25 pp. 143 156 ASYMPTOTIC SPEED OF SPREAD AND TRAVELING WAVES FOR A NONLOCAL EPIDEMIC MODEL
More informationAnalysis and Monte Carlo simulations of a model for the spread of infectious diseases in heterogeneous metapopulations
PHYSICAL REVIEW E 8, 492 29 Analysis and Monte Carlo simulations of a model for the spread of infectious diseases in heterogeneous metapopulations David Juher,* Jordi Ripoll, and Joan Saldaña Departament
More informationSOME SIMPLE EPIDEMIC MODELS. Fred Brauer. In honor of Professor Zhien Ma s 70th birthday
MATHEMATICAL BIOSCIENCES http://www.mbejournal.org/ AND ENGINEERING Volume 3, Number 1, January 26 pp. 1 15 SOME SIMPLE EPIDEMIC MODELS Fred Brauer Department of Mathematics University of British Columbia
More information