rmat3: Sample from Matern-III marked point process

In guiludwig/mat3c: Marked spatial point processes with inhibition based on Matern-III clustered marks

Description

Generates a sample from Matern-III marked spatial point process

Usage

rmat3(beta, phi, gamma, sigma, kappa, win = owin(c(0, 1), c(0, 1)),
  R_clusters = 0.005, R_centers = 0.02, time = 30, death = 1,
  returnTimes = FALSE, returnThinned = FALSE, formatData = NULL, ...)

Arguments

beta	Intensity of the dominant Poisson process, stricitly greater than 0.
phi	Centers' inhibition parameter; rate in which centers are deleted based on how many fingers from other centers surround it. We only allow inhibition, so phi >= 1.
gamma	Average number of fingers generated, must be strictly greater than 0.
sigma	Length (in standard deviations, according to half-Normal distribution) of fingers.
kappa	Concentration parameter from von Mises distribution.
win	Square window in which the process is to be sampled. See owin.
R_clusters	Hardcore radius of inhibition for finger process.
R_centers	Radius of inhibition for center process.
time	Length of time for which a birth-and-death process will be sampled. Defaults to 30.
death	Rate in which the birth-and-death process centers will die. Defaults to 1.
returnTimes	Whether the data should be returned with birth times, or just the sampled spatial locations. Defaults to FALSE. Birth times are assumed unknown in sampled data, so the birth times should be used for debugging only.
formatData	Saves the data as text files in the current directory, with the corresponding file names. Set to NULL to override.
...	further arguments passed to or from other methods.

Value

rmat3 returns an object of class "mat3" containing at least the following components

centers	The spatial point pattern.
fingers	The marks.

Author(s)

Guilherme Ludwig and Nancy Garcia

References

Garcia, N., Guttorp, P. and Ludwig, G. (2020) "Interacting cluster point process model for epidermal nerve fibers", to appear.

See Also

fitmat3

Examples

set.seed(1234)
x <- rmat3(70, 2, 5, 0.05, 3)
plot(x)
z <- sapply(x, function(x) x$centers)
plot(Kest(ppp(z[1,],z[2,])))
plot(envelope(ppp(z[1,],z[2,])))
# Saves data to  files "centers.txt" and "fingers.txt"
## Not run: set.seed(1234)
x <- rmat3(70, 2, 5, 0.05, 3, formatData = c("centers.txt", "fingers.txt"))
# When phi < 1, uses MCMC
x <- rmat3(70, 1/2, 5, 0.05, 3)

## End(Not run)

guiludwig/mat3c documentation built on Dec. 2, 2019, 1:32 a.m.