Compute Expansion Window for Simulation

Defines the default expansion window for simulation of a fitted point process model.

spatial, datagen
default.expand(object, m=2, epsilon=1e-6)
A fitted point process model (object of class "ppp").
A single numeric value. The window will be expanded by a distance m * reach(object) along each side.
Threshold argument passed to reach.ppm to determine reach(object).

This function computes a default value for the expansion window (the argument expand in rmhcontrol) given a fitted point process model object. This default is used by envelope, qqplot.ppm and other functions.

Suppose we wish to generate simulated realisations of a fitted point process model inside a window w. It is advisable to first simulate the pattern on a larger window, and then clip it to the original window w. This avoids edge effects in the simulation. It is called expansion of the simulation window. Accordingly, for the Metropolis-Hastings simulation algorithm rmh, the algorithm control parameters specified by rmhcontrol include an argument expand that determines the expansion of the simulation window.

The function default.expand determines the default expansion window for a fitted point process model object. If the model depends on external covariates (i.e. covariates other than the Cartesian covariates x and y and the marks) then no expansion is performed, and default.expand returns the original window w = as.owin(object).

Otherwise, the original window w = as.owin(object) is expanded by a distance m * rr, where rr is the interaction range of the model, computed by reach. If w is a rectangle then each edge of w is displaced outward by distance m * rr. If w is not a rectangle then w is dilated by distance m * rr using dilation.


  • A window (object of class "owin").

See Also

rmhcontrol, rmh, envelope, qqplot.ppm

  • default.expand
  fit <- ppm(cells, ~1, Strauss(0.07))
Documentation reproduced from package spatstat, version 1.23-3, License: GPL (>= 2)

Community examples

Looks like there are no examples yet.