[R-es] Evaluación del numero de modas en una distribución
Carlos J. Gil Bellosta
cgb en datanalytics.com
Mar Dic 8 15:48:04 CET 2009
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
Más información sobre la lista de distribución R-help-es