これはlme()
構文に関する質問です。私の応答変数は「応答」です。私の固定変数は「年」です。私は2つの確率変数を持っています:「学校」内にネストされている「学生」。
year*school
インタラクションを含めたいのですが、含めたくありませんyear*student
。
これは私がこれまでに持っている構文ですが、これには1つだけが必要な2つのランダムな相互作用が含まれているようです。
lme1 = lme(response ~ 1 + year, random = ~ year | school/student,
method = "REML", data = data)
これはlmer
、lme4
パッケージから行うのがはるかに簡単です。各学生が一意の識別子を持っていると仮定すると、次のようになります(テストされていません)。
lmer(response ~ year + (year | school) + (1 | student), data = data)
この記事はインターネットから収集されたものであり、転載の際にはソースを示してください。
侵害の場合は、連絡してください[email protected]
コメントを追加