Wie werden zwei zufällige Effekte getrennt in lme eingefügt?
Ich mache ein lineares Modell mit gemischten Effekten, das von REML in das nlme-Paket eingepasst wird. Und das sind Codes, die für mich funktionieren:
# Linear mixed-effects model fit by REML (intercept and not slope)
x <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker)
summary(x)
# Linear mixed-effects model fit by REML (slope and no intercept)
x1 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~IV3-1|speaker)
summary(x1)
# Linear mixed-effects model fit by REML (slope and intercept)
x2 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~IV3|speaker)
summary(x2)
#nested random effect
x5 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker/item)
summary(x5)
Was ich wirklich gerne machen würde, ist ein Modell, bei dem sowohl Sprecher als auch Gegenstand als zufällige Effekte getrennt sind. Ich hatte versucht, diese Formel zu verwenden:
x4 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker + 1|item)
Diese Formel gibt mir jedoch die folgende Warnmeldung:
Warning message:
In Ops.factor(speaker, 1) : + not meaningful for factors
Haben Sie eine Ahnung, was das bedeutet? Und wie kann ich Lautsprecher und Objekt als zufällige Effekte separat anpassen?