dens.par: Non-separable adaptive spatio-temporal intensity estimator

In kernstadapt: Adaptive Kernel Estimators for Point Process Intensities on Linear Networks

Non-separable adaptive spatio-temporal intensity estimator

Description

Provides an adaptive-bandwidth kernel estimate for spatio-temporal point patterns in a non-separable fashion by using binning of the bandwidth values.

Usage

dens.par(
  X,
  t = NULL,
  dimyx = 128,
  dimt = 128,
  bw.xy = NULL,
  bw.t = NULL,
  ngroups.xy = NULL,
  ngroups.t = NULL,
  at = c("bins", "points")
)

Arguments

X	A spatial point pattern (an object of class ppp) with the spatial coordinates of the observations. It may contain marks representing times.
t	A numeric vector of temporal coordinates with equal length to the number of points in X. This gives the time associated with each spatial point. This argument is not necessary if time marks are provided to the point pattern X.
dimyx	Spatial pixel resolution. The default is 128 for each axes.
dimt	Temporal bin vector dimension. The default is 128.
bw.xy	Numeric vector of spatial smoothing bandwidths for each point in X. By default this is computed using bw.abram.
bw.t	Numeric vector of temporal smoothing bandwidths for each point in t. By default this is computed using bw.abram.temp.
ngroups.xy	Number of groups in which the spatial bandwidths should be partitioned. If this number is 1, then a classical non-adaptive estimator will be used for the spatial part with a bandwidth selected as the median of the bw.xy vector.
ngroups.t	Number of groups in which the temporal bandwidths should be partitioned. If this number is 1, then a classical non-adaptive estimator will be used for the temporal part with a bandwidth selected as the median of the bw.t vector.
at	String specifying whether to estimate the intensity at a mesh (at = "bins") or only at the points of X (at = "points").

Details

This function computes a non-separable spatio-temporal adaptive kernel estimate of the intensity. It starts from a planar point pattern X and a vector of times t and partition (cells) the spatial and temporal components to apply a non-separable kernel estimator within each cell. The arguments bw.xy and bw.t specify the smoothing bandwidth vectors to be applied to each of the points in X and t. They should be a numeric vectors of bandwidths. The method partition the range of bandwidths into intervals, subdividing the points of the pattern X and t into sub-patterns according to the bandwidths, and applying fixed-bandwidth smoothing to each sub-pattern. Specifying ngroups.xy = 1 is the same as fixed-bandwidth smoothing with bandwidth sigma = median(bw.xy) in the spatial case and ngroups.t = 1 is the same as fixed-bandwidth smoothing with bandwidth sigma = median(bw.xy).

Value

If at = "points" (the default), the result is a numeric vector with one entry for each data point in X. if at = "bins" is a list named (by time-point) list of pixel images (im objects) corresponding to the joint spatio-temporal intensity over space at each discretised time bin.

Author(s)

Jonatan A. González

References

González J.A. and Moraga P. (2018) An adaptive kernel estimator for the intensity function of spatio-temporal point processes http://arxiv.org/pdf/2208.12026

Examples

data(lGCpp)
stIntensity <- dens.par(lGCpp, dimt = 16)
plot(spatstat.geom::as.solist(stIntensity[13:16]), ncols = 4,
     main = 'Non-separable Example', equal.ribbon = TRUE)

kernstadapt documentation built on Sept. 30, 2024, 9:44 a.m.