Je voudrais exécuter des tests de Kolmogorov-Smironov bidimensionnels pour déterminer si une distribution bidimensionnelle correspond à une référence.
Existe-t-il un package ou une application que je pourrais utiliser de manière relativement simple? Ou existe-t-il un algorithme différent préférable? J'ai juste une connaissance statistique de base.
kolmogorov-smirnov
Manuel J Gomez
la source
la source
Réponses:
Une extension bidimensionnelle du test de Kolmogorov-Smirnov a été décrite par Justel, Pena et Zamar dans un «test multivarié de Komogorov-Smirnov de qualité de l'ajustement» . @ Les commentaires du procrastinateur suggèrent qu'il pourrait y avoir d'autres propositions de ce type.
Cependant, je n'ai pas vu de package avec une implémentation simple.
Selon ce que vous voulez faire, kde.test () dans le paquet ks de Tarn Duong pour R pourrait être plus utile.
la source
Implémentation de Python
J'ai écrit une implémentation python en utilisant numpy. Vous pouvez trouver le code ici , vous pouvez trouver plus d'informations dans la docstring du code.
Et voici un autre (pas par moi). Ce Notebook fournit une implémentation Python pour le test 2D KS avec 2 échantillons. Le
.py
fichier peut être téléchargé ici . Le code semble être une traduction directe duC
code, l'efficacité peut être un problème si la taille de l'échantillon est grande.Cependant, vous feriez mieux de vérifier les codes (peu importe lequel) avec les papiers / livres originaux avant de les utiliser. Les implémentations python du test 2d KS sont beaucoup moins vérifiées que celles de R.
Plus d'informations
L'algorithme est d'abord développé dans deux articles (comme je le vois)
Une belle introduction et l'
C
implémentation se trouvent dansAppuyez sur, WH et al. 1992, Recettes numériques en C , section 14.7, p645.
Vous pouvez trouver l'
C++/Fortran
implémentation dans d'autres versions du livre.Voici un article intitulé Attention, le test de Kolmogorov-Smirnov est également lié au sujet, vous voudrez peut-être y jeter un œil. Il encourage l'utilisation de la méthode de rééchantillonnage pour évaluer la valeur de p avec une distance KS donnée.
la source
vous pouvez trouver ce code Matlab utile.
http://www.mathworks.com/matlabcentral/fileexchange/38617-two-dimensional-2d-paired-kolmogorov-smirnov-test
la source