Как совместить два случайных эффекта отдельно в lme?

Question

Я делаю линейную модель смешанных эффектов REML в пакете nlme. И вот коды, которые работают для меня:

# 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)

То, что я действительно хотел бы сделать, это иметь модель с динамиком и предметом в качестве случайных эффектов отдельно. Я пытался использовать эту формулу:

x4 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker + 1|item)

Однако эта формула дает мне следующее предупреждающее сообщение:

Warning message:
In Ops.factor(speaker, 1) : + not meaningful for factors

Ты хоть представляешь, что это значит? И как я могу подобрать как динамик, так и элемент как отдельные эффекты?

Answer 1

Я думаю, что можно включить два случайных эффекта отдельно (один для динамика и один для времени), используя lme() по следующему коду:

x4 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~ speaker + item -1 | id),

с id переменная более высокого уровня, в которой оба speaker а также item вложены Если у вас нет такой переменной, вы можете ввести ее как новую переменную со значением 1 для всех наблюдений.