Einfaktorielle ANOVA

Einfaktorielle ANOVA ist ein Test auf Differenzen innerhalb der Gruppenmittelwerte

Eine einfaktorielle ANOVA ist eine statistische Methode zum Testen der Nullhypothese (H₀), dass drei oder mehr Populationsmittelwerte gleich sind gegen die alternative Hypothese (H_a), dass mindestens ein Mittelwert sich vom Rest unterscheidet. Unter Verwendung der formellen Notation für statistische Hypothesen schreiben wir für k Mittelwerte:

$ H_0:\mu_1=\mu_2=\cdots=\mu_k $

$ H_a:\mathrm{Nicht\mathrm{\ }alle\ Mittelwerte\ sind\ gleich} $

wobei $\mu_i$ der Mittelwert der unabhängigen Variablen für die i-te Stufe des Faktors ist.

Sie denken sich vielleicht: In Ordnung, aber in welchen Situationen müsste ich feststellen, ob die Mittelwerte mehrerer Populationen gleich oder unterschiedlich sind? Ein häufiges Szenario ist, dass Sie vermuten, dass eine bestimmte unabhängige Prozessvariable sich auf ein wichtiges Ergebnis dieses Prozesses auswirkt. Zum Beispiel haben Sie vielleicht den Verdacht, dass unterschiedliche Fertigungslose, Mitarbeiter oder Rohmaterialchargen die Leistung (ein Qualitäts-Maßstab) eines Produktionsprozesses beeinflussen.

Um Ihren Verdacht zu bestätigen, könnten Sie den Prozess unter Verwendung von drei oder mehr Variationen (Stufen) dieser unabhängigen Variablen (Faktor) ablaufen lassen und anschließend eine Stichprobe an Beobachtungen aus den Ergebnissen jedes Laufs nehmen. Falls Sie beim Vergleich der Mittelwerte aus jeder Gruppe von Beobachtungen mittels einer ANOVA Differenzen feststellen, dann (vorausgesetzt, Sie haben alles richtig gemacht!) haben Sie den Beweis, dass Ihr Verdacht stimmte – der Faktor, den Sie untersucht haben, scheint eine Rolle für das Ergebnis zu spielen!

Beispiel für eine einfaktorielle ANOVA

Sehen wir uns ein Beispiel für eine einfaktorielle ANOVA genauer an. Stellen Sie sich vor, Sie arbeiten für ein Unternehmen, das ein Klebegel herstellt, das in kleinen Gefäßen verkauft wird. Die Viskosität des Gels ist wichtig: ist es zu dick, ist es schwer auftragbar; ist es zu dünn, gibt es Einbußen bei der Haftfestigkeit. Sie haben Beschwerden von einigen unzufriedenen Kunden erhalten, dass die Viskosität Ihres Klebstoffes nicht so gleichbleibend ist wie zuvor. Ihr Vorgesetzter hat Sie gebeten, dies zu untersuchen.

Sie entscheiden, dass es als erster Schritt sinnvoll wäre, die durchschnittliche Viskosität der letzten fünf Fertigungslose zu untersuchen. Falls Sie Differenzen zwischen den Losen finden, würde dies vermutlich bestätigen, dass tatsächlich ein Problem besteht. Es könnte Ihnen ebenfalls helfen, Hypothesen über Faktoren zu entwickeln, die Inkonsistenzen zwischen Losen verursachen könnten.

Sie messen die Viskosität mithilfe eines Instruments, das eine in das Klebstoffgefäß getauchte Spindel rotieren lässt. Dieser Test ergibt eine Messgröße namens Drehmomentwiderstand. Sie testen je fünf Gefäße, die zufällig aus jedem der fünf letzten Lose ausgewählt wurden. Sie erhalten den Drehmomentwiderstands-Messwert für jedes Gefäß und stellen die Daten in einem Diagramm dar.

In dem Diagramm beobachten Sie, dass die Drehmoment-Messwerte der Gefäße aus Los 3 niedriger sind als die Drehmoment-Messwerte der Stichproben aus den anderen Losen. Nachdem Sie die Mittelwerte aus allen Ihren Messwerten berechnet haben, sehen Sie, dass der Mittelwert des Drehmoments für Los 3 bei 26,77 liegt – wesentlich geringer als die anderen vier Lose, deren Mittelwert etwa 30 beträgt.

Tabelle 1: Mittelwert der Drehmoments-Messwerte von fünf getesteten Losen Klebstoff

Die ANOVA-Tabelle

ANOVA-Ergebnisse werden normalerweise in einer ANOVA-Tabelle angezeigt. Eine ANOVA-Tabelle beinhaltet:

• Quelle: die Variationsquellen einschließlich des zu untersuchenden Faktors (in unserem Fall, das Los), Fehler und Summe.
• FG: Freiheitsgrade für jede Variationsquelle.
• Quadratsumme: Summe der Quadrate (SQ) für jede Variationsquelle sowie die Summe aus allen Quellen.
• Mittleres Quadrat: Summe der Quadrate geteilt durch die zugehörigen Freiheitsgrade.
• F-Wert: das mittlere Quadrat des Faktors (Los) geteilt durch das mittlere Quadrat des Fehlers.
• Wahrsch. > F: der p-Wert.

Tabelle 2: ANOVA-Tabelle mit Ergebnissen aus unseren Drehmohment-Messwerten

Wir erklären weiter unten, wie die Komponenten dieser Tabelle abgeleitet werden. Ein wesentliches Element dieser Tabelle, auf das wir uns zuerst konzentrieren, ist der p-Wert. Der p-Wert wird verwendet, um auszuwerten, ob die Nullhypothese zutrifft, dass alle Mittelwerte gleich sind. In unserem Beispiel beträgt der p-Wert (Wahrsch. > F) 0,0012. Dieser geringe p-Wert kann als Beweis dienen, dass die Mittelwerte nicht alle gleich sind. Unsere Stichproben liefern den Beweis dafür, dass zwischen den Drehmomentwiderstands-Mittelwerten von einem oder mehreren der fünf Lose eine Differenz besteht.

Auf Grundlage dieser Ergebnisse entscheiden Sie sich, Los 3 für weitere Tests zurückzubehalten. In Ihrem Bericht würden Sie beispielsweise schreiben: Es wurde das Drehmoment von je fünf Gefäßen des Produkts aus jedem der letzten fünf Fertigungslose gemessen. Eine ANOVA-Analyse ergab, dass die Beobachtungen eine Differenz der Drehmoment-Mittelwerte zwischen den Losen unterstützen (p = 0,0012). Ein Diagramm zeigt, dass Los 3 im Vergleich zu den anderen vier Losen einen geringeren Drehmoment-Mittelwert (26,77) aufwies. Wir werden Los 3 zur weiteren Auswertung zurückbehalten.

Berechnung der einfaktoriellen ANOVA

Betrachten wir nun unser Beispiel für Drehmoment-Messwerte genauer. Erinnern Sie sich, dass wir fünf Lose Material hatten. Aus jedem Los wählten wir zufällig fünf Gefäße zum Test aus. Dies nennt man einfaktorielles Design. Der eine Faktor, das Los, hat fünf Stufen. Jede Stufe wird fünf Mal wiederholt (getestet). Die Ergebnisse des Tests sind unten aufgeführt.

Tabelle 3: Drehmoment-Messwerte nach Los

Um die Berechnungen nachzuvollziehen, die der oben stehenden ANOVA-Tabelle (Tabelle 2) zugrunde liegen, legen wir zunächst die folgenden Definitionen fest:

$n_i$ = Anzahl der Beobachtungen für Behandlung $i$ (in unserem Beispiel, Los $i$)

$N$ = Gesamtanzahl der Beobachtungen

$Y_{ij}$ = Die j-te Beobachtung an der i-ten Behandlung

$\overline{Y}_i$ = Der Stichprobenmittelwert für die i-te Behandlung

$\overline{\overline{Y}}$ = Der Mittelwert aller Beobachtungen (Gesamtmittelwert)

Summe Quadrate

Unter Berücksichtigung dieser Definitionen beginnen wir mit der Spalte „Summe Quadrate“ aus der ANOVA-Tabelle. Die Summe der Quadrate bietet uns einen Weg, die Variabilität innerhalb eines Datensatzes zu quantifizieren, indem wir uns auf die Differenz zwischen den einzelnen Datenpunkten und dem Mittelwert aller Datenpunkte innerhalb dieses Datensatzes konzentrieren. Die unten stehende Formel teilt die Gesamtvariabilität in zwei Teile: die Variabilität aufgrund des Modells oder der Faktorstufen und die Variabilität aufgrund von zufälligen Fehlern.  

$$ \sum_{i=1}^{a}\sum_{j=1}^{n_i}(Y_{ij}-\overline{\overline{Y}})^2\;=\;\sum_{i=1}^{a}n_i(\overline{Y}_i-\overline{\overline{Y}})^2+\sum_{i=1}^{a}\sum_{j=1}^{n_i}(Y_{ij}-\overline{Y}_i)^2 $$

$$ SS(Total)\;     =     \;SS(Factor)\;     +     \;SS(Error) $$

Diese Gleichung mag zwar kompliziert erscheinen, aber sie ist leichter verständlich, wenn wir uns auf die einzelnen Elemente konzentrieren. Die unten stehende Tabelle 4 führt die einzelnen Komponenten der Formel auf und entwickelt sie zu den quadrierten Termen, die die Summe der Quadrate ergeben. Die erste Datenspalte ($Y_{ij}$) enthält die Drehmoment-Messwerte, die wir oben in Tabelle 3 zusammengetragen hatten.

Tabelle 4: Berechnung der Summe der Quadrate

Freiheitsgrade (FG)

Zu jeder Summe der Quadrate gehört eine Größe namens Freiheitsgrade (FG). Die Freiheitsgrade zeigen die Anzahl der unabhängigen Informationen an, die zur Berechnung der einzelnen Summen der Quadrate verwendet werden. Bei einem einfaktoriellen Design mit einem Faktor mit k Stufen (fünf Lose in unserem Beispiel) und insgesamt N Beobachtungen (fünf Gefäße pro Los und insgesamt 25) sind die Freiheitsgrade folgendermaßen:

Tabelle 5: Ermittlung der Freiheitsgrade

Mittlere Quadrate (MQ) und F-Wert

Wir teilen jede Summe der Quadrate durch die entsprechenden Freiheitsgrade, um mittlere Quadrate zu erhalten. Wenn die Nullhypothese wahr ist (d. h. die Mittelwerte sind gleich), sind MQ (Faktor) und MQ (Fehler) beide Schätzwerte der Fehlervarianz und wären etwa gleich groß. Ihr Verhältnis, oder der F-Wert, wäre nahe Eins. Wenn die Nullhypothese nicht wahr ist, dann ist das MQ (Faktor) größer als das MQ (Fehler) und ihr Verhältnis ist größer als 1. In unserem Beispiel des Klebstofftests stellt der berechnete F-Wert von 6,90 einen signifikanten Einwand gegen die Nullhypothese dar, dass die Mittelwerte gleich sind.

Tabelle 6: Berechnung der mittleren Quadrate und des F-Wertes