Kernel Estimate Intensity at non-Data locations

Usage

intensity_somewhere(x, loc, bw, kernel = "gaussian", border = "local")

Arguments

x: point pattern
loc: coordinate matrix for non-data locations at which to estimate the intensity.
bw: bandwidth. Gaussian sd=bw, Epanechnicov domain [-bw, bw].
kernel: Either 'gaussian' or 'epanenchnikov'
border: Border correction to apply. One of 'none', 'local', 'global', 'toroidal'.
normalise: renormalise so that inverse sum = volume
loo: leave-one-out -estimate?

Details

For border correction, the bounding box of the coordinates of x will be used, so at the moment it works only for axis-aligned cuboids.

The kernel will be a product of one dimensional kernels, so Epanechnikov in many dimensions will not be truly isotropic but prefers diagonal directions.