From Counts to Functionals
A Poisson random measure describes random points through the numbers that fall in measurable sets. Many quantities of interest, however, depend on the points through a function: a total weight, a spatial cost, or a product of contributions. Two formulas connect these point-based quantities to the intensity measure. The Laplace functional gives the transform of a sum over points, while Campbell’s formula gives its expectation.
Work on a probability space \((\Omega,\mathcal{F},\mathbb{P})\), and let \((E,\mathcal{E})\) be a measurable space with a sigma-finite measure \(\nu\). A Poisson random measure \(N\) on \(E\) with intensity measure \(\nu\) has the following properties: for every \(A\in\mathcal{E}\) with \(\nu(A)<\infty\), \(N(A)\) has a Poisson distribution with mean \(\nu(A)\); and counts on any finite collection of pairwise disjoint measurable sets of finite measure are independent. Here \(N(A)\) counts points in \(A\), including multiplicities if points coincide.
For a nonnegative measurable function \(f:E\to[0,\infty]\), its integral against \(N\), written \(\int_E f(x)\,N(dx)\), is the sum of \(f\) over the points, with multiplicities. It is defined first for indicator functions, then for nonnegative simple functions, and then by increasing approximation for all nonnegative measurable functions. The integral may be infinite. An integral against \(\nu\) is defined in the same way.
The Laplace functional describes the distributional transform of a point sum, not just its mean. Campbell’s formula extracts the mean directly. Both apply to nonnegative measurable functions without requiring boundedness or finite support; when a right-hand side is infinite, the corresponding expectation is infinite.
Why the Formulas Hold
Begin with a measurable set \(A\) of finite intensity \(m=\nu(A)\). By the Poisson distribution, for \(a\ge0\),
This calculation also covers \(m=0\): then \(N(A)=0\) almost surely, and both sides equal \(1\). For the expectation, the Poisson mean gives \(\mathbb{E}N(A)=\nu(A)\), again including the zero-measure case.
Proof. First take a nonnegative simple function \(f=\sum_{\ell=1}^r a_\ell\mathbf{1}_{A_\ell}\), where the \(A_\ell\) are pairwise disjoint measurable sets of finite \(\nu\)-measure and \(a_\ell\ge0\). The random integral equals \(\sum_{\ell=1}^r a_\ell N(A_\ell)\). Independence of these counts factors the expectation of its exponential into \(r\) terms. Applying the Poisson calculation above to each term gives $$ \mathbb{E}\exp\!\left(-\int_E f\,dN\right) =\prod_{\ell=1}^r\exp\!\bigl(-\nu(A_\ell)(1-e^{-a_\ell})\bigr) =\exp\!\left(-\int_E(1-e^{-f})\,d\nu\right). $$ For a simple nonnegative \(g=\sum_{\ell=1}^r b_\ell\mathbf{1}_{A_\ell}\) on such sets, linearity of expectation and \(\mathbb{E}N(A_\ell)=\nu(A_\ell)\) give $$ \mathbb{E}\int_E g\,dN =\sum_{\ell=1}^r b_\ell\nu(A_\ell) =\int_E g\,d\nu. $$ To pass to general nonnegative measurable functions, use sigma-finiteness to choose measurable sets \(E_m\) of finite measure increasing to \(E\). Each \(f\) can be approximated from below by nonnegative simple functions supported on these sets; choose \(f_m\uparrow f\), and likewise \(g_m\uparrow g\). Then \(\int f_m\,dN\uparrow\int f\,dN\) and \(\int g_m\,dN\uparrow\int g\,dN\). The exponentials \(e^{-\int f_m\,dN}\) decrease to \(e^{-\int f\,dN}\), and are bounded by \(1\), so bounded convergence passes their expectations to the limit. Also \(1-e^{-f_m}\uparrow1-e^{-f}\), so monotone convergence passes the intensity integrals to the limit. Monotone convergence applied to \(g_m\) passes both sides of Campbell’s formula to the limit. Thus both identities hold, including when the nonnegative integrals are infinite. \(\square\)
A useful consequence is that if \(\int_E g\,d\nu<\infty\), then \(\int_E g\,dN\) is finite almost surely. Indeed, a nonnegative random variable that equals infinity with positive probability has infinite expectation. The formulas also specify what to do when integrability fails: Campbell’s identity remains valid, but may give an infinite expectation.
Worked Calculations
Worked Example: A Transform from Two Disjoint Regions
Suppose \(A_1,A_2\) are disjoint measurable sets with \(\nu(A_1)=2\) and \(\nu(A_2)=1\). Define \(Y=(\log 2)N(A_1)+(\log 3)N(A_2)\). Find \(\mathbb{E}e^{-Y}\) and \(\mathbb{E}Y\).
Use the Laplace functional with \(f=(\log 2)\mathbf{1}_{A_1}+(\log 3)\mathbf{1}_{A_2}\). On \(A_1\), \(1-e^{-f}=1-\tfrac12=\tfrac12\); on \(A_2\), it is \(1-\tfrac13=\tfrac23\); elsewhere it is \(0\). Therefore
For the mean, Campbell’s formula gives
Check: the independent-count calculation gives \(\mathbb{E}2^{-N(A_1)}=\exp[-2(1-\tfrac12)]=e^{-1}\) and \(\mathbb{E}3^{-N(A_2)}=\exp[-1(1-\tfrac13)]=e^{-2/3}\). Their product is \(e^{-5/3}\), as obtained from the functional. Also \(\mathbb{E}N(A_1)=2\) and \(\mathbb{E}N(A_2)=1\), so the weighted mean is \(2\log2+\log3=\log12\).
Worked Example: A Weighted Sum of Poisson Counts
Let \(A\) and \(B\) be disjoint, with \(\nu(A)=\tfrac32\) and \(\nu(B)=\tfrac12\). Find the Laplace transform and mean of \(Z=2N(A)+5N(B)\).
For \(\theta\ge0\), choose \(f=2\theta\mathbf{1}_A+5\theta\mathbf{1}_B\). The Laplace functional gives
Campbell’s formula, applied to \(g=2\mathbf{1}_A+5\mathbf{1}_B\), gives
Check: at \(\theta=0\), each term \(1-e^{-c\theta}\) is zero, so the transform equals \(1\), as any Laplace transform must at zero. Differentiating its displayed expression at \(\theta=0\) gives \(-(\tfrac32\cdot2+\tfrac12\cdot5)=-\tfrac{11}{2}\), the negative of the mean just computed.
Worked Example: Radial Weights in a Spatial Poisson Model
Let \(E=\mathbb{R}^2\), and let \(N\) be a homogeneous Poisson random measure with intensity \(\nu(dx)=\tfrac1\pi\,dx\), where \(dx\) is area. In the disk \(D=\{x:\|x\|\le1\}\), let \(S_D=\int_D\|x\|\,N(dx)\), the sum of distances of points in the disk from its centre. Find its mean and Laplace transform.
Campbell’s formula and polar coordinates give
For \(\theta>0\), apply the Laplace functional to \(f(x)=\theta\|x\|\mathbf{1}_D(x)\). The intensity integral is
For the last equality, integration by parts gives \(\int_0^1re^{-\theta r}\,dr=[1-(1+\theta)e^{-\theta}]/\theta^2\), and \(\int_0^1r\,dr=\tfrac12\). Hence
At \(\theta=0\), the transform is \(1\), either by its defining expectation or by taking the limit in the formula. Check at \(\theta=1\): the intensity integral is \(1-2(1-2e^{-1})=4e^{-1}-1\), so \(\mathbb{E}e^{-S_D}=\exp(1-4e^{-1})\). This is at most \(1\), since the intensity integral \(4e^{-1}-1\) is positive; for example, \(e<4\) implies \(4e^{-1}>1\).
Where the Assumptions Matter
The independence condition concerns disjoint sets. Counts on overlapping sets are generally dependent, so multiplying their separate Laplace transforms can give a wrong answer. For a concrete calculation, suppose \(\nu(A\cap B)=\tfrac12\) and \(\nu(A\setminus B)=\nu(B\setminus A)=\tfrac12\). Then \(\nu(A)=\nu(B)=1\). For \(Y=(\log2)\{N(A)+N(B)\}\), incorrectly treating \(N(A)\) and \(N(B)\) as independent would give \(\mathbb{E}e^{-Y}=e^{-1}\cdot e^{-1}=e^{-2}\). The correction is to split the union into the three disjoint regions. On \(A\cap B\), the function value is \(2\log2\), and on each exclusive region it is \(\log2\). Thus the correct exponent in the Laplace functional is
Therefore the correct transform is \(e^{-7/8}\), not \(e^{-2}\). The overlap contributes to both counts, which is why the independent-count shortcut fails.
Another point to keep distinct is what each formula returns. The Laplace functional requires a nonnegative \(f\); Campbell’s formula requires a nonnegative \(g\). If \(\int g\,d\nu\) is finite, Campbell’s formula gives a finite expected sum. If it is infinite, the formula still holds in the extended sense and does not promise a finite sum. For functions that take both positive and negative values, one must separately control their positive and negative parts before using an expectation identity; the nonnegative version alone does not justify subtracting two infinite quantities.
These formulas are useful whenever a model records random points but the quantity of interest is built by adding contributions point by point. A simple function turns the calculation into a finite collection of Poisson counts. General nonnegative functions follow by increasing approximation, so spatial weights and other measurable contributions can be handled without requiring the points to lie in a finite set.
Check Your Understanding
For each question, identify the relevant function and check whether the sets used in a count calculation are disjoint.
- If \(N(A)\) is Poisson with mean \(m\), derive \(\mathbb{E}e^{-aN(A)}\) for \(a\ge0\).
- For disjoint \(A,B\) with \(\nu(A)=2\) and \(\nu(B)=1\), find \(\mathbb{E}e^{-\{(\log2)N(A)+(\log4)N(B)\}}\).
- If \(g=3\mathbf{1}_A+2\mathbf{1}_B\), where \(A,B\) are disjoint and have finite intensity, express \(\mathbb{E}\int g\,dN\) in terms of \(\nu(A)\) and \(\nu(B)\).
- What does Campbell’s formula imply about the almost-sure finiteness of \(\int g\,dN\) when \(\int g\,d\nu<\infty\)?
- Why can the Laplace transforms of \(N(A)\) and \(N(B)\) not generally be multiplied when \(A\) and \(B\) overlap?