Kolmogorov-Smirnov bidimensional

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.

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.

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 del Ccó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)

• Peacock, JA 1983, Pruebas bidimensionales de bondad de ajuste en astronomía
• Fasano, G. y Franceschini, A. 1987, una versión multidimensional de la prueba de Kolmogorov-Smirnov .

Una buena introducción y la Cimplementación se pueden encontrar en

• Presione, 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.

Es posible que este código de Matlab sea útil.

http://www.mathworks.com/matlabcentral/fileexchange/38617-two-dimensional-2d-paired-kolmogorov-smirnov-test