Der Run-Test (auch Runs-Test, Wald-Wolfowitz-Test nach Abraham Wald und Jacob Wolfowitz, Iterationstest oder Geary-Test) ist ein nichtparametrischer Test auf Zufälligkeit einer Folge. Ausgegangen wird von einem Urnenmodell mit zwei Sorten Kugeln (dichotome Grundgesamtheit). Es werden n Kugeln entnommen und es soll die Hypothese geprüft werden, dass die Entnahme zufällig erfolgt ist.

Vorgehensweise

Bearbeiten

Es wurden einer dichotomen Grundgesamtheit   Kugeln entnommen. Die Ergebnisse liegen in ihrer chronologischen Abfolge vor. Es werden nun alle benachbarten Ergebnisse gleicher Ausprägung zu einem Lauf oder Run zusammengefasst. Wenn die Folge tatsächlich zufällig ist, sollten nicht zu wenige Runs vorliegen, aber auch nicht zu viele.

Es wird die Nullhypothese aufgestellt: Die Entnahme erfolgte zufällig.

Für die Festlegung der Zahl der Runs, bei der die Hypothese abgelehnt wird, wird die Verteilung der Runs benötigt: Es seien   die Zahl der Kugeln erster Sorte und   der zweiten Sorte; es sei   die Zahl der Runs. Nach dem Symmetrieprinzip ist die Wahrscheinlichkeit für jede beliebige Folge der Kugeln bei zufälliger Entnahme gleich groß. Es gibt insgesamt

 

Möglichkeiten der Entnahme.

Testverteilung

Bearbeiten

Die Testverteilung ist die Verteilung der Zahl der Iterationen (Runs) bei Richtigkeit der Nullhypothese.

Bezüglich der Verteilung der Zahl der Runs unterscheidet man die Fälle:

1. Die Anzahl   der Runs ist geradzahlig:

Es liegen   Runs der Kugeln der ersten Sorte und auch   Runs der Kugeln der zweiten Sorte vor. Die Wahrscheinlichkeit, dass genau   Runs eingetreten sind, ist dann
 

2. Die Anzahl   der Runs ist ungeradzahlig:

Es liegen   Runs der Kugeln der ersten Sorte und   Runs der Kugeln der zweiten Sorte vor oder der umgekehrte Fall. Die Wahrscheinlichkeit, dass genau   Runs eingetreten sind, berechnet sich dann als Summe aus diesen beiden Möglichkeiten
 

Dabei bezeichnet   den Binomialkoeffizienten. Zu berücksichtigen ist  , falls  .

Testdurchführung

Bearbeiten

Als Prüfgröße wird die Zufallsvariable   verwendet. Die Testverteilung, d. h. die Verteilung der Prüfgröße bei Richtigkeit der Nullhypothese, wurde zuvor beschrieben. Die Nullhypothese wird bei zweiseitigem Testen dann abgelehnt, wenn die beobachtete Anzahl von Iterationen   zu klein oder zu groß ist. Bei einem vorgegebenen Signifikanzniveau   wird die Nullhypothese dann abgelehnt, wenn

  oder  

gilt. Dabei bezeichnet   das  -Quantil der Testverteilung. Da diese Verteilung diskret ist, kann das vorgegebene Signifikanzniveau   durch einen (nicht randomisierten) Test im Allgemeinen nicht exakt erreicht werden. Es werden daher die beiden kritischen Wert so gewählt, dass die Fehlerwahrscheinlichkeit 1. Art höchstens so groß wie das vorgegebene Signifikanzniveau   ist. Dieses Vorgehen wird auch als konservatives Testen bezeichnet. Da die Berechnung der kritischen Werte umständlich ist, bedient man sich häufig einer Tabelle.[1]

Einfaches Beispiel

Bearbeiten

Für eine Podiumsdiskussion mit zwei politischen Parteien wurden die Sprecher angeblich zufällig ermittelt. Es wurde ausgelost, dass von der Partei Supi 4 Vertreter und von der Partei Toll 5 Vertreter in der folgenden Reihe sprechen dürfen:

S S  T  S  T T T  S  T

Ein Vertreter von Toll beschwerte sich, dass S bevorzugt würde. Es wurde ein Run-Test vorgenommen:

Es ist n1 = 4 und n2 = 5. Man erhielt r = 6 Runs.

Deutlich ist, dass im Falle vieler Runs kein Verdacht besteht auf Bevorzugung einer der Parteien. Die Nullhypothese wird also abgelehnt, wenn es zu wenig Runs gibt. Nach der Tabelle des Run-Testes wird H0 abgelehnt, wenn r ≤ 2. Also liegt die Prüfgröße r = 6 nicht im Ablehnungsbereich; man kann nach den Kriterien des Run-Testes nicht darauf schließen, dass die Reihenfolge der Sprecher nicht zufällig ist.

Übrigens wird auch im nächsten Fall:

S S S  T  S  T T T T

mit r = 4 Runs, die Nullhypothese nicht abgelehnt, obwohl fast jeder einen Verdacht haben wird, dass Supi vorgezogen wurde. Man kann aber wegen der relativ geringen Anzahl der Beobachtungen nicht ausschließen, dass das Ergebnis auf Zufall beruht.

Ergänzungen

Bearbeiten

Parameter der Verteilung von R

Bearbeiten

Der Erwartungswert von R ist

 

und die Varianz

 .

Grundgesamtheit mit mehr als zwei Ausprägungen des Merkmals

Bearbeiten

Liegt eine endliche Folge reeller Zahlen   eines metrischen Merkmals vor, wird die Folge dichotomisiert: Man bestimmt zunächst den Median z der Folge. Werte   werden dann als Kugeln der ersten Sorte, Werte   als Kugeln der zweiten Sorte interpretiert. Die entstandene dichotome Folge kann dann wieder auf Zufälligkeit getestet werden (siehe Beispiel unten).

Liegt eine nichtnumerische Symbolsequenz mit mehr als zwei Ausprägungen vor, muss zunächst eine numerische Reihe erzeugt werden, wobei hier das Problem bestehen kann, dass die Symbole nicht geordnet werden können.

Normalapproximation

Bearbeiten

Für Stichprobenumfänge n1,n2 > 20 ist die Zahl der Runs R annähernd normalverteilt mit Erwartungswert und Varianz wie oben. Man erhält die standardisierte Prüfgröße

 

Die Hypothese wird abgelehnt, wenn

  oder  

mit   als Quantil der Standardnormalverteilung für die Wahrscheinlichkeit   .

Anwendungen

Bearbeiten

Der Runtest kann angewendet werden, um Stationarität bzw. Nicht-Korrelation in einer Zeitreihe oder anderen Sequenz zu überprüfen, vor allem wenn die Verteilung des Merkmals unbekannt ist. Die Nullhypothese ist hier, dass aufeinanderfolgende Werte unkorreliert sind.

Der Run-Test kann mit dem Chi-Quadrat-Test kombiniert werden, da beide Prüfgrößen asymptotisch unabhängig voneinander sind.

Beispiel für ein metrisches Merkmal

Bearbeiten

Es liegt die Folge

13	 3	14	14	1	14	3	8	14	17	9	14	13	2	16	1	3	12	13	14

vor. Sie wird mit dem Median z = 13 dichotomisiert. Für die erste Ausprägung wird + gesetzt, für die zweite Ausprägung -.

0	-10	1	1	-12	1	-10	-5	1	4	-4	1	0	-11	3	-12	-10	-1	0	1
+	-	+	+	-	+	-	-	+	+	-	+	+	-	+	-	-	-	+	+

Man erhält bei n1 = 11 (+) und n2 = 9 (-) r = 13 Runs. R ist annähernd normalverteilt mit dem Erwartungswert

 

und der Varianz

 .

Die Prüfgröße z errechnet sich dann als

 

Bei einem Signifikanzniveau von 0,05 wird H0 abgelehnt, wenn |z| > 1,96. Dies ist nicht der Fall.

Entscheidung: Die Hypothese wird nicht abgelehnt. Die Elemente der Stichprobe sind vermutlich zufällig entnommen worden.

Da der Run-Test aber kein parametrischer Test ist, ist das Resultat mit Vorsicht zu genießen. Bei einem Konfidenzniveau von 90 % könnte z. B. die Nullhypothese abgelehnt werden. Der parametrische Shapiro-Wilk-Test zeigt nämlich, dass bei der vorliegenden Zahlenreihe die Normalverteilung nicht gegeben ist!

Literatur

Bearbeiten

Siehe auch

Bearbeiten

Autokorrelation, Zufallszahlengenerator, Pseudozufallszahlen, Trend

Einzelnachweise

Bearbeiten
  1. Eine Tabelle für   findet sich hier: Horst Rinne: Taschenbuch der Statistik. 4. Auflage. Harri Deutsch, Frankfurt am Main 2008, ISBN 978-3-8171-1827-4, S. 564. Eine weitere Quelle ist: Herbert Büning, Götz Trenkler: Nichtparametrische Statistische Methoden. 2. völlig neu bearbeitete Auflage. Walter de Gruyter, Berlin / New York 1994, ISBN 3-11-013860-3, S. 393, doi:10.1515/9783110902990.
Bearbeiten