Density estimation using Locfit
This function provides an interface to Locfit, in the syntax of
(a now old version of) the S-Plus density
function. This can reproduce
density
results, but allows additional
locfit.raw
arguments, such as the degree of fit, to be given.
It also works in double precision, whereas density
only works
in single precision.
density.lf(x, n = 50, window = "gaussian", width, from, to, cut = if(iwindow == 4.) 0.75 else 0.5, ev = lfgrid(mg = n, ll = from, ur = to), deg = 0, family = "density", link = "ident", ...)
x |
numeric vector of observations whose density is to be estimated. |
n |
number of evaluation points.
Equivalent to the |
window |
Window type to use for estimation.
Equivalent to the |
width |
Window width. Following |
from |
Lower limit for estimation domain. |
to |
Upper limit for estimation domain. |
cut |
Controls default expansion of the domain. |
ev |
Locfit evaluation structure – default |
deg |
Fitting degree – default 0 for kernel estimation. |
family |
Fitting family – default is |
link |
Link function – default is the |
... |
Additional arguments to |
A list with components x
(evaluation points) and y
(estimated density).
density
,
locfit
,
locfit.raw
data(geyser) density.lf(geyser, window="tria") # the same result with density, except less precision. density(geyser, window="tria")
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.