Levene-Test

Verteilung des Nettoeinkommens in Deutschland 2008 (ALLBUS) nach Geschlecht und Geburtsmonats des Befragten.

Der Levene-Test[1][2] bezeichnet in der Statistik einen Signifikanztest, der auf Gleichheit der Varianzen (Homoskedastizität) von zwei oder mehr Grundgesamtheiten (Gruppen) prüft. Der Brown–Forsythe Test ist aus dem Levene-Test abgeleitet. Er stammt von Howard Levene.

Ähnlich dem Bartlett-Test prüft der Levene-Test die Nullhypothese darauf, dass alle Gruppenvarianzen gleich sind. Die Alternativhypothese lautet demnach, dass mindestens ein Gruppenpaar ungleiche Varianzen besitzt (Heteroskedastizität):

Nullhypothese: H 0 : σ 1 2 = σ 2 2 = = σ k 2 {\displaystyle H_{0}\colon \sigma _{1}^{2}=\sigma _{2}^{2}=\ldots =\sigma _{k}^{2}}
Alternativhypothese: H 1 : σ i 2 σ j 2 {\displaystyle H_{1}\colon \sigma _{i}^{2}\neq \sigma _{j}^{2}}   für mindestens ein Gruppenpaar i , j {\displaystyle i,j} mit i j {\displaystyle i\neq j}

Befindet sich der p-Wert des Tests unter einem zuvor bestimmten Niveau, so sind die Unterschiede in den Varianzen der Stichproben überzufällig (signifikant) und die Nullhypothese der Varianzgleichheit kann abgelehnt werden.[3]

Beispiel

Die Grafik oben zeigt die Verteilung des Nettoeinkommens nach Geschlecht und Geburtsmonat. Die Ausgabe von car::leveneTest in R:

  • Der Levene-Test nach Geschlecht ergibt einen p-Wert kleiner als 2 , 2 × 10 16 {\displaystyle 2{,}2\times 10^{-16}} und ist damit hochsignifikant:
Levene’s Test for Homogeneity of Variance
        Df F value Pr(>F)
group 1  106.09 < 2.2e-16 ***
      2404
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Bei einem solchen p-Wert kann davon ausgegangen werden, dass die Varianzen in der Population unterschiedlich sind. Die Nullhypothese gleicher Varianzen wird entsprechend verworfen.

  • Der Levene-Test nach Geburtsmonat ergibt einen p-Wert von 0,076 {\displaystyle 0{,}076} und ist bei einem vorgegebenen Signifikanzniveau von 5 % nicht signifikant:
Levene’s Test for Homogeneity of Variance
        Df F value Pr(>F)
group 11  1.6621  0.076.
      2384
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Teststatistik

Sind X j i {\displaystyle X_{ji}} ( j = 1 , , k {\displaystyle j=1,\dots ,k} und i = 1 , , n j {\displaystyle i=1,\dots ,n_{j}} ) die Stichprobenvariablen und

Y j i = | X j i X ¯ j | {\displaystyle Y_{ji}=|X_{ji}-{\bar {X}}_{j}|}

mit k {\displaystyle k} Anzahl der Gruppen (Stichproben), n j {\displaystyle n_{j}} die Anzahl der Beobachtungen in Gruppe j {\displaystyle j} und X ¯ j {\displaystyle {\bar {X}}_{j}} der Stichprobenmittelwert der Gruppe j {\displaystyle j} . Dann ist die Teststatistik

L = 1 k 1 j = 1 k n j ( Y ¯ j Y ¯ ) 2 1 n k j = 1 k i = 1 n j ( Y j i Y ¯ j ) 2 {\displaystyle L={\frac {{\frac {1}{k-1}}\sum _{j=1}^{k}n_{j}({\bar {Y}}_{j}-{\bar {Y}})^{2}}{{\frac {1}{n-k}}\sum _{j=1}^{k}\sum _{i=1}^{n_{j}}(Y_{ji}-{\bar {Y}}_{j})^{2}}}}

annähernd F ( k 1 , n k ) {\displaystyle F(k-1,n-k)} -verteilt mit n {\displaystyle n} die Anzahl aller Beobachtungen:

n = j = 1 k n j {\displaystyle n=\sum _{j=1}^{k}n_{j}} ,

Y ¯ {\displaystyle {\bar {Y}}} der Stichprobenmittelwert über alle Gruppen und Y ¯ j {\displaystyle {\bar {Y}}_{j}} der Stichprobenmittelwert über Gruppe j {\displaystyle j} .

Die Teststatistik bzgl. Y j i {\displaystyle Y_{ji}} ist identisch mit der Teststatistik der einfachen Varianzanalyse (Test auf Gleichheit von k {\displaystyle k} Gruppenmittelwerten). Durch die Transformation von X j i {\displaystyle X_{ji}} auf Y j i {\displaystyle Y_{ji}} sind die Gruppenmittelwerte

Y ¯ j = 1 n j i = 1 n j Y j i = 1 n j i = 1 n j | X j i X ¯ j | {\displaystyle {\bar {Y}}_{j}={\frac {1}{n_{j}}}\sum _{i=1}^{n_{j}}Y_{ji}={\frac {1}{n_{j}}}\sum _{i=1}^{n_{j}}|X_{ji}-{\bar {X}}_{j}|}

robuste Schätzfunktionen der Gruppenvarianzen. Die Normalverteilungsannahme für die Varianzanalyse gilt zwar nicht, jedoch haben die Y j i {\displaystyle Y_{ji}} oft eine rechtsschiefe Verteilung für die die Varianzanalyse angewandt werden kann.[4]

Brown–Forsythe-Test

Im Brown–Forsythe-Test wird bei Berechnung von Y j i {\displaystyle Y_{ji}} statt des Gruppenmittelwertes der Gruppenmedian benutzt.[5] Um eine gute Teststärke zu erhalten, muss der Lageparameter in Abhängigkeit von der zugrunde liegenden Verteilung gewählt werden. Brown und Forsythe zeigten in Simulationsstudien, dass der Mittelwert eine gute Wahl ist, wenn die Verteilung symmetrisch und „normale“ Verteilungsenden (Exzess {\displaystyle \approx } 0) hat, z. B. einer Normalverteilung ähnlich ist. Der Median sollte benutzt werden, wenn die Verteilungen stark schief sind, und der getrimmte Mittelwert, wenn die Verteilung schwere Verteilungsenden hat (Exzess<0).

Einzelnachweise

  1. Howard Levene: Robust tests for equality of variances. In: Ingram Olkin, Harold Hotelling et al. (Hrsg.): Contributions to Probability and Statistics: Essays in Honor of Harold Hotelling. Stanford University Press, 1960, ISBN 0-8047-0596-8, S. 278–292}. 
  2. Joseph L. Gastwirth, Yulia R. Gel, Weiwen Miao: The impact of Levene´s test of equality of variances on statistical theory and practice. In: Statistical Science. Band 24, Nr. 3, S. 343–360, doi:10.1214/09-STS301. 
  3. Jürgen Janssen, Wilfried Laatz: Statistische Datenanalyse mit SPSS für Windows. 8. Auflage. Springer Verlag, 2007, S. 246. 
  4. Maxwell J. Roberts, Riccardo Russo: Student’s Guide to Analysis of Variance. Routledge Chapman & Hall, 1999, ISBN 978-0-415-16565-5, S. 71. 
  5. Morton B. Brown, Alan B. Forsythe: Robust tests for equality of variances. In: Journal of the American Statistical Association. Band 69, 1974, S. 364–367, doi:10.1080/01621459.1974.10482955. 

Literatur

  • Biostatistik: Eine Einführung für Biowissenschaftler. (2008). München: Pearson Studium. S. 150–154.