Poisson-type random measure

{{Short description|Family of three random counting measures}}

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.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 The PT family of distributions is also known as the Katz family of distributions,Katz L.. Classical and Contagious Discrete Distributions ch. Unified treatment of a broad class of discrete probability distributions, :175-182. Pergamon Press, Oxford 1965. the Panjer or (a,b,0) class of distributionsPanjer Harry H.. Recursive Evaluation of a Family of Compound Distributions. 1981;12(1):22-26 and may be retrieved through the Conway–Maxwell–Poisson distribution.Conway R. W., Maxwell W. L.. A Queuing Model with State Dependent Service Rates. Journal of Industrial Engineering. 1962;12.

Throwing stones

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

The random counting measure N on (E,\mathcal{E}) depends on the pair of deterministic probability measures (\kappa,\nu) through the stone throwing construction (STC) Cinlar Erhan. Probability and Stochastics. Springer-Verlag New York; 2011

\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 K has law \kappa and iid X_1,X_2,\dotsb have law \nu. N is a mixed binomial processKallenberg Olav. Random Measures, Theory and Applications. Springer; 2017

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

\quad\mathbb{E} e^{-N f} =\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 \psi(\cdot) is the generating function of K. The mean and variance are given by

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

and

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

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

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

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

:

\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 N=(\kappa,\nu) to the case when the collection \mathbf{X} is expanded to (\mathbf{X},\mathbf{Y})=\{(X_i,Y_i)\} where Y_i is a random transformation of X_i. Heuristically, Y_i represents some properties (marks) of X_i. We assume that the conditional law of Y follows some transition kernel according to \mathbb{P}(Y\in B|X=x)=Q(x,B).

= Theorem: Marked STC =

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

The following corollary is an immediate consequence.

= Corollary: Restricted STC =

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

Note \psi(\nu e^{-f}\mathbb{I}_A+1-a)=\psi_A(\nu_A e^{-f}) where we use \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 K_A = N\mathbb{I}_A be the counting variable of K restricted to A\subset E. When \{N\mathbb{I}_A: A\subset E\} and K=N\mathbb{I}_E share the same family of laws subject to a rescaling h_a(\theta) of the parameter \theta, then K is a called a bone distribution. The bone condition for the pgf is given by

\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 K\sim \kappa_\theta with pgf \psi_\theta belongs to the canonical non-negative power series (NNPS) family of distributions and \{0,1\}\subset\text{supp}(K). Consider the random measure N=(\kappa_\theta,\nu) on the space (E,\mathcal{E}) and assume that \nu is diffuse. Then for any A\subset E with \nu(A)=a>0 there exists a mapping h_a:\Theta\rightarrow\Theta such that the restricted random measure is N_A=(\kappa_{h_a(\theta)},\nu_A), that is,

\quad \mathbb{E} e^{-N_A f} = \psi_{h_a(\theta)}(\nu_A e^{-f})\quad \text{for}\quad f\in\mathcal{E}_+

iff 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 N\mathbb{I}_A share the same family of distribution as 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 p\rightarrow 1, n\rightarrow c.

Distributional self-similarity applications

The "bone" condition on the pgf \psi_\theta of K encodes a distributional self-similarity property whereby all counts in restrictions (thinnings) to subspaces (encoded by pgf \psi_A) are in the same family as \psi_\theta of K through rescaling of the canonical parameter. These ideas appear closely connected to those of self-decomposability and stability of discrete random variables.Steutel FW, Van Harn K. Discrete analogues of self-decomposability and stability. The Annals of Probability. 1979;:893–899. Binomial thinning is a foundational model to count time-series.Al-Osh M. A., Alzaid A. A.. First-order integer-valued autoregressive (INAR(1)) process. Journal of Time Series Analysis. 1987;8(3):261–275.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. 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 Lévy processes, the jumps of Kolmogorov equations, and the excursions of Brownian motion.Cinlar Erhan. Probability and Stochastics. Springer-Verlag New York; 2011. 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