Birth–death process

From Wikipedia, the free encyclopedia

The birth–death process (or birth-and-death process) is a special case of continuous-time Markov process where the state transitions are of only two types: "births", which increase the state variable by one and "deaths", which decrease the state by one. The model's name comes from a common application, the use of such models to represent the current size of a population where the transitions are literal births and deaths. Birth–death processes have many applications in demography, queueing theory, performance engineering, epidemiology, biology and other areas. They may be used, for example, to study the evolution of bacteria, the number of people with a disease within a population, or the number of customers in line at the supermarket.

When a birth occurs, the process goes from state n to n + 1. When a death occurs, the process goes from state n to state n − 1. The process is specified by birth rates and death rates .

State diagram of a birth-death process

Recurrence and transience[]

For recurrence and transience in Markov processes see Section 5.3 from Markov chain.

Conditions for recurrence and transience[]

Conditions for recurrence and transience were established by Samuel Karlin and James McGregor.[1]

A birth-and-death process is recurrent if and only if
A birth-and-death process is ergodic if and only if
A birth-and-death process is null-recurrent if and only if

By using Extended Bertrand's test (see Section 4.1.4 from Ratio test) the conditions for recurrence, transience, ergodicity and null-recurrence can be derived in a more explicit form.[2]

For integer let denote the th iterate of natural logarithm, i.e. and for any , .

Then, the conditions for recurrence and transience of a birth-and-death process are as follows.

The birth-and-death process is transient if there exist and such that for all

where the empty sum for is assumed to be 0.

The birth-and-death process is recurrent if there exist and such that for all

Application[]

Consider one-dimensional random walk that is defined as follows. Let , and where takes values , and the distribution of is defined by the following conditions:

where satisfy the condition .

The random walk described here is a discrete time analogue of the birth-and-death process (see Markov chain) with the birth rates

and the death rates

.

So, recurrence or transience of the random walk is associated with recurrence or transience of the birth-and-death process.[2]

The random walk is transient if there exist , and such that for all

where the empty sum for is assumed to be zero.

The random walk is recurrent if there exist and such that for all

Stationary solution[]

If a birth-and-death process is ergodic, then there exists steady-state probabilities where is the probability that the birth-and-death process is in state at time The limit exists, independent of the initial values and is calculated by the relations:

These limiting probabilities are obtained from the infinite system of differential equations for

and the initial condition

In turn, the last system of differential equations is derived from the system of that describes the dynamic of the system in a small time . During this small time only three types of transitions are considered as one death, or one birth, or no birth nor death. The probability of the first two of these transitions has the order of . Other transitions during this small interval such as more than one birth, or more than one death, or at least one birth and at least one death have the probabilities that are of smaller order than , and hence are negligible in derivations. If the system is in state k, then the probability of birth during an interval is , the probability of death is , and the probability of no birth and no death is . For a population process, "birth" is the transition towards increasing the population size by 1 while "death" is the transition towards decreasing the population size by 1.

Examples of birth-death processes[]

A pure birth process is a birth–death process where for all .

A pure death process is a birth–death process where for all .

M/M/1 model and M/M/c model, both used in queueing theory, are birth–death processes used to describe customers in an infinite queue.

Use in phylodynamics[]

Birth-death processes are in phylodynamics as a generative model of phylogenies, ie a binary tree in which birth events correspond to branches of the tree and death events correspond to leaf nodes.[3] Notably, they are used in viral phylodynamics[4] to understand the transmission process.

Use in queueing theory[]

In queueing theory the birth–death process is the most fundamental example of a queueing model, the M/M/C/K//FIFO (in complete Kendall's notation) queue. This is a queue with Poisson arrivals, drawn from an infinite population, and C servers with exponentially distributed service times with K places in the queue. Despite the assumption of an infinite population this model is a good model for various telecommunication systems.

M/M/1 queue[]

The M/M/1 is a single server queue with an infinite buffer size. In a non-random environment the birth–death process in queueing models tend to be long-term averages, so the average rate of arrival is given as and the average service time as . The birth and death process is an M/M/1 queue when,

The differential equations for the probability that the system is in state k at time t are

Pure birth process associated with an M/M/1 queue[]

Pure birth process with is a particular case of the M/M/1 queueing process. We have the following system of differential equations:

Under the initial condition and , the solution of the system is

That is, a (homogeneous) Poisson process is a pure birth process.

M/M/c queue[]

The M/M/C is a multi-server queue with C servers and an infinite buffer. It characterizes by the following birth and death parameters:

and

with

The system of differential equations in this case has the form:

Pure death process associated with an M/M/C queue[]

Pure death process with is a particular case of the M/M/C queueing process. We have the following system of differential equations:

Under the initial condition and we obtain the solution

that presents the version of binomial distribution depending on time parameter (see Binomial process).

M/M/1/K queue[]

The M/M/1/K queue is a single server queue with a buffer of size K. This queue has applications in telecommunications, as well as in biology when a population has a capacity limit. In telecommunication we again use the parameters from the M/M/1 queue with,

In biology, particularly the growth of bacteria, when the population is zero there is no ability to grow so,

Additionally if the capacity represents a limit where the individual dies from over population,

The differential equations for the probability that the system is in state k at time t are

Equilibrium[]

A queue is said to be in equilibrium if the steady state probabilities exist. The condition for the existence of these steady-state probabilities in the case of M/M/1 queue is and in the case of M/M/C queue is . The parameter is usually called load parameter or utilization parameter. Sometimes it is also called traffic intensity.

Using the M/M/1 queue as an example, the steady state equations are

This can be reduced to

So, taking into account that , we obtain

See also[]

Notes[]

  1. ^ Karlin, Samuel; McGregor, James (1957). "The classification of birth and death processes" (PDF). Transactions of the American Mathematical Society. 86 (2): 366–400.
  2. ^ a b Abramov, Vyacheslav M. (2020). "Extension of the Bertrand–De Morgan test and its application". The American Mathematical Monthly. 127 (5): 444--448. arXiv:1901.05843. doi:10.1080/00029890.2020.1722551.
  3. ^ Stadler T (December 2010). "Sampling-through-time in birth-death trees". Journal of Theoretical Biology. 267 (3): 396–404. doi:10.1016/j.jtbi.2010.09.010. PMID 20851708.
  4. ^ Kühnert D, Wu CH, Drummond AJ (December 2011). "Phylogenetic and epidemic modeling of rapidly evolving infectious diseases". Infection, Genetics and Evolution. 11 (8): 1825–41. doi:10.1016/j.meegid.2011.08.005. PMC 7106223. PMID 21906695.

References[]

  • Latouche, G.; Ramaswami, V. (1999). "Quasi-Birth-and-Death Processes". Introduction to Matrix Analytic Methods in Stochastic Modelling (1st ed.). ASA SIAM. ISBN 0-89871-425-7.
  • Nowak, M. A. (2006). Evolutionary Dynamics: Exploring the Equations of Life. Harvard University Press. ISBN 0-674-02338-2.
  • Virtamo, J. "Birth-death processes" (PDF). 38.3143 Queueing Theory. Retrieved 2 December 2019.
Retrieved from ""