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

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

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

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

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

r interaction mixed-models
person It Figures    schedule 03.06.2014    source источник
comment
Я не вижу, как ученик соотносится со школой: не могли бы вы показать, как выглядят ваши данные?   -  person Vincent    schedule 03.06.2014
comment
Чтобы быть более точным, я бы сказал, что у вас есть две группирующие переменные (школа и ученик [в школе]) и либо два случайных эффекта, либо три в зависимости от того, как вы ставите задачу (перехват по школе, год за школой, перехват учеником, где первые два могут быть объединены в векторнозначный случайный эффект)   -  person Ben Bolker    schedule 03.06.2014

Ответы (1)


arrow_upward
1
arrow_downward

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

lmer(response ~ year + (year | school) + (1 | student), data = data)
person Aaron left Stack Overflow    schedule 03.06.2014