Poisson-type random measure

Poisson-type random measures are a family of three random counting measures which are closed under restriction to a subspace, i.e. closed under thinning. They are the only distributions in the canonical non-negative power series family of distributions to possess this property and include the Poisson distribution, negative binomial distribution, and binomial distribution. ^[1] The PT family of distributions is also known as the Katz family of distributions, ^[2] the Panjer or (a,b,0) class of distributions ^[3] and may be retrieved through the Conway–Maxwell–Poisson distribution. ^[4]

Throwing stones

Let ${\displaystyle K}$ be a non-negative integer-valued random variable ${\displaystyle K\in \mathbb {N} _{\geq 0}=\mathbb {N} _{>0}\cup \{0\}}$) with law ${\displaystyle \kappa }$, mean ${\displaystyle c\in (0,\infty )}$ and when it exists variance ${\displaystyle \delta ^{2}>0}$. Let ${\displaystyle \nu }$ be a probability measure on the measurable space ${\displaystyle (E,{\mathcal {E}})}$. Let ${\displaystyle \mathbf {X} =\{X_{i}\}}$ be a collection of iid random variables (stones) taking values in ${\displaystyle (E,{\mathcal {E}})}$ with law ${\displaystyle \nu }$.

The random counting measure ${\displaystyle N}$ on ${\displaystyle (E,{\mathcal {E}})}$ depends on the pair of deterministic probability measures ${\displaystyle (\kappa ,\nu )}$ through the stone throwing construction (STC) ^[5]

${\displaystyle \quad N_{\omega }(A)=N(\omega ,A)=\sum _{i=1}^{K(\omega )}\mathbb {I} _{A}(X_{i}(\omega ))\quad {\text{for}}\quad \omega \in \Omega ,\,\,\,A\in {\mathcal {E}}}$

where ${\displaystyle K}$ has law ${\displaystyle \kappa }$ and iid ${\displaystyle X_{1},X_{2},\dotsb }$ have law ${\displaystyle \nu }$. ${\displaystyle N}$ is a mixed binomial process ^[6]

Let ${\displaystyle {\mathcal {E}}_{+}=\{f:E\mapsto \mathbb {R} _{+}\}}$ be the collection of positive ${\displaystyle {\mathcal {E}}}$-measurable functions. The probability law of ${\displaystyle N}$ is encoded in the Laplace functional

${\displaystyle \quad \mathbb {E} e^{-Nf}=\mathbb {E} (\mathbb {E} e^{-f(X)})^{K}=\mathbb {E} (\nu e^{-f})^{K}=\psi (\nu e^{-f})\quad {\text{for}}\quad f\in {\mathcal {E}}_{+}}$

where ${\displaystyle \psi (\cdot )}$ is the generating function of ${\displaystyle K}$. The mean and variance are given by

${\displaystyle \quad \mathbb {E} Nf=c\nu f}$

and

${\displaystyle \quad \mathbb {V} {\text{ar}}Nf=c\nu f^{2}+(\delta ^{2}-c)(\nu f)^{2}}$

The covariance for arbitrary ${\displaystyle f,g\in {\mathcal {E}}_{+}}$ is given by

${\displaystyle \quad \mathbb {C} {\text{ov}}(Nf,Ng)=c\nu (fg)+(\delta ^{2}-c)\nu f\nu g}$

When ${\displaystyle K}$ is Poisson, negative binomial, or binomial, it is said to be Poisson-type (PT). The joint distribution of the collection ${\displaystyle N(A),\ldots ,N(B)}$ is for ${\displaystyle i,\ldots ,j\in \mathbb {N} }$ and ${\displaystyle i+\cdots +j=k}$

${\displaystyle \mathbb {P} (N(A)=i,\ldots ,N(B)=j)=\mathbb {P} (N(A)=i,\ldots ,N(B)=j|K=k)\,\mathbb {P} (K=k)={\frac {k!}{i!\cdots j!}}\,\nu (A)^{i}\cdots \nu (B)^{j}\,\mathbb {P} (K=k)}$

The following result extends construction of a random measure ${\displaystyle N=(\kappa ,\nu )}$ to the case when the collection ${\displaystyle \mathbf {X} }$ is expanded to ${\displaystyle (\mathbf {X} ,\mathbf {Y} )=\{(X_{i},Y_{i})\}}$ where ${\displaystyle Y_{i}}$ is a random transformation of ${\displaystyle X_{i}}$. Heuristically, ${\displaystyle Y_{i}}$ represents some properties (marks) of ${\displaystyle X_{i}}$. We assume that the conditional law of ${\displaystyle Y}$ follows some transition kernel according to ${\displaystyle \mathbb {P} (Y\in B|X=x)=Q(x,B)}$.

Theorem: Marked STC

Consider random measure ${\displaystyle N=(\kappa ,\nu )}$ and the transition probability kernel ${\displaystyle Q}$ from ${\displaystyle (E,{\cal {E)}}}$ into ${\displaystyle (F,{\cal {F)}}}$. Assume that given the collection ${\displaystyle \mathbf {X} }$ the variables ${\displaystyle \mathbf {Y} =\{Y_{i}\}}$ are conditionally independent with ${\displaystyle Y_{i}\sim Q(X_{i},\cdot )}$. Then ${\displaystyle M=(\kappa ,\nu \times Q)}$ is a random measure on ${\displaystyle (E\times F,{\cal {E\otimes F)}}}$. Here ${\displaystyle \mu =\nu \times Q}$ is understood as ${\displaystyle \mu (dx,dy)=\nu (dx)Q(x,dy)}$. Moreover, for any ${\displaystyle f\in ({\cal {E}}\otimes {\cal {F}})_{+}}$ we have that ${\displaystyle \mathbb {E} e^{-Mf}=\psi (\nu e^{-g})}$ where ${\displaystyle \psi (\cdot )}$ is pgf of ${\displaystyle K}$ and ${\displaystyle g\in {\mathcal {E}}_{+}}$ is defined as ${\displaystyle e^{-g(x)}=\int _{F}Q(x,dy)e^{-f(x,y)}.}$

The following corollary is an immediate consequence.

Corollary: Restricted STC

The quantity ${\displaystyle N_{A}=(N\mathbb {I} _{A},\nu _{A})}$ is a well-defined random measure on the measurable subspace ${\displaystyle (E\cap A,{\mathcal {E}}_{A})}$ where ${\displaystyle {\mathcal {E}}_{A}=\{A\cap B:B\in {\mathcal {E}}\}}$ and ${\displaystyle \nu _{A}(B)=\nu (A\cap B)/\nu (A)}$. Moreover, for any ${\displaystyle f\in {\mathcal {E}}_{+}}$, we have that ${\displaystyle \mathbb {E} e^{-N_{A}f}=\psi (\nu e^{-f}\mathbb {I} _{A}+b)}$ where ${\displaystyle b=1-\nu (A)}$.

Note ${\displaystyle \psi (\nu e^{-f}\mathbb {I} _{A}+1-a)=\psi _{A}(\nu _{A}e^{-f})}$ where we use ${\displaystyle \nu e^{-f}\mathbb {I} _{A}=a\nu _{A}e^{-f}}$.

Collecting Bones

The probability law of the random measure is determined by its Laplace functional and hence generating function.

Definition: Bone

Let ${\displaystyle K_{A}=N\mathbb {I} _{A}}$ be the counting variable of ${\displaystyle K}$ restricted to ${\displaystyle A\subset E}$. When ${\displaystyle \{N\mathbb {I} _{A}:A\subset E\}}$ and ${\displaystyle K=N\mathbb {I} _{E}}$ share the same family of laws subject to a rescaling ${\displaystyle h_{a}(\theta )}$ of the parameter ${\displaystyle \theta }$, then ${\displaystyle K}$ is a called a bone distribution. The bone condition for the pgf is given by ${\displaystyle \psi _{\theta }(at+1-a)=\psi _{h_{a}(\theta )}(t)}$.

Equipped with the notion of a bone distribution and condition, the main result for the existence and uniqueness of Poisson-type (PT) random counting measures is given as follows.

Theorem: existence and uniqueness of PT random measures

Assume that ${\displaystyle K\sim \kappa _{\theta }}$ with pgf ${\displaystyle \psi _{\theta }}$ belongs to the canonical non-negative power series (NNPS) family of distributions and ${\displaystyle \{0,1\}\subset {\text{supp}}(K)}$. Consider the random measure ${\displaystyle N=(\kappa _{\theta },\nu )}$ on the space ${\displaystyle (E,{\mathcal {E}})}$ and assume that ${\displaystyle \nu }$ is diffuse. Then for any ${\displaystyle A\subset E}$ with ${\displaystyle \nu (A)=a>0}$ there exists a mapping ${\displaystyle h_{a}:\Theta \rightarrow \Theta }$ such that the restricted random measure is ${\displaystyle N_{A}=(\kappa _{h_{a}(\theta )},\nu _{A})}$, that is,

${\displaystyle \quad \mathbb {E} e^{-N_{A}f}=\psi _{h_{a}(\theta )}(\nu _{A}e^{-f})\quad {\text{for}}\quad f\in {\mathcal {E}}_{+}}$

iff ${\displaystyle K}$ is Poisson, negative binomial, or binomial (Poisson-type).

The proof for this theorem is based on a generalized additive Cauchy equation and its solutions. The theorem states that out of all NNPS distributions, only PT have the property that their restrictions ${\displaystyle N\mathbb {I} _{A}}$ share the same family of distribution as ${\displaystyle K}$, that is, they are closed under thinning. The PT random measures are the Poisson random measure, negative binomial random measure, and binomial random measure. Poisson is additive with independence on disjoint sets, whereas negative binomial has positive covariance and binomial has negative covariance. The binomial process is a limiting case of binomial random measure where ${\displaystyle p\rightarrow 1,n\rightarrow c}$.

Distributional self-similarity applications

The "bone" condition on the pgf ${\displaystyle \psi _{\theta }}$ of ${\displaystyle K}$ encodes a distributional self-similarity property whereby all counts in restrictions (thinnings) to subspaces (encoded by pgf ${\displaystyle \psi _{A}}$) are in the same family as ${\displaystyle \psi _{\theta }}$ of ${\displaystyle K}$ through rescaling of the canonical parameter. These ideas appear closely connected to those of self-decomposability and stability of discrete random variables. ^[7] Binomial thinning is a foundational model to count time-series. ^{[8] [9]} The Poisson random measure has the well-known splitting property, is prototypical to the class of additive (completely random) random measures, and is related to the structure of Levy processes, the jumps of Kolmogorov equations (Markov jump process), and the excursions of Brownian motion. ^[10] Hence the self-similarity property of the PT family is fundamental to multiple areas. The PT family members are "primitives" or prototypical random measures by which many random measures and processes can be constructed.

References

1. Caleb Bastian, Gregory Rempala. Throwing stones and collecting bones: Looking for Poisson-like random measures, Mathematical Methods in the Applied Sciences, 2020. doi:10.1002/mma.6224
2. Katz L.. Classical and Contagious Discrete Distributions ch. Unified treatment of a broad class of discrete probability distributions, :175-182. Pergamon Press, Oxford 1965.
3. Panjer Harry H.. Recursive Evaluation of a Family of Compound Distributions. 1981;12(1):22-26
4. Conway R. W., Maxwell W. L.. A Queuing Model with State Dependent Service Rates. Journal of Industrial Engineering. 1962;12.
5. Cinlar Erhan. Probability and Stochastics. Springer-Verlag New York; 2011
6. Kallenberg Olav. Random Measures, Theory and Applications. Springer; 2017
7. Steutel FW, Van Harn K. Discrete analogues of self-decomposability and stability. The Annals of Probability. 1979;:893–899.
8. Al-Osh M. A., Alzaid A. A.. First-order integer-valued autogressive (INAR(1)) process. Journal of Time Series Analysis. 1987;8(3):261–275.
9. Scotto Manuel G., Weiß Christian H., Gouveia Sónia. Thinning models in the analysis of integer-valued time series: a review. Statistical Modelling. 2015;15(6):590–618.
10. Cinlar Erhan. Probability and Stochastics. Springer-Verlag New York; 2011.