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

Question

Это вопрос о lme() синтаксис.Моя переменная ответа — «ответ».Моя фиксированная переменная — «год».У меня есть 2 случайные переменные:«ученик», который вложен в «школу».

Я хочу включить year*school взаимодействие, но я НЕ хочу включать year*student один.

Это синтаксис, который у меня есть до сих пор, но, похоже, он включает в себя 2 случайных взаимодействия, а мне нужно только 1.

lme1 = lme(response ~ 1 + year, random = ~ year | school/student, 
            method = "REML", data = data)

Answer 1

Это было бы гораздо проще сделать в lmer из lme4 упаковка.Если предположить, что у каждого учащегося есть уникальный идентификатор, это будет выглядеть так (непроверено).

lmer(response ~ year + (year | school) + (1 | student), data = data)