Simulate Matern Model I
Generate a random point pattern, a simulated realisation of the Matern Model I inhibition process model.
rMaternI(kappa, r, win = owin(c(0,1),c(0,1)), stationary=TRUE, ..., nsim=1, drop=TRUE)
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 |
stationary |
Logical. Whether to start with a stationary process of proposal points
( |
... |
Ignored. |
nsim |
Number of simulated realisations to be generated. |
drop |
Logical. If |
This algorithm generates one or more realisations
of Matern'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 Matern's Model I.
A point pattern
if nsim=1
, or a list of point patterns if nsim > 1
.
Each point pattern is normally an object of class "ppp"
,
but may be of class "pp3"
or "ppx"
depending on the window.
Adrian Baddeley Adrian.Baddeley@curtin.edu.au
, Ute Hahn, Rolf Turner r.turner@auckland.ac.nz
and Ege Rubak rubak@math.aau.dk
X <- rMaternI(20, 0.05) Y <- rMaternI(20, 0.05, stationary=FALSE)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.