rgtp function

Simulation of generalized Thomas process

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 )

Arguments

  • 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.

Returns

A list(X, C), where X is Generalized Thomas process, and C is Process of cluster centers for Generalized Thomas process.

Examples

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)