ConceptComplete

The Poisson Process

The Poisson process is the most fundamental continuous-time stochastic process, modeling the occurrence of random events over time. It is characterized by independent, stationary increments and plays a central role in queueing theory, insurance mathematics, and the theory of point processes.


Definition and basic properties

Definition2.1Poisson process

A stochastic process (Nt)t0(N_t)_{t \geq 0} taking values in {0,1,2,}\{0, 1, 2, \ldots\} is a Poisson process with rate (or intensity) λ>0\lambda > 0 if:

  1. N0=0N_0 = 0.
  2. Independent increments: For any 0t1<t2<<tn0 \leq t_1 < t_2 < \cdots < t_n, the increments Nt2Nt1,Nt3Nt2,,NtnNtn1N_{t_2} - N_{t_1}, N_{t_3} - N_{t_2}, \ldots, N_{t_n} - N_{t_{n-1}} are independent.
  3. Stationary increments: For any s,t0s, t \geq 0, the distribution of Nt+sNsN_{t+s} - N_s depends only on tt, not on ss.
  4. Poisson distribution: For all t0t \geq 0, NtPoisson(λt)N_t \sim \text{Poisson}(\lambda t), i.e.,

P(Nt=k)=eλt(λt)kk!,k=0,1,2,\mathbb{P}(N_t = k) = e^{-\lambda t} \frac{(\lambda t)^k}{k!}, \quad k = 0, 1, 2, \ldots

Intuitively, NtN_t counts the number of "arrivals" (or "events") in the time interval [0,t][0, t]. The rate λ\lambda is the expected number of arrivals per unit time: E[Nt]=λt\mathbb{E}[N_t] = \lambda t.

RemarkEquivalent characterizations

The Poisson process can also be defined via:

  • Infinitesimal properties: For small h>0h > 0:

    • P(Nh=1)=λh+o(h)\mathbb{P}(N_h = 1) = \lambda h + o(h),
    • P(Nh2)=o(h)\mathbb{P}(N_h \geq 2) = o(h),
    • P(Nh=0)=1λh+o(h)\mathbb{P}(N_h = 0) = 1 - \lambda h + o(h).
  • Inter-arrival times: The times between successive events are i.i.d. exponential random variables with rate λ\lambda.

These characterizations are all equivalent and provide different intuitions for the Poisson process.


Inter-arrival times

Theorem2.1Inter-arrival times are exponential

Let T1,T2,T_1, T_2, \ldots be the times of the successive arrivals in a Poisson process with rate λ\lambda, and let τn=TnTn1\tau_n = T_n - T_{n-1} (with T0=0T_0 = 0) be the inter-arrival times. Then τ1,τ2,\tau_1, \tau_2, \ldots are i.i.d. Exp(λ)\text{Exp}(\lambda) random variables.

Conversely, if τ1,τ2,\tau_1, \tau_2, \ldots are i.i.d. Exp(λ)\text{Exp}(\lambda), and we define Tn=k=1nτkT_n = \sum_{k=1}^n \tau_k and Nt=max{n:Tnt}N_t = \max\{n : T_n \leq t\}, then (Nt)(N_t) is a Poisson process with rate λ\lambda.

ExampleBus arrivals

Suppose buses arrive at a stop according to a Poisson process with rate λ=5\lambda = 5 per hour. Then:

  • The expected number of buses in 2 hours is E[N2]=52=10\mathbb{E}[N_2] = 5 \cdot 2 = 10.
  • The probability of exactly 3 buses arriving in 30 minutes is

P(N0.5=3)=e2.5(2.5)33!0.2138.\mathbb{P}(N_{0.5} = 3) = e^{-2.5} \frac{(2.5)^3}{3!} \approx 0.2138.

  • The time until the first bus arrives is Exp(5)\text{Exp}(5), with mean 1/51/5 hour = 12 minutes.

Superposition and thinning

Theorem2.2Superposition

If (Nt(1))(N_t^{(1)}) and (Nt(2))(N_t^{(2)}) are independent Poisson processes with rates λ1\lambda_1 and λ2\lambda_2, then their superposition

Nt=Nt(1)+Nt(2)N_t = N_t^{(1)} + N_t^{(2)}

is a Poisson process with rate λ1+λ2\lambda_1 + \lambda_2.

This is the additive property of Poisson processes: merging two independent Poisson streams yields another Poisson stream.

Theorem2.3Thinning

If (Nt)(N_t) is a Poisson process with rate λ\lambda, and each arrival is independently retained with probability pp (and discarded with probability 1p1-p), then the retained arrivals form a Poisson process with rate pλp\lambda.

Thinning is the converse of superposition: splitting a Poisson stream into independent substreams.

ExampleCustomer types in a queueing system

Customers arrive at a service center according to a Poisson process with rate λ=10\lambda = 10 per hour. Each customer is type A with probability 0.60.6 and type B with probability 0.40.4, independently. Then:

  • Type A arrivals form a Poisson process with rate 100.6=610 \cdot 0.6 = 6 per hour.
  • Type B arrivals form a Poisson process with rate 100.4=410 \cdot 0.4 = 4 per hour.
  • The two streams are independent.

Conditional properties

Theorem2.4Uniform order statistics

Conditionally on Nt=nN_t = n, the arrival times T1,,TnT_1, \ldots, T_n are distributed as the order statistics of nn i.i.d. Uniform(0,t)\text{Uniform}(0, t) random variables.

This beautiful result says that, given nn arrivals in [0,t][0, t], the arrival times are "uniformly scattered" throughout the interval.

ExampleConditional probability

Suppose N1=2N_1 = 2 (two arrivals in [0,1][0,1]). What is the probability that both arrivals occurred in the first half [0,0.5][0, 0.5]?

By the uniform order statistics property, each arrival independently falls in [0,0.5][0, 0.5] with probability 0.50.5, so

P(T1,T20.5N1=2)=(0.5)2=0.25.\mathbb{P}(T_1, T_2 \leq 0.5 \mid N_1 = 2) = (0.5)^2 = 0.25.


Compound Poisson processes

Definition2.2Compound Poisson process

A compound Poisson process is a process of the form

Xt=i=1NtYi,X_t = \sum_{i=1}^{N_t} Y_i,

where (Nt)(N_t) is a Poisson process with rate λ\lambda and {Yi}\{Y_i\} are i.i.d. random variables (independent of NtN_t), called the jump sizes.

If YiY_i represents the "size" or "value" of the ii-th arrival, then XtX_t is the cumulative value up to time tt.

ExampleInsurance claims

An insurance company receives claims according to a Poisson process with rate λ=20\lambda = 20 per year. Each claim amount is i.i.d. with mean μ=1000\mu = 1000 and variance σ2=500,000\sigma^2 = 500{,}000. The total claims up to time tt is

Xt=i=1NtYi.X_t = \sum_{i=1}^{N_t} Y_i.

Then:

  • E[Xt]=E[Nt]E[Y1]=20t1000=20,000t\mathbb{E}[X_t] = \mathbb{E}[N_t] \mathbb{E}[Y_1] = 20t \cdot 1000 = 20{,}000 t.
  • Var(Xt)=E[Nt]E[Y12]=20t(1,000,000+500,000)=30,000,000t\text{Var}(X_t) = \mathbb{E}[N_t] \mathbb{E}[Y_1^2] = 20t \cdot (1{,}000{,}000 + 500{,}000) = 30{,}000{,}000 t.

Waiting time paradox

ExampleWaiting time paradox (renewal paradox)

You arrive at a bus stop at a random time. The buses arrive according to a Poisson process with rate λ\lambda. What is the expected time until the next bus?

Answer: 1/λ1/\lambda (the same as the mean inter-arrival time). This is a special feature of the memoryless property of the exponential distribution: the remaining time until the next arrival is always Exp(λ)\text{Exp}(\lambda), regardless of when you arrive.

In contrast, for a non-Poisson renewal process (e.g., deterministic inter-arrival times), the expected waiting time for a random arrival is typically greater than half the inter-arrival time. This is the inspection paradox or waiting time paradox.


Summary

The Poisson process is characterized by:

  • Independent, stationary increments: Memoryless evolution.
  • Exponential inter-arrival times: τiExp(λ)\tau_i \sim \text{Exp}(\lambda), i.i.d.
  • Superposition and thinning: Merging and splitting preserve the Poisson property.
  • Uniform order statistics: Conditional arrivals are uniformly distributed.

These properties make the Poisson process tractable for analysis and ubiquitous in applications: queueing, insurance, telecommunications, reliability theory, and more.