2つのランダム効果がネストされている場合のLME（）単一ランダム相互作用項目のための適切な構文

Question

lme()構文に関する質問です。私の応答変数は '応答'です。私の固定変数は「年」です。私は2つのランダム変数を持っています： '学校'の中に入れ子になっています。

year*schoolの相互作用を含みたいが、year*studentを含めたくない。

これは私がこれまでの構文ですが、これは私が望んでいる2つのランダムな相互作用を含むようです。

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)

.