smacpod (version 1.2.3)

spdensity: Kernel smoothed spatial density of point pattern

Description

spdensity computes a kernel smoothed spatial density function from a point pattern. This is essentially a slight modification of the density.ppp function from the spatstat package, which computes the spatial intensity of a point pattern.

Usage

spdensity(x, sigma = NULL, ..., weights = NULL, edge = TRUE,
  varcov = NULL, at = "pixels", leaveoneout = TRUE, adjust = 1,
  diggle = FALSE)

Arguments

x

Point pattern (object of class "ppp").

sigma

Standard deviation of isotropic Gaussian smoothing kernel. Either a numerical value, or a function that computes an appropriate value of sigma.

...

Additional arguments passed to pixellate.ppp and as.mask to determine the pixel resolution, or passed to sigma if it is a function.

weights

Optional weights to be attached to the points. A numeric vector, numeric matrix, or an expression.

edge

Logical flag: if TRUE, apply edge correction.

varcov

Variance-covariance matrix of anisotropic Gaussian kernel. Incompatible with sigma.

at

String specifying whether to compute the intensity values at a grid of pixel locations (at="pixels") or only at the points of x (at="points").

leaveoneout

Logical value indicating whether to compute a leave-one-out estimator. Applicable only when at="points".

adjust

Optional. Adjustment factor for the smoothing parameter.

diggle

Logical. If TRUE, use Diggle's edge correction, which is more accurate but slower to compute than the correction described under Details.

Value

• This function produces an object of the same type as the density.ppp function from the spatstat package, with the intensity values scaled so that a true spatial density function is produced.

References

Waller, L.A. and Gotway, C.A. (2005). Applied Spatial Statistics for Public Health Data. Hoboken, NJ: Wiley.

Examples

data(grave)
contour(spdensity(grave))