Sintaxis adecuada para el término de interacción aleatoria única lme() cuando se anidan 2 efectos aleatorios

Question

Esta es una pregunta sobre lme() sintaxis.Mi variable de respuesta es 'respuesta'.Mi variable fija es 'año'.Tengo 2 variables aleatorias:'estudiante' que está anidado dentro de 'escuela'.

quiero incluir un year*school interacción, pero NO quiero incluir una year*student uno.

Esta es la sintaxis que tengo hasta ahora, pero parece incluir 2 interacciones aleatorias donde solo quiero 1.

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

Answer 1

Esto sería mucho más fácil de hacer en lmer desde el lme4 paquete.Suponiendo que cada estudiante tenga un identificador único, se vería así (no probado).

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