Part 8 - Poisson Processes II

Introduction

In this part we will introduce more properties of Poisson processes and Spatial Poisson processes.

Superposition and Thinning of Poisson Processes

Lemma 1 (Superposition of Poisson Processes)

Let {Nt(1)}t0,,{Nt(n)}t0\{N_t^{(1)}\}_{t \geq 0}, \ldots, \{N_t^{(n)}\}_{t \geq 0} be independent Poisson processes with parameters λp1,,λpn\lambda p_1, \ldots, \lambda p_n, respectively, where p=(p1,,pn)p = (p_1, \ldots, p_n) is a probability vector. If c=(c1,,cn)c = (c_1, \ldots, c_n) are the counts after time tt (i.e., ci=Nt(i)c_i = N_t^{(i)}), then, the conditional distribution of (Nt(1),,Nt(n))(N_t^{(1)}, \ldots, N_t^{(n)}), an equivalent model is,

cMultinomial(Nt,p),c \sim \mathrm{Multinomial}(N_t, p),

where {Nt}t0\{N_t\}_{t \geq 0} is a Poisson process with parameter λ\lambda.

Proof

Using the model with independent Poisson processes, the probability of observing the count vector cc after time tt is (denoting N=c1++cnN = c_1 + \ldots + c_n),

i=1nPoisson(ci;λpit)=i=1neλpit(λpit)cici!=eλt(λt)Ni=1npicici!=eλt(λt)NN!N!c1!c2!cn!p1c1p2c2pncn=Poisson(N;λt)Multinomial(c;N,p)\begin{align*} \prod_{i = 1}^n \mathrm{Poisson}(c_i; \lambda p_i t) & = \prod_{i = 1}^n e^{-\lambda p_i t} \frac{(\lambda p_i t)^{c_i}}{c_i!} \newline & = e^{-\lambda t} (\lambda t)^N \prod_{i = 1}^n \frac{p_i^{c_i}}{c_i!} \newline & = e^{-\lambda t} \frac{(\lambda t)^N}{N!} \cdot \frac{N!}{c_1! c_2! \ldots c_n!} p_1^{c_1} p_2^{c_2} \ldots p_n^{c_n} \newline & = \mathrm{Poisson}(N; \lambda t) \cdot \mathrm{Multinomial}(c; N, p) \newline \end{align*}

The process for NN inherits independent and stationary increments from the sub-processes, so it follows it also a Poisson process with parameter λ\lambda.

Uniformly Distributed Arrival Times

Lemma 2 (Uniformly Distributed Arrival Times)

Let {Nt}t0\{N_t\}_{t \geq 0} be a Poisson process with parameter λ\lambda. If we fix that Nt=kN_t = k, and we select uniformly randomly one of the kk arrivals, then its arrival time is uniformly distributed on the interval [0,t][0, t].

Proof

Let S1,S2,,SkS_1, S_2, \ldots, S_k be the arrival times given that Nt=kN_t = k.

P(Sksk uniformly random in 1,,n,Nt=n)=1nk=1nP(SksNt=n)=1nk=1nj=0k1P(Ns=jNt=n)=1nk=1nj=0k1P(Ns=j)P(Nts=nj)P(Nt=n)=1nk=1nj=0k1eλs(λs)jj!eλ(ts)(λ(ts))nj(nj)!eλt(λt)nn!=1nj=0n1(nj)n!j!(nj)!(st)j(1st)nj=[j=0n1n!j!(nj1)!(st)j(1st)nj1](1st)=1st\begin{align*} P(S_k \geq s \mid k \text{ uniformly random in } {1, \ldots, n}, N_t = n) & = \frac{1}{n} \sum_{k = 1}^{n} P(S_k \geq s \mid N_t = n) \newline & = \frac{1}{n} \sum_{k = 1}^{n} \sum_{j = 0}^{k - 1} P(N_s = j \mid N_t = n) \newline & = \frac{1}{n} \sum_{k = 1}^{n} \sum_{j = 0}^{k - 1} \frac{P(N_s = j) P(N_{t - s} = n - j)}{P(N_t = n)} \newline & = \frac{1}{n} \sum_{k = 1}^{n} \sum_{j = 0}^{k - 1} \frac{e^{-\lambda s} \frac{(\lambda s)^j}{j!} e^{-\lambda (t - s)} \frac{(\lambda (t - s))^{n - j}}{(n - j)!}}{e^{-\lambda t} \frac{(\lambda t)^n}{n!}} \newline & = \frac{1}{n} \sum_{j = 0}^{n - 1} (n - j) \frac{n!}{j! (n - j)!} \left(\frac{s}{t}\right)^j \left(1 - \frac{s}{t}\right)^{n - j} \newline & = \left[\sum_{j = 0}^{n - 1} \frac{n!}{j! (n - j - 1)!} \left(\frac{s}{t}\right)^j \left(1 - \frac{s}{t}\right)^{n - j - 1}\right] \cdot \left(1 - \frac{s}{t}\right) \newline & = 1 - \frac{s}{t} \end{align*}

Spatial Poisson Processes and Non-Homogeneous Poisson Processes

Definition 1 (Spatial Poisson Process)

A collection of random variables {NA}ARd\{N_A\}_{A \subseteq \mathbb{R}^d} is a spatial Poisson process with parameter λ\lambda if,

  • For each bounded set ARdA \subseteq \mathbb{R}^d, NAN_A has a Poisson distribution with parameter λA\lambda |A|.
  • Whenever AB,NANBA \subseteq B, N_A \leq N_B (i.e., spatial counting process).
  • Whenever AA and BB are disjoint, NAN_A and NBN_B are independent.
Definition 2 (Non-Homogeneous Poisson Process)

A counting process {Nt}t0\{N_t\}_{t \geq 0} is a non-homogeneous Poisson process with intensity function λ(t)\lambda(t) if,

  • N0=0N_0 = 0.
  • For 0<s<t0 < s < t,
NtNsPoisson(stλ(x) dx).N_t - N_s \sim \mathrm{Poisson}\left(\int_{s}^{t} \lambda(x) \ dx\right).
  • It has independent increments.