PairPiece

0th

Percentile

The Piecewise Constant Pairwise Interaction Point Process Model

Creates an instance of a pairwise interaction point process model with piecewise constant potential function. The model can then be fitted to point pattern data.

Keywords
spatial
Usage
PairPiece(r)
Arguments
r
vector of discontinuity points for the potential
Details

A pairwise interaction point process in a bounded region is a stochastic point process with probability density of the form $$f(x_1,\ldots,x_n) = \alpha \prod_i b(x_i) \prod_{i < j} h(x_i, x_j)$$ where $x_1,\ldots,x_n$ represent the points of the pattern. The first product on the right hand side is over all points of the pattern; the second product is over all unordered pairs of points of the pattern.

Thus each point $x_i$ of the pattern contributes a factor $b(x_i)$ to the probability density, and each pair of points $x_i, x_j$ contributes a factor $h(x_i,x_j)$ to the density.

The pairwise interaction term $h(u, v)$ is called piecewise constant if it depends only on the distance between $u$ and $v$, say $h(u,v) = H(||u-v||)$, and $H$ is a piecewise constant function (a function which is constant except for jumps at a finite number of places). The function mpl(), which fits point process models to point pattern data, requires an argument of class "interact" describing the interpoint interaction structure of the model to be fitted. The appropriate description of the piecewise constant pairwise interaction is yielded by the function PairPiece(). See the examples below.

The entries of r must be strictly increasing, positive numbers. They are interpreted as the points of discontinuity of $H$. It is assumed that $h(s) =1$ for all $s > r_{max}$ where $r_{max}$ is the maximum value in r. Thus the model has as many regular parameters (see mpl) as there are entries in r. The $i$-th regular parameter $\theta_i$ is the logarithm of the value of the interaction function $H$ on the interval $(r_{i-1},r_i)$.

If r is a single number, this reduces to the Strauss process, see Strauss.

The analogue of this model for multitype point processes has not yet been implemented.

Value

• An object of class "interact" describing the interpoint interaction structure of a point process. The process is a pairwise interaction process, whose interaction potential is piecewise constant, with jumps at the distances given in the vector $r$.

mpl, pairwise.family, ppm.object

• PairPiece
Examples
library(spatstat)
PairPiece(c(0.1,0.2))
# prints a sensible description of itself
data(cells)
mpl(cells, ~1, PairPiece(r = c(0.05, 0.1, 0.2)), rbord=0.2)
# fit a stationary piecewise constant pairwise interaction process
mpl(cells, ~polynom(x,y,3), PairPiece(c(0.05, 0.1)), rbord=0.1)
# nonstationary process with log-cubic polynomial trend
Documentation reproduced from package spatstat, version 1.3-2, License: GPL version 2 or newer

Community examples

Looks like there are no examples yet.