NScluster (version 1.1.0)

PalmThomas: Non-Parametric Estimate of The Palm Intensity of The Thomas Model

Description

Calculate the non-parametric Palm intensity function of the Thomas model estimated directly from a set of point pattern data.

Usage

PalmThomas(offspring, pa, delta, ty, plot=TRUE)

Arguments

offspring

the list of coordinates $(x,y)$ of offspring points.

pa

the parameters $(mu_i, nu_i, sigma_i)$, $i=1,2,...,I

delta

a width for the non-parametric Palm intensity function.

ty

the variable $Ty$ for the standardized coordinates of points in the rectangular region $[0,1]*[0,Ty]$ without loss of generality except for the scaling.

plot

logical. If TRUE (default) the non-parametric estimate and the curves of the true and MPLE (the maximum Palm likelihood estimator) parameters are shown.

Value

References

U. Tanaka, Y. Ogata and K. Katsura, Simulation and estimation of the Neyman-Scott type spatial cluster models, Computer Science Monographs No.34, 2008, 1-44. The Institute of Statistical Mathematics.

Examples

  ## simulation
  seeds <- c(822, 913, 905)
  pa1 <- c(50, 30, 0.03)
  ty <- 1
  z <- SimulateThomas(seeds, pa1, ty)

  ## estimation
  delta <- 0.001
  pa2 <- c(0.51311e+02, 0.26104e+02, 0.29071e-01)
  pa <- matrix(c(pa1,pa2), 2, 3, byrow=TRUE) 
  PalmThomas(z$offspring, pa, delta, ty)