Re: [R] using kernel density estimates to infer mode of distribution

From: Adelchi Azzalini <azzalini_at_stat.unipd.it>
Date: Thu 16 Feb 2006 - 20:28:41 EST

On Wed, 15 Feb 2006 18:28:25 -0500, Dan Rabosky wrote:

DR> 
DR>  Is it possible to use "density" or another kernel density
DR>  estimator to  identify the mode of a distribution?  When I use
DR>  'density', the resulting 

a simple option is of the form

   fit$eval[fit$estimate==max(fit$estimate)] assuming that fit$eval is the vector of evaluation points, and fit$estimate the corrisponding density estimates (this is the sort of output produced by sm.density)

Here I have assumed there is single mode and we are in the scalar case, for simplicity. Some variant required in the more general case.

best regards,

Adelchi Azzalini

-- 
Adelchi Azzalini  <azzalini@stat.unipd.it>
Dipart.Scienze Statistiche, UniversitÓ di Padova, Italia
tel. +39 049 8274147,  http://azzalini.stat.unipd.it/

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Thu Feb 16 20:36:25 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:31 EST