Sintaxis adecuada para el término de interacción aleatoria única lme() cuando se anidan 2 efectos aleatorios
-
21-12-2019 - |
Pregunta
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)
Solución
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)
Licenciado bajo: CC-BY-SA con atribución
No afiliado a StackOverflow