R help es - Dec 2009 - Evaluación del numero de modas en una distribución

Hola a todos,

Estoy intentando encontrar algún paquete que me permita chequear el número
de modas de una distribución, aplicando por ejemplo el test de Li(1996). En
la actualidad utilizo el paquete sm para testear la igualdad de densidades
kernel, y salvo error no permite lo que necesito.

Alguien puede ayudarme.

Muchas gracias y saludos
Justo

	[[alternative HTML version deleted]]

Hola, ¿qué tal?

Puedes usar la función "peaks" de

http://finzi.psych.upenn.edu/R/Rhelp02a/archive/33097.html

Por ejemplo,

d <- density(faithful$eruptions, bw = "sj")
sum( peaks( d$y ) )

El resultado es 2 porque la distribución "es" bimodal.

Entrecomillo "es" porque el número de "modas" depende del
parámetro "bw"
de "density": si éste es muy pequeño, el número de modas será igual al
de valores distintos de la distribución original...

Un saludo,

Carlos J. Gil Bellosta
http://www.datanalytics.com


Justo de Jorge Moreno wrote:
> Hola a todos,
> 
> Estoy intentando encontrar algún paquete que me permita chequear el número
> de modas de una distribución, aplicando por ejemplo el test de Li(1996). En
> la actualidad utilizo el paquete sm para testear la igualdad de densidades
> kernel, y salvo error no permite lo que necesito.
> 
> Alguien puede ayudarme.
> 
> Muchas gracias y saludos
> Justo
> 
> 	[[alternative HTML version deleted]]
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es