Birth process

A birth process with birth rates ${\lambda }_0,{\lambda }_1,{\lambda }_2,...$.

In probability theory, a birth process or a pure birth processTemplate:Sfnp is a special case of a continuous-time Markov process and a generalisation of a Poisson process. It defines a continuous process which takes values in the natural numbers and can only increase by one (a "birth") or remain unchanged. This is a type of birth–death process with no deaths. The rate at which births occur is given by an exponential random variable whose parameter depends only on the current value of the process

Definition

Birth rates definition

A birth process with birth rates $({\lambda }_n,n\in \mathbb{N})$ and initial value $k\in \mathbb{N}$ is a minimal right-continuous process $(X_t,t\ge 0)$ such that $X_0=k$ and the interarrival times $T_i=\inf \{t\ge 0:X_t=i+1\}-\inf \{t\ge 0:X_t=i\}$ are independent exponential random variables with parameter ${\lambda }_i$.Template:Sfnp

Infinitesimal definition

A birth process with rates $({\lambda }_n,n\in \mathbb{N})$ and initial value $k\in \mathbb{N}$ is a process $(X_t,t\ge 0)$ such that:

• $X_0=k$
• $\mathrm{\forall }s,t\ge 0:s<t⟹X_s\le X_t$
• $\mathbb{P}(X_{t+h}=X_t+1)={\lambda }_{X_t}h+o(h)$
• $\mathbb{P}(X_{t+h}=X_t)=o(h)$
• $\mathrm{\forall }s,t\ge 0:s<t⟹X_t-X_s$ is independent of $(X_u,u<s)$

(The third and fourth conditions use little o notation.)

These conditions ensure that the process starts at $i$, is non-decreasing and has independent single births continuously at rate ${\lambda }_n$, when the process has value $n$.Template:Sfnp

Continuous-time Markov chain definition

A birth process can be defined as a continuous-time Markov process (CTMC) $(X_t,t\ge 0)$ with the non-zero Q-matrix entries $q_{n,n+1}={\lambda }_n=-q_{n,n}$ and initial distribution $i$ (the random variable which takes value $i$ with probability 1).Template:Sfnp

$Q=\left(\begin{array}{ccccc}-{\lambda }_0& {\lambda }_0& 0& 0& \cdots \\ 0& -{\lambda }_1& {\lambda }_1& 0& \cdots \\ 0& 0& -{\lambda }_2& {\lambda }_2& \cdots \\ ⋮& ⋮& ⋮& & ⋮\ddots \end{array}\right)$

Variations

Some authors require that a birth process start from 0 i.e. that $X_0=0$,Template:Sfnp while others allow the initial value to be given by a probability distribution on the natural numbers.Template:Sfnp The state space can include infinity, in the case of an explosive birth process.Template:Sfnp The birth rates are also called intensities.Template:Sfnp

Properties

As for CTMCs, a birth process has the Markov property. The CTMC definitions for communicating classes, irreducibility and so on apply to birth processes. By the conditions for recurrence and transience of a birth–death process,Template:Sfnp any birth process is transient. The transition matrices $((p_{i,j}(t){)}_{i,j\in \mathbb{N}}),t\ge 0)$ of a birth process satisfy the Kolmogorov forward and backward equations.

The backwards equations are:Template:Sfnp

$p{\text{'}}_{i,j}(t)={\lambda }_i(p_{i+1,j}(t)-p_{i,j}(t))$ (for $i,j\in \mathbb{N}$)

The forward equations are:Template:Sfnp

$p{\text{'}}_{i,i}(t)=-{\lambda }_i{p}_{i,i}(t)$ (for $i\in \mathbb{N}$)

$p{\text{'}}_{i,j}(t)={\lambda }_{j-1}{p}_{i,j-1}(t)-{\lambda }_j{p}_{i,j}(t)$ (for $j\ge i+1$)

From the forward equations it follows that:Template:Sfnp

$p_{i,i}(t)=e^{-{\lambda }_{i}t}$ (for $i\in \mathbb{N}$)

$p_{i,j}(t)={\lambda }_{j-1}{e}^{-{\lambda }_{j}t}\underset{0}{\overset{t}{\int }}{e}^{{\lambda }_{j}s}{p}_{i,j-1}(s)\text{d}s$ (for $j\ge i+1$)

Unlike a Poisson process, a birth process may have infinitely many births in a finite amount of time. We define $T_{\mathrm{\infty }}=\sup \{T_n:n\in \mathbb{N}\}$ and say that a birth process explodes if $T_{\mathrm{\infty }}$ is finite. If $\sum _{n=0}^{\mathrm{\infty }}\frac{1}{{\lambda }_n}<\mathrm{\infty }$ then the process is explosive with probability 1; otherwise, it is non-explosive with probability 1 ("honest").Template:SfnpTemplate:Sfnp

Examples

A Poisson process is a special case of a birth process.

A Poisson process is a birth process where the birth rates are constant i.e. ${\lambda }_n=\lambda$ for some $\lambda >0$.Template:Sfnp

Simple birth process

A simple birth process, where birth rates are equal to the size of the current population.

A simple birth process is a birth process with rates ${\lambda }_n=n\lambda$.Template:Sfnp It models a population in which each individual gives birth repeatedly and independently at rate $\lambda$. Udny Yule studied the processes, so they may be known as Yule processes.Template:Sfnp

The number of births in time $t$ from a simple birth process of population $n$ is given by:Template:Sfnp

$p_{n,n+m}(t)=(\genfrac{}{}{0ex}{}{n}{m})(\lambda t{)}^m(1-\lambda t{)}^{n-m}+o(h)$

In exact form, the number of births is the negative binomial distribution with parameters $n$ and $e^{-\lambda t}$. For the special case $n=1$, this is the geometric distribution with success rate $e^{-\lambda t}$.Template:Sfnp

The expectation of the process grows exponentially; specifically, if $X_0=1$ then $𝔼(X_t)=e^{\lambda t}$.Template:Sfnp

A simple birth process with immigration is a modification of this process with rates ${\lambda }_n=n\lambda +\nu$. This models a population with births by each population member in addition to a constant rate of immigration into the system.Template:Sfnp

References

• Template:Cite book
• Template:Cite journal
• Template:Cite book
• Template:Cite book
• Template:Cite book