provenance (version 1.3)

KDE: Create a kernel density estimate

Description

Turns a vector of numbers into an object of class KDE using a combination of the Botev (2010) bandwidth selector and the Abramson (1982) adaptive kernel bandwidth modifier.

Usage

KDE(x, from = NA, to = NA, bw = NA, adaptive = TRUE, log = FALSE, ...)

Arguments

x

a vector of numbers

from

minimum age of the time axis. If NULL, this is set automatically

to

maximum age of the time axis. If NULL, this is set automatically

bw

the bandwidth of the KDE. If NULL, bw will be calculated automatically using botev()

adaptive

boolean flag controlling if the adaptive KDE modifier of Abramson (1982) is used

log

transform the ages to a log scale if TRUE

...

optional arguments to be passed on to density

Value

• an object of class KDE, i.e. a list containing the following items:

  x: horizontal plot coordinates

  y: vertical plot coordinates

  bw: the base bandwidth of the density estimate

  ages: the data values from the input to the KDE function

See Also

KDEs

Examples

data(Namib)
samp <- Namib$DZ$x[['N1']]
dens <- KDE(samp,0,3000,kernel="epanechnikov")
plot(dens)