asselin {modeest} | R Documentation |
This mode estimator is based on the algorithm described in Asselin de Beauville (1978).
asselin(x, bw = 1, ...)
x |
numeric. Vector of observations. |
bw |
numeric. A number in (0, 1] . If bw = 1 , the selected 'modal chain' may be too long. |
... |
further arguments to be passed to the quantile function. |
If bw
is missing, then bw = (1:length(x))^(-1/7)
, which is the default value advised by Djeddour et al (2003).
If a
is missing, then a = (1:length(x))^(-alpha)
(with alpha = 0.9
is alpha
is missing), which is the default value advised by Djeddour et al (2003).
A numeric value is returned, the mode estimate.
The user should preferentially call asselin
through
mlv(x, method = "asselin", ...)
.
This returns an object of class mlv
.
Paul Poncet paulponcet@yahoo.fr
mlv
for general mode estimation
x <- rbeta(1000, shape1 = 2, shape2 = 5) ## True mode: betaMode(shape1 = 2, shape2 = 5) ## Estimation: asselin(x, bw = 1) asselin(x, bw = 1/2) M <- mlv(x, method = "asselin") print(M) plot(M)