asselin {modeest}R Documentation

The Asselin de Beauville Mode Estimator

Description

This mode estimator is based on the algorithm described in Asselin de Beauville (1978).

Usage

  asselin(x, 
          bw = 1, 
          ...)

Arguments

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.

Details

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).

Value

A numeric value is returned, the mode estimate.

Note

The user should preferentially call asselin through mlv(x, method = "asselin", ...). This returns an object of class mlv.

Author(s)

Paul Poncet paulponcet@yahoo.fr

References

See Also

mlv for general mode estimation

Examples

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)

[Package modeest version 1.09 Index]