Spherical contact distribution function

From HandWiki

In probability and statistics, a spherical contact distribution function, first contact distribution function,[1] or empty space function[2] is a mathematical function 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.[1][3] More specifically, a spherical contact distribution function is defined as probability distribution of the radius of a sphere when it first encounters or makes contact with a point in a point process. This function can be contrasted with the nearest neighbour function, which is defined in relation to some point in the point process as being the probability distribution of the distance from that point to its nearest neighbouring point in the same point process. The spherical contact function is also referred to as the contact distribution function,[2] but some authors[1] define the contact distribution function in relation to a more general set, and not simply a sphere as in the case of the spherical contact distribution function.

Spherical contact distribution functions are used in the study of point processes[2][3][4] as well as the related fields of stochastic geometry[1] and spatial statistics,[2][5] which are applied in various scientific and engineering disciplines such as biology, geology, physics, and telecommunications.[1][3][6][7]

Point process notation

Main page: 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 [math]\displaystyle{ \textstyle \textbf{R}^{ d} }[/math], but they can be defined on more abstract mathematical spaces.[4]

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

[math]\displaystyle{ \textstyle x\in {N}, }[/math]

and represents the point process being interpreted as a random set. Alternatively, the number of points of [math]\displaystyle{ \textstyle {N} }[/math] located in some Borel set [math]\displaystyle{ \textstyle B }[/math] is often written as:[1][5][6]

[math]\displaystyle{ \textstyle {N}(B), }[/math]

which reflects a random measure interpretation for point processes. These two notations are often used in parallel or interchangeably.[1][5][6]

Definitions

Spherical contact distribution function

The spherical contact distribution function is defined as:

[math]\displaystyle{ H_s(r)=1-P({N}(b(o,r))=0). }[/math]

where b(o,r) is a ball with radius r centered at the origin o. In other words, spherical contact distribution function is the probability there are no points from the point process located in a hyper-sphere of radius r.

Contact distribution function

The spherical contact distribution function can be generalized for sets other than the (hyper-)sphere in [math]\displaystyle{ \textstyle \textbf{R}^{ d} }[/math]. For some Borel set [math]\displaystyle{ \textstyle B }[/math] with positive volume (or more specifically, Lebesgue measure), the contact distribution function (with respect to [math]\displaystyle{ \textstyle B }[/math]) for [math]\displaystyle{ \textstyle r\geq0 }[/math] is defined by the equation:[1]

[math]\displaystyle{ H_B(r)=P({N}(rB)=0). }[/math]

Examples

Poisson point process

For a Poisson point process [math]\displaystyle{ \textstyle {N} }[/math] on [math]\displaystyle{ \textstyle \textbf{R}^d }[/math] with intensity measure [math]\displaystyle{ \textstyle \Lambda }[/math] this becomes

[math]\displaystyle{ H_s(r)=1-e^{-\Lambda(b(o,r))}, }[/math]

which for the homogeneous case becomes

[math]\displaystyle{ H_s(r)=1-e^{-\lambda |b(o,r)|}, }[/math]

where [math]\displaystyle{ \textstyle |b(o,r)| }[/math] denotes the volume (or more specifically, the Lebesgue measure) of the ball of radius [math]\displaystyle{ \textstyle r }[/math]. In the plane [math]\displaystyle{ \textstyle \textbf{R}^2 }[/math], this expression simplifies to

[math]\displaystyle{ H_s(r)=1-e^{-\lambda \pi r^2}. }[/math]

Relationship to other functions

Nearest neighbour function

In general, the spherical contact distribution function and the corresponding nearest neighbour function are not equal. However, these two functions are identical for Poisson point processes.[1] In fact, this characteristic is due to a unique property of Poisson processes and their Palm distributions, which forms part of the result known as the Slivnyak-Mecke[6] or Slivnyak's theorem.[2]

J-function

The fact that the spherical distribution function Hs(r) and nearest neighbour function Do(r) are identical for the Poisson point process can be used to statistically test if point process data appears to be that of a Poisson point process. For example, in spatial statistics the J-function is defined for all r ≥ 0 as:[1]

[math]\displaystyle{ J(r)=\frac{1-D_o(r)}{1-H_s(r)} }[/math]

For a Poisson point process, the J function is simply J(r)=1, hence why it is used as a non-parametric test for whether data behaves as though it were from a Poisson process. It is, however, thought possible to construct non-Poisson point processes for which J(r)=1,[8] but such counterexamples are viewed as somewhat 'artificial' by some and exist for other statistical tests.[9]

More generally, J-function serves as one way (others include using factorial moment measures[2]) to measure the interaction between points in a point process.[1]

See also

References

  1. 1.00 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.10 1.11 1.12 D. Stoyan, W. S. Kendall, J. Mecke, and L. Ruschendorf. Stochastic geometry and its applications, edition 2. Wiley Chichester, 1995.
  2. 2.0 2.1 2.2 2.3 2.4 2.5 A. Baddeley, I. Bárá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.
  3. 3.0 3.1 3.2 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.
  4. 4.0 4.1 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.
  5. 5.0 5.1 5.2 J. Moller and R. P. Waagepetersen. Statistical inference and simulation for spatial point processes. CRC Press, 2003.
  6. 6.0 6.1 6.2 6.3 F. Baccelli and B. Błaszczyszyn. Stochastic Geometry and Wireless Networks, Volume I – Theory, volume 3, No 3-4 of Foundations and Trends in Networking. NoW Publishers, 2009.
  7. 7.0 7.1 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.
  8. Bedford, T, Van den Berg, J (1997). "A remark on the Van Lieshout and Baddeley J-function for point processes". Advances in Applied Probability (JSTOR) 29 (1): 19–25. doi:10.2307/1427858. https://ir.cwi.nl/pub/1373. 
  9. Foxall, Rob, Baddeley, Adrian (2002). "Nonparametric measures of association between a spatial point process and a random set, with geological applications". Journal of the Royal Statistical Society, Series C (Wiley Online Library) 51 (2): 165–182. doi:10.1111/1467-9876.00261.