Compound Poisson process

A compound Poisson process is a continuous-time (random) stochastic process with jumps. The jumps arrive randomly according to a Poisson process and the size of the jumps is also random, with a specified probability distribution. A compound Poisson process, parameterised by a rate ${\displaystyle \lambda >0}$ and jump size distribution G, is a process ${\displaystyle \{\,Y(t):t\geq 0\,\}}$ given by

${\displaystyle Y(t)=\sum _{i=1}^{N(t)}D_{i}}$

where, ${\displaystyle \{\,N(t):t\geq 0\,\}}$ is a Poisson process with rate ${\displaystyle \lambda }$, and ${\displaystyle \{\,D_{i}:i\geq 1\,\}}$ are independent and identically distributed random variables, with distribution function G, which are also independent of ${\displaystyle \{\,N(t):t\geq 0\,\}.\,}$

When ${\displaystyle Y(t)}$ are non-negative discrete random variable,then this compound Poisson process is named stuttering Poisson process which has the feature that two or more events occur in a very short time .

Properties of the compound Poisson process

Using conditional expectation, the expected value of a compound Poisson process can be calculated as:

${\displaystyle \,E(Y(t))=E(E(Y(t)|N(t)))=E(N(t)E(D))=E(N(t))E(D)=\lambda tE(D).}$ 

Making similar use of the law of total variance, the variance can be calculated as:

${\displaystyle {\begin{aligned}\operatorname {var} (Y(t))&=E(\operatorname {var} (Y(t)|N(t)))+\operatorname {var} (E(Y(t)|N(t)))\\&=E(N(t)\operatorname {var} (D))+\operatorname {var} (N(t)E(D))\\&=\operatorname {var} (D)E(N(t))+E(D)^{2}\operatorname {var} (N(t))\\&=\operatorname {var} (D)\lambda t+E(D)^{2}\lambda t\\&=\lambda t(\operatorname {var} (D)+E(D)^{2})\\&=\lambda tE(D^{2}).\end{aligned}}}$ 

Lastly, using the law of total probability, the moment generating function can be given as follows:

${\displaystyle \,\Pr(Y(t)=i)=\sum _{n}\Pr(Y(t)=i|N(t)=n)\Pr(N(t)=n)}$ 

${\displaystyle {\begin{aligned}E(e^{sY})&=\sum _{i}e^{si}\Pr(Y(t)=i)\\&=\sum _{i}e^{si}\sum _{n}\Pr(Y(t)=i|N(t)=n)\Pr(N(t)=n)\\&=\sum _{n}\Pr(N(t)=n)\sum _{i}e^{si}\Pr(Y(t)=i|N(t)=n)\\&=\sum _{n}\Pr(N(t)=n)\sum _{i}e^{si}\Pr(D_{1}+D_{2}+\cdots +D_{n}=i)\\&=\sum _{n}\Pr(N(t)=n)M_{D}(s)^{n}\\&=\sum _{n}\Pr(N(t)=n)e^{n\ln(M_{D}(s))}\\&=M_{N(t)}(\ln(M_{D}(s))\\&=e^{\lambda t\left(M_{D}(s)-1\right)}.\end{aligned}}}$ 

Exponentiation of measures

Let N, Y, and D be as above. Let μ be the probability measure according to which D is distributed, i.e.

${\displaystyle \mu (A)=\Pr(D\in A).\,}$ 

Let δ₀ be the trivial probability distribution putting all of the mass at zero. Then the probability distribution of Y(t) is the measure

${\displaystyle \exp(\lambda t(\mu -\delta _{0}))\,}$ 

where the exponential exp(ν) of a finite measure ν on Borel subsets of the real line is defined by

${\displaystyle \exp(\nu )=\sum _{n=0}^{\infty }{\nu ^{*n} \over n!}}$ 

and

${\displaystyle \nu ^{*n}=\underbrace {\nu *\cdots *\nu } _{n{\text{ factors}}}}$ 

is a convolution of measures, and the series converges weakly.

Fitting a compound Poisson process

The parameters for independent observations of a compound Poisson process can be chosen using a maximum likelihood estimator using Simar's algorithm,^[1] which has been shown to converge.^[2]

See also

• Poisson process
• Poisson distribution
• Non-homogeneous Poisson process

References

1. Attention: This template ({{cite doi}}) is deprecated. To cite the publication identified by doi:10.1214/aos/1176343651, please use {{cite journal}} (if it was published in a bona fide academic journal, otherwise {{cite report}} with |doi=10.1214/aos/1176343651 instead.
2. Attention: This template ({{cite doi}}) is deprecated. To cite the publication identified by doi:10.1214/aos/1176345890, please use {{cite journal}} (if it was published in a bona fide academic journal, otherwise {{cite report}} with |doi=10.1214/aos/1176345890 instead.