NScluster (version 1.1.0)

PalmTypeB: Non-Parametric Estimate of The Palm Intensity of Type B Model

Description

Calculate the non-parametric Palm intensity function of Type B Model estimated directly from a set of point pattern data.

Usage

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

Arguments

offspring

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

pa

the parameters $(mu_i, nu_i, a_i, sigma_{1,i}, sigma_{2,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)
  mu1 <- 10; mu2 <- 40; nu <- 30; sig1 <- 0.01; sig2 <- 0.03
  ty <- 1
  z <- SimulateTypeB(seeds, c(mu1,mu2,nu,sig1,sig2 ), ty, 100, 150)

  ## estimation
  delta <- 0.001
  pa1 <- c(50, 30, 0.2, 0.01, 0.03)
  pa2 <- c(0.44100e+02, 0.30691e+02, 0.15119, 0.94494e-02, 0.26822e-01)
  pa <- matrix(c(pa1,pa2), 2, 5, byrow=TRUE) 
  PalmTypeB(z$offspring, pa, delta, ty)