Syntaxe appropriée pour le terme d'interaction aléatoire unique lme() lorsque 2 effets aléatoires sont imbriqués

Question

Il s'agit d'une question sur lme() syntaxe.Ma variable de réponse est « réponse ».Ma variable fixe est « année ».J'ai 2 variables aléatoires :« étudiant » qui est imbriqué dans « école ».

Je veux inclure un year*school interaction, mais je ne veux PAS inclure un year*student un.

C'est la syntaxe que j'ai jusqu'à présent, mais cela semble inclure 2 interactions aléatoires là où je n'en veux qu'une.

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

Answer 1

Ce serait beaucoup plus facile à faire dans lmer du lme4 emballer.En supposant que chaque étudiant possède un identifiant unique, cela ressemblerait à ceci (non testé).

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