p548 lecture2v2

Information about p548 lecture2v2

Published on March 3, 2008

Author: Ulisse

Source: authorstream.com

Content

Population Modeling:  Population Modeling Mathematical Biology Lecture 2 James A. Glazier (Partially Based on Brittain Chapter 1) Population Models:  Population Models Simple and a Good Introduction to Methods Two Types Continuum [Britton, Chapter 1] Discrete [Britton, Chapter 2] Continuum Population Models:  Continuum Population Models Given a Population, N0 of an animal, cell, bacterium,… at time t=0, What is the Population N(t) at time t? Assume that the population is large so treat N as a continuous variable. Naively: Continuum Models are Generally More Stable than discrete models (no chaos or oscillations) Malthusian Model (Exponential Growth):  Malthusian Model (Exponential Growth) For a Fertility Rate, b, a Death Rate, d, and no Migration:  In reality have a saturation: limited food, disease, predation, reduced birth-rate from crowding… Density Dependent Effects:  Density Dependent Effects How to Introduce Density Dependent Effects? Decide on Essential Characteristics of Data. Write Simplest form of f(N) which Gives these Characteristics. Choose Model parameters to Fit Data Generally, Growth is Sigmoidal, i.e. small for small and large populations  f(0) = f(K) = 0, where K is the Carrying Capacity and f(N) has a unique maximum for some value of N, Nmax The Simplest Possible Solution is the Verhulst or Logistic Equation Verhulst or Logistic Equation:  Verhulst or Logistic Equation A Key Equation—Will Use Repeatedly Assume Death Rate  N, or that Birth Rate Declines with Increasing N, Reaching 0 at the Carrying Capacity, K: The Logistic Equation has a Closed-Form Solution: No Chaos in Continuum Logistic Equation Click for Solution Details Solving the Logistic Equation:  Solving the Logistic Equation 2) Now let: 3) Substitute: 1) Start with Logistic Equation: 4) Solve for N(t): General Issues in Modeling:  General Issues in Modeling Not a model unless we can explain why the death rate d~N/K. Can always improve fit using more parameters. Meaningless unless we can justify them. Logistic Map has only three parameters N0, K, r – doesn't fit real populations, very well. But we are not just curve fitting. Don't introduce parameters unless we know they describe a real mechanism in biology. Fitting changes in response to different parameters is much more useful than fitting a curve with a single set of parameters. Idea: Steady State or Fixed Point:  Idea: Steady State or Fixed Point For a Differential Equation of Form is a Fixed Point  So the Logistic Equation has Two Fixed Points, N=0 and N=K Fixed Points are also often designated x* Idea: Stability:  Idea: Stability Is the Fixed Point Stable? I.e. if you move a small distance e away from x0 does x(t) return to x0? If so x0 is Stable, if not, x0 is Unstable. Calculating Stability: Linear Stability Analysis:  Calculating Stability: Linear Stability Analysis Consider a Fixed Point x0 and a Perturbation e. Assume that: Taylor Expand f around x0: So, if Response Timescale, t, for disturbance to grow or shrink by a factor of e is: Example Logistic Equation:  Example Logistic Equation Start with the Logistic Eqn. Fixed Points at N0=0 and N0=K For N0=0 For N0=K stable unstable Phase Portraits :  Phase Portraits Idea: Describe Stability Behavior Graphically Arrows show direction of Flow Generally: Solution of the Logistic Equation:  Solution of the Logistic Equation Solution of the Logistic Equation: For N0>K, N(t) decreases exponentially to K For N0<2, K/N(t) increases sigmoidally to K For K/2<N0<K, N(t) increases exponentially to K Example Stability in Population Competition:  Example Stability in Population Competition Consider two species, N1 and N2, with growth rates r1 and r2 and carrying capacities K1 and K2, competing for the same resource. Both obey Logistic Equation. If one species has both bigger carrying capacity and faster growth rate, it will displace the other. What if one species has faster growth rate and the other a greater carrying capacity? An example of a serious evolutionary/ecological question answerable with simple mathematics. Population Competition—Contd.:  Population Competition—Contd. Start with all N1 and no N2. Represent population as a vector (N1, N2) Steady state is (K1,0). What if we introduce a few N2? In two dimensions we need to look at the eigenvalues of the Jacobian Matrix evaluated at the fixed point. Evaluate at (K1,0). Stability in Two Dimensions:  Stability in Two Dimensions 1) Both Eigenvalues Positive—Unstable Cases: 2) One Eigenvalue Positive, One Negative—Unstable 1) Both Eigenvalues Negative—Stable Population Competition—Contd.:  Population Competition—Contd. -r1 always < 0 so fixed point is stable  r2(1-K1/K2)<0 i.e. if K1>K2. Fixed Point Unstable (i.e. species 2 Invades Successfully)  K2>K1 Independent of r2! So high carrying capacity wins out over high fertility (called K-selection in evolutionary biology). A surprising result. The opposite of what is generally observed in nature. Eigenvalues are solutions of

Related presentations


Other presentations created by Ulisse

wtoprimitivo
08. 05. 2008
0 views

wtoprimitivo

M Naoi IPES2007
07. 05. 2008
0 views

M Naoi IPES2007

atos maxwell jamescfh
02. 05. 2008
0 views

atos maxwell jamescfh

06ESADE
30. 04. 2008
0 views

06ESADE

1 4 Eiropas Kopienas programmas
22. 04. 2008
0 views

1 4 Eiropas Kopienas programmas

5 11
21. 04. 2008
0 views

5 11

phan
18. 04. 2008
0 views

phan

Brown Velocity
17. 04. 2008
0 views

Brown Velocity

SME in Pakistan
10. 01. 2008
0 views

SME in Pakistan

visual perception 1
14. 01. 2008
0 views

visual perception 1

Ross Thompson Presentation
15. 01. 2008
0 views

Ross Thompson Presentation

ArmorofGod ShieldofFaithNew
17. 01. 2008
0 views

ArmorofGod ShieldofFaithNew

Hospital Employee Immunization 1
17. 01. 2008
0 views

Hospital Employee Immunization 1

EnhancingCoreWeb
10. 01. 2008
0 views

EnhancingCoreWeb

A C 8 US CHAPTER
10. 01. 2008
0 views

A C 8 US CHAPTER

6 SMS
18. 01. 2008
0 views

6 SMS

textiles grad 2006
22. 01. 2008
0 views

textiles grad 2006

SessIIBYeboahEtAl
22. 01. 2008
0 views

SessIIBYeboahEtAl

hfcourse frzn
19. 01. 2008
0 views

hfcourse frzn

Heezen overview of cases
22. 01. 2008
0 views

Heezen overview of cases

2c
23. 01. 2008
0 views

2c

language
14. 01. 2008
0 views

language

10 Cash Control
04. 02. 2008
0 views

10 Cash Control

gsoverview
04. 02. 2008
0 views

gsoverview

Cuisine
22. 01. 2008
0 views

Cuisine

schap011
07. 02. 2008
0 views

schap011

Frontier1
11. 02. 2008
0 views

Frontier1

resiliency
09. 01. 2008
0 views

resiliency

library orientation
28. 01. 2008
0 views

library orientation

WeddingPlanning
29. 01. 2008
0 views

WeddingPlanning

harpring
21. 01. 2008
0 views

harpring

EnabGrid
07. 02. 2008
0 views

EnabGrid

owen holland
12. 01. 2008
0 views

owen holland

breakout 2ARavishankara
20. 02. 2008
0 views

breakout 2ARavishankara

Michelozzo and Bernardo
10. 01. 2008
0 views

Michelozzo and Bernardo

absent
15. 01. 2008
0 views

absent

p1 calyon de dumast
28. 02. 2008
0 views

p1 calyon de dumast

SB 18 Productive Consultation
08. 03. 2008
0 views

SB 18 Productive Consultation

tortola
20. 01. 2008
0 views

tortola

what
20. 01. 2008
0 views

what

ONEPowerPointPresent ation
19. 03. 2008
0 views

ONEPowerPointPresent ation

hall peter
25. 03. 2008
0 views

hall peter

EID Intro
31. 03. 2008
0 views

EID Intro

cottrell icfa korea may05
03. 04. 2008
0 views

cottrell icfa korea may05

Maps
07. 04. 2008
0 views

Maps

ErgonomicsPresentati on 000
07. 03. 2008
0 views

ErgonomicsPresentati on 000

London2012
14. 04. 2008
0 views

London2012

Samurai Hoedown
07. 02. 2008
0 views

Samurai Hoedown

m2f2
04. 02. 2008
0 views

m2f2

Olympic Mascots
14. 04. 2008
0 views

Olympic Mascots

stdslides2003
12. 03. 2008
0 views

stdslides2003

Packaging ru
14. 01. 2008
0 views

Packaging ru

phpulsefeb
14. 02. 2008
0 views

phpulsefeb

ALADI May2 2005
25. 01. 2008
0 views

ALADI May2 2005

041014 TEAHA
10. 04. 2008
0 views

041014 TEAHA

tds gap
22. 01. 2008
0 views

tds gap

AIDS f06
16. 01. 2008
0 views

AIDS f06

UniversityForumFeb23 2005
16. 01. 2008
0 views

UniversityForumFeb23 2005

Presentation to LGA 23Jan07 V3
24. 01. 2008
0 views

Presentation to LGA 23Jan07 V3

NW web presentation7 03
29. 02. 2008
0 views

NW web presentation7 03

build ssp system
21. 01. 2008
0 views

build ssp system

laighside
29. 01. 2008
0 views

laighside

landis
15. 01. 2008
0 views

landis

Service design for active ageing
21. 11. 2016
0 views

Service design for active ageing