rMaternI {spatstat}R Documentation

Simulate Matern Model I

Description

Generate a random point pattern, a simulated realisation of the Mat\'ern Model I inhibition process model.

Usage

 rMaternI(kappa, r, win = owin(c(0,1),c(0,1)), stationary=TRUE)

Arguments

kappa

Intensity of the Poisson process of proposal points. A single positive number.

r

Inhibition distance.

win

Window in which to simulate the pattern. An object of class "owin" or something acceptable to as.owin.

stationary

Logical. Whether to start with a stationary process of proposal points (stationary=TRUE) or to generate the proposal points only inside the window (stationary=FALSE).

Details

This algorithm generates a realisation of Mat\'ern's Model I inhibition process inside the window win.

The process is constructed by first generating a uniform Poisson point process of “proposal” points with intensity kappa. If stationary = TRUE (the default), the proposal points are generated in a window larger than win that effectively means the proposals are stationary. If stationary=FALSE then the proposal points are only generated inside the window win.

A proposal point is then deleted if it lies within r units' distance of another proposal point. Otherwise it is retained.

The retained points constitute Mat\'ern's Model I.

Value

The simulated point pattern (an object of class "ppp").

Author(s)

Adrian Baddeley Adrian.Baddeley@csiro.au http://www.maths.uwa.edu.au/~adrian/, Ute Hahn, and Rolf Turner r.turner@auckland.ac.nz

See Also

rpoispp, rMatClust

Examples

 X <- rMaternI(20, 0.05)
 Y <- rMaternI(20, 0.05, stationary=FALSE)

[Package spatstat version 1.25-3 Index]