MA 777: Topics in Mathematical Biology

Similar documents
MA 138: Calculus II for the Life Sciences

A&S 320: Mathematical Modeling in Biology

Mathematical Models of Biological Systems

1 (t + 4)(t 1) dt. Solution: The denominator of the integrand is already factored with the factors being distinct, so 1 (t + 4)(t 1) = A

Stochastic Simulation.

MA 137 Calculus 1 with Life Science Application A First Look at Differential Equations (Section 4.1.2)

Lecture 4 The stochastic ingredient

Today. Introduction to Differential Equations. Linear DE ( y = ky ) Nonlinear DE (e.g. y = y (1-y) ) Qualitative analysis (phase line)

Chapter 2 Lecture. Density dependent growth and intraspecific competition ~ The Good, The Bad and The Ugly. Spring 2013

8 Ecosystem stability

Mathematical and Computational Methods for the Life Sciences

Math 3B: Lecture 14. Noah White. February 13, 2016

6 Continuous-Time Birth and Death Chains

Mathematics 22: Lecture 4

3 Single species models. Reading: Otto & Day (2007) section 4.2.1

Dynamical Systems. 1.0 Ordinary Differential Equations. 2.0 Dynamical Systems

1 Mathematical models of the population growth

Math 266: Autonomous equation and population dynamics

MATH3203 Lecture 1 Mathematical Modelling and ODEs

Dynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325

Lecture 1 Modeling in Biology: an introduction

The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor

Introduction to course: BSCI 462 of BIOL 708 R

2 One-dimensional models in discrete time

Population modeling of marine mammal populations

4 Insect outbreak model

Figure 5: Bifurcation diagram for equation 4 as a function of K. n(t) << 1 then substituting into f(n) we get (using Taylor s theorem)

Lecture 3. Dynamical Systems in Continuous Time

BIOL 410 Population and Community Ecology. Density-dependent growth 1

Transmission in finite populations

Difference Equations

Lecture 5: Modelling Population Dynamics N(t)

Math 142-2, Homework 2

Essential Ideas of Mathematical Modeling in Population Dynamics

SARAH P. OTTO and TROY DAY

Lecture 11: Logistic growth models #3

Assume closed population (no I or E). NB: why? Because it makes it easier.

1. The growth of a cancerous tumor can be modeled by the Gompertz Law: dn. = an ln ( )

Density Dependence, Stability & Allee effects ECOL 4000/6000

COMPETITION OF FAST AND SLOW MOVERS FOR RENEWABLE AND DIFFUSIVE RESOURCE

NOTES ON CHAOS. Chaotic dynamics stem from deterministic mechanisms, but they look very similar to random fluctuations in appearance.

Modeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool

A&S 320: Mathematical Modeling in Biology

population size at time t, then in continuous time this assumption translates into the equation for exponential growth dn dt = rn N(0)

Continuous-Time Dynamical Systems: Sketching solution curves and bifurcations

Autonomous Equations and Stability Sections

UNIVERSITY OF CALIFORNIA SANTA BARBARA DEPARTMENT OF CHEMICAL ENGINEERING. CHE 154: Engineering Approaches to Systems Biology Spring Quarter 2004

Objective. Single population growth models

Student Background Readings for Bug Dynamics: The Logistic Map

Reading: Most of the material in today s lecture comes from pages 1 6 of Chapter 1in

Qualitative analysis of differential equations: Part I

Final Project Descriptions Introduction to Mathematical Biology Professor: Paul J. Atzberger. Project I: Predator-Prey Equations

1.2. Introduction to Modeling

Some Generalizations of Logistic Model

Population Games and Evolutionary Dynamics

MODELS ONE ORDINARY DIFFERENTIAL EQUATION:

Calculus and Differential Equations II

Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations

An Introduction to Evolutionary Game Theory: Lecture 2

Analysis of bacterial population growth using extended logistic Growth model with distributed delay. Abstract INTRODUCTION

Derivation of Itô SDE and Relationship to ODE and CTMC Models

MA : Introductory Probability

PERSISTENCE AND EXTINCTION OF SINGLE POPULATION IN A POLLUTED ENVIRONMENT

Science One Math. October 23, 2018

Math Assignment 2

1.2. Introduction to Modeling. P (t) = r P (t) (b) When r > 0 this is the exponential growth equation.

MTH210 DIFFERENTIAL EQUATIONS. Dr. Gizem SEYHAN ÖZTEPE

Math 31S. Rumbos Fall Solutions to Exam 1

Scenarios for the transition to chaos

Review Problems for Exam 2

8. Qualitative analysis of autonomous equations on the line/population dynamics models, phase line, and stability of equilibrium points (corresponds

1 Simple one-dimensional dynamical systems birth/death and migration processes, logistic

Math 256: Applied Differential Equations: Final Review

History and meaning of the word Ecology A. Definition 1. Oikos, ology - the study of the house - the place we live

The logistic difference equation and the route to chaotic behaviour

perfect mixing paradox

3.5 Competition Models: Principle of Competitive Exclusion

1 This book is a translation of an amended and enlarged version of Part 1 of the following monograph

Stochastic Models. John M. Drake & Pejman Rohani

1.Introduction: 2. The Model. Key words: Prey, Predator, Seasonality, Stability, Bifurcations, Chaos.

Use separation of variables to solve the following differential equations with given initial conditions. y 1 1 y ). y(y 1) = 1

Mathematical Ecology

Introduction to SEIR Models

SIR Epidemic Model with total Population size

Computer Simulation and Applications in Life Sciences. Dr. Michael Emmerich & Dr. Andre Deutz LIACS

Introduction to stochastic multiscale modelling of tumour growth

Figure 1: Ca2+ wave in a Xenopus oocyte following fertilization. Time goes from top left to bottom right. From Fall et al., 2002.

An Application of Perturbation Methods in Evolutionary Ecology

Simulating stochastic epidemics

Mathematical Modelling Lecture 10 Difference Equations

Chapter 2: Growth & Decay

Lecture 1. Scott Pauls 1 3/28/07. Dartmouth College. Math 23, Spring Scott Pauls. Administrivia. Today s material.

From cell biology to Petri nets. Rainer Breitling, Groningen, NL David Gilbert, London, UK Monika Heiner, Cottbus, DE

Differential Equations ( DEs) MA 137 Calculus 1 with Life Science Application A First Look at Differential Equations (Section 4.1.

Multiple choice 2 pts each): x 2 = 18) Essay (pre-prepared) / 15 points. 19) Short Answer: / 2 points. 20) Short Answer / 5 points

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad

BIOS 3010: ECOLOGY. Dr Stephen Malcolm. Laboratory 6: Lotka-Volterra, the logistic. equation & Isle Royale

then the substitution z = ax + by + c reduces this equation to the separable one.

Dynamical Systems and Deep Learning: Overview. Abbas Edalat

Predator-Prey Population Dynamics

Transcription:

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 1: Mathematical Models MA 777 Spring 2018 1 / 15

Figure: from the Society for Industrial and Applied Mathematics (SIAM). David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 2 / 15

What is a Mathematical Model? A mathematical model is a description of a system using mathematical concepts and language. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 3 / 15

What is a Mathematical Model? A mathematical model is a description of a system using mathematical concepts and language. Modeling is a process that uses math to represent, analyze, make predictions, or otherwise provide insight into real world phenomena. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 3 / 15

What is a Mathematical Model? A mathematical model is a description of a system using mathematical concepts and language. Modeling is a process that uses math to represent, analyze, make predictions, or otherwise provide insight into real world phenomena. The process of developing a mathematical model is termed mathematical modeling. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 3 / 15

What is a Mathematical Model? A mathematical model is a description of a system using mathematical concepts and language. Modeling is a process that uses math to represent, analyze, make predictions, or otherwise provide insight into real world phenomena. The process of developing a mathematical model is termed mathematical modeling. Mathematical modeling is used in Physics Biology Economy etc. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 3 / 15

Types of Mathematical Models Many dynamic modeling approaches have been used over the last 6-7 decades for modeling biological systems. As a result, a large variety of models exists today. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 4 / 15

Types of Mathematical Models Many dynamic modeling approaches have been used over the last 6-7 decades for modeling biological systems. As a result, a large variety of models exists today. Generally, dynamic models can be classified according to how the time and the population of gene products are treated. Thus, there exist models based on continuous populations and continuous time such as systems of ordinary differential equations, David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 4 / 15

Types of Mathematical Models Many dynamic modeling approaches have been used over the last 6-7 decades for modeling biological systems. As a result, a large variety of models exists today. Generally, dynamic models can be classified according to how the time and the population of gene products are treated. Thus, there exist models based on continuous populations and continuous time such as systems of ordinary differential equations, models based on discrete populations and continuous time such as models based on the Gillespie formulation and their generalizations, and David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 4 / 15

Types of Mathematical Models Many dynamic modeling approaches have been used over the last 6-7 decades for modeling biological systems. As a result, a large variety of models exists today. Generally, dynamic models can be classified according to how the time and the population of gene products are treated. Thus, there exist models based on continuous populations and continuous time such as systems of ordinary differential equations, models based on discrete populations and continuous time such as models based on the Gillespie formulation and their generalizations, and models based on discrete population and discrete time frameworks such as Boolean Networks (BNs) and their stochastic variants. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 4 / 15

Other types of models 1 Statistical models. 2 Agent based models. 3 Animal models. 4 Logical models. 5 Petri nets. 6 PDEs and stochastic PDEs. 7 etc. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 5 / 15

Deterministic vs Stochastic { 1. All parameters are fixed, Deterministic 2. Dynamics is fixed. 1. Even when all parameters are know, the system may follow different trajectories. Stochastic 2. Although the final state might not be determined exactly, there are some parameters that govern the outputs. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 6 / 15

Variability during phage lambda infection Figure: A. Arkin, J. Ross, and H. H. McAdams. Genetics, 149(4), 1998. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 7 / 15

Modeling framework trade-off Desirable Requirements Accuracy Feasibility Fine model high high low Coarse model low "medium" high Requirements Accuracy Feasibility Model medium high high David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 8 / 15

Simulating Dynamics of Chemically Reacting Systems Data requirement Master Equa+on Intractable either numerically or analy+cally. Accurate but very slow. Con+nuous Gillespie Algorithm Sta+s+cally correct but s+ll slow Simple procedure Discrete Stochas+c ODE Determinis+c plus a stochas+c part Kine+cs needed Stochas+c DDS Stochas+c No kine+cs needed but some probability parameters are needed ODEs Determinis+c Kine+cs needed Discrete Dynamical Systems Determinis+c No kine+cs needed David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 9 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. Bifurcation analysis. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. Bifurcation analysis. Sensitivity analysis. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. Bifurcation analysis. Sensitivity analysis. Poincare-Bendixson theorem. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. Bifurcation analysis. Sensitivity analysis. Poincare-Bendixson theorem. Gillespie simulations. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. Bifurcation analysis. Sensitivity analysis. Poincare-Bendixson theorem. Gillespie simulations. Boolean networks. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. Bifurcation analysis. Sensitivity analysis. Poincare-Bendixson theorem. Gillespie simulations. Boolean networks. Markov chains. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. Bifurcation analysis. Sensitivity analysis. Poincare-Bendixson theorem. Gillespie simulations. Boolean networks. Markov chains. Perron-Frobenius theorem. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

How much you know about math modeling? Have you seen the following topics before? Dynamical systems. Bifurcation analysis. Sensitivity analysis. Poincare-Bendixson theorem. Gillespie simulations. Boolean networks. Markov chains. Perron-Frobenius theorem. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 10 / 15

Modeling population growth for single species Let N(t) be the population of a species at time t, then the rate of change can be modeled by dn dt = births deaths + migration David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 11 / 15

Modeling population growth for single species Let N(t) be the population of a species at time t, then the rate of change can be modeled by dn dt = births deaths + migration Example (Exponential growth) The simplest model has no migration and the birth and death are proportional to N, dn dt = bn dn N(t) = N 0 e (b d)t where b and d are positive constant and the initial population N(0) = N 0. This approach, due to Malthus in 1798, is fairly unrealistic. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 11 / 15

Modeling population growth for single species Example (Logistic growth) Verhulst (1838-1845) proposed that a self-limiting process should operate when a population becomes too large. He proposed ( dn = rn 1 N ) ( 1 dn or = r 1 N ) dt K N dt K where r and K are positive constants. In this model, the per capita growth rate decreases linearly with population size, r(1 N/K ). David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 12 / 15

Modeling population growth for single species Example (Logistic growth) Verhulst (1838-1845) proposed that a self-limiting process should operate when a population becomes too large. He proposed ( dn = rn 1 N ) ( 1 dn or = r 1 N ) dt K N dt K where r and K are positive constants. In this model, the per capita growth rate decreases linearly with population size, r(1 N/K ). If N(0) = N 0, the solution is N(t) = N 0 Ke rt K + N 0 (e rt 1) K as t The constant K is the carrying capacity. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 12 / 15

Modeling population growth for single species Figure: Logistic growth. Notice the qualitative difference when N 0 < K /2 and K /2 < N 0 < K. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 13 / 15

Modeling population growth for single species Example In general, we will consider dn dt = f (N) where f (N) is a nonlinear function. The equilibrium points N are the solutions of f (N) = 0. N is locally stable if f (N ) < 0. N is unstable if f (N ) > 0. Proof. Let N(t) = N + n(t). Show that n(t) = n(0)e rt, where r = f (N ). David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 14 / 15

The year of Mathematical Biology 2018 is a joint venture of the European Mathematical Society (EMS) and the European Society for Mathematical and Theoretical Biology (ESMTB). For more information, see http://euro-math-soc.eu/year-mathematical-biology-2018. David Murrugarra (University of Kentucky) Lecture 1: Mathematical Models MA 777 Spring 2018 15 / 15