Czy istnieje test dobroci dopasowania Andersona dla dwóch zestawów danych?

10

Wiem, że ad.test () może być używany do testowania normalności.

Czy można uzyskać ad.test w celu porównania dystrybucji z dwóch próbek danych?

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

Jak mogę wykonać test Andersona-Darlinga na 2 próbkach?

evdstat
źródło
2
Artykuł w Wikipedii na temat testu AD wspomina o tym pod nagłówkiem „Nieparametryczne testy k-próbki”. Jego odniesienie, artykuł JASA z 1987 roku autorstwa Sholza i Stephena, jest dostępny bezpłatnie na stronie cithep.caltech.edu/~fcp/statistics/hypothesisTest/ ...
whuber
Jeśli pytanie brzmi: jak mogę to zrobić w R (jak sugeruje tag): dobre pytanie (+1) (a odpowiedź brzmi prawdopodobnie: sam to sfałszuj), choć nieco źle tutaj ( StackOverflow jest lepszym miejscem dla tego rodzaju pytania).
Nick Sabbe
@ Nick Znalezienie lub wdrożenie testu GoF, zarówno w języku R, jak i innym języku, leży w naszym interesie statystycznym.
whuber
1
@ whuber: Poprawiłem się: Właśnie przeczytałem odpowiednią część FAQ. Jest to jednak cienka granica między miłością a nienawiścią. Ale nie głosowałem na migrację :-)
Nick Sabbe,
2
@Nick Zgadzam się co do cienkiej linii. Kiedy pytanie koncentruje się wyłącznie na mechanice programowania, jego trafność tutaj staje się wątpliwa. Możesz znaleźć okresowe dyskusje na ten temat na stronie meta.
whuber

Odpowiedzi:

7

Pakiet adkzostał zastąpiony przez pakiet kSamples:

Próbować:

install.packages("kSamples")  
library(kSamples)
ad.test(runif(50), rnorm(30))
Fritz Scholz
źródło
kSamples::ad.testfunkcja jest raczej powolne. Czy istnieje bardziej skuteczna alternatywa?
Nemesi
5

Robi to pakiet adk dla R. http://cran.r-project.org/web/packages/adk/

install.packages("adk")
library(adk)
adk.test(runif(50), rnorm(30))
momeara
źródło
2
Pakiet adk został usunięty z CRAN. Czy jest alternatywa?
JE42