Simulation of generalized Thomas process
Simulation of generalized Thomas process.
rgtp( kappa, omega, lambda, theta, win = owin(c(0, 1), c(0, 1)), nsim = 1, expand = 4 * omega )
kappa
: intensity of cluster centers.
omega
: standard deviation of normal distribution specifying the clusters spread.
lambda
: parameter of generalised Poisson distribution controlling over or under dispersion.
theta
: parameter of generalised Poisson distribution controlling the mean number of points in a cluster.
win
: window in which to simulate the pattern. An object in the
spatstat.geom::owin()
format of the spatstat
package.
nsim
: number of simulations.
expand
: the size of expansion of window to simulate the centers of clusters.
A list(X, C), where X is Generalized Thomas process, and C is Process of cluster centers for Generalized Thomas process.
library(spatstat) kappa = 10 omega = .1 lambda= .5 theta = 10 X = rgtp(kappa, omega, lambda, theta, win = owin(c(0, 1), c(0, 1))) plot(X$X) plot(X$C)