Existe-t-il un test d'adéquation d'Anderson pour deux ensembles de données?

10

Je sais que ad.test () peut être utilisé pour tester la normalité.

Est-il possible d'obtenir ad.test pour comparer les distributions de deux échantillons de données?

x <- rnorm(1000)
y <- rgev(2000)
ad.test(x,y)

Comment puis-je effectuer le test Anderson-Darling sur 2 échantillons?

evdstat
la source
2
L' article de Wikipédia sur le test AD le mentionne sous le titre "Tests non paramétriques à échantillon k". Sa référence, un article JASA de 1987 par Sholz et Stephens, est disponible gratuitement sur cithep.caltech.edu/~fcp/statistics/hypothesisTest/… .
whuber
Si la question est: comment puis-je le faire en R (comme le suggère le tag): bonne question (+1) (et la réponse est probablement: montez-le vous-même), bien que quelque peu déplacé ici ( StackOverflow est un meilleur endroit pour ce genre question).
Nick Sabbe
@Nick Trouver ou implémenter un test du GoF, que ce soit en R ou dans tout autre langage, correspond parfaitement à notre intérêt pour tout ce qui concerne les statistiques.
whuber
1
@whuber: Je suis corrigé: je viens de lire la partie pertinente de la FAQ. Pourtant, c'est une fine ligne entre l'amour et la haine. Mais je n'ai pas voté pour migrer :-)
Nick Sabbe
2
@ Nick, je suis d'accord sur la ligne mince. Lorsqu'une question se concentre uniquement sur la mécanique de la programmation, sa pertinence devient ici douteuse. Vous pouvez trouver des discussions périodiques à ce sujet sur meta.
whuber

Réponses:

7

Le package a adkété remplacé par le package kSamples:

Essayer:

install.packages("kSamples")  
library(kSamples)
ad.test(runif(50), rnorm(30))
Fritz Scholz
la source
la kSamples::ad.testfonction est plutôt lente. Existe-t-il une alternative plus efficace?
Nemesi
5

Le package adk pour R fait cela. http://cran.r-project.org/web/packages/adk/

install.packages("adk")
library(adk)
adk.test(runif(50), rnorm(30))
momeara
la source
2
Le package adk a été supprimé de CRAN. Y a-t-il une alternative?
JE42