TransWikia.com

On stochastic Lotka-Volterra predator-prey model

Biology Asked by Chaos on January 4, 2021

I am currently working (as a mathematician) on some estimations involving an stochastic predator-prey type model in which some of the coefficients have been perturbed by a Brownian Motion yielding to a system of stochastic differential equations.

The basic formulation of this problem can be found for instance in this article.

We have the following model

begin{equation}label{0}
begin{cases}
dX(t)=X(t)left(a-bX(t)-frac{sY(t)}{beta+Y(t)}right)dt
X(0)=x>0,
dY(t)=Y(t)left(frac{h X(t)}{beta+Y(t)}-c-fY(t)right)dt
Y(0)=y>0.
end{cases}
end{equation}

where we use Holling II response function.

In the literature (mostly mathematics literature) the authors propose to perturb the birth rate of the preys $amapsto a+dot{B}_1(t)$ and the mortality rate of predators $cmapsto c+dot{B}_2(t)$.

This yields to the following (stochastic) system

begin{equation}label{1}
begin{cases}
dX(t)=X(t)left(a-bX(t)-frac{sY(t)}{beta+Y(t)}right)dt+sigma_1 X(t)dB_1(t),
X(0)=x>0,
dY(t)=Y(t)left(frac{h X(t)}{beta+Y(t)}-c-fY(t)right)dt+sigma_2 Y(t)dB_2(t),
Y(0)=y>0
end{cases}
end{equation}

I know nothing about biology but I am concerned about a couple of things regarding this particular formulation:

  1. If we assume that the Brownian motions are perturbing the birth/mortality rates, how can I interpret the fact that this new “perturbed” rates can become negative (due to the effect of the BM). Does it have any sense (in the framework of this particular model) to talk about “negative birth rate” or “negative mortality rate”?
  2. What if the “perturbed” mortality rate goes beyond $100 %$? Mathematically nothing “wrong” will happen, but what about the interpretation?

Notice that the “noise” and the model as a whole can be interpreted differently (ignoring the fact that we are actually modelling the dynamics of two species) but my main issue is that if I(together with many authors) am stating that we perturb a certain parameter I believe that we must respect the basic assumptions of the model!

I hope everything’s clear and I thank you all in advance, any opinion or suggestion will be welcome!

One Answer

  1. To interpret $a$ and $c$ as a birth and mortality rates is somewhat inaccurate, as $a$ and $c$, as written are density-independent rates that grow or shrink populations exponentially. If $a$ were negative, the prey population would shrink and most people would find a problem like that unexciting unless it was, say, a sink population maintained by immigration from another population. This is true of $c$ being negative, where the predator population would grow independent of the prey population. More accurately, $a$ can be interpreted as the difference between birth, $b$, and deaths, $delta$; i.e., $a = b - delta$. In this case, depending on the kinds of stochasticity being added to the model, negative $a$ is not just normal, but expected. The same is true for $c$. Ultimately, a perturbations like the white noise in this equation is simultaneously perturbing the difference between birth and death rates, which explains the negative difference.
  2. The value of the rates seem to have the same interpretation if they are constant or are perturbed to an arbitrary magnitude.

I hope this helps!

Correct answer by Chris Moore on January 4, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP