asselin function

The Asselin de Beauville mode estimator

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

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

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.

Returns

A numeric value is returned, the mode estimate.

Note

The user may call asselin through mlv(x, method = "asselin", ...).

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)
mlv(x, method = "asselin")

References

• Asselin de Beauville J.-P. (1978). Estimation non parametrique de la densite et du mode, exemple de la distribution Gamma. Revue de Statistique Appliquee, 26 (3):47-70.

See Also

mlv for general mode estimation.