Since point processes were historically developed by different communities, there are different mathematical interpretations of a point process, such as a random counting measure or a random set,910 and different notations. The notations are described in detail on the point process notation page.
Some authors regard a point process and stochastic process as two different objects such that a point process is a random object that arises from or is associated with a stochastic process,1112 though it has been remarked that the difference between point processes and stochastic processes is not clear.13 Others consider a point process as a stochastic process, where the process is indexed by sets of the underlying space14 on which it is defined, such as the real line or n {\displaystyle n} -dimensional Euclidean space.1516 Other stochastic processes such as renewal and counting processes are studied in the theory of point processes.1718 Sometimes the term "point process" is not preferred, as historically the word "process" denoted an evolution of some system in time, so point process is also called a random point field.19
In mathematics, a point process is a random element whose values are "point patterns" on a set S. While in the exact mathematical definition a point pattern is specified as a locally finite counting measure, it is sufficient for more applied purposes to think of a point pattern as a countable subset of S that has no limit points.
To define general point processes, we start with a probability space ( Ω , F , P ) {\displaystyle (\Omega ,{\mathcal {F}},P)} , and a measurable space ( S , S ) {\displaystyle (S,{\mathcal {S}})} where S {\displaystyle S} is a locally compact second countable Hausdorff space and S {\displaystyle {\mathcal {S}}} is its Borel σ-algebra. Consider now an integer-valued locally finite kernel ξ {\displaystyle \xi } from ( Ω , F ) {\displaystyle (\Omega ,{\mathcal {F}})} into ( S , S ) {\displaystyle (S,{\mathcal {S}})} , that is, a mapping Ω × S ↦ Z + {\displaystyle \Omega \times {\mathcal {S}}\mapsto \mathbb {Z} _{+}} such that:
This kernel defines a random measure in the following way. We would like to think of ξ {\displaystyle \xi } as defining a mapping which maps ω ∈ Ω {\displaystyle \omega \in \Omega } to a measure ξ ω ∈ M ( S ) {\displaystyle \xi _{\omega }\in {\mathcal {M}}({\mathcal {S}})} (namely, Ω ↦ M ( S ) {\displaystyle \Omega \mapsto {\mathcal {M}}({\mathcal {S}})} ), where M ( S ) {\displaystyle {\mathcal {M}}({\mathcal {S}})} is the set of all locally finite measures on S {\displaystyle S} . Now, to make this mapping measurable, we need to define a σ {\displaystyle \sigma } -field over M ( S ) {\displaystyle {\mathcal {M}}({\mathcal {S}})} . This σ {\displaystyle \sigma } -field is constructed as the minimal algebra so that all evaluation maps of the form π B : μ ↦ μ ( B ) {\displaystyle \pi _{B}:\mu \mapsto \mu (B)} , where B ∈ S {\displaystyle B\in {\mathcal {S}}} is relatively compact, are measurable. Equipped with this σ {\displaystyle \sigma } -field, then ξ {\displaystyle \xi } is a random element, where for every ω ∈ Ω {\displaystyle \omega \in \Omega } , ξ ω {\displaystyle \xi _{\omega }} is a locally finite measure over S {\displaystyle S} .
Now, by a point process on S {\displaystyle S} we simply mean an integer-valued random measure (or equivalently, integer-valued kernel) ξ {\displaystyle \xi } constructed as above. The most common example for the state space S is the Euclidean space Rn or a subset thereof, where a particularly interesting special case is given by the real half-line [0,∞). However, point processes are not limited to these examples and may among other things also be used if the points are themselves compact subsets of Rn, in which case ξ is usually referred to as a particle process.
Despite the name point process since S might not be a subset of the real line, as it might suggest that ξ is a stochastic process.
Every instance (or event) of a point process ξ can be represented as
where δ {\displaystyle \delta } denotes the Dirac measure, n is an integer-valued random variable and X i {\displaystyle X_{i}} are random elements of S. If X i {\displaystyle X_{i}} 's are almost surely distinct (or equivalently, almost surely ξ ( x ) ≤ 1 {\displaystyle \xi (x)\leq 1} for all x ∈ R d {\displaystyle x\in \mathbb {R} ^{d}} ), then the point process is known as simple.
Another different but useful representation of an event (an event in the event space, i.e. a series of points) is the counting notation, where each instance is represented as an N ( t ) {\displaystyle N(t)} function, a continuous function which takes integer values: N : R → Z 0 + {\displaystyle N:{\mathbb {R} }\rightarrow {\mathbb {Z} _{0}^{+}}} :
which is the number of events in the observation interval ( t 1 , t 2 ] {\displaystyle (t_{1},t_{2}]} . It is sometimes denoted by N t 1 , t 2 {\displaystyle N_{t_{1},t_{2}}} , and N T {\displaystyle N_{T}} or N ( T ) {\displaystyle N(T)} mean N 0 , T {\displaystyle N_{0,T}} .
Main article: Intensity measure
The expectation measure Eξ (also known as mean measure) of a point process ξ is a measure on S that assigns to every Borel subset B of S the expected number of points of ξ in B. That is,
The Laplace functional Ψ N ( f ) {\displaystyle \Psi _{N}(f)} of a point process N is a map from the set of all positive valued functions f on the state space of N, to [ 0 , ∞ ) {\displaystyle [0,\infty )} defined as follows:
They play a similar role as the characteristic functions for random variable. One important theorem says that: two point processes have the same law if their Laplace functionals are equal.
Main article: Moment measure
The n {\displaystyle n} th power of a point process, ξ n , {\displaystyle \xi ^{n},} is defined on the product space S n {\displaystyle S^{n}} as follows :
By monotone class theorem, this uniquely defines the product measure on ( S n , B ( S n ) ) . {\displaystyle (S^{n},B(S^{n})).} The expectation E ξ n ( ⋅ ) {\displaystyle E\xi ^{n}(\cdot )} is called the n {\displaystyle n} th moment measure. The first moment measure is the mean measure.
Let S = R d {\displaystyle S=\mathbb {R} ^{d}} . The joint intensities of a point process ξ {\displaystyle \xi } w.r.t. the Lebesgue measure are functions ρ ( k ) : ( R d ) k → [ 0 , ∞ ) {\displaystyle \rho ^{(k)}:(\mathbb {R} ^{d})^{k}\to [0,\infty )} such that for any disjoint bounded Borel subsets B 1 , … , B k {\displaystyle B_{1},\ldots ,B_{k}}
Joint intensities do not always exist for point processes. Given that moments of a random variable determine the random variable in many cases, a similar result is to be expected for joint intensities. Indeed, this has been shown in many cases.20
A point process ξ ⊂ R d {\displaystyle \xi \subset \mathbb {R} ^{d}} is said to be stationary if ξ + x := ∑ i = 1 N δ X i + x {\displaystyle \xi +x:=\sum _{i=1}^{N}\delta _{X_{i}+x}} has the same distribution as ξ {\displaystyle \xi } for all x ∈ R d . {\displaystyle x\in \mathbb {R} ^{d}.} For a stationary point process, the mean measure E ξ ( ⋅ ) = λ ‖ ⋅ ‖ {\displaystyle E\xi (\cdot )=\lambda \|\cdot \|} for some constant λ ≥ 0 {\displaystyle \lambda \geq 0} and where ‖ ⋅ ‖ {\displaystyle \|\cdot \|} stands for the Lebesgue measure. This λ {\displaystyle \lambda } is called the intensity of the point process. A stationary point process on R d {\displaystyle \mathbb {R} ^{d}} has almost surely either 0 or an infinite number of points in total. For more on stationary point processes and random measure, refer to Chapter 12 of Daley & Vere-Jones.21 Stationarity has been defined and studied for point processes in more general spaces than R d {\displaystyle \mathbb {R} ^{d}} .
Main article: Point process operation
A point process transformation is a function that maps a point process to another point process.
We shall see some examples of point processes in R d . {\displaystyle \mathbb {R} ^{d}.}
Main article: Poisson point process
The simplest and most ubiquitous example of a point process is the Poisson point process, which is a spatial generalisation of the Poisson process. A Poisson (counting) process on the line can be characterised by two properties : the number of points (or events) in disjoint intervals are independent and have a Poisson distribution. A Poisson point process can also be defined using these two properties. Namely, we say that a point process ξ {\displaystyle \xi } is a Poisson point process if the following two conditions hold
1) ξ ( B 1 ) , … , ξ ( B n ) {\displaystyle \xi (B_{1}),\ldots ,\xi (B_{n})} are independent for disjoint subsets B 1 , … , B n . {\displaystyle B_{1},\ldots ,B_{n}.}
2) For any bounded subset B {\displaystyle B} , ξ ( B ) {\displaystyle \xi (B)} has a Poisson distribution with parameter λ ‖ B ‖ , {\displaystyle \lambda \|B\|,} where ‖ ⋅ ‖ {\displaystyle \|\cdot \|} denotes the Lebesgue measure.
The two conditions can be combined and written as follows : For any disjoint bounded subsets B 1 , … , B n {\displaystyle B_{1},\ldots ,B_{n}} and non-negative integers k 1 , … , k n {\displaystyle k_{1},\ldots ,k_{n}} we have that
The constant λ {\displaystyle \lambda } is called the intensity of the Poisson point process. Note that the Poisson point process is characterised by the single parameter λ . {\displaystyle \lambda .} It is a simple, stationary point process. To be more specific one calls the above point process a homogeneous Poisson point process. An inhomogeneous Poisson process is defined as above but by replacing λ ‖ B ‖ {\displaystyle \lambda \|B\|} with ∫ B λ ( x ) d x {\displaystyle \int _{B}\lambda (x)\,dx} where λ {\displaystyle \lambda } is a non-negative function on R d . {\displaystyle \mathbb {R} ^{d}.}
A Cox process (named after Sir David Cox) is a generalisation of the Poisson point process, in that we use random measures in place of λ ‖ B ‖ {\displaystyle \lambda \|B\|} . More formally, let Λ {\displaystyle \Lambda } be a random measure. A Cox point process driven by the random measure Λ {\displaystyle \Lambda } is the point process ξ {\displaystyle \xi } with the following two properties :
It is easy to see that Poisson point process (homogeneous and inhomogeneous) follow as special cases of Cox point processes. The mean measure of a Cox point process is E ξ ( ⋅ ) = E Λ ( ⋅ ) {\displaystyle E\xi (\cdot )=E\Lambda (\cdot )} and thus in the special case of a Poisson point process, it is λ ‖ ⋅ ‖ . {\displaystyle \lambda \|\cdot \|.}
For a Cox point process, Λ ( ⋅ ) {\displaystyle \Lambda (\cdot )} is called the intensity measure. Further, if Λ ( ⋅ ) {\displaystyle \Lambda (\cdot )} has a (random) density (Radon–Nikodym derivative) λ ( ⋅ ) {\displaystyle \lambda (\cdot )} i.e.,
then λ ( ⋅ ) {\displaystyle \lambda (\cdot )} is called the intensity field of the Cox point process. Stationarity of the intensity measures or intensity fields imply the stationarity of the corresponding Cox point processes.
There have been many specific classes of Cox point processes that have been studied in detail such as:
By Jensen's inequality, one can verify that Cox point processes satisfy the following inequality: for all bounded Borel subsets B {\displaystyle B} ,
where ξ α {\displaystyle \xi _{\alpha }} stands for a Poisson point process with intensity measure α ( ⋅ ) := E ξ ( ⋅ ) = E Λ ( ⋅ ) . {\displaystyle \alpha (\cdot ):=E\xi (\cdot )=E\Lambda (\cdot ).} Thus points are distributed with greater variability in a Cox point process compared to a Poisson point process. This is sometimes called clustering or attractive property of the Cox point process.
An important class of point processes, with applications to physics, random matrix theory, and combinatorics, is that of determinantal point processes.28
Main article: Hawkes process
A Hawkes process N t {\displaystyle N_{t}} , also known as a self-exciting counting process, is a simple point process whose conditional intensity can be expressed as
where ν : R + → R + {\displaystyle \nu :\mathbb {R} ^{+}\rightarrow \mathbb {R} ^{+}} is a kernel function which expresses the positive influence of past events T i {\displaystyle T_{i}} on the current value of the intensity process λ ( t ) {\displaystyle \lambda (t)} , μ ( t ) {\displaystyle \mu (t)} is a possibly non-stationary function representing the expected, predictable, or deterministic part of the intensity, and { T i : T i < T i + 1 } ∈ R {\displaystyle \{T_{i}:T_{i}<T_{i+1}\}\in \mathbb {R} } is the time of occurrence of the i-th event of the process.29
Given a sequence of non-negative random variables { X k , k = 1 , 2 , … } {\textstyle \{X_{k},k=1,2,\dots \}} , if they are independent and the cdf of X k {\displaystyle X_{k}} is given by F ( a k − 1 x ) {\displaystyle F(a^{k-1}x)} for k = 1 , 2 , … {\displaystyle k=1,2,\dots } , where a {\displaystyle a} is a positive constant, then { X k , k = 1 , 2 , … } {\displaystyle \{X_{k},k=1,2,\ldots \}} is called a geometric process (GP).30
The geometric process has several extensions, including the α- series process31 and the doubly geometric process.32
Historically the first point processes that were studied had the real half line R+ = [0,∞) as their state space, which in this context is usually interpreted as time. These studies were motivated by the wish to model telecommunication systems,33 in which the points represented events in time, such as calls to a telephone exchange.
Point processes on R+ are typically described by giving the sequence of their (random) inter-event times (T1, T2, ...), from which the actual sequence (X1, X2, ...) of event times can be obtained as
If the inter-event times are independent and identically distributed, the point process obtained is called a renewal process.
The intensity λ(t | Ht) of a point process on the real half-line with respect to a filtration Ht is defined as
Ht can denote the history of event-point times preceding time t but can also correspond to other filtrations (for example in the case of a Cox process).
In the N ( t ) {\displaystyle N(t)} -notation, this can be written in a more compact form:
The compensator of a point process, also known as the dual-predictable projection, is the integrated conditional intensity function defined by
The Papangelou intensity function of a point process N {\displaystyle N} in the n {\displaystyle n} -dimensional Euclidean space R n {\displaystyle \mathbb {R} ^{n}} is defined as
where B δ ( x ) {\displaystyle B_{\delta }(x)} is the ball centered at x {\displaystyle x} of a radius δ {\displaystyle \delta } , and σ [ N ( R n ∖ B δ ( x ) ) ] {\displaystyle \sigma [N(\mathbb {R} ^{n}\setminus B_{\delta }(x))]} denotes the information of the point process N {\displaystyle N} outside B δ ( x ) {\displaystyle B_{\delta }(x)} .
The logarithmic likelihood of a parameterized simple point process conditional upon some observed data is written as
The analysis of point pattern data in a compact subset S of Rn is a major object of study within spatial statistics. Such data appear in a broad range of disciplines,35 amongst which are
The need to use point processes to model these kinds of data lies in their inherent spatial structure. Accordingly, a first question of interest is often whether the given data exhibit complete spatial randomness (i.e. are a realization of a spatial Poisson process) as opposed to exhibiting either spatial aggregation or spatial inhibition.
In contrast, many datasets considered in classical multivariate statistics consist of independently generated datapoints that may be governed by one or several covariates (typically non-spatial).
Apart from the applications in spatial statistics, point processes are one of the fundamental objects in stochastic geometry. Research has also focussed extensively on various models built on point processes such as Voronoi tessellations, random geometric graphs, and Boolean models.
Kallenberg, O. (1986). Random Measures, 4th edition. Academic Press, New York, London; Akademie-Verlag, Berlin. ISBN 0-12-394960-2, MR854102. /wiki/Olav_Kallenberg ↩
Daley, D.J, Vere-Jones, D. (1988). An Introduction to the Theory of Point Processes. Springer, New York. ISBN 0-387-96666-8, MR950166. /wiki/ISBN_(identifier) ↩
Last, G., Brandt, A. (1995).Marked point processes on the real line: The dynamic approach. Probability and its Applications. Springer, New York. ISBN 0-387-94547-4, MR1353912 /wiki/ISBN_(identifier) ↩
Gilbert E.N. (1961). "Random plane networks". Journal of the Society for Industrial and Applied Mathematics. 9 (4): 533–543. doi:10.1137/0109045. /wiki/Edgar_N._Gilbert ↩
Diggle, P. (2003). Statistical Analysis of Spatial Point Patterns, 2nd edition. Arnold, London. ISBN 0-340-74070-1. /wiki/ISBN_(identifier) ↩
Baddeley, A. (2006). Spatial point processes and their applications. In A. Baddeley, I. Bárány, R. Schneider, and W. Weil, editors, Stochastic Geometry: Lectures given at the C.I.M.E. Summer School held in Martina Franca, Italy, September 13–18, 2004, Lecture Notes in Mathematics 1892, Springer. ISBN 3-540-38174-0, pp. 1–75 /wiki/ISBN_(identifier) ↩
Brown E. N., Kass R. E., Mitra P. P. (2004). "Multiple neural spike train data analysis: state-of-the-art and future challenges". Nature Neuroscience. 7 (5): 456–461. doi:10.1038/nn1228. PMID 15114358. S2CID 562815.{{cite journal}}: CS1 maint: multiple names: authors list (link) /wiki/Doi_(identifier) ↩
Engle Robert F., Lunde Asger (2003). "Trades and Quotes: A Bivariate Point Process" (PDF). Journal of Financial Econometrics. 1 (2): 159–188. doi:10.1093/jjfinec/nbg011. https://escholarship.org/content/qt8bh079sq/qt8bh079sq.pdf?t=li5awc ↩
Sung Nok Chiu; Dietrich Stoyan; Wilfrid S. Kendall; Joseph Mecke (27 June 2013). Stochastic Geometry and Its Applications. John Wiley & Sons. p. 108. ISBN 978-1-118-65825-3. 978-1-118-65825-3 ↩
Martin Haenggi (2013). Stochastic Geometry for Wireless Networks. Cambridge University Press. p. 10. ISBN 978-1-107-01469-5. 978-1-107-01469-5 ↩
D.J. Daley; D. Vere-Jones (10 April 2006). An Introduction to the Theory of Point Processes: Volume I: Elementary Theory and Methods. Springer Science & Business Media. p. 194. ISBN 978-0-387-21564-8. 978-0-387-21564-8 ↩
Cox, D. R.; Isham, Valerie (1980). Point Processes. CRC Press. p. 3. ISBN 978-0-412-21910-8. 978-0-412-21910-8 ↩
In the context of point processes, the term "state space" can mean the space on which the point process is defined such as the real line,[13][14] which corresponds to the index set in stochastic process terminology. ↩
Samuel Karlin; Howard E. Taylor (2 December 2012). A First Course in Stochastic Processes. Academic Press. p. 31. ISBN 978-0-08-057041-9. 978-0-08-057041-9 ↩
Volker Schmidt (24 October 2014). Stochastic Geometry, Spatial Statistics and Random Fields: Models and Algorithms. Springer. p. 99. ISBN 978-3-319-10064-7. 978-3-319-10064-7 ↩
D.J. Daley; D. Vere-Jones (10 April 2006). An Introduction to the Theory of Point Processes: Volume I: Elementary Theory and Methods. Springer Science & Business Media. ISBN 978-0-387-21564-8. 978-0-387-21564-8 ↩
Sung Nok Chiu; Dietrich Stoyan; Wilfrid S. Kendall; Joseph Mecke (27 June 2013). Stochastic Geometry and Its Applications. John Wiley & Sons. p. 109. ISBN 978-1-118-65825-3. 978-1-118-65825-3 ↩
Moller, J.; Syversveen, A. R.; Waagepetersen, R. P. (1998). "Log Gaussian Cox Processes". Scandinavian Journal of Statistics. 25 (3): 451. CiteSeerX 10.1.1.71.6732. doi:10.1111/1467-9469.00115. S2CID 120543073. /wiki/CiteSeerX_(identifier) ↩
Moller, J. (2003) Shot noise Cox processes, Adv. Appl. Prob., 35.[page needed] /wiki/Wikipedia:Citing_sources ↩
Moller, J. and Torrisi, G.L. (2005) "Generalised Shot noise Cox processes", Adv. Appl. Prob., 37. ↩
Hellmund, G., Prokesova, M. and Vedel Jensen, E.B. (2008) "Lévy-based Cox point processes", Adv. Appl. Prob., 40. [page needed] /wiki/Eva_Vedel_Jensen ↩
Mccullagh,P. and Moller, J. (2006) "The permanental processes", Adv. Appl. Prob., 38.[page needed] /wiki/Wikipedia:Citing_sources ↩
Adams, R. P., Murray, I. MacKay, D. J. C. (2009) "Tractable inference in Poisson processes with Gaussian process intensities", Proceedings of the 26th International Conference on Machine Learning doi:10.1145/1553374.1553376 /wiki/Doi_(identifier) ↩
Hough, J. B., Krishnapur, M., Peres, Y., and Virág, B., Zeros of Gaussian analytic functions and determinantal point processes. University Lecture Series, 51. American Mathematical Society, Providence, RI, 2009. ↩
Patrick J. Laub, Young Lee, Thomas Taimre, The Elements of Hawkes Processes, Springer, 2022. ↩
Lin, Ye (Lam Yeh) (1988). "Geometric processes and replacement problem". Acta Mathematicae Applicatae Sinica. 4 (4): 366–377. doi:10.1007/BF02007241. S2CID 123338120. /wiki/Doi_(identifier) ↩
Braun, W. John; Li, Wei; Zhao, Yiqiang Q. (2005). "Properties of the geometric and related processes". Naval Research Logistics. 52 (7): 607–616. CiteSeerX 10.1.1.113.9550. doi:10.1002/nav.20099. S2CID 7745023. /wiki/CiteSeerX_(identifier) ↩
Wu, Shaomin (2018). "Doubly geometric processes and applications" (PDF). Journal of the Operational Research Society. 69: 66–77. doi:10.1057/s41274-017-0217-4. S2CID 51889022. https://kar.kent.ac.uk/60730/1/JORS_Wu.pdf ↩
Palm, C. (1943). Intensitätsschwankungen im Fernsprechverkehr (German). Ericsson Technics no. 44, (1943). MR11402 /wiki/MR_(identifier) ↩
Rubin, I. (Sep 1972). "Regular point processes and their detection". IEEE Transactions on Information Theory. 18 (5): 547–557. doi:10.1109/tit.1972.1054897. /wiki/Doi_(identifier) ↩
Baddeley, A., Gregori, P., Mateu, J., Stoica, R., and Stoyan, D., editors (2006). Case Studies in Spatial Point Pattern Modelling, Lecture Notes in Statistics No. 185. Springer, New York. ISBN 0-387-28311-0. /wiki/ISBN_(identifier) ↩