Beispiel für Modell mit gemischten Effekten anpassen

Ein Wissenschaftler untersucht die Ausbeute von sechs Luzernesorten auf vier zufällig ausgewählten Feldern. Der Ertrag auf den einzelnen Feldern wurde für die einzelnen Sorten aufgezeichnet.

Der Wissenschaftler möchte feststellen, ob sich die Sorte der Luzernepflanzen auf die mittlere Ausbeute auswirkt. Er verfügt über 4 Felder, auf denen er Daten erfassen kann. Der Wissenschaftler möchte jedoch modellieren können, wie die Luzernepflanzen auf Feldern wachsen, die in der Untersuchung nicht eingeschlossen sind. Daher legt er das Feld, auf dem die Luzernepflanzen wachsen, als Zufallsfaktor fest. Der Wissenschaftler verwendet ein Modell mit gemischten Effekten, um sowohl feste als auch zufällige Effekte zu untersuchen.

  1. Öffnen Sie die Beispieldaten Luzerne.MTW.
  2. Wählen Sie Statistik > Varianzanalyse (ANOVA) > Modell mit gemischten Effekten > Modell mit gemischten Effekten anpassen aus.
  3. Geben Sie im Feld Antworten die Spalte Ausbeute ein.
  4. Geben Sie im Feld Zufallsfaktoren (erforderlich) die Spalte Feld ein.
  5. Geben Sie im Feld Feste Faktoren die Spalte Sorte ein.
  6. Klicken Sie auf Grafiken.
  7. Wählen Sie im Feld Residuen für Diagramme die Option Bedingt standardisiert aus.
  8. Wählen Sie im Feld Residuendiagramme die Option Vier-in-Eins aus.
  9. Klicken Sie in den einzelnen Dialogfeldern auf OK.

Interpretieren der Ergebnisse

In der Varianzkomponententabelle ist der p-Wert für „Feld“ gleich 0,124. Im Hypothesentest werden keine Hinweise ersichtlich, die darauf hindeuten, dass sich die Varianzkomponente von 0 unterscheidet. Der p-Wert für die Varianzkomponente für Fehler beträgt 0,003. Da der p-Wert kleiner als das Signifikanzniveau 0,05 ist, kann der Wissenschaftler schlussfolgern, dass die Varianzkomponente für Fehler ungleich 0 ist.

Der p-Wert von annähernd 0 für den Term mit festem Faktor „Sorte“ zeigt, dass sich der Effekt von mindestens einer Luzernesorte signifikant von dem der anderen fünf Sorten unterscheidet.

Die Koeffizienten für die Haupteffekte stellen die Differenz zwischen den einzelnen Stufenmittelwerten und dem Gesamtmittelwert dar. Sorte 1 ist beispielsweise einer Luzerneausbeute zugeordnet, die um ca. 0,385 Einheiten über dem Gesamtmittelwert liegt. Der p-Wert von annähernd 0 für diesen Koeffizienten zeigt, dass sich der Effekt von Sorte 1 signifikant vom Effekt einer anderen Stufe des Terms „Stufe“ unterscheidet. Um zu bestimmen, welche Stufeneffekte statistisch gleich sind und welche Stufeneffekte sich statistisch unterscheiden, beabsichtigt der Wissenschaftler, eine Analyse mit Mehrfachvergleichen für den Term durchzuführen.

Der R2-Wert zeigt, dass das Modell ca. 92 % der Streuung in der Ausbeute erklärt. Der Wert von R-Qd(kor) ist mit einem Wert von annähernd 90,2 % ebenfalls hoch. Der Wissenschaftler vergleicht anhand dieses Werts Modelle mit unterschiedlicher Anzahl von Prädiktoren.

Beobachtungen 1 und 5 sind ungewöhnliche Beobachtungen, da ihre standardisierten Residuen größer als 2 sind. Der Wissenschaftler untersucht die Daten, um sich zu vergewissern, dass die Werte der Antwortvariablen für die betreffenden Beobachtungen korrekt sind.

Die Residuen im Wahrscheinlichkeitsnetz für Normalverteilung folgen annähernd einer Geraden, und die Punkte sind im Diagramm der Residuen vs. Anpassung anscheinend zufällig um 0 gestreut.

Methode

VarianzschätzungEingeschränkte Maximum-Likelihood
DF für feste EffekteKenward-Roger

Faktorinformationen

FaktorTypStufenWerte
FeldZufällig41; 2; 3; 4
SorteFest61; 2; 3; 4; 5; 6

Varianzkomponenten

QuelleVar% von GesamtSE Varz-Wertp-Wert
Feld0,07791972,93%0,0675801,1529960,124
Fehler0,02892427,07%0,0105622,7386130,003
Gesamt0,106843       
-2 Log-Likelihood = 7,736012

Tests auf feste Effekte

TermDF des
Zählers
DF des
Nenners
F-Wertp-Wert
Sorte5,0015,0026,290,000

Zusammenfassung des Modells

SR-QdR-Qd(kor)AICcBIC
0,17007192,33%90,20%12,5413,52

Koeffizienten

TermKoefSE KoefDFt-Wertp-Wert
Konstante3,0945830,1438223,0021,5166920,000
Sorte         
  10,3854170,07762615,004,9650160,000
  20,1454170,07762615,001,8732870,081
  30,1079170,07762615,001,3902050,185
  4-0,3195830,07762615,00-4,1169380,001
  50,3954170,07762615,005,0938380,000

Randanpassungen und Bewertung für ungewöhnliche Beobachtungen

BeobAusbeuteAnpassungResidStd. Resid
14,1000003,4800000,6200002,190221R
54,2200003,4900000,7300002,578808R
R  Großes Residuum

Bedingte Anpassungen und Bewertung für ungewöhnliche Beobachtungen

BeobAusbeuteAnpassungResidStd. Resid
54,2200003,8953390,3246612,400733R
R  Großes Residuum