Stochastic Viral Dynamics with Beddington-DeAngelis Functional Response
|
|
- Tamsyn Stewart
- 5 years ago
- Views:
Transcription
1 Stochastic Viral Dynamics with Beddington-DeAngelis Functional Response Junyi Tu, Yuncheng You University of South Florida, USA IMA Workshop in Memory of George R. Sell June 016
2 Outline 1. Stochastic SIR model with the functional response. Positive invariant region 3. Stationary distribution 4. Pathwise and moment bounds 5. Moment Lyapunov exponent 6. Persistence and extinction
3 1. Stochastic Model Equations ( βxy dx = λ δx 1 + ax + by ( βxy dy = 1 + ax + by qy dz = (ky γz)dt + σ 3 z db 3 (t), ) dt + σ 1 x db 1 (t), ) dt + σ y db (t), (1) where the susceptible cells with density x(t) are generated at a constant rate λ, die at a density-proportional rate δx, and become infective cells y(t) at a rate βxy/(1 + ax + by) called Beddington-DeAngelis functional response. The density of recovered cells is z(t).
4 Beddington-DeAngelis Response Classical infectious and epidemic disease models use bilinear incidence or interaction terms such as βxy. The nonlinear incidence rate with this Holling type II response βxy/(1 + ax + by) was introduced separately by Beddington (1975) and DeAngelis (1975) in ecology. It was also used in the study of the interaction dynamics of HIV-1 virus and uninfected CD4 + T cells, cf. Huang et al. (009).
5 Deterministic Dynamics and Stability If the comprehensive reproduction ratio of the virus, R 0 = kβλ δqγ + aqγλ 1, the disease-free equilibrium E 0 = (λ/δ, 0, 0) is asymptotically stable. Disease Extinction. If R 0 > 1, the endemic equilibrium E 1 = (x 0, y 0, z 0 ) in R 3 + is globally asymptotically stable Persistence, where cf. Huang-Ma-Takeuch (009) qγ + kbλ x 0 = kβ aqγ + bkδ, kβλ aλqγ δqγ y 0 = q(kβ aqγ + bkδ), k(kβλ aλqγ δqγ) z 0 = qγ(kβ aqγ + bkδ).
6 Related Models and Researches Quite a few research work on stochastic population dynamics, in particular Lotka-Volterra type equations for predator-prey models whose growth/death rates perturbed by noise and adding various functional responses: bilibear response, ratio-dependent response, Leslie-Gower response, and Beddington-DeAngelis response on many topics: Khasminskii and Klebaner (001), Rudnicki (003), Rudnicki-Pichor (007), Deng et al (008), Ji-Jiang-Li (011), Mao (011), Ton-Yagi (011), Li-Shuang (013), Zou-Fan-Wang (013), etc. Diffusive predator-prey Model with Bedington-DeAngelis response: Chen-Wang (005).
7 Findings by Shown Results With stochastic perturbation, there will be no positive equilibrium points. Instead one can try to prove that (open problem) there exists a stationary distribution on the positive invariant region R n +, which can be viewed as a counterpart in the stochastic dynamics arena. If intensity of noise is relatively small, hopefully (not always) one can show that stochastic asymptotic dynamics imitates the corresponding deterministic dynamics in some way, almost surely, in moment, or in time average. If noise is sufficiently large, something dramatic and/or unexpected in dynamics may occur: extinction of all species, suppress or express exponential growth or explosion. Large random environment, weather, or epidemic diseases effect seems to be the decisive responsible for the extinction of some species in the nature.
8 Preliminaries Assume that B i (t) are independent standard Brownian motion defined on the canonical probability space (Ω, F, P). R n + = {x R n : x i > 0 for all 1 i n}. Itô formula: For V C,1 (R n [0, ), R + ), dv (X (t), t) = LV (X, t) dt + V x (X, t)g(x, t) db, for a SDE dx = f (X, t) dt + g(x, t) db, and LV (X (t), t) = V t (X, t) + V x (X, t)f (X, t) + 1 Trace [ g T (X, t)v xx (X, t)g(x, t) ].
9 . Global Existence of Positive Solutions Theorem Under the condition q k, for any positive initial data (x 0, y 0, z 0 ), there is a unique positive solution (x(t), y(t), z(t)) of the model SDE (1) such that the solution will remain in R 3 + for all t 0 with probability one. It shows that R 3 + is a positive invariant region for the model system. The solution X (t) = (x(t), y(t), z(t)) is a time-homogeneous Markov process. Let the maximal existence interval of a unique local pathwise positive solution be τ e. We show τ e = a.s. V (x, y, z) = (x 1 log x) + (y 1 log y) + (z 1 log z).
10 Sketch of Proof dv Kdt + σ 1 (x 1) db 1 + σ (y 1) db + σ 3 (z 1) db 3, K = 1 (σ 1 + σ + σ 3) + λ + δ + q + γ + β b. It suffices to prove τ = lim m τ m = a.s. { ( ) } 3 1 τ m = inf t [0, τ e ) : (x(t), y(t), z(t)) m, m. Otherwise, if τ <, then P{τ T } > ε. It leads to [ ε (m 1 log m) > V (x 0, y 0, z 0 ) + KT ( 1 m 1 + log m )].
11 3. Stationary Distribution Definition Let P X0,t( ) be the probability measure induced by a stochastic process {X (t)} t 0 in R n + (or R n ) over (Ω, F, P) with X (0) = X 0, P X0,t(S) = P{ω Ω : X (t, ω) S}, S B(R n +), If there is a probability measure µ( ) on the measurable space (R n +, B(R n +)) such that P X0,t( ) µ( ) as t, in distribution for any X 0 R n +, then we say that X (t) has a stationary distribution µ( ).
12 Khasminskii Assumption Assumption: There is a bounded open set U R n + with regular boundary and the following properties: 1. In a neighborhood of U, the smallest eigenvalue of the diffusion matrix g(x) is uniformly bounded away from zero.. For any x R n +\U, the mean time τ at which a path from x reaches the set U satisfies sup x K E x [τ] <, for every compact set K R n +.
13 Khasminskii Lemma Lemma Under the Khasminskii Assumption, the Markov process X (t) has a stationary distribution µ( ), µ(r n \ R n +) = 0. For any bounded continuous function Q(ξ), it holds that E[Q(X (t, ξ)] dµ(ξ) = Q(ξ) dµ(ξ), t 0. R n + For any integrable function F (X ) with respect to µ, the R n + ergodic property holds, { 1 T P lim F (X (t, ξ)) dt = T T 0 R n + F (ξ) dµ(ξ) } = 1.
14 The Existence of Stationary Distribution Theorem If there is a positive equilibrium point (x 0, y 0, z 0 ) R 3 + for the corresponding deterministic system of ODE, then there exists a stationary distribution with respect to the Markov process generated by the positive solutions (x(t), y(t), z(t)) of the system of SDE (1). Proof. Construct another V -function ) 1 + by 0 V (x, y, z) = (x x 0 x 0 log xx0 1 + ax 0 + by 0 +y y 0 y 0 log y + q ) (z z 0 z 0 log zz0. y 0 k
15 Sketch of Proof By the equilibrium conditiom we have where and LV (x(t), y(t), z(t)) = Q(x(t), x(t), z(t)) + σ σ = (1 + by 0)x 0 σ1 (1 + ax 0 + by 0 ) + y 0σ + qz 0σ3 > 0, k Q(x, y, z) = δ(1 + by 0) x(1 + ax 0 + by 0 ) (x x 0) qy 0 b(1 + ax)(y y 0 ) z 0 (1 + ax + by)(1 + ax + by 0 ) + qy 0R(x, y, z) 0, R(x, y, z) = 4 (1 + ax + by 0)x 0 (1 + ax 0 + by 0 )x (1 + ax 0 + by 0 )xy 0 z (1 + ax + by)x 0 yz 0 yz 0 y 0 z (1 + ax + by) (1 + ax + by 0 ) 0.
16 Continued Q(x, y, z) = 0, iff (x, y, z) = (x 0, y 0, z 0 ). Q(x, y, z) > 0, if otherwise in R 3 +. Q(x, y, z) as (x, y, z). The diffusion matrix g(x, y, z) = diag (σ 1x, σ y, σ 3z ) has the smallest eigenvalue uniformly bounded from zero, on a nbhd of any compact set in R 3 +\U. The two conditions in Khasminskii Assumption are satisfied for the bounded open set U = {(x, y, z) : Q(x, y, z) < σ} R 3 +. By Khasminskii Lemma, the system of SDE (1) has a unique stationary distribution µ( ) on R 3 +.
17 4. Moment and Pathwise Estimates Lemma For p > 1, the solution of the Bernoulli equation [ ) ] dv dt = p v(t) (δ σ (p 1) + λv 1p (t), v(0) = v 0, is given by where v(t) = [ v 1 p 0 e θpt + λ θ p ( 1 e θ pt )] p, θ p = δ σ (p 1).
18 Uniform bounds of the p-th Moments Theorem Suppose the following condition is satisfied, δ σ 1 (p 1) > 0, q σ (p 1) > 0, γ σ 3 (p 1) > 0. Then for p > 1 and all positive solution of the stochastic viral equations (1) with (x 0, y 0, z 0 ) R 3 +, lim sup t lim sup t lim sup t E[x p (t)] L 1 (p), E[y p (t)] L (p), E[z p (t)] L 3 (p). Remark. For p = 1, bounds are given thru linear equations.
19 Uniform Moment Bounds ( ) p λ L 1 (p) = δ σ 1 (p 1) ( ) p 4βbλ L (p) = b(δ σ1 (p 1))(q σ (p 1)) ( 8kβλ L 3 (p) = b(δ σ1 (p 1))(q σ (p 1))(γ σ 3 (p 1)) where δ, q, γ are death rates, λ, β, k are production or growth rates, and σ 1, σ, σ 3 are noise intensities. ) p
20 Sketch of Proof For positive solutions, d(x p ) px p 1 ( λ δx + σ 1 ) (p 1)x dt + px p σ 1 db 1 (t). de[x p (t)] dt ( ) p δ σ 1 (p 1) E[x p ] + pλ E[x p ] 1 1 p. Use Lemma for Bernoulli equations and comparison theorem. For the y-equation, the nonlinear term of functional response βxy 1 + ax + by βx b.
21 Pathwise Upper/Lower Bounds Theorem Every positive solution (x(t), y(t), z(t)) of the system (1) with (x 0, y 0, z 0 ) R 3 + satisfies Φ l (t) x(t) Φ u (t), Ψ l (t) y(t) Ψ u (t), Γ l (t) z(t) Γ u (t), t 0, a.s. where { ( } t δ+ β b Φ u (t) = λ exp + σ 1 )(t s)+σ 1 (B 1 (t) B 1 (s)) ds 0 +x 0 exp { ( δ+ β b + σ 1 )t+σ 1 B 1 (t) }
22 Continued Ψ u (t) = β b Γ u (t) = k t 0 + y 0 e t 0 + z 0 e ( ) Φ u (s)e q+ σ (t s)+σ (B (t) B (s)) ds ( ) q+ σ t+σ B (t), ( ) Ψ u (s)e γ+ σ 3 (t s)+σ 3 (B 3 (t) B 3 (s)) ds ( ) γ+ σ 3 t+σ 3 B 3 (t).
23 Absorbing Property in Time Average Theorem Under the same condition that the noise intensities are relatively small, δ σ 1 (p 1) > 0, q σ (p 1) > 0, γ σ 3 (p 1) > 0, all the positive solutions has the property a.s. 1 t lim x p (s) ds = ξ p dµ(ξ, η, ζ) L 1 (p), t t 0 R+ 3 1 t lim y p (s) ds = η p dµ(ξ, η, ζ) L (p), t t 0 R+ 3 1 t lim z p (s) ds = ζ p dµ(ξ, η, ζ) L 3 (p). t t 0 R 3 +
24 5. Moment Lyapunov Exponent Definition The p-th moment Lyapunov exponent of a pathwise solution X (t, X 0 ) of the SDE (1) is defined by Λ(p) = lim sup t 1 t log E X (t, X 0) p, p 1. Theorem Under the same assumption as above, for p 1, Λ(p) = lim sup t 1 t log E X (t, X 0) p 0, R 3 + for any positive solution X (t, X 0 ) = (x(t, x 0 ), y(t, y 0 ), z(t, z 0 )) of the stochastic viral equations (1) with X 0 R 3 +.
25 Sketch of Proof The p-th moment of the geometric Brownian motion for p > 1 is given by ds = δs(t) dt + σs(t) db E S(t) p = S(0) p exp [ ( ) ] (p 1) p σ δ t. Directly estimate: in view of log(1 + x) x for x > 0, lim sup t lim sup t + lim sup t 1 t log (E x(t, x 0) 1 ) lim sup t t log (E Φ u(t) ) ( 1 t log 1 + x ) 0 (δ σ1) exp{(σ λ 1 δ)t} 1 t log λ δ σ 1 = 0. when p =. Then bootstrap for the y and z components.
26 6. Persistence and Extinction For a deterministic or stochastic model in population dynamics for ecology or viral and epidemic dynamics for disease control, the two most important questions are the analysis and prediction of persistence and extinction. Definition A component x(t) of the system is said to be persistent in mean, if lim inf E[x(t)] > 0, t or persistent a.s. if lim inf t x(t) > 0. Other types of persistence and extinction can be defined.
27 Useful Lemma for Stochastic Predator-Prey Eqns Lemma Consider a one-dimensional SDE dx = X (t)(a bx (t)) dt + σx (t) db. If a > σ, then for any solution X (t) with X 0 > 0, one has 1 lim log X (t) = 0, a.s. t t and 1 t lim t t X (s) ds = a σ, a.s. 0 b Thus the solutions are persistent in time mean. But this does not work for the stochastic SIR model with the Beddington-DeAngelis functional response.
28 x(t) Φ u (t) Lemma lim 1 t log x(t) = 0 +x 0 exp λ δ + β + σ 1 b [ ( exp [σ 1 ( δ + β b + σ 1 B 1 (t) min B 1(s) 0 s t ) ] t + σ 1 B 1 (t) )] The running max of BM, B 1 (t) min 0 s t B 1 (s) = B 1 (t) in distribution, and log(1 + x) < x lead to lim sup t + lim sup t 1 t log Φ u(t) lim sup 1 t log ( 1 [σ 1 B 1 (t) min t t B 1(s) 0 s t [ 1 + x ] 0(δ + β/b + σ1/) exp (Ratio) λ )] = 0.
29 The Exponential Transformation Let x(t) = e u(t), y(t) = u v(t), z(t) = e w(t). The SDE of this stochastic viral model is converted to ( ( ) du = δ + σ 1 + λe u βe v 1 + a e u + be ( ( ) v dv = q + σ βe u a e u + be ( ( )) v dw = ke v w γ + σ 3 dt + σ 3 db 3. Here we see dv ) dt + σ db, ( ( )) β a q + σ dt + σ db. ) dt + σ 1 db 1,
30 Extinction of the Infectives Proposition If β/a < q + σ /, then lim t y(t) = 0. The infective cells tends to extinction. By the stochastic comparison theorem, cf. Ikeda and Watanabe (1981), p. 35, ( ( )) β v(t) v 0 + a q + σ t + σ B (t) ( ( )) β = v 0 + a q + σ σ B (t) t 1 + ( ( )), β q + σ a t by SLLN for Brownian motion. Then lim y(t) = lim e v(t) = 0.
31 Persistence of the Susceptible From the u-equation, ( ( ) du = δ + σ 1 + λe u βe v 1 + a e u + be ( ( v λe u δ + σ 1 + β )) dt + σ 1 db 1, b ) dt + σ 1 db 1, and dũ = (λe u θ) dt + σ 1 db 1, θ = δ + σ 1/ + β/b, has a unique stationary distribution whose probability density is a solution of the stationary Fokker-Planck equation The solution is 1 d p σ 1 dξ d dξ [( λe ξ θ ) p] = 0. p(ξ) = C exp where C is the normalizing const. { θ ξ + λ } e ξ, σ1 σ1
32 Persistence of the Susceptible Accordingly, cf. Rudnicki (003), Skorokhod (1987), u(t) ũ(t) which converges in distribution to a stationary solution ũ whoe probability density is p(ξ), as t. Proposition The susceptible cells will be persistent in mean and = C lim inf t 0 E[x(t)] e ξ p(ξ) dξ 0 {( exp 1 θ ) ξ + λ } e ξ dξ σ1 σ1 Cσ 1 σ 1 θ = Cσ 1 (δ + β/b) > 0.
33 Persistence of the Recovered Proposition Under the condition β/a < q + σ /, The recovered cells will be persistent in mean. From the w-equation ( ( )) dw = ke v w γ + σ 3 dt + σ 3 db 3 with y(t) = e v(t) 0 under the condition, the proof is similar. Many questions still unanswered: stability, bifurcation, stochastic pattern formation for SPDE models, etc.
34 Conclusion The study shows that under the condition q k the positive cone R 3 + is invariant with probability one for the pathwise solutions of this model of stochastic viral dynamics. Construction of Lyapunov functions shows that there is a unique stationary distribution whose probability density function is a steady state of the Fokker-Planck equation. Moments and pathwise upper/lower bounds and some asymptotic estimates are obtained, which make possible to discuss the persistence and extinction of the virus diseases. All these can be used for better reaveal of features of models of stochastic SIR equations with different functional responses and for better understanding of medical/epidemic incident events.
35 THANKS
A Stochastic Viral Infection Model with General Functional Response
Nonlinear Analysis and Differential Equations, Vol. 4, 16, no. 9, 435-445 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.1988/nade.16.664 A Stochastic Viral Infection Model with General Functional Response
More informationStationary distribution and pathwise estimation of n-species mutualism system with stochastic perturbation
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 6), 936 93 Research Article Stationary distribution and pathwise estimation of n-species mutualism system with stochastic perturbation Weiwei
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationApplied Mathematics Letters. Stationary distribution, ergodicity and extinction of a stochastic generalized logistic system
Applied Mathematics Letters 5 (1) 198 1985 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Stationary distribution, ergodicity
More informationA Delayed HIV Infection Model with Specific Nonlinear Incidence Rate and Cure of Infected Cells in Eclipse Stage
Applied Mathematical Sciences, Vol. 1, 216, no. 43, 2121-213 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.216.63128 A Delayed HIV Infection Model with Specific Nonlinear Incidence Rate and
More informationAsymptotic behaviour of the stochastic Lotka Volterra model
J. Math. Anal. Appl. 87 3 56 www.elsevier.com/locate/jmaa Asymptotic behaviour of the stochastic Lotka Volterra model Xuerong Mao, Sotirios Sabanis, and Eric Renshaw Department of Statistics and Modelling
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 informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationOn a non-autonomous stochastic Lotka-Volterra competitive system
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 7), 399 38 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa On a non-autonomous stochastic Lotka-Volterra
More informationResearch Article On the Stability Property of the Infection-Free Equilibrium of a Viral Infection Model
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume, Article ID 644, 9 pages doi:.55//644 Research Article On the Stability Property of the Infection-Free Equilibrium of a Viral
More informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationOptimal harvesting policy of a stochastic delay predator-prey model with Lévy jumps
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 1 (217), 4222 423 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa Optimal harvesting policy of
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationNonlinear representation, backward SDEs, and application to the Principal-Agent problem
Nonlinear representation, backward SDEs, and application to the Principal-Agent problem Ecole Polytechnique, France April 4, 218 Outline The Principal-Agent problem Formulation 1 The Principal-Agent problem
More informationTheoretical Tutorial Session 2
1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations
More informationOn continuous time contract theory
Ecole Polytechnique, France Journée de rentrée du CMAP, 3 octobre, 218 Outline 1 2 Semimartingale measures on the canonical space Random horizon 2nd order backward SDEs (Static) Principal-Agent Problem
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 informationI forgot to mention last time: in the Ito formula for two standard processes, putting
I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
More informationSIS epidemics on Networks
SIS epidemics on etworks Piet Van Mieghem in collaboration with Eric Cator, Ruud van de Bovenkamp, Cong Li, Stojan Trajanovski, Dongchao Guo, Annalisa Socievole and Huijuan Wang 1 EURADOM 30 June-2 July,
More informationInteractions. Yuan Gao. Spring Applied Mathematics University of Washington
Interactions Yuan Gao Applied Mathematics University of Washington yuangao@uw.edu Spring 2015 1 / 27 Nonlinear System Consider the following coupled ODEs: dx = f (x, y). dt dy = g(x, y). dt In general,
More informationBrownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion
Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of
More informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationRandom and Deterministic perturbations of dynamical systems. Leonid Koralov
Random and Deterministic perturbations of dynamical systems Leonid Koralov - M. Freidlin, L. Koralov Metastability for Nonlinear Random Perturbations of Dynamical Systems, Stochastic Processes and Applications
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationSome Properties of NSFDEs
Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationNew ideas in the non-equilibrium statistical physics and the micro approach to transportation flows
New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows Plenary talk on the conference Stochastic and Analytic Methods in Mathematical Physics, Yerevan, Armenia,
More informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationIntroduction to numerical simulations for Stochastic ODEs
Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
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 informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationGlobal Threshold Dynamics of a Stochastic Differential Equation SIS Model
arxiv:156.232v2 [math.ds] 22 Aug 216 Global Threshold Dynamics of a Stochastic Differential Equation SIS Model Chuang Xu Department of Mathematics, Harbin Institute of Technology, Harbin, Heilongjiang
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 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 informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationA new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009
A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance
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 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 informationStochastic differential equations in neuroscience
Stochastic differential equations in neuroscience Nils Berglund MAPMO, Orléans (CNRS, UMR 6628) http://www.univ-orleans.fr/mapmo/membres/berglund/ Barbara Gentz, Universität Bielefeld Damien Landon, MAPMO-Orléans
More informationPartial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula
Partial Differential Equations with Applications to Finance Seminar 1: Proving and applying Dynkin s formula Group 4: Bertan Yilmaz, Richard Oti-Aboagye and Di Liu May, 15 Chapter 1 Proving Dynkin s formula
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationMean-field SDE driven by a fractional BM. A related stochastic control problem
Mean-field SDE driven by a fractional BM. A related stochastic control problem Rainer Buckdahn, Université de Bretagne Occidentale, Brest Durham Symposium on Stochastic Analysis, July 1th to July 2th,
More informationLecture 17 Brownian motion as a Markov process
Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is
More informationStrong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term
1 Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term Enrico Priola Torino (Italy) Joint work with G. Da Prato, F. Flandoli and M. Röckner Stochastic Processes
More informationConvergence to equilibrium for rough differential equations
Convergence to equilibrium for rough differential equations Samy Tindel Purdue University Barcelona GSE Summer Forum 2017 Joint work with Aurélien Deya (Nancy) and Fabien Panloup (Angers) Samy T. (Purdue)
More informationA numerical method for solving uncertain differential equations
Journal of Intelligent & Fuzzy Systems 25 (213 825 832 DOI:1.3233/IFS-12688 IOS Press 825 A numerical method for solving uncertain differential equations Kai Yao a and Xiaowei Chen b, a Department of Mathematical
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 informationThe Effect of Stochastic Migration on an SIR Model for the Transmission of HIV. Jan P. Medlock
The Effect of Stochastic Migration on an SIR Model for the Transmission of HIV A Thesis Presented to The Faculty of the Division of Graduate Studies by Jan P. Medlock In Partial Fulfillment of the Requirements
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More information9.2 Branching random walk and branching Brownian motions
168 CHAPTER 9. SPATIALLY STRUCTURED MODELS 9.2 Branching random walk and branching Brownian motions Branching random walks and branching diffusions have a long history. A general theory of branching Markov
More informationSensitivity and Stability Analysis of Hepatitis B Virus Model with Non-Cytolytic Cure Process and Logistic Hepatocyte Growth
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 3 2016), pp. 2297 2312 Research India Publications http://www.ripublication.com/gjpam.htm Sensitivity and Stability Analysis
More informationExistence, Uniqueness and Stability of Invariant Distributions in Continuous-Time Stochastic Models
Existence, Uniqueness and Stability of Invariant Distributions in Continuous-Time Stochastic Models Christian Bayer and Klaus Wälde Weierstrass Institute for Applied Analysis and Stochastics and University
More informationAn Introduction to Malliavin calculus and its applications
An Introduction to Malliavin calculus and its applications Lecture 3: Clark-Ocone formula David Nualart Department of Mathematics Kansas University University of Wyoming Summer School 214 David Nualart
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 information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More information1 Brownian Local Time
1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =
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 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 informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationSTUDY OF THE DYNAMICAL MODEL OF HIV
STUDY OF THE DYNAMICAL MODEL OF HIV M.A. Lapshova, E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper is devoted to the study of the dynamical model of HIV. An application
More informationDISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS. Ed Allen. NIMBioS Tutorial: Stochastic Models With Biological Applications
DISCRETE-TIME STOCHASTIC MODELS, SDEs, AND NUMERICAL METHODS Ed Allen NIMBioS Tutorial: Stochastic Models With Biological Applications University of Tennessee, Knoxville March, 2011 ACKNOWLEDGEMENT I thank
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationLYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR
LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu
More informationStability of epidemic model with time delays influenced by stochastic perturbations 1
Mathematics and Computers in Simulation 45 (1998) 269±277 Stability of epidemic model with time delays influenced by stochastic perturbations 1 Edoardo Beretta a,*, Vladimir Kolmanovskii b, Leonid Shaikhet
More informationOn the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem
On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family
More informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationBayesian Inference and the Symbolic Dynamics of Deterministic Chaos. Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2
How Random Bayesian Inference and the Symbolic Dynamics of Deterministic Chaos Christopher C. Strelioff 1,2 Dr. James P. Crutchfield 2 1 Center for Complex Systems Research and Department of Physics University
More informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationRough Burgers-like equations with multiplicative noise
Rough Burgers-like equations with multiplicative noise Martin Hairer Hendrik Weber Mathematics Institute University of Warwick Bielefeld, 3.11.21 Burgers-like equation Aim: Existence/Uniqueness for du
More informationUNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.
More informationCOMPETITIVE OR WEAK COOPERATIVE STOCHASTIC LOTKA-VOLTERRA SYSTEMS CONDITIONED ON NON-EXTINCTION. Université de Toulouse
COMPETITIVE OR WEAK COOPERATIVE STOCHASTIC LOTKA-VOLTERRA SYSTEMS CONITIONE ON NON-EXTINCTION PATRICK CATTIAUX AN SYLVIE MÉLÉAR Ecole Polytechnique Université de Toulouse Abstract. We are interested in
More informationMean-field dual of cooperative reproduction
The mean-field dual of systems with cooperative reproduction joint with Tibor Mach (Prague) A. Sturm (Göttingen) Friday, July 6th, 2018 Poisson construction of Markov processes Let (X t ) t 0 be a continuous-time
More informationExact Simulation of Diffusions and Jump Diffusions
Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction
More informationStrong Markov property of determinatal processes
Strong Markov property of determinatal processes Hideki Tanemura Chiba university (Chiba, Japan) (August 2, 2013) Hideki Tanemura (Chiba univ.) () Markov process (August 2, 2013) 1 / 27 Introduction The
More informationAn Introduction to Stochastic Epidemic Models
An Introduction to Stochastic Epidemic Models Linda J. S. Allen Department of Mathematics and Statistics Texas Tech University Lubbock, Texas 79409-1042, U.S.A. linda.j.allen@ttu.edu 1 Introduction The
More informationMATH 220 solution to homework 5
MATH 220 solution to homework 5 Problem. (i Define E(t = k(t + p(t = then E (t = 2 = 2 = 2 u t u tt + u x u xt dx u 2 t + u 2 x dx, u t u xx + u x u xt dx x [u tu x ] dx. Because f and g are compactly
More informationSolution of Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
More informationRandom attractors and the preservation of synchronization in the presence of noise
Random attractors and the preservation of synchronization in the presence of noise Peter Kloeden Institut für Mathematik Johann Wolfgang Goethe Universität Frankfurt am Main Deterministic case Consider
More informationKai Lai Chung
First Prev Next Go To Go Back Full Screen Close Quit 1 Kai Lai Chung 1917-29 Mathematicians are more inclined to build fire stations than to put out fires. Courses from Chung First Prev Next Go To Go Back
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationStochastic Homogenization for Reaction-Diffusion Equations
Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires
More informationIntroduction to Stochastic SIR Model
Introduction to Stochastic R Model Chiu- Yu Yang (Alex), Yi Yang R model is used to model the infection of diseases. It is short for Susceptible- Infected- Recovered. It is important to address that R
More informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationAn Uncertain Control Model with Application to. Production-Inventory System
An Uncertain Control Model with Application to Production-Inventory System Kai Yao 1, Zhongfeng Qin 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 2 School of Economics
More information