Geyer(r,sat)"interact"
describing the interpoint interaction
structure of Geyer's saturation point process
with interaction radius $r$ and saturation threshold sat.Strauss)
in which the total contribution
to the potential from each point (from its pairwise interaction with all
other points) is trimmed to a maximum value $s$.
This model is implemented in the function Geyer(). The saturation point process with interaction radius $r$,
saturation threshold $s$, and
parameters $\beta$ and $\gamma$,
is the point process
in which each point
$x_i$ in the pattern $X$
contributes a factor
$$\beta \gamma^{\min(s, t(x_i, X))}$$
to the probability density of the point pattern,
where $t(x_i, X)$ denotes the
number of
If the saturation threshold $s$ is set to infinity,
this model reduces to the Strauss process (see Strauss)
with interaction parameter $\gamma^2$.
If $s = 0$, the model reduces to the Poisson point process.
If $s$ is a finite positive number, then the interaction parameter
$\gamma$ may take any positive value (unlike the case
of the Strauss process), with
values $\gamma < 1$
describing an ppm(), 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 saturation process interaction is
yielded by Geyer(r, sat) where the
arguments r and sat specify
the Strauss interaction radius $r$ and the saturation threshold
$s$, respectively. See the examples below.
Note the only arguments are the interaction radius r
and the saturation threshold sat.
When r and sat are fixed,
the model becomes an exponential family.
The canonical parameters $\log(\beta)$
and $\log(\gamma)$
are estimated by ppm(), not fixed in
Geyer().
ppm,
pairwise.family,
ppm.object,
Strauss,
SatPiecedata(cells)
ppm(cells, ~1, Geyer(r=0.07, sat=2))
# fit the stationary saturation process to `cells'Run the code above in your browser using DataLab