[R-es] Evaluación del numero de modas en una distribución

[Carlos J. Gil Bellosta]

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