Factorial moment measure

In probability and statistics, a factorial moment measure is a mathematical quantity, function or, more precisely, measure that is defined in relation to mathematical objects known as point processes, which are types of stochastic processes often used as mathematical models of physical phenomena representable as randomly positioned points in time, space or both. Moment measures generalize the idea of factorial moments, which are useful for studying non-negative integer-valued random variables.[1]

The first factorial moment measure of a point process coincides with its first moment measure or intensity measure,[2] which gives the expected or average number of points of the point process located in some region of space. In general, if the number of points in some region is considered as a random variable, then the moment factorial measure of this region is the factorial moment of this random variable.[3] Factorial moment measures completely characterize a wide class of point processes, which means they can be used to uniquely identify a point process.

If a factorial moment measure is absolutely continuous, then with respect to the Lebesgue measure it is said to have a density (which is a generalized form of a derivative), and this density is known by a number of names such as factorial moment density and product density, as well as coincidence density,[1] joint intensity[4] , correlation function or multivariate frequency spectrum[5] The first and second factorial moment densities of a point process are used in the definition of the pair correlation function, which gives a way to statistically quantify the strength of interaction or correlation between points of a point process.[6]

Factorial moment measures serve as useful tools in the study of point processes[1][6][7] as well as the related fields of stochastic geometry[3] and spatial statistics,[6][8] which are applied in various scientific and engineering disciplines such as biology, geology, physics, and telecommunications.[1][3][9]

Point process notation

Point processes are mathematical objects that are defined on some underlying mathematical space. Since these processes are often used to represent collections of points randomly scattered in space, time or both, the underlying space is usually d-dimensional Euclidean space denoted here by R^d, but they can be defined on more abstract mathematical spaces.[7]

Point processes have a number of interpretations, which is reflected by the various types of point process notation.[3][9] For example, if a point ${\displaystyle \textstyle x}$ belongs to or is a member of a point process, denoted by N, then this can be written as:[3]

${\displaystyle \textstyle x\in {N},}$

and represents the point process being interpreted as a random set. Alternatively, the number of points of N located in some Borel set B is often written as:[2][3][8]

${\displaystyle \textstyle {N}(B),}$

which reflects a random measure interpretation for point processes. These two notations are often used in parallel or interchangeably.[3][8][2]

Definitions

n th factorial power of a point process

For some positive integer ${\displaystyle \textstyle n=1,2,\ldots }$, the ${\displaystyle \textstyle n}$-th factorial power of a point process ${\displaystyle \textstyle {N}}$ on ${\displaystyle \textstyle {\textbf {R}}^{d}}$ is defined as:[2]

${\displaystyle {N}^{(n)}(B_{1}\times \cdots \times B_{n})=\sum _{(x_{1}\neq \dots \neq x_{n})\in {N}}\prod _{i=1}^{n}\mathbf {1} _{B_{i}}(x_{i})}$

where ${\displaystyle \textstyle B_{1},...,B_{n}}$ is a collection of not necessarily disjoint Borel sets in ${\displaystyle \textstyle {\textbf {R}}^{d}}$, which form an ${\displaystyle \textstyle n}$-fold Cartesian product of sets denoted by:

${\displaystyle B_{1}\times \cdots \times B_{n}.}$

The symbol ${\displaystyle \textstyle \mathbf {1} }$ denotes an indicator function such that ${\displaystyle \textstyle \mathbf {1} _{B_{1}}}$is a Dirac measure for the set ${\displaystyle \textstyle B_{n}}$. The summation in the above expression is performed over all ${\displaystyle \textstyle n}$-tuples of distinct points, including permutations, which can be contrasted with the definition of the n-th power of a point process. The symbol ${\displaystyle \textstyle \Pi }$ denotes multiplication while the existence of various point process notation means that the n-th factorial power of a point process is sometimes defined using other notation.[2]

n th factorial moment measure

The n th factorial moment measure or n th order factorial moment measure is defined as:

${\displaystyle M^{(n)}(B_{1}\times \cdots \times B_{n})=E[{N}^{(n)}(B_{1}\times \cdots \times B_{n})],}$

where the E denotes the expectation (operator) of the point process N. In other words, the n-th factorial moment measure is the expectation of the n th factorial power of some point process.

The n th factorial moment measure of a point process N is equivalently defined[3] by:

${\displaystyle \int _{{\textbf {R}}^{nd}}f(x_{1},\ldots ,x_{n})M^{(n)}(dx_{1},\ldots ,dx_{n})=E\left[\sum _{(x_{1}\neq \cdots \neq x_{n})\in {N}}f(x_{1},\ldots ,x_{n})\right],}$

where ${\displaystyle f}$ is any non-negative measurable function on ${\displaystyle {\textbf {R}}^{nd}}$, and the above summation is performed over all ${\displaystyle n}$ tuples of distinct points, including permutations. Consequently, the factorial moment measure is defined such that there are no points repeating in the product set, as opposed to the moment measure.[7]

First factorial moment measure

The first factorial moment measure ${\displaystyle \textstyle M^{1}}$ coincides with the first moment measure:[2]

${\displaystyle M^{(1)}(B)=M^{1}(B)=E[{N}(B)],}$

where ${\displaystyle \textstyle M^{1}}$ is known, among other terms, as the intensity measure[3] or mean measure,[10] and is interpreted as the expected number of points of ${\displaystyle \textstyle {N}}$ found or located in the set ${\displaystyle \textstyle B}$

Second factorial moment measure

The second factorial moment measure for two Borel sets ${\displaystyle \textstyle A}$ and ${\displaystyle \textstyle B}$ is:

${\displaystyle M^{(2)}(A\times B)=M^{2}(A\times B)-M^{1}(A\cap B).}$

Name explanation

For some Borel set ${\displaystyle \textstyle B}$, the namesake of this measure is revealed when the ${\displaystyle \textstyle n\,}$th factorial moment measure reduces to:

${\displaystyle M^{(n)}(B\times \cdots \times B)=E[{N}(B)({N}(B)-1)\cdots ({N}(B)-n+1)],}$

which is the ${\displaystyle \textstyle n\,}$-th factorial moment of the random variable ${\displaystyle \textstyle {N}(B)}$.[3]

Factorial moment density

If a factorial moment measure is absolutely continuous, then it has a density (or more precisely, a Radon–Nikodym derivative or density) with respect to the Lebesgue measure and this density is known as the factorial moment density or product density, joint intensity, correlation function, or multivariate frequency spectrum. Denoting the ${\displaystyle \textstyle n}$-th factorial moment density by ${\displaystyle \textstyle \mu ^{(n)}(x_{1},\ldots ,x_{n})}$, it is defined in respect to the equation:[3]

${\displaystyle M^{(n)}(B_{1}\times \ldots \times B_{n})=\int _{B_{1}}\cdots \int _{B_{n}}\mu ^{(n)}(x_{1},\ldots ,x_{n})dx_{1}\cdots dx_{n}.}$

Furthermore, this means the following expression

${\displaystyle E\left[\sum _{(x_{1}\neq \cdots \neq x_{n})\in {N}}f(x_{1},\ldots ,x_{n})\right]=\int _{{\textbf {R}}^{nd}}f(x_{1},\ldots ,x_{n})\mu ^{(n)}(x_{1},\ldots ,x_{n})dx_{1}\cdots dx_{n},}$

where ${\displaystyle \textstyle f}$ is any non-negative bounded measurable function defined on ${\displaystyle \textstyle {\textbf {R}}^{n}}$.

Pair correlation function

In spatial statistics and stochastic geometry, to measure the statistical correlation relationship between points of a point process, the pair correlation function of a point process ${\displaystyle {N}}$ is defined as:[3][6]

${\displaystyle \rho (x_{1},x_{2})={\frac {\mu ^{(2)}(x_{1},x_{2})}{\mu ^{(1)}(x_{1})\mu ^{(1)}(x_{2})}},}$

where the points ${\displaystyle x_{1},x_{2}\in R^{d}}$. In general, ${\displaystyle \rho (x_{1},x_{2})\geq 0}$ whereas ${\displaystyle \rho (x_{1},x_{2})=1}$ corresponds to no correlation (between points) in the typical statistical sense.[6]

Examples

Poisson point process

For a general Poisson point process with intensity measure ${\displaystyle \textstyle \Lambda }$ the ${\displaystyle \textstyle n}$-th factorial moment measure is given by the expression:[3]

${\displaystyle M^{(n)}(B_{1}\times \cdots \times B_{n})=\prod _{i=1}^{n}[\Lambda (B_{i})],}$

where ${\displaystyle \textstyle \Lambda }$ is the intensity measure or first moment measure of ${\displaystyle \textstyle {N}}$, which for some Borel set ${\displaystyle \textstyle B}$ is given by:

${\displaystyle \Lambda (B)=M^{1}(B)=E[{N}(B)].}$

For a homogeneous Poisson point process the ${\displaystyle \textstyle n}$-th factorial moment measure is simply:[2]

${\displaystyle M^{(n)}(B_{1}\times \cdots \times B_{n})=\lambda ^{n}\prod _{i=1}^{n}|B_{i}|,}$

where ${\displaystyle \textstyle |B_{i}|}$ is the length, area, or volume (or more generally, the Lebesgue measure) of ${\displaystyle \textstyle B_{i}}$. Furthermore, the ${\displaystyle \textstyle n}$-th factorial moment density is:[3]

${\displaystyle \mu ^{(n)}(x_{1},\ldots ,x_{n})=\lambda ^{n}.}$

The pair-correlation function of the homogeneous Poisson point process is simply

${\displaystyle \rho (x_{1},x_{2})=1,}$

which reflects the lack of interaction between points of this point process.

Factorial moment expansion

The expectations of general functionals of simple point processes, provided some certain mathematical conditions, have (possibly infinite) expansions or series consisting of the corresponding factorial moment measures.[11][12] In comparison to the Taylor series, which consists of a series of derivatives of some function, the nth factorial moment measure plays the roll as that of the n th derivative the Taylor series. In other words, given a general functional f of some simple point process, then this Taylor-like theorem for non-Poisson point processes means an expansion exists for the expectation of the function E, provided some mathematical condition is satisfied, which ensures convergence of the expansion.

See also

• Factorial moment
• Moment
• Moment measure

References

1. D. J. Daley and D. Vere-Jones. An introduction to the theory of point processes. Vol. I. Probability and its Applications (New York). Springer, New York, second edition, 2003.
2. Baccelli, François (2009). "Stochastic Geometry and Wireless Networks: Volume I Theory" (PDF). Foundations and Trends in Networking. 3 (3–4): 249–449. doi:10.1561/1300000006. ISSN 1554-057X.
3. D. Stoyan, W. S. Kendall, J. Mecke, and L. Ruschendorf. Stochastic geometry and its applications, volume 2. Wiley Chichester, 1995.
4. Hough, J Ben, Krishnapur, Manjunath, Peres, Yuval, Vir{\'a}g, B{\'a}lint (2006). "Determinantal processes and independence". Probability Surveys. 3: 206–229. arXiv:math/0503110. doi:10.1214/154957806000000078.CS1 maint: multiple names: authors list (link)
5. K. Handa. The two-parameter {Poisson-Dirichlet} point process. Bernoulli, 15(4):1082–1116, 2009.
6. A. Baddeley, I. B{\'a}r{\'a}ny, and R. Schneider. Spatial point processes and their applications. Stochastic Geometry: Lectures given at the CIME Summer School held in Martina Franca, Italy, September 13–18, 2004, pages 1–75, 2007.
7. D. J. Daley and D. Vere-Jones. An introduction to the theory of point processes. Vol. {II}. Probability and its Applications (New York). Springer, New York, second edition, 2008
8. Møller, Jesper; Waagepetersen, Rasmus Plenge (2003). Statistical Inference and Simulation for Spatial Point Processes. C&H/CRC Monographs on Statistics & Applied Probability. 100. CiteSeerX 10.1.1.124.1275. doi:10.1201/9780203496930. ISBN 978-1-58488-265-7.
9. F. Baccelli and B. Błaszczyszyn. Stochastic Geometry and Wireless Networks, Volume II – Applications, volume 4, No 1–2 of Foundations and Trends in Networking. NoW Publishers, 2009.
10. J. F. C. Kingman. Poisson processes, volume 3. Oxford university press, 1992.
11. B. Blaszczyszyn. Factorial-moment expansion for stochastic systems. Stoch. Proc. Appl., 56:321–335, 1995.
12. D. P. Kroese and V. Schmidt. Light-traffic analysis for queues with spatially distributed arrivals. Mathematics of Operations Research, 21(1):pp. 135–157, 1996.