[R-es] ayuda con la función hist

[Olivier Nuñez]

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