Kernel Smoothed Intensity of Point Pattern
Compute a kernel smoothed intensity function from a point pattern.
## S3 method for class 'ppp': density(x, sigma, \dots, weights, edge=TRUE, varcov=NULL)
- Point pattern (object of class
"ppp") to be smoothed.
- Standard deviation of isotropic Gaussian smoothing kernel.
- Optional vector of weights to be attached to the points. May include negative values.
- Arguments passed to
as.maskto determine the pixel resolution.
- Logical flag: if
TRUE, apply edge correction.
- Variance-covariance matrix of anisotropic Gaussian kernel.
This is a method for the generic function
A kernel estimate of the intensity function of the point pattern
is computed. The result is
the convolution of the isotropic Gaussian kernel of
sigma with point masses at each of the data
points. The default is to assign
a unit weight to each point.
weights is present, the point masses have these
weights (which may be signed real numbers).
edge=TRUE, the intensity estimate is corrected for
edge effect bias by dividing it by the convolution of the
Gaussian kernel with the window of observation.
Instead of the isotropic Gaussian kernel with standard deviation
sigma, the smoothing kernel may be chosen to be any Gaussian
kernel, by giving the variance-covariance matrix
varcov are incompatible.
sigma may be a vector of length 2 giving the
standard deviations of two independent Gaussian coordinates,
thus equivalent to
varcov = diag(sigma^2).
Computation is performed using the Fast Fourier Transform.
Accuracy depends on the pixel resolution, controlled by the arguments
... passed to
- A pixel image (object of class
data(cells) Z <- density.ppp(cells, 0.05) plot(Z)