Me gustaría realizar algunas pruebas bidimensionales de Kolmogorov-Smironov para determinar si una distribución bidimensional se ajusta a una referencia.
¿Hay algún paquete o aplicación que pueda usar de una manera relativamente sencilla? ¿O hay un algoritmo diferente que es preferible? Solo tengo un conocimiento estadístico básico.
kolmogorov-smirnov
Manuel J Gomez
fuente
fuente
Respuestas:
Justel, Pena y Zamar han descrito una extensión bidimensional de la prueba de Kolmogorov-Smirnov en una "Prueba de bondad de ajuste Komogorov-Smirnov multivariante" . Los comentarios de @ Procrastinator sugieren que puede haber otras propuestas de este tipo.
Sin embargo, no he visto un paquete con una implementación sencilla.
Dependiendo de lo que quiera hacer, kde.test () en el paquete ks de Tarn Duong para R podría ser más útil.
fuente
Implementación de Python
He escrito una implementación de python usando numpy. Puede encontrar el código aquí , puede encontrar más información en la cadena de documentación del código.
Y aquí hay otro (no por mí). Este cuaderno proporciona una implementación de Python para la prueba 2D KS con 2 muestras. El
.pyarchivo se puede descargar aquí . El código parece ser una traducción directa delCcódigo, la eficiencia podría ser un problema si el tamaño de la muestra es grande.Sin embargo, es mejor que compruebe los códigos (sin importar cuál) con los documentos / libros originales antes de usarlos. Las implementaciones de Python de la prueba 2d KS están mucho menos controladas que las de R.
Más información
El algoritmo se desarrolló por primera vez en dos documentos (como veo)
Una buena introducción y la
Cimplementación se pueden encontrar enPresione, WH et al. 1992, Recetas numéricas en C , Sección 14.7, p645.
Puede encontrar la
C++/Fortranimplementación en otras versiones del libro.Aquí hay una publicación titulada Cuidado con que la prueba de Kolmogorov-Smirnov también está relacionada con el tema, es posible que desee echar un vistazo. Alienta el uso del método de remuestreo para evaluar el valor p con una distancia KS dada.
fuente
Es posible que este código de Matlab sea útil.
http://www.mathworks.com/matlabcentral/fileexchange/38617-two-dimensional-2d-paired-kolmogorov-smirnov-test
fuente