bw.abram.temp: Abramson's adaptive temporal bandwidths

bw.abram.tempR Documentation

Abramson's adaptive temporal bandwidths

Description

Computes adaptive smoothing bandwidth in the temporal case according to the inverse-square-root rule of Abramson (1982).

Usage

bw.abram.temp(t, h0 = NULL, nt = 128, trim = NULL, at = "points")

Arguments

t

A vector (a temporal point pattern) from which the bandwidths should be computed.

h0

The global smoothing bandwidth. The default is Silverman's rule of thumb (bw.nrd0).

nt

The number of equally spaced points at which the temporal density is to be estimated.

trim

A trimming value to cut extreme large bandwidths.

at

Character string specifying whether to compute bandwidths at the points (at = "points", the default) or to compute bandwidths at every bin in a bin grid (at = "bins").

Details

This function returns a set of temporal adaptive smoothing bandwidths driven by the methods of Abramson (1982) and Hall and Marron (1988). The bandwidth at location v is given by

h0 * pmin((pilotvalues ^ (-0.5)) / gamma, trim)

where λ^{\mbox{t}}(v) is a pilot estimate of the temporally varying intensity and γ is a scaling constant depending on the pilot estimate.

Value

If at = "points" (the default), the result is a numeric vector with one bandwidth for each data point in X. If at = "bins", the output is an object with class "density" where y component is a vector with the estimated intensity values (see density).

Author(s)

Jonatan A. González

References

Abramson, I. (1982) On bandwidth variation in kernel estimates — a square root law. Annals of Statistics, 10(4), 1217-1223.

Davies, T.M. and Baddeley, A. (2018) Fast computation of spatially adaptive kernel estimates. Statistics and Computing, 28(4), 937-956.

Davies, T.M., Marshall, J.C., and Hazelton, M.L. (2018) Tutorial on kernel estimation of continuous spatial and spatiotemporal relative risk. Statistics in Medicine, 37(7), 1191-1221.

Hall, P. and Marron, J.S. (1988) Variable window width kernel density estimates of probability densities. Probability Theory and Related Fields, 80, 37-49.

Silverman, B.W. (1986) Density Estimation for Statistics and Data Analysis. Chapman and Hall, New York.

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

Examples

t <- 2 * rbeta(100, 1.5, 5.5, 0.2) #Simulated temporal point pattern
bw.abram.temp(t)


kernstadapt documentation built on Nov. 16, 2022, 1:12 a.m.