Buenas tardes, Estoy utilizando la función hist para calcular los histogramas de seis grupos de datos de la siguiente forma HIST <- hist (misdatos, main = "MG2008", ylab="", xlab ="", cex.axis=2, breaks=9), y quisiera que los seis histogramas tuvieran el mismo número de $counts, es decir, el mismo número de clases, pero no lo logro. Especifico breaks o nclass en la función, pero R termina haciendo el gráfico con el número de clases que él decide. Agradezco mucho su ayuda, Angela C. ---- *Angela Andrea Camargo Sanabria* Estudiante Doctorado en Ciencias Biológicas Laboratorio de Ecología de poblaciones y comunidades tropicales Centro de Investigaciones en Ecosistemas (CIEco) UNAM, campus Morelia Antigua Carretera a Pátzcuaro # 8701 Col. Ex-Hacienda de San José de la Huerta, CP 58190 Morelia, Michoacán, México Tel.: 443-3222706 ext. 42511 e-mail: aacamargo@cieco.unam.mx skype: angela.camargo26 [[alternative HTML version deleted]]
Hola Angela, Aunque no tengo claro el por que quieres algo asi, ni tampoco tengo una respuesta a tu pregunta, presumo que no quieres el _mismo_ numero de counts sino el mismo numero de clases. En principio, la opcion plot = FALSE en hist(), e.g. myh <- hist(x, plot = FALSE) str(myh) parece proporcionar informacion relevante para resolver lo que finalmente quieres hacer. Otras opciones incluyen [1] y [2], pero debes manipular los datos previamente. En [3] hay un ejemplo que puedes adaptar facilmente. Si no hay inconveniente, podrias por favor darnos un poco mas de detalles acerca del por que quieres algo asi? Que intentas mostrar? Quizas estos detalles nos ayuden a ayudarte un poco mas. Un saludo, Jorge Ivan Velez [1] http://had.co.nz/ggplot2/geom_histogram.html [2] http://lmdvr.r-forge.r-project.org/figures/figures.html [3] http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=55 2011/6/14 Angela Andrea Camargo Sanabria <angela.andrea.camargo@gmail.com>> Buenas tardes, > > Estoy utilizando la función hist para calcular los histogramas de seis > grupos de datos de la siguiente forma > HIST <- hist (misdatos, main = "MG2008", ylab="", xlab ="", cex.axis=2, > breaks=9), > y quisiera que los seis histogramas tuvieran el mismo número de $counts, es > decir, el mismo número de clases, pero no lo logro. Especifico breaks o > nclass en la función, pero R termina haciendo el gráfico con el número de > clases que él decide. > > Agradezco mucho su ayuda, > > Angela C. > ---- > *Angela Andrea Camargo Sanabria* > Estudiante Doctorado en Ciencias Biológicas > Laboratorio de Ecología de poblaciones y comunidades tropicales > Centro de Investigaciones en Ecosistemas (CIEco) > UNAM, campus Morelia > Antigua Carretera a Pátzcuaro # 8701 > Col. Ex-Hacienda de San José de la Huerta, CP 58190 > Morelia, Michoacán, México > Tel.: 443-3222706 ext. 42511 > e-mail: aacamargo@cieco.unam.mx > skype: angela.camargo26 > > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >[[alternative HTML version deleted]]
Hola Jorge, Por supuesto no hay inconveniente en explicar, por el contrario mucho mejor! Trataré de hacerme entender. Básicamente quiero comparar si una serie de datos tienen o no la misma distribución. Quiero utilizar el test de KS y según entiendo, requiero la frecuencia o los "counts" para hacer la comparación. Mi idea era calcular los histogramas y guardar el $counts del objeto HIST en otro objeto y hacer las comparaciones posteriormente. Entiendo también que requiere el mismo número de clases para realizar la comparación con este test, de ahí mi inquietud. Gracias a todos, ACS ---- *Angela Andrea Camargo Sanabria* Estudiante Doctorado en Ciencias Biológicas Laboratorio de Ecología de poblaciones y comunidades tropicales Centro de Investigaciones en Ecosistemas (CIEco) UNAM, campus Morelia Antigua Carretera a Pátzcuaro # 8701 Col. Ex-Hacienda de San José de la Huerta, CP 58190 Morelia, Michoacán, México Tel.: 443-3222706 ext. 42511 e-mail: aacamargo@cieco.unam.mx skype: angela.camargo26 2011/6/14 Jorge Ivan Velez <jorgeivanvelez@gmail.com>> Hola Angela, > > Aunque no tengo claro el por que quieres algo asi, ni tampoco tengo una > respuesta a tu pregunta, presumo que no quieres el _mismo_ numero de counts > sino el mismo numero de clases. > > En principio, la opcion plot = FALSE en hist(), e.g. > > myh <- hist(x, plot = FALSE) > str(myh) > > parece proporcionar informacion relevante para resolver lo que finalmente > quieres hacer. Otras opciones incluyen [1] y [2], pero debes manipular los > datos previamente. En [3] hay un ejemplo que puedes adaptar facilmente. > > Si no hay inconveniente, podrias por favor darnos un poco mas de detalles > acerca del por que quieres algo asi? Que intentas mostrar? Quizas estos > detalles nos ayuden a ayudarte un poco mas. > > Un saludo, > Jorge Ivan Velez > > [1] http://had.co.nz/ggplot2/geom_histogram.html > [2] http://lmdvr.r-forge.r-project.org/figures/figures.html > [3] http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=55 > > > > 2011/6/14 Angela Andrea Camargo Sanabria <angela.andrea.camargo@gmail.com> > >> Buenas tardes, >> >> Estoy utilizando la función hist para calcular los histogramas de seis >> grupos de datos de la siguiente forma >> HIST <- hist (misdatos, main = "MG2008", ylab="", xlab ="", cex.axis=2, >> breaks=9), >> y quisiera que los seis histogramas tuvieran el mismo número de $counts, >> es >> decir, el mismo número de clases, pero no lo logro. Especifico breaks o >> nclass en la función, pero R termina haciendo el gráfico con el número de >> clases que él decide. >> >> Agradezco mucho su ayuda, >> >> Angela C. >> ---- >> *Angela Andrea Camargo Sanabria* >> Estudiante Doctorado en Ciencias Biológicas >> Laboratorio de Ecología de poblaciones y comunidades tropicales >> Centro de Investigaciones en Ecosistemas (CIEco) >> UNAM, campus Morelia >> Antigua Carretera a Pátzcuaro # 8701 >> Col. Ex-Hacienda de San José de la Huerta, CP 58190 >> Morelia, Michoacán, México >> Tel.: 443-3222706 ext. 42511 >> e-mail: aacamargo@cieco.unam.mx >> skype: angela.camargo26 >> >> [[alternative HTML version deleted]] >> >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es >> >> >[[alternative HTML version deleted]]
Hola Angela, Muchas gracias por los detalles. Creo que no son necesarios los "counts" si tus datos y la distribucion de referencia son continuos. Ahora, para el caso discreto, bastaria con hacer un test chi^2 apropiadamente (ya sea construyendo tu propia funcion o utilizando la que tiene R). Para mas detalles dale una mirada a ?ks.test ?chisq.test y http://cran.r-project.org/doc/contrib/Ricci-distributions-en.pdf Un saludo, Jorge 2011/6/14 Angela Andrea Camargo Sanabria <>> Hola Jorge, > > Por supuesto no hay inconveniente en explicar, por el contrario mucho > mejor! Trataré de hacerme entender. > > Básicamente quiero comparar si una serie de datos tienen o no la misma > distribución. Quiero utilizar el test de KS y según entiendo, requiero la > frecuencia o los "counts" para hacer la comparación. Mi idea era calcular > los histogramas y guardar el $counts del objeto HIST en otro objeto y hacer > las comparaciones posteriormente. Entiendo también que requiere el mismo > número de clases para realizar la comparación con este test, de ahí mi > inquietud. > > Gracias a todos, > > ACS > > > > > > ---- > > *Angela Andrea Camargo Sanabria* > Estudiante Doctorado en Ciencias Biológicas > Laboratorio de Ecología de poblaciones y comunidades tropicales > Centro de Investigaciones en Ecosistemas (CIEco) > UNAM, campus Morelia > Antigua Carretera a Pátzcuaro # 8701 > Col. Ex-Hacienda de San José de la Huerta, CP 58190 > Morelia, Michoacán, México > Tel.: 443-3222706 ext. 42511 > e-mail: aacamargo@cieco.unam.mx > skype: angela.camargo26 > > > > > 2011/6/14 Jorge Ivan Velez <> > >> Hola Angela, >> >> Aunque no tengo claro el por que quieres algo asi, ni tampoco tengo una >> respuesta a tu pregunta, presumo que no quieres el _mismo_ numero de counts >> sino el mismo numero de clases. >> >> En principio, la opcion plot = FALSE en hist(), e.g. >> >> myh <- hist(x, plot = FALSE) >> str(myh) >> >> parece proporcionar informacion relevante para resolver lo que finalmente >> quieres hacer. Otras opciones incluyen [1] y [2], pero debes manipular los >> datos previamente. En [3] hay un ejemplo que puedes adaptar facilmente. >> >> Si no hay inconveniente, podrias por favor darnos un poco mas de detalles >> acerca del por que quieres algo asi? Que intentas mostrar? Quizas estos >> detalles nos ayuden a ayudarte un poco mas. >> >> Un saludo, >> Jorge Ivan Velez >> >> [1] http://had.co.nz/ggplot2/geom_histogram.html >> [2] http://lmdvr.r-forge.r-project.org/figures/figures.html >> [3] http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=55 >> >> >> >> 2011/6/14 Angela Andrea Camargo Sanabria <angela.andrea.camargo@gmail.com >> > >> >>> Buenas tardes, >>> >>> Estoy utilizando la función hist para calcular los histogramas de seis >>> grupos de datos de la siguiente forma >>> HIST <- hist (misdatos, main = "MG2008", ylab="", xlab ="", cex.axis=2, >>> breaks=9), >>> y quisiera que los seis histogramas tuvieran el mismo número de $counts, >>> es >>> decir, el mismo número de clases, pero no lo logro. Especifico breaks o >>> nclass en la función, pero R termina haciendo el gráfico con el número de >>> clases que él decide. >>> >>> Agradezco mucho su ayuda, >>> >>> Angela C. >>> ---- >>> *Angela Andrea Camargo Sanabria* >>> Estudiante Doctorado en Ciencias Biológicas >>> Laboratorio de Ecología de poblaciones y comunidades tropicales >>> Centro de Investigaciones en Ecosistemas (CIEco) >>> UNAM, campus Morelia >>> Antigua Carretera a Pátzcuaro # 8701 >>> Col. Ex-Hacienda de San José de la Huerta, CP 58190 >>> Morelia, Michoacán, México >>> Tel.: 443-3222706 ext. 42511 >>> e-mail: aacamargo@cieco.unam.mx >>> skype: angela.camargo26 >>> >>> [[alternative HTML version deleted]] >>> >>> >>> _______________________________________________ >>> R-help-es mailing list >>> R-help-es@r-project.org >>> https://stat.ethz.ch/mailman/listinfo/r-help-es >>> >>> >> >[[alternative HTML version deleted]]
Una solución un poco chapucera puede ser poner breaks=seq(min(misdatos)-1,max(misdatos)+1,length.out=9). De esa forma le expecificas claramente los 9 puntos de corte. Lo de restarle 1 al mínimo y sumarle 1 al máximo puedes adaptarlo a tus datos: se trata sólo de que no los dejes en el extremo del intervalo. De hecho, si los datos son más o menos equiparables, ponles a todos los mismos puntos de corte.> > Message: 3 > Date: Tue, 14 Jun 2011 13:26:20 -0500 > From: Angela Andrea Camargo Sanabria <angela.andrea.camargo en gmail.com> > To: r-help-es en r-project.org > Subject: [R-es] ayuda con la función hist > Message-ID: <BANLkTinV=v=uWVfQNc5yS4mXXKGsoUtLvg en mail.gmail.com> > Content-Type: text/plain > > Buenas tardes, > > Estoy utilizando la función hist para calcular los histogramas de seis > grupos de datos de la siguiente forma > HIST <- hist (misdatos, main = "MG2008", ylab="", xlab ="", cex.axis=2, > breaks=9), > y quisiera que los seis histogramas tuvieran el mismo número de $counts, > es > decir, el mismo número de clases, pero no lo logro. Especifico breaks o > nclass en la función, pero R termina haciendo el gráfico con el número de > clases que él decide. > > Agradezco mucho su ayuda, > > Angela C. > ---- > *Angela Andrea Camargo Sanabria* > Estudiante Doctorado en Ciencias Biológicas > Laboratorio de Ecología de poblaciones y comunidades tropicales > Centro de Investigaciones en Ecosistemas (CIEco) > UNAM, campus Morelia > Antigua Carretera a Pátzcuaro # 8701 > Col. Ex-Hacienda de San José de la Huerta, CP 58190 > Morelia, Michoacán, México > Tel.: 443-3222706 ext. 42511 > e-mail: aacamargo en cieco.unam.mx > skype: angela.camargo26 > > [[alternative HTML version deleted]] > > >-- Dr. Antonio José Sáez Castillo Departamento de Estadística e Investigación Operativa Escuela Politécnica Superior de Linares Universidad de Jaén c/ Alfonso X El Sabio 28. Dep. A229 Tfno y FAX 953648578
Angela, el test de Kolmogorov Smirnov no está basado en una comparación de los histogramas (aproximación empírica de la densidad), sino más bien en la comparación de las la funciones de probabilidades acumuladas (empíricas). De hecho está implementado en R. Así, para contrastar que x e y tienen la misma distribución, basta con escribir: ks.test(x, y) Un saludo. Olivier -- ____________________________________ Olivier G. Nuñez Email: onunez en iberstat.es Tel : +34 663 03 69 09 Web: http://www.iberstat.es ____________________________________ El 15/06/2011, a las 0:41, Angela Andrea Camargo Sanabria escribió:> Hola Jorge, > > Por supuesto no hay inconveniente en explicar, por el contrario > mucho mejor! > Trataré de hacerme entender. > > Básicamente quiero comparar si una serie de datos tienen o no la misma > distribución. Quiero utilizar el test de KS y según entiendo, > requiero la > frecuencia o los "counts" para hacer la comparación. Mi idea era > calcular > los histogramas y guardar el $counts del objeto HIST en otro > objeto y hacer > las comparaciones posteriormente. Entiendo también que requiere el > mismo > número de clases para realizar la comparación con este test, de ahí mi > inquietud. > > Gracias a todos, > > ACS > > > > > > ---- > *Angela Andrea Camargo Sanabria* > Estudiante Doctorado en Ciencias Biológicas > Laboratorio de Ecología de poblaciones y comunidades tropicales > Centro de Investigaciones en Ecosistemas (CIEco) > UNAM, campus Morelia > Antigua Carretera a Pátzcuaro # 8701 > Col. Ex-Hacienda de San José de la Huerta, CP 58190 > Morelia, Michoacán, México > Tel.: 443-3222706 ext. 42511 > e-mail: aacamargo en cieco.unam.mx > skype: angela.camargo26 > > > > > 2011/6/14 Jorge Ivan Velez <jorgeivanvelez en gmail.com> > >> Hola Angela, >> >> Aunque no tengo claro el por que quieres algo asi, ni tampoco >> tengo una >> respuesta a tu pregunta, presumo que no quieres el _mismo_ numero >> de counts >> sino el mismo numero de clases. >> >> En principio, la opcion plot = FALSE en hist(), e.g. >> >> myh <- hist(x, plot = FALSE) >> str(myh) >> >> parece proporcionar informacion relevante para resolver lo que >> finalmente >> quieres hacer. Otras opciones incluyen [1] y [2], pero debes >> manipular los >> datos previamente. En [3] hay un ejemplo que puedes adaptar >> facilmente. >> >> Si no hay inconveniente, podrias por favor darnos un poco mas de >> detalles >> acerca del por que quieres algo asi? Que intentas mostrar? >> Quizas estos >> detalles nos ayuden a ayudarte un poco mas. >> >> Un saludo, >> Jorge Ivan Velez >> >> [1] http://had.co.nz/ggplot2/geom_histogram.html >> [2] http://lmdvr.r-forge.r-project.org/figures/figures.html >> [3] http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=55 >> >> >> >> 2011/6/14 Angela Andrea Camargo Sanabria >> <angela.andrea.camargo en gmail.com> >> >>> Buenas tardes, >>> >>> Estoy utilizando la función hist para calcular los histogramas de >>> seis >>> grupos de datos de la siguiente forma >>> HIST <- hist (misdatos, main = "MG2008", ylab="", xlab ="", >>> cex.axis=2, >>> breaks=9), >>> y quisiera que los seis histogramas tuvieran el mismo número de >>> $counts, >>> es >>> decir, el mismo número de clases, pero no lo logro. Especifico >>> breaks o >>> nclass en la función, pero R termina haciendo el gráfico con el >>> número de >>> clases que él decide. >>> >>> Agradezco mucho su ayuda, >>> >>> Angela C. >>> ---- >>> *Angela Andrea Camargo Sanabria* >>> Estudiante Doctorado en Ciencias Biológicas >>> Laboratorio de Ecología de poblaciones y comunidades tropicales >>> Centro de Investigaciones en Ecosistemas (CIEco) >>> UNAM, campus Morelia >>> Antigua Carretera a Pátzcuaro # 8701 >>> Col. Ex-Hacienda de San José de la Huerta, CP 58190 >>> Morelia, Michoacán, México >>> Tel.: 443-3222706 ext. 42511 >>> e-mail: aacamargo en cieco.unam.mx >>> skype: angela.camargo26 >>> >>> [[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 >>> >>> >> > > [[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
Una alternativa es: # Creamos unos datos de prueba x = 2 * rnorm(100) + 1 y = rnorm(100) # Y hacemos los gráficos plot(density(x), col="red") lines(density(y), col="blue") # Y esto es todo y con esto te aparecen curvas suavizadas de las distribuciones de x e y en el mismo gráfico. -----Mensaje original----- De: r-help-es-bounces en r-project.org [mailto:r-help-es-bounces en r-project.org] En nombre de Angela Andrea Camargo Sanabria Enviado el: Martes, 14 de Junio de 2011 06:42 p.m. Para: Jorge Ivan Velez CC: R-help-es Asunto: Re: [R-es] ayuda con la función hist Hola Jorge, Por supuesto no hay inconveniente en explicar, por el contrario mucho mejor! Trataré de hacerme entender. Básicamente quiero comparar si una serie de datos tienen o no la misma distribución. Quiero utilizar el test de KS y según entiendo, requiero la frecuencia o los "counts" para hacer la comparación. Mi idea era calcular los histogramas y guardar el $counts del objeto HIST en otro objeto y hacer las comparaciones posteriormente. Entiendo también que requiere el mismo número de clases para realizar la comparación con este test, de ahí mi inquietud. Gracias a todos, ACS ---- *Angela Andrea Camargo Sanabria* Estudiante Doctorado en Ciencias Biológicas Laboratorio de Ecología de poblaciones y comunidades tropicales Centro de Investigaciones en Ecosistemas (CIEco) UNAM, campus Morelia Antigua Carretera a Pátzcuaro # 8701 Col. Ex-Hacienda de San José de la Huerta, CP 58190 Morelia, Michoacán, México Tel.: 443-3222706 ext. 42511 e-mail: aacamargo en cieco.unam.mx skype: angela.camargo26 2011/6/14 Jorge Ivan Velez <jorgeivanvelez en gmail.com>> Hola Angela, > > Aunque no tengo claro el por que quieres algo asi, ni tampoco tengo > una respuesta a tu pregunta, presumo que no quieres el _mismo_ numero > de counts sino el mismo numero de clases. > > En principio, la opcion plot = FALSE en hist(), e.g. > > myh <- hist(x, plot = FALSE) > str(myh) > > parece proporcionar informacion relevante para resolver lo que > finalmente quieres hacer. Otras opciones incluyen [1] y [2], pero > debes manipular los datos previamente. En [3] hay un ejemplo que puedesadaptar facilmente.> > Si no hay inconveniente, podrias por favor darnos un poco mas de > detalles acerca del por que quieres algo asi? Que intentas mostrar? > Quizas estos detalles nos ayuden a ayudarte un poco mas. > > Un saludo, > Jorge Ivan Velez > > [1] http://had.co.nz/ggplot2/geom_histogram.html > [2] http://lmdvr.r-forge.r-project.org/figures/figures.html > [3] http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=55 > > > > 2011/6/14 Angela Andrea Camargo Sanabria > <angela.andrea.camargo en gmail.com> > >> Buenas tardes, >> >> Estoy utilizando la función hist para calcular los histogramas de >> seis grupos de datos de la siguiente forma HIST <- hist (misdatos, >> main = "MG2008", ylab="", xlab ="", cex.axis=2, breaks=9), y quisiera >> que los seis histogramas tuvieran el mismo número de $counts, es >> decir, el mismo número de clases, pero no lo logro. Especifico breaks >> o nclass en la función, pero R termina haciendo el gráfico con el >> número de clases que él decide. >> >> Agradezco mucho su ayuda, >> >> Angela C. >> ---- >> *Angela Andrea Camargo Sanabria* >> Estudiante Doctorado en Ciencias Biológicas Laboratorio de Ecología >> de poblaciones y comunidades tropicales Centro de Investigaciones en >> Ecosistemas (CIEco) UNAM, campus Morelia Antigua Carretera a >> Pátzcuaro # 8701 Col. Ex-Hacienda de San José de la Huerta, CP 58190 >> Morelia, Michoacán, México >> Tel.: 443-3222706 ext. 42511 >> e-mail: aacamargo en cieco.unam.mx >> skype: angela.camargo26 >> >> [[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 >> >> >[[alternative HTML version deleted]] *************************************************************************************************** 1 de Agosto de 1961, creación del Ministerio del Transporte.
Angela, efectivamente, si quieres comparar dos distribuciones empíricas, la versión del test KS adecuada es el test de Lilliefors ( función lillie.test() de R) Este test lo encontraras en el paquete nortest. http://www.stat.ucl.ac.be/ISdidactique/Rhelp/library/nortest/html/ lillie.test.html Un saludo. Olivier -- ____________________________________ Olivier G. Nuñez Email: onunez en iberstat.es Tel : +34 663 03 69 09 Web: http://www.iberstat.es ____________________________________ El 15/06/2011, a las 13:14, Angela Andrea Camargo Sanabria escribió:> Hola Oliver, > > Es que al utilizar este test me salía el siguiente warning > Mensajes de aviso perdidos > In ks.test(c2010[, 4], tm2010[, 4]) : > cannot compute correct p-values with ties > > y por esto pensé que una razón sería que no tienen igual número de > clases... > > Pero gracias, estoy probando varias de las soluciones q me han dado. > > ACS > ---- > Angela Andrea Camargo Sanabria > Estudiante Doctorado en Ciencias Biológicas > Laboratorio de Ecología de poblaciones y comunidades tropicales > Centro de Investigaciones en Ecosistemas (CIEco) > UNAM, campus Morelia > Antigua Carretera a Pátzcuaro # 8701 > Col. Ex-Hacienda de San José de la Huerta, CP 58190 > Morelia, Michoacán, México > Tel.: 443-3222706 ext. 42511 > e-mail: aacamargo en cieco.unam.mx > skype: angela.camargo26 > >
Hola a tod en s. Tengo que modificar unas rutinas en FoxPro y estoy considerando el migrarlas a R. Sabéis si existe alguna manera de que R las lea para evitar reescribirlas? NOTA: Son ficheros .PRG con código de las versiones antiguas de FoxPro (no de Visual FoxPro). Un Saludo, Miguel. Nota: A información contida nesta mensaxe e os seus posibles documentos adxuntos é privada e confidencial e está dirixida únicamente ó seu destinatario/a. Se vostede non é o/a destinatario/a orixinal desta mensaxe, por favor elimínea. A distribución ou copia desta mensaxe non está autorizada. Nota: La información contenida en este mensaje y sus posibles documentos adjuntos es privada y confidencial y está dirigida únicamente a su destinatario/a. Si usted no es el/la destinatario/a original de este mensaje, por favor elimínelo. La distribución o copia de este mensaje no está autorizada. See more languages: http://www.sergas.es/aviso_confidencialidad.htm
Perdona Angela, mis recuerdos sobre el test de Lilliefors eran erróneos. Lilliefors es para contrastar la normalidad cuando la media y la varianza son desconocidas. No toma en cuenta mi ultimo mensaje. Un saludo. Olivier El 15/06/2011, a las 13:40, Olivier Nuñez escribió:> Angela, > > efectivamente, si quieres comparar dos distribuciones empíricas, > la versión del test KS adecuada es el test de Lilliefors ( función > lillie.test() de R) > Este test lo encontraras en el paquete nortest. > > http://www.stat.ucl.ac.be/ISdidactique/Rhelp/library/nortest/html/ > lillie.test.html > > Un saludo. Olivier > -- ____________________________________ > > Olivier G. Nuñez > Email: onunez@iberstat.es > Tel : +34 663 03 69 09 > Web: http://www.iberstat.es > > ____________________________________ > > > > > El 15/06/2011, a las 13:14, Angela Andrea Camargo Sanabria escribió: > >> Hola Oliver, >> >> Es que al utilizar este test me salía el siguiente warning >> Mensajes de aviso perdidos >> In ks.test(c2010[, 4], tm2010[, 4]) : >> cannot compute correct p-values with ties >> >> y por esto pensé que una razón sería que no tienen igual número de >> clases... >> >> Pero gracias, estoy probando varias de las soluciones q me han dado. >> >> ACS >> ---- >> Angela Andrea Camargo Sanabria >> Estudiante Doctorado en Ciencias Biológicas >> Laboratorio de Ecología de poblaciones y comunidades tropicales >> Centro de Investigaciones en Ecosistemas (CIEco) >> UNAM, campus Morelia >> Antigua Carretera a Pátzcuaro # 8701 >> Col. Ex-Hacienda de San José de la Huerta, CP 58190 >> Morelia, Michoacán, México >> Tel.: 443-3222706 ext. 42511 >> e-mail: aacamargo@cieco.unam.mx >> skype: angela.camargo26 >> >> > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es[[alternative HTML version deleted]]
Angela, Una respuesta más adecuada a tu duda sería sugerirte el uso de la función ks.boot() del paquete Matching, que proporciona p-valores correctos aún si las distribuciones que compararas son parcialmente continuas. http://sekhon.berkeley.edu/matching/ks.boot.html Un saludo. Olivier El 15/06/2011, a las 15:50, Olivier Nuñez escribió:> Perdona Angela, > mis recuerdos sobre el test de Lilliefors eran erróneos. > Lilliefors es para contrastar la normalidad cuando la media y la > varianza son desconocidas. > No toma en cuenta mi ultimo mensaje. > Un saludo. Olivier > > > El 15/06/2011, a las 13:40, Olivier Nuñez escribió: > >> Angela, >> >> efectivamente, si quieres comparar dos distribuciones empíricas, >> la versión del test KS adecuada es el test de Lilliefors ( función >> lillie.test() de R) >> Este test lo encontraras en el paquete nortest. >> >> http://www.stat.ucl.ac.be/ISdidactique/Rhelp/library/nortest/html/ >> lillie.test.html >> >> Un saludo. Olivier >> -- ____________________________________ >> >> Olivier G. Nuñez >> Email: onunez en iberstat.es >> Tel : +34 663 03 69 09 >> Web: http://www.iberstat.es >> >> ____________________________________ >> >> >> >> >> El 15/06/2011, a las 13:14, Angela Andrea Camargo Sanabria escribió: >> >>> Hola Oliver, >>> >>> Es que al utilizar este test me salía el siguiente warning >>> Mensajes de aviso perdidos >>> In ks.test(c2010[, 4], tm2010[, 4]) : >>> cannot compute correct p-values with ties >>> >>> y por esto pensé que una razón sería que no tienen igual número >>> de clases... >>> >>> Pero gracias, estoy probando varias de las soluciones q me han dado. >>> >>> ACS >>> ---- >>> Angela Andrea Camargo Sanabria >>> Estudiante Doctorado en Ciencias Biológicas >>> Laboratorio de Ecología de poblaciones y comunidades tropicales >>> Centro de Investigaciones en Ecosistemas (CIEco) >>> UNAM, campus Morelia >>> Antigua Carretera a Pátzcuaro # 8701 >>> Col. Ex-Hacienda de San José de la Huerta, CP 58190 >>> Morelia, Michoacán, México >>> Tel.: 443-3222706 ext. 42511 >>> e-mail: aacamargo en cieco.unam.mx >>> skype: angela.camargo26 >>> >>> >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es en r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es >