Essential Ideas of Mathematical Modeling in Population Dynamics
|
|
- Claribel Franklin
- 5 years ago
- Views:
Transcription
1 Essential Ideas of Mathematical Modeling in Population Dynamics Toward the application for the disease transmission dynamics Hiromi SENO Research Center for Pure and Applied Mathematics Department of Computer and Mathematical Sciences Graduate School of Information Sciences Tohoku University, Japan FOR: Department of Mathematics, Universitas Indonesia, August 23 26, 2017
2 Outline of This Course Prologue Introduction to the mathematical modeling of population dynamics Modeling for malthusian and logistic growths Law of mass-action Extension of logistic equation model Time-discrete model Mathematical modeling of disease transmission dynamics Kermack McKendrick model Modeling of a discrete model for the disease transmission dynamics Mathematical nature of the discrete model Extension with the introduction of time step size Limit to zero time step size Relation to Kermack McKendrick model Basic reproduction number R 0 Key factors and idea for the more sophisticated modeling Epilogue
3 Prologue
4 Hiromi SENO Research Center for Pure and Applied Mathematics, Department of Computer and Mathematical Sciences Graduate School of Information Sciences Tohoku University mathematical expressions mathematical translations BIOLOGICAL PROBLEMS biological hypotheses & assumptions design of mathematical analyses mathematical analyses modeling hypotheses & assumptions biological researches mathematical results mathematical discussions MATHEMATICAL PROBLEMS biological translations biological discussions
5 Mathematical model for biological phenomenon Mathematical Model for Application Biology Advanced Mathematical Model Basic Mathematical Model Mathematics
6 Mathematical model for biological phenomenon Mathematical model for the systemization of knowledges for biological phenomena Mathematical model for the quantitative understanding of biological phenomena Mathematical model for the qualitative understanding of biological phenomena Mathematical model for the research on a specific biological phenomena
7 Mathematical model for biological phenomenon explanation experiment description systemization prediction model development understanding mathematical interest
8 Mathematical modeling
9 Mathematical modeling of population dynamics is the nature of the spatio-temporal variation of biological population size (i.e. density etc.).
10 Mathematical modeling of population dynamics is the nature of the spatio-temporal variation of biological population size (i.e. density etc.).
11 Mathematical modeling of population dynamics
12 Mathematical modeling of population dynamics
13 Hiromi SENO Research Center for Pure and Applied Mathematics, Department of Computer and Mathematical Sciences Graduate School of Information Sciences Tohoku University What mathematical model is reasonable from the biological viewpoint?
14 Hiromi SENO Research Center for Pure and Applied Mathematics, Department of Computer and Mathematical Sciences Graduate School of Information Sciences Tohoku University What mathematical structure is appropriate for the reasonable modeling?
15 Hiromi SENO Research Center for Pure and Applied Mathematics, Department of Computer and Mathematical Sciences Graduate School of Information Sciences Tohoku University Reasonability of modeling depends on i) purpose of theoretical research; ii) available data/knowledge/hypothesis; iii) design of mathematical analysis.
16 Introduction to the mathematical modeling of population dynamics
17 Modeling of malthusian growth
18 Modeling of malthusian growth dn(t) = rn(t) N(t) = N(0)e rt Population size at time t Intrinsic growth rate
19 Malthusian growth N t
20 Malthusian growth
21 Modeling of malthusian growth dn(t) = rn(t) N(t) = N(0)e rt
22 Modeling of malthusian growth dn(t) = rn(t) t 0 lim N(t + t) N(t) t = rn(t)
23 Modeling of malthusian growth dn(t) = rn(t) N(t + t) N(t) t rn(t)
24 Modeling of malthusian growth dn(t) = rn(t) N(t + t) N(t) N(t) t r Time-independent and density-independent constant
25 Modeling of malthusian growth { N(t + t) N(t) } N(t) t r Mean per capita increment of population size in t
26 Modeling of malthusian growth { N(t + t) N(t) } N(t) t r Mean per capita increment velocity of population size in t i.e. Mean per capita growth rate Physiological condition of individual
27 Modeling of malthusian growth dn(t) = rn(t) N(t + t) N(t) N(t) t r Time-independent and density-independent constant Mean per capita growth rate Physiological condition of individual
28 Modeling of logistic growth
29 Logistic equations [L-1] [L-2] [L-3] [L-4] dn(t) dn(t) dn(t) dn(t) = { r 0 βn(t) } N(t) { = r 0 1 N(t) } N(t) K = r 0 N(t) b{n(t)} 2 = { r 0 βn(t) } N(t) b{n(t)} 2 r 0 : intrinsic growth rate
30 Logistic equation [L-1] dn(t) = { r 0 βn(t) } N(t)
31 Logistic growth N(t) = r 0 /β 1 + { r 0 /β N(0) 1} e r 0t
32 Logistic growth N t
33 Logistic growth
34 Logistic growth
35 Logistic growth Drosophila
36 Logistic growth? Tasmanian sheep
37 Logistic equation [L-1] dn(t) = { r 0 βn(t) } N(t)
38 Modeling of logistic growth dn(t) = { r 0 N(t) }N(t) N(t + t) N(t) lim t 0 t = { r 0 N(t) }N(t)
39 Modeling of logistic growth dn(t) = { r 0 N(t) }N(t) N(t + t) N(t) t { N(t) }N(t) r 0
40 Modeling of logistic growth dn(t) = { r 0 N(t) }N(t) N(t + t) N(t) N(t) t r 0 N(t) Density-dependent Mean per capita growth rate Physiological condition of individual
41 Modeling of logistic growth N(t + t) N(t) N(t) t r 0 N(t) Density Effect
42 Modeling of logistic growth Per capita growth rate r 0 r 0 N
43 Logistic growth (!?) Per capita growth rate azuki bean beetle Callosobruchus chinensis (Linnaeus) Bruchidae Population density
44 Logistic growth (!?) Population Density azuki bean beetle Callosobruchus chinensis (Linnaeus) Bruchidae Generation
45 Logistic growth (!??) The great tit Parus major in Holland
46 Number of independent offsprings per reproductive individual Logistic growth (!?) Number of reproductive individuals
47 Logistic equations L-1 and L-2 [L-1] [L-2] dn(t) dn(t) = { r 0 βn(t) } N(t) { = r 0 1 N(t) } N(t) K
48 Logistic growth N carrying capacity t
49 Logistic equations L-1 and L-2 Carrying capacity [L-1] [L-2] dn(t) dn(t) = { r 0 βn(t) } N(t) { = r 0 1 N(t) K } N(t) [L-1] r 0 β [L-2] K Mathematical equivalence between L-1 and L-2 is not that of modeling!
50 Logistic equations L-1 and L-2 Carrying capacity [L-1] [L-2] dn(t) dn(t) = { r 0 βn(t) } N(t) { = r 0 1 N(t) K } N(t) [L-1] r 0 β [L-2] K They are significantly different from each other in the dependence of the carrying capacity on the intrinsic growth rate!
51 Logistic equation L-2 dn(t) { = r 0 1 N(t) K } N(t)
52 Logistic equations [L-1] [L-2] [L-3] [L-4] dn(t) dn(t) dn(t) dn(t) = { r 0 βn(t) } N(t) { = r 0 1 N(t) } N(t) K = r 0 N(t) b{n(t)} 2 = { r 0 βn(t) } N(t) b{n(t)} 2 r 0 : intrinsic growth rate
53 Law of mass-action
54 Law of mass-action αa + βb γc
55 Law of mass-action αa + βb γc Reaction velocity V = 1 α d[a] = 1 β d[b] = 1 γ d[c]
56 Law of mass-action in chemical reaction αa + βb γc Reaction velocity: law of kinetic mass-action V = V([A], [B], [C]) = κ [A] n A [B] n B [C] n C
57 Law of mass-action Lotka Volterra type of mass-action assumption V = V([A], [B]) = κ [A][B]
58 Logistic equation L-3 [L-3] Malthusian growth Decrease due to intra-specific reaction
59 Logistic equation L-4 [L-4] = [L-1] + [L-3] Logistic growth Decrease due to intra-specific reaction
60 Logistic equations [L-1] [L-2] [L-3] [L-4] dn(t) dn(t) dn(t) dn(t) = { r 0 βn(t) } N(t) { = r 0 1 N(t) } N(t) K = r 0 N(t) b{n(t)} 2 = { r 0 βn(t) } N(t) b{n(t)} 2 r 0 : intrinsic growth rate
61 Logistic equation with the mass-action assumption for resource consumption dn(t) dr(t) = cρr(t)n(t) = ρr(t)n(t) cf. bacteria population in culture
62 Logistic equation with the mass-action assumption for resource consumption dr/ dn/ = dr dn = ρrn cρrn = 1 c (= const.)
63 Logistic equation with the mass-action assumption for resource consumption R(t) = 1 c N(t)+C = 1 c N(t)+R(0)+ 1 c N(0)
64 Logistic equation with the mass-action assumption for resource consumption dn(t) dr(t) = cρr(t)n(t) = ρr(t)n(t) cf. bacteria population in culture
65 Logistic equation with the mass-action assumption for resource consumption dn(t) = ρr(t)n(t) = ρ { cr(0)+n(0) N(t) } N(t) Carrying capacity: cr(0)+n(0)
66 Logistic growth
67 Logistic growth
68 Logistic growth Drosophila
69 Extension of logistic equation model
70 Extended logistic equation with the generalized density effect function dn(t) = [ ] r 0 β{n(t)} α N(t) }{{} Density-dependent per capita growth rate Gilpin & Ayala (1973)
71 Extended logistic equation with the generalized density effect function dn(t) = [ r 0 β{n(t)} α] N(t)
72 Logistic growth (!?) Per capita growth rate azuki bean beetle Callosobruchus chinensis (Linnaeus) Bruchidae Population density
73 Number of independent offsprings per reproductive individual Logistic growth (!?) Number of reproductive individuals
74 Extended logistic equation with the generalized mass-action assumption dn(t) = r 0 N(t) γ{n(t)} θ }{{} Decrease due to the intraspecific reaction
75 Extended logistic equation with the generalized mass-action assumption dn(t) = r 0 N(t) γ{n(t)} θ (θ < 1)
76 Extended logistic equation dn = (r 0 βn α ) N γn }{{}}{{} θ Per capita growth rate Intraspecific reaction
77 Extended logistic equation dn(t) =(r 0 βn α )N γn θ (θ < 1) Allee effect
78 Time-discrete model of malthusian growth
79 Time-discrete model of malthusian growth dn(t) = r 0 N(t) N(t) =N(0) e rt
80 Time-discrete model of malthusian growth dn(t) = r 0 N(t) N(t) =N(0) e rt N(t) = N(0) e rt N(t + h) = N(0) e r(t+h) = e rh N(0) e rt h time step size for the discrete dynamics
81 Time-discrete model of malthusian growth dn(t) = r 0 N(t) N(t) =N(0) e rt N(t + h) =e rh N(t) h time step size for the discrete dynamics
82 Time-discrete model of malthusian growth dn(t) = r 0 N(t) N(t) =N(0) e rt N k+1 = RN k with time step size h and R = e rh.
83 Malthusian growth N t
84 Time-discrete model of logistic growth
85 Time-discrete model of logistic growth dn(t) = { r 0 βn(t) } N(t) N(t) = r 0 /β 1 + { r 0 /β N(0) 1} e r 0t
86 Time-discrete model of logistic growth dn(t) = { r 0 βn(t) } N(t) N(t) = r 0 /β 1 + { r 0 /β N(0) 1} e r 0t N(t + h) = e rh N(t) 1 + er 0 h 1 r 0 βn(t) h time step size for the discrete dynamics
87 Time-discrete model of logistic growth dn(t) = { r 0 βn(t) } N(t) N(t) = r 0 /β 1 + { r 0 /β N(0) 1} e r 0t Verhulst Beverton Holt model N k+1 = RN k 1 + bn k with time step size h, R = e rh and b = er 0 h 1 r 0 β.
88 Time-discrete model of logistic growth 4 3 N time
89 Time-discrete model of logistic growth Discretization with the Euler scheme dn(t) = { r 0 βn(t) } N(t) Ñ(t + h) Ñ(t) h = { r 0 βñ(t) } Ñ(t) h time step size for the discretization
90 Time-discrete model of logistic growth Discretization with the Euler scheme dn(t) = { r 0 βn(t) } N(t) Ñ(t + h) = { 1 + r 0 h βhñ(t) } Ñ(t) h time step size for the discretization
91 Time-discrete model of logistic growth logistic map N k+1 = r(1 N k ) N k with time step size h, r = 1 + r 0 h and N k = βhñ(kh) 1+r 0 h.
92 Time-discrete model of logistic growth logistic map ^ N k ^ N k
93 Time-discrete model of logistic growth Bifurcation diagram N^ ^ r
94 Time-discrete model of logistic growth 4 Simple Euler Scheme 3 N time
95 Time-discrete model of logistic growth The other discretization with the Euler scheme dn(t) = { r 0 βn(t) } N(t) d log N(t) = r 0 βn(t) log Ñ(t + h) log Ñ(t) h = r 0 βñ(t) h time step size for the discretization
96 Time-discrete model of logistic growth The other discretization with the Euler scheme dn(t) = { r 0 βn(t) } N(t) Ñ(t + h) =e r 0h βhñ(t) Ñ(t) h time step size for the discretization
97 Time-discrete model of logistic growth Ricker model N k+1 = R e N k N k with time step size h, R = e r 0h and N k = βhñ(kh).
98 Logistic growth (!?) Population Density azuki bean beetle Callosobruchus chinensis (Linnaeus) Bruchidae Generation
99 Mathematical modeling of disease transmission dynamics
100 Kermack McKendrick model
101 Kermack McKendrick model Susceptible S β I (1 m)q Infective I R Removed/Recovered mq
102 Kermack McKendrick model ds(t) di(t) dr(t) = βi(t)s(t)+(1 m)qi(t)+θr(t) = βi(t)s(t) qi(t) = mqi(t) θr(t) S(t) : Susceptible population size at time t; I(t) : Infective population size at time t; R(t) : Recovered immune population size at time t; N : Total population size (time-independent constant; = S(t)+I(t)+R(t))
103 Kermack McKendrick model ds(t) di(t) dr(t) = β I(t) N = β I(t) N cn S(t)+(1 m)qi(t)+θr(t) cn S(t) qi(t) = mqi(t) θr(t) ( β = βc) S(t) : Susceptible population size at time t; I(t) : Infective population size at time t; R(t) : Recovered immune population size at time t; N : Total population size (time-independent constant; = S(t)+I(t)+R(t))
104 Kermack McKendrick model Susceptible S β I N cn (1 m)q Infective I R Removed/Recovered mq
105 Kermack McKendrick SIR model (θ = 0 and m = 1) Susceptible S β I N cn Infective I R Removed/Recovered q
106 Kermack McKendrick SIR model (θ = 0 and m = 1) ds(t) di(t) dr(t) = βi(t)s(t)+(1 m)qi(t)+θr(t) = βi(t)s(t) qi(t) = mqi(t) θr(t)
107 Kermack McKendrick SIR model (θ = 0 and m = 1) S R I I time q β S
108 Kermack McKendrick SIR model (θ = 0 and m = 1) ds(t) = βi(t)s(t) di(t) = βi(t)s(t) qi(t)
109 Kermack McKendrick SIR model (θ = 0 and m = 1) 1 S(t) ds(t) = βi(t) ds(t) + di(t) = qi(t)
110 Kermack McKendrick SIR model (θ = 0 and m = 1) ds(t) + di(t) q β 1 S(t) ds(t) = 0
111 Kermack McKendrick SIR model (θ = 0 and m = 1) S(t)+I(t) q β log S(t) =const.
112 Kermack McKendrick SIR model (θ = 0 and m = 1) Conserved quantity of Kermack McKendrick SIR model S(t)+I(t) q β log S(t) =S(0)+I(0) q β log S(0)
113 Kermack McKendrick SIR model (θ = 0 and m = 1) Conserved quantity of Kermack McKendrick SIR model S(t)+I(t) q β log S(t) =S(0)+I(0) q β log S(0) I 0 q β S
114 Modeling of a discrete model for the disease transmission dynamics
115 Assumptions A population with an infectious disease of negligible fatality; Disease transmission dynamics in such a time scale that the temporal variation of total population size with birth, death and migration could be negligible; Approximation of complete mixing of individuals in terms of disease transmission process.
116 Assumptions Probability that the cumulative number of contacts of an individual with others per day, P(i); Expected cumulative number of contacts of an individual with others per day, π = i=0 ip(i); Probability that a susceptible individual is infected per contact, β; Probability that a susceptible individual with j times contacts to infectives can successfully escape the disease transmission, (1 β) j.
117 Modeling: probability about the number of contacted infectives Probability that l of j contacts are with infective at the k th day ( j l )( ) l ( Ik 1 I ) j l k N N ( ) j := l l! l!(j l)! I k N : probability that an encounter is with an infective
118 Modeling: probability about successful escape from disease transmission Probability that a susceptible can escape from the infection when he/she contacts j times with others at the k th day j l=0(1 β) l ( j l = j l=0 ( j l )( ) l ( Ik 1 I ) j l k N N ){ } Ik /N l ( (1 β) 1 I k /N } j ( { Ik /N = (1 β)+1 1 I k /N ( = 1 β I ) j k N 1 I k N 1 I k N ) j ) j
119 Discrete model for the disease transmission dynamics S k+1 = I k+1 = ( 1 β I ) j k P(j)S k +(1 m)qi k + θr k N j=0 { ( 1 1 β I ) } j k P(j)S k +(1 q)i k N j=0 R k+1 = mqi k +(1 θ)r k q : Probability of recovery with loss of infectivity; m : Probability of successful establishment of immunity; θ : Probability of loss of immunity.
120 Discrete model for the disease transmission dynamics
121 Discrete model for the disease transmission dynamics In case of Poisson distribution of contact numbers with others per day P(j) = γj e γ j! ( π = γ ) ( 1 β I ) j k P(j) =e βγi k/n j=0 N
122 Discrete model for the disease transmission dynamics S k+1 = S k e βγi k/n +(1 m)qi k + θr k I k+1 = S k (1 e βγi k/n )+(1 q)i k R k+1 = mqi k +(1 θ)r k
123 Discrete model for the disease transmission dynamics Susceptible S 1 e (1 m)q βγ I k N Infective I R mq Removed/Recovered θ=0 and m=1 SIR model
124 Discrete SIR model for the disease transmission dynamics θ=0 and m=1 Susceptible S 1 e βγ I k N Infective I R Removed/Recovered q
125 Discrete SIR model for the disease transmission dynamics θ=0 and m=1 S k+1 = S k e βγi k/n +(1 m)qi k + θr k I k+1 = S k (1 e βγi k/n )+(1 q)i k R k+1 = mqi k + (1 θ)r k
126 Discrete SIR model for the disease transmission dynamics I k S k Trajectories of our time-discrete SIR model. S(0) =0.99; I(0) =0.01; R(0) =0.0; q = 0.2; β = 0.01; γ = 50.0.
127 Mathematical nature of the discrete model
128 Discrete SIR model for the disease transmission dynamics S k+1 = S k e βγi k/n I k+1 = S k (1 e βγi k/n )+(1 q)i k
129 Discrete SIR model for the disease transmission dynamics log S k+1 S k = βγ I k N S k+1 + I k+1 = S k + I k qi k
130 Discrete SIR model for the disease transmission dynamics S k+1 + I k+1 qn βγ log S k+1 = S k + I k qn βγ log S k
131 Discrete SIR model for the disease transmission dynamics Conserved quantity of discrete SIR model S k + I k qn βγ log S k = S 0 + I 0 qn βγ log S 0
132 Extension of the discrete model with the introduction of time step size
133 Introduction of time step size h Discrete model for the disease transmission dynamics S k+1 = S k e βγi k/n +(1 m)qi k + θr k I k+1 = S k (1 e βγi k/n )+(1 q)i k R k+1 = mqi k +(1 θ)r k (S k, I k, R k ) ( S(t), I(t), R(t) ) ; (S k+1, I k+1, R k+1 ) ( S(t + h), I(t + h), R(t + h) ) ; P(j) = γj e γ j! (γh)j e γh ; γ γh; j! q qh (0 qh 1); θ θh (0 θh < 1).
134 Introduction of time step size h Discrete model wth the time step size h S(t + h) =S(t)e βγh I(t)/N +(1 m)qh I(t)+θh R(t) I(t + h) =S(t){1 e βγh I(t)/N } +(1 qh)i(t) R(t + h) =mqh I(t)+(1 θh)r(t) (S k, I k, R k ) ( S(t), I(t), R(t) ) ; (S k+1, I k+1, R k+1 ) ( S(t + h), I(t + h), R(t + h) ) ; P(j) = γj e γ j! (γh)j e γh ; γ γh; j! q qh (0 qh 1); θ θh (0 θh < 1).
135 Limit to zero time step size
136 Limit to zero time step size Discrete model with the time step size h S(t + h) =S(t)e βγh I(t)/N +(1 m)qh I(t)+θh R(t) I(t + h) =S(t){1 e βγh I(t)/N } +(1 qh)i(t) R(t + h) =mqh I(t)+(1 θh)r(t)
137 Limit to zero time step size Discrete model with the time step size h S(t + h) S(t) h = S(t) 1 e βγh I(t)/N h +(1 m)qi(t)+θr(t) I(t + h) I(t) h = S(t) 1 e βγh I(t)/N h qi(t) R(t + h) R(t) h = mqi(t) θr(t) h 0
138 Limit to zero time step size h 0 ODE model at the limit of h 0 ds(t) di(t) dr(t) = β I(t) γs(t)+(1 m)qi(t)+θr(t) N = β I(t) γs(t) qi(t) N = mqi(t) θr(t)
139 Relation to Kermack McKendrick model
140 Relation to Kermack McKendrick model ODE model at the limit of h 0 ds(t) di(t) dr(t) = β I(t) γs(t)+(1 m)qi(t)+θr(t) N = β I(t) γs(t) qi(t) N = mqi(t) θr(t)
141 Relation to Kermack McKendrick model Kermack McKendrick model ds(t) di(t) dr(t) = β I(t) N = β I(t) N cn S(t)+(1 m)qi(t)+θr(t) cn S(t) qi(t) = mqi(t) θr(t)
142 Relation to Kermack McKendrick model ODE model at the limit of h 0 ds(t) di(t) = β I(t) γs(t)+(1 m)qi(t)+θr(t) N = β I(t) γs(t) qi(t) N Kermack McKendrick model ds(t) di(t) = β I(t) N = β I(t) N cn S(t)+(1 m)qi(t)+θr(t) cn S(t) qi(t)
143 Relation to Kermack McKendrick model γ = π = cn
144 Relation to Kermack McKendrick model Correspondence of modeling to Kermack McKendrick model With Poisson distribution {P(j)} and γ = π = cn, our time-discrete SIR model coincides with Kermack McKendrick model at the limit of zero time step size. This result was expected, because of the coincidence in their modelings. Further for the SIR models,...
145 Dynamical consistency in SIR model Dynamical consistency in SIR model With Poisson distribution {P(j)} and γ = π = cn, our time-discrete SIR model has the quantitative dynamical consistency with Kermack McKendrick SIR model, independently of time step size h.
146 Dynamical consistency in SIR model Kermack McKendrick SIR model ( β = cβ) S(t)+I(t) q q log S(t) =S(0)+I(0) log S(0) cβ cβ Discrete SIR model (γ = π = cn) S k + I k q cβ log S k = S 0 + I 0 q cβ log S 0 Remark: This is independent of time step size h.
147 SIR models Relation to Kermack McKendrick model I k S k Trajectories of our time-discrete SIR model. S(0) =0.99; I(0) =0.01; R(0) =0.0; N = 1; q = 0.2; β = 0.01; γ = cn = 50.0.
148 Relation to Kermack McKendrick model Trajectories of our time-discrete model and Kermack McKendrick model. (a) m = 0.5; (b)m = 1.0. S(0) =0.999; I(0) =0.001; R(0) =0.0; q = 0.02; β = 0.1; γ = cn = 0.5; θ = 0.0; N = 1.0; h = In(b),thetrajectoryof time-discrete model is always on the solution curve of Kermack McKendrick model.
149 Basic Basic reproduction reproduction number number R R00
150 Basic reproductive number R 0 Definition of the basic reproductive number R 0 The expected number of new cases of an infection caused by an infected individual, in a population consisting of susceptible contacts only. Since the frequency of susceptibles decreases as the epidemic process goes on, the basic reproductive number defines, rigorously saying, the supremum of the expected number of new cases of an infection caused by an infected individual at the initial stage of epidemic process.
151 Basic reproductive number R 0 Definition of the basic reproductive number R 0 The expected number of new cases of an infection caused by an infected individual, in a population consisting of susceptible contacts only. R 0 > 1 R 0 < 1 Initial increase of infective frequency Monotonic disappearance of infection
152 Basic reproductive number R 0 Definition of the basic reproductive number R 0 The expected number of new cases of an infection caused by an infected individual, in a population consisting of susceptible contacts only. R 0 > 1 R 0 < 1 Invasion success of epidemic disease Invasion failure of epidemic disease
153 Discrete model for the disease transmission dynamics General formula of discrete model ( S k+1 = 1 β I ) j k P(j)S N k +(1 m)qi k + θr k j=0 { ( I k+1 = 1 1 β I ) } j k P(j)S N k +(1 q)i k j=0 R k+1 = mqi k +(1 θ)r k
154 Basic reproductive number R 0 For the general discrete model with I 0 N, R 0 = 0, ands 0 N I 1 = = { ( 1 1 β I ) } j 0 P(j)S 0 +(1 q)i 0 N j=0 j=0 { β jβi 0 P(j)+(1 q)i 0 } jp(j)+(1 q) I 0 j=0 = {β π +(1 q)}i 0
155 Basic reproductive number R 0 At the invasion stage in the discrete model Invasion success of infectious disease: I 1 > I 0 I 1 I 0 β π +(1 q) > 1
156 Basic reproductive number R 0 At the invasion stage in the discrete model Invasion success of infectious disease: I 1 > I 0 R 0 = β π q > 1
157 Kermack McKendrick model ds(t) di(t) dr(t) = βi(t)s(t)+(1 m)qi(t)+θr(t) = βi(t)s(t) qi(t) = mqi(t) θr(t) S(t) : Susceptible population size at time t; I(t) : Infective population size at time t; R(t) : Recovered immune population size at time t; N : Total population size (time-independent constant; = S(t)+I(t)+R(t))
158 Basic reproductive number R 0 Kermack McKendrick model with I(0) N, R(0) =0, ands(0) N di(t) = β I(0) cns(0) qi(0) t=0 N β I(0) cn N qi(0) N =(βcn q)i(0)
159 Basic reproductive number R 0 At the invasion stage in Kermack McKendrick model Invasion success of infectious disease: di(t) > 0 t=0 di(t) (βcn q)i(0) > 0 t=0
160 Basic reproductive number R 0 At the invasion stage in Kermack McKendrick model Invasion success of infectious disease: di(t) > 0 t=0 R 0 = βcn q > 1
161 Basic reproductive number R 0 For Kermack McKendrick model R 0 = βcn q For the discrete model R 0 = β π q
162 Consistency of R 0 The basic reproduction number R 0 is coincident between the Kermack McKendrick model and the present discrete model, with the correspondence π = cn. This result is reasonable from the coincidence in their modelings!
163 Discrete model and Kermack McKendrick model Trajectories of our time-discrete model and Kermack McKendrick model. Numerical calculation for the case that the system approaches an endemic state. S(0) =0.999; I(0) =0.001; R(0) =0.0; q = 0.02; m = 0.5; θ = 0.001; β = 0.1; γ = cn = 0.5; N = 1.0; h = 10.0.
164 Key factors and idea for the more sophisticated modeling
165 Key factors and idea for the more sophisticated modeling Spatio-temporal scales mean-field approximation, quasi-stationary state approximation, pair approximation, singular perturbation, metapopulation, renormalization Spatio-temporal heterogeneity reaction diffusion equations, lattice/cellular space, network space, non-standard mean-field approximation Human activities and community structure optimal theory, control theory, game theory, mathematical population genetics
166 Epilogue
167 Epilogue Discrete model has flexibility wider than the ODE model does with respect to the reasonable modeling, although it is in general less tractable for the mathematical analysis. Despite such potential difficulty in mathematics, it is primarily important for the theoretical research in mathematical biology that the structure of mathematical model must be reasonably determined/chosen/designed depending on the biological focus of the theoretical research with it. Such reasonable modeling necessarily requires both the biological/medical knowledge about the focused phenomenon and the mathematical knowledge/sense about the connection between the natures of phenomenon and the structures of mathematical factor applied for the modeling.
(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 informationpopulation size at time t, then in continuous time this assumption translates into the equation for exponential growth dn dt = rn N(0)
Appendix S1: Classic models of population dynamics in ecology and fisheries science Populations do not grow indefinitely. No concept is more fundamental to ecology and evolution. Malthus hypothesized that
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 informationMODELS ONE ORDINARY DIFFERENTIAL EQUATION:
MODELS ONE ORDINARY DIFFERENTIAL EQUATION: opulation Dynamics (e.g. Malthusian, Verhulstian, Gompertz, Logistic with Harvesting) Harmonic Oscillator (e.g. pendulum) A modified normal distribution curve
More informationFinal Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger. Project I: Predator-Prey Equations
Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger Project I: Predator-Prey Equations The Lotka-Volterra Predator-Prey Model is given by: du dv = αu βuv = ρβuv
More informationAnalysis of Numerical and Exact solutions of certain SIR and SIS Epidemic models
Journal of Mathematical Modelling and Application 2011, Vol. 1, No. 4, 51-56 ISSN: 2178-2423 Analysis of Numerical and Exact solutions of certain SIR and SIS Epidemic models S O Maliki Department of Industrial
More informationEpidemics in Networks Part 2 Compartmental Disease Models
Epidemics in Networks Part 2 Compartmental Disease Models Joel C. Miller & Tom Hladish 18 20 July 2018 1 / 35 Introduction to Compartmental Models Dynamics R 0 Epidemic Probability Epidemic size Review
More informationMA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, Exam Scores. Question Score Total
MA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, 2016 Exam Scores Question Score Total 1 10 Name: Section: Last 4 digits of student ID #: No books or notes may be used. Turn off all
More informationMA 777: Topics in Mathematical Biology
MA 777: Topics in Mathematical Biology David Murrugarra Department of Mathematics, University of Kentucky http://www.math.uky.edu/~dmu228/ma777/ Spring 2018 David Murrugarra (University of Kentucky) Lecture
More informationWe have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v),
Let us apply the approach presented above to the analysis of population dynamics models. 9. Lotka-Volterra predator-prey model: phase plane analysis. Earlier we introduced the system of equations for prey
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 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 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 informationMathematical Models of Biological Systems
Mathematical Models of Biological Systems CH924 November 25, 2012 Course outline Syllabus: Week 6: Stability of first-order autonomous differential equation, genetic switch/clock, introduction to enzyme
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 informationIntroduction to Epidemic Modeling
Chapter 2 Introduction to Epidemic Modeling 2.1 Kermack McKendrick SIR Epidemic Model Introduction to epidemic modeling is usually made through one of the first epidemic models proposed by Kermack and
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 informationMathematical and Computational Methods for the Life Sciences
Mathematical and Computational Methods for the Life Sciences Preliminary Lecture Notes Adolfo J. Rumbos November 24, 2007 2 Contents 1 Preface 5 2 Introduction: An Example from Microbial Genetics 7 I Deterministic
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 informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
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 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 informationThree Disguises of 1 x = e λx
Three Disguises of 1 x = e λx Chathuri Karunarathna Mudiyanselage Rabi K.C. Winfried Just Department of Mathematics, Ohio University Mathematical Biology and Dynamical Systems Seminar Ohio University November
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 informationIstván Faragó, Róbert Horváth. Numerical Analysis of Evolution Equations, Innsbruck, Austria October, 2014.
s István, Róbert 1 Eötvös Loránd University, MTA-ELTE NumNet Research Group 2 University of Technology Budapest, MTA-ELTE NumNet Research Group Analysis of Evolution Equations, Innsbruck, Austria 14-17
More informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 6: Population dynamics Ilya Potapov Mathematics Department, TUT Room TD325 Living things are dynamical systems Dynamical systems theory
More informationOn epidemic models with nonlinear cross-diffusion
20th International Congress on Modelling and Simulation, Adelaide, Australia, 1 6 December 2013 www.mssanz.org.au/modsim2013 On epidemic models with nonlinear cross-diffusion S. Berres a and J. Gonzalez-Marin
More informationStability analysis of delayed SIR epidemic models with a class of nonlinear incidence rates
Published in Applied Mathematics and Computation 218 (2012 5327-5336 DOI: 10.1016/j.amc.2011.11.016 Stability analysis of delayed SIR epidemic models with a class of nonlinear incidence rates Yoichi Enatsu
More informationC2 Differential Equations : Computational Modeling and Simulation Instructor: Linwei Wang
C2 Differential Equations 4040-849-03: Computational Modeling and Simulation Instructor: Linwei Wang Part IV Dynamic Systems Equilibrium: Stable or Unstable? Equilibrium is a state of a system which does
More informationEpidemic model for influenza A (H1N1)
Epidemic model for influenza A (H1N1) Modeling the outbreak of the pandemic in Kolkata, West Bengal, India, 2010 Sandipan Dey India July 8, 2017 2 Summary In this report, the spread of the pandemic influenza
More informationOn the Spread of Epidemics in a Closed Heterogeneous Population
On the Spread of Epidemics in a Closed Heterogeneous Population Artem Novozhilov Applied Mathematics 1 Moscow State University of Railway Engineering (MIIT) the 3d Workshop on Mathematical Models and Numerical
More information2 One-dimensional models in discrete time
2 One-dimensional models in discrete time So far, we have assumed that demographic events happen continuously over time and can thus be written as rates. For many biological species with overlapping generations
More informationON THE GLOBAL STABILITY OF AN SIRS EPIDEMIC MODEL WITH DISTRIBUTED DELAYS. Yukihiko Nakata. Yoichi Enatsu. Yoshiaki Muroya
Manuscript submitted to AIMS Journals Volume X, Number X, XX 2X Website: ttp://aimsciences.org pp. X XX ON THE GLOBAL STABILITY OF AN SIRS EPIDEMIC MODEL WITH DISTRIBUTED DELAYS Yukiiko Nakata Basque Center
More informationMathematical Modelling in Biology Lecture Notes
Mathematical Modelling in Biology Lecture Notes Ruth Baker Trinity Term 2016 Contents 1 Discrete-time models for a single species 1 1.1 Examples........................................ 1 1.2 Dynamic behaviour...................................
More informationBiological Models With Time Delay
Biological Models With Time Delay Project Number: MA-RYL-2016 A Major Qualifying Project submitted to the Faculty of the WORCESTER POLYTECHNIC INSTITUTE in partial fulfillment of the requirements for the
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 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 informationA FRACTIONAL ORDER SEIR MODEL WITH DENSITY DEPENDENT DEATH RATE
Hacettepe Journal of Mathematics and Statistics Volume 4(2) (211), 287 295 A FRACTIONAL ORDER SEIR MODEL WITH DENSITY DEPENDENT DEATH RATE Elif Demirci, Arzu Unal and Nuri Özalp Received 21:6 :21 : Accepted
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for COMPUTATIONAL BIOLOGY A. COURSE CODES: FFR 110, FIM740GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for COMPUTATIONAL BIOLOGY A COURSE CODES: FFR 110, FIM740GU, PhD Time: Place: Teachers: Allowed material: Not allowed: June 8, 2018, at 08 30 12 30 Johanneberg Kristian
More informationMathematical modelling and controlling the dynamics of infectious diseases
Mathematical modelling and controlling the dynamics of infectious diseases Musa Mammadov Centre for Informatics and Applied Optimisation Federation University Australia 25 August 2017, School of Science,
More informationAn intrinsic connection between Richards model and SIR model
An intrinsic connection between Richards model and SIR model validation by and application to pandemic influenza data Xiang-Sheng Wang Mprime Centre for Disease Modelling York University, Toronto (joint
More informationPARAMETER ESTIMATION IN EPIDEMIC MODELS: SIMPLIFIED FORMULAS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 19, Number 4, Winter 211 PARAMETER ESTIMATION IN EPIDEMIC MODELS: SIMPLIFIED FORMULAS Dedicated to Herb Freedman on the occasion of his seventieth birthday
More informationIntroduction: What one must do to analyze any model Prove the positivity and boundedness of the solutions Determine the disease free equilibrium
Introduction: What one must do to analyze any model Prove the positivity and boundedness of the solutions Determine the disease free equilibrium point and the model reproduction number Prove the stability
More informationThe SIR Disease Model Trajectories and MatLab
The SIR Disease Model Trajectories and MatLab James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University November 17, 2013 Outline Reviewing the SIR
More informationHepatitis C Mathematical Model
Hepatitis C Mathematical Model Syed Ali Raza May 18, 2012 1 Introduction Hepatitis C is an infectious disease that really harms the liver. It is caused by the hepatitis C virus. The infection leads to
More informationA New Mathematical Approach for. Rabies Endemy
Applied Mathematical Sciences, Vol. 8, 2014, no. 2, 59-67 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.39525 A New Mathematical Approach for Rabies Endemy Elif Demirci Ankara University
More informationNon-Linear Models Cont d: Infectious Diseases. Non-Linear Models Cont d: Infectious Diseases
Cont d: Infectious Diseases Infectious Diseases Can be classified into 2 broad categories: 1 those caused by viruses & bacteria (microparasitic diseases e.g. smallpox, measles), 2 those due to vectors
More informationMarkov-modulated interactions in SIR epidemics
Markov-modulated interactions in SIR epidemics E. Almaraz 1, A. Gómez-Corral 2 (1)Departamento de Estadística e Investigación Operativa, Facultad de Ciencias Matemáticas (UCM), (2)Instituto de Ciencias
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 informationDiscrete time mathematical models in ecology. Andrew Whittle University of Tennessee Department of Mathematics
Discrete time mathematical models in ecology Andrew Whittle University of Tennessee Department of Mathematics 1 Outline Introduction - Why use discrete-time models? Single species models Geometric model,
More informationAnalysis of bacterial population growth using extended logistic Growth model with distributed delay. Abstract INTRODUCTION
Analysis of bacterial population growth using extended logistic Growth model with distributed delay Tahani Ali Omer Department of Mathematics and Statistics University of Missouri-ansas City ansas City,
More informationCopyright is owned by the Author of the thesis. Permission is given for a copy to be downloaded by an individual for the purpose of research and
Copyright is owned by the Author of the thesis. Permission is given for a copy to be downloaded by an individual for the purpose of research and private study only. The thesis may not be reproduced elsewhere
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 informationSYMBIOTIC MODELS WITH AN SIR DISEASE
SYMBIOTIC MODELS WITH AN SIR DISEASE A Thesis Presented to the Faculty of California State Polytechnic University, Pomona In Partial Fulfillment Of the Requirements for the Degree Master of Science In
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 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 informationEstimating the Exponential Growth Rate and R 0
Junling Ma Department of Mathematics and Statistics, University of Victoria May 23, 2012 Introduction Daily pneumonia and influenza (P&I) deaths of 1918 pandemic influenza in Philadelphia. 900 800 700
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 informationKasetsart University Workshop. Mathematical modeling using calculus & differential equations concepts
Kasetsart University Workshop Mathematical modeling using calculus & differential equations concepts Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu
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 informationDensity Dependence, Stability & Allee effects ECOL 4000/6000
Density Dependence, Stability & Allee effects ECOL 4000/6000 Quiz1: For each statement, write True or False e.g. A: True, B: True A: intrinsic growth rate for gray populanon > intrinsic growth rate for
More informationMath 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations
Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Junping Shi College of William and Mary, USA Equilibrium Model: x n+1 = f (x n ), here f is a nonlinear function
More informationMath 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA Molecular biology and Biochemical kinetics Molecular biology is one of
More informationComputer Simulation and Data Analysis in Molecular Biology and Biophysics
Victor Bloomfield Computer Simulation and Data Analysis in Molecular Biology and Biophysics An Introduction Using R Springer Contents Part I The Basics of R 1 Calculating with R 3 1.1 Installing R 3 1.1.1
More informationSection 8.1 Def. and Examp. Systems
Section 8.1 Def. and Examp. Systems Key Terms: SIR Model of an epidemic o Nonlinear o Autonomous Vector functions o Derivative of vector functions Order of a DE system Planar systems Dimension of a system
More informationModels of Infectious Disease Formal Demography Stanford Summer Short Course James Holland Jones, Instructor. August 15, 2005
Models of Infectious Disease Formal Demography Stanford Summer Short Course James Holland Jones, Instructor August 15, 2005 1 Outline 1. Compartmental Thinking 2. Simple Epidemic (a) Epidemic Curve 1:
More informationChapter 2: Growth & Decay
Chapter 2: Growth & Decay 107/226 Introduction In this chapter, model biological systems whose population is so large or where growth is so fine-grained that continuity can be assumed. These continuous
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 informationObjective. Single population growth models
Objective Single population growth models We are given a table with the population census at different time intervals between a date a and a date b, and want to get an expression for the population. This
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More informationSTABILITY IN A MODEL FOR GENETICALLY ALTERED MOSQUITOS WITH PERIODIC PARAMETERS
STABILITY IN A MODEL FOR GENETICALLY ALTERED MOSQUITOS WITH PERIODIC PARAMETERS Hubertus F. von Bremen and Robert J. Sacker Department of Mathematics and Statistics, California State Polytechnic University,
More informationAge (x) nx lx. Population dynamics Population size through time should be predictable N t+1 = N t + B + I - D - E
Population dynamics Population size through time should be predictable N t+1 = N t + B + I - D - E Time 1 N = 100 20 births 25 deaths 10 immigrants 15 emmigrants Time 2 100 + 20 +10 25 15 = 90 Life History
More informationNumerical Modeling of the Transmission Dynamics of Influenza
The First International Symposium on Optimization and Systems Biology (OSB 07) Beijing, China, August 8 0, 2007 Copyright 2007 ORSC & APORC pp. 52 59 Numerical Modeling of the Transmission Dynamics of
More informationFixed Point Analysis of Kermack Mckendrick SIR Model
Kalpa Publications in Computing Volume, 17, Pages 13 19 ICRISET17. International Conference on Research and Innovations in Science, Engineering &Technology. Selected Papers in Computing Fixed Point Analysis
More informationBackward bifurcation and oscillations in a nested immuno-eco-epidemiological model
Journal of Biological Dynamics ISSN: 75-3758 Print) 75-3766 Online) Journal homepage: http://www.tandfonline.com/loi/tjbd Backward bifurcation and oscillations in a nested immuno-eco-epidemiological model
More informationResilience and stability of harvested predator-prey systems to infectious diseases in the predator
Resilience and stability of harvested predator-prey systems to infectious diseases in the predator Morgane Chevé Ronan Congar Papa A. Diop November 1, 2010 Abstract In the context of global change, emerging
More informationME 406 S-I-R Model of Epidemics Part 2 Vital Dynamics Included
ME 406 S-I-R Model of Epidemics Part 2 Vital Dynamics Included sysid Mathematica 6.0.3, DynPac 11.01, 1ê13ê9 1. Introduction Description of the Model In this notebook, we include births and deaths in the
More informationStochastic Differential Equations in Population Dynamics
Stochastic Differential Equations in Population Dynamics Numerical Analysis, Stability and Theoretical Perspectives Bhaskar Ramasubramanian Abstract Population dynamics in the presence of noise in the
More informationApplications in Biology
11 Applications in Biology In this chapter we make use of the techniques developed in the previous few chapters to examine some nonlinear systems that have been used as mathematical models for a variety
More informationMathematical Ecology
Mathematical Ecology Joachim Hermisson, Claus Rueffler & Meike Wittmann January 14, 2016 Literature and Software Sarah P. Otto, Troy Day: A Biologist s Guide to Mathematical Modeling in Ecology and Evolution,
More informationSlow and fast dynamics model of a Malaria with Sickle-Cell genetic disease with multi-stage infections of the mosquitoes population
Journal of Physics: Conference Series PAPER OPEN ACCESS Slow and fast dynamics model of a Malaria with Sickle-Cell genetic disease with multi-stage infections of the mosquitoes population To cite this
More informationDifference Equations
Difference Equations Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2017 M. Macauley (Clemson) Difference equations Math
More informationPredation is.. The eating of live organisms, regardless of their identity
Predation Predation Predation is.. The eating of live organisms, regardless of their identity Predation 1)Moves energy and nutrients through the ecosystem 2)Regulates populations 3)Weeds the unfit from
More informationIntermediate Differential Equations. John A. Burns
Intermediate Differential Equations Delay Differential Equations John A. Burns jaburns@vt.edu Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute and State University Blacksburg,
More informationLecture 3. Dynamical Systems in Continuous Time
Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver Yue-Xian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a
More informationBIOL 410 Population and Community Ecology. Density-dependent growth 1
BIOL 410 Population and Community Ecology Density-dependent growth 1 Exponential growth N t = N 0 e rt # Model 6, Exponential growth t
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 informationChapter 2 Lecture. Density dependent growth and intraspecific competition ~ The Good, The Bad and The Ugly. Spring 2013
Chapter 2 Lecture Density dependent growth and intraspecific competition ~ The Good, The Bad and The Ugly Spring 2013 2.1 Density dependence, logistic equation and carrying capacity dn = rn K-N Dt K Where
More informationM469, Fall 2010, Practice Problems for the Final
M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More informationEcology Regulation, Fluctuations and Metapopulations
Ecology Regulation, Fluctuations and Metapopulations The Influence of Density on Population Growth and Consideration of Geographic Structure in Populations Predictions of Logistic Growth The reality of
More informationMathematical Analysis of Epidemiological Models: Introduction
Mathematical Analysis of Epidemiological Models: Introduction Jan Medlock Clemson University Department of Mathematical Sciences 8 February 2010 1. Introduction. The effectiveness of improved sanitation,
More informationMANAGEMENT AND ANALYSIS OF BIOLOGICAL POPULATIONS
',' Developments in Agricultural and Managed-Forest Ecology, 8 MANAGEMENT AND ANALYSIS OF BIOLOGICAL POPULATIONS by BEAN-SAN GOH Department ofmathematics, University of Western Australia, Nedlands, W.A.
More informationThe Spreading of Epidemics in Complex Networks
The Spreading of Epidemics in Complex Networks Xiangyu Song PHY 563 Term Paper, Department of Physics, UIUC May 8, 2017 Abstract The spreading of epidemics in complex networks has been extensively studied
More informationPartially Observed Stochastic Epidemics
Institut de Mathématiques Ecole Polytechnique Fédérale de Lausanne victor.panaretos@epfl.ch http://smat.epfl.ch Stochastic Epidemics Stochastic Epidemic = stochastic model for evolution of epidemic { stochastic
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 informationToday. Introduction to Differential Equations. Linear DE ( y = ky ) Nonlinear DE (e.g. y = y (1-y) ) Qualitative analysis (phase line)
Today Introduction to Differential Equations Linear DE ( y = ky ) Nonlinear DE (e.g. y = y (1-y) ) Qualitative analysis (phase line) Differential equations (DE) Carbon dating: The amount of Carbon-14 in
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 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 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 informationFigure The Threshold Theorem of epidemiology
K/a Figure 3 6. Assurne that K 1/ a < K 2 and K 2 / ß < K 1 (a) Show that the equilibrium solution N 1 =0, N 2 =0 of (*) is unstable. (b) Show that the equilibrium solutions N 2 =0 and N 1 =0, N 2 =K 2
More information