Itô versus Stratonovich calculus in random population growth

Math Biosci. 2007 Mar;206(1):81-107. doi: 10.1016/j.mbs.2004.09.002. Epub 2005 Oct 6.

Abstract

The context is the general stochastic differential equation (SDE) model dN/dt=N(g(N)+sigmaepsilon(t)) for population growth in a randomly fluctuating environment. Here, N=N(t) is the population size at time t, g(N) is the 'average' per capita growth rate (we work with a general almost arbitrary function g), and sigmaepsilon(t) is the effect of environmental fluctuations (sigma>0, epsilon(t) standard white noise). There are two main stochastic calculus used to interpret the SDE, Itô calculus and Stratonovich calculus. They yield different solutions and even qualitatively different predictions (on extinction, for example). So, there is a controversy on which calculus one should use. We will resolve the controversy and show that the real issue is merely semantic. It is due to the informal interpretation of g(x) as being an (unspecified) 'average' per capita growth rate (when population size is x). The implicit assumption usually made in the literature is that the 'average' growth rate is the same for both calculi, when indeed this rate should be defined in terms of the observed process. We prove that, when using Itô calculus, g(N) is indeed the arithmetic average growth rate R(a)(x) and, when using Stratonovich calculus, g(N) is indeed the geometric average growth rate R(g)(x). Writing the solutions of the SDE in terms of a well-defined average, R(a)(x) or R(g)(x), instead of an undefined 'average' g(x), we prove that the two calculi yield exactly the same solution. The apparent difference was due to the semantic confusion of taking the informal term 'average growth rate' as meaning the same average.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Algorithms*
  • Animals
  • Biomass
  • Extinction, Biological
  • Models, Biological*
  • Population Dynamics
  • Population Growth
  • Probability
  • Stochastic Processes
  • Time Factors