Generate a random set by taking a random selection of tiles of a given tessellation.
rMosaicSet(X, p=0.5)A tessellation (object of class "tess").
Probability of including a given tile. A number strictly between 0 and 1.
A window (object of class "owin").
Given a tessellation X, this function
randomly selects some of the tiles of X,
including each tile with probability \(p\) independently of
the other tiles. The selected tiles are then combined to form a
set in the plane.
One application of this is Switzer's (1965) example of a random set
which has a Markov property. It is constructed by generating X
according to a Poisson line tessellation (see rpoislinetess).
Switzer, P. A random set process in the plane with a Markovian property. Annals of Mathematical Statistics 36 (1965) 1859--1863.
# NOT RUN {
# Switzer's random set
X <- rpoislinetess(3)
plot(rMosaicSet(X, 0.5), col="green", border=NA)
# another example
plot(rMosaicSet(dirichlet(runifpoint(30)), 0.4))
# }
Run the code above in your browser using DataCamp Workspace