rmat3: Sample from Matern-III marked point process

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/rmat3.R

Description

Generates a sample from Matern-III marked spatial point process

Usage

1
2
3
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
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.