The SPP Procedure

Nonparametric Intensity Estimation

The KERNEL option in the PROCESS statement enables you to perform nonparametric intensity estimation. You can use five different kernel types: Epanechnikov, Gaussian, uniform, triangular, and quartic (Silverman 1986), whose kernel functions are as follows, where , , are the grid point coordinates, x and y are the point coordinates, and h is the bandwidth parameter:

Epanechnikov
Gaussian
uniform
triangular
quartic

Given the preceding kernel definitions, the nonparametric intensity estimate can be computed as

where h is the fixed bandwidth. In practice, nonparametric intensity estimation also involves an edge correction. By default, PROC SPP divides the nonparametric estimate by an edge correction factor

where A is the study area. The choice of the bandwidth parameter that nonparametric intensity estimation requires is more important than the choice of the kernel type itself (Silverman 1986). The bandwidth can be spatially fixed or spatial varying. If the bandwidth is spatially varying, it is called adaptive kernel estimation. For adaptive kernel estimation, the SPP procedure uses the technique suggested in Silverman (1986, p. 101) and Diggle, Rowlingson, and Su (2005, p. 426), which is computed in two steps:

Use an initial bandwidth h to compute pilot estimates of the first-order intensity as
where K(.) is a kernel.
Compute bandwidth factors as
where is the geometric mean of the pilot estimates .

Based on the computed bandwidth estimates, , the nonparametric intensity estimates are computed as

In PROC SPP, adaptive kernel estimation does not incorporate edge correction.