Justo de Jorge Moreno
2009-Dec-08 13:42 UTC
[R-es] 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]]
Carlos J. Gil Bellosta
2009-Dec-08 14:48 UTC
[R-es] Evaluación del numero de modas en una distribución
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