3

120人の被験者を6回測定した縦断研究を分析するために、Rのlmer()関数（パッケージ）を使用しています。lme4最初に、次のようなモデルを指定しました。

library(lme4)
model1 = lmer(DV ~ 1 + X1*X2  + (1+X1|SubjectID), REML="false")

X1は時変変数 (レベル 1) でX2あり、被験者レベルの変数 (レベル 2) です。

これらの被験者は複数のチーム内にネストされているため、チーム レベル (レベル 3) でランダム インターセプトを含めるようにアドバイスされました。ただし、ランダム切片と勾配の両方を含める方法しか見つけられません。

model2 = lmer(DV ~ 1 + X1*X2  + (1+X1|TeamID/SubjectID), REML="false")

モデル1にレベル3のランダム切片のみを追加する方法を知っている人はいますか?

4

1 に答える 1

1

という用語(1|SubjectID)を使用することで、 の異なるインスタンスに対してのみ異なるベースラインを期待するようにモデルに伝えていますSubjectID。固定効果 に対する this の異なる応答を期待するようにモデルに指示するにはX1、 を使用します(1+X1|SubjectID)。したがって、必要なのは条件だけです

(1|TeamID) + (1+X1|SubjectID)

あなたのモデルで。

ちなみに、 Cross Validatedにはこれに関する良い情報がたくさんあります。

于 2016-11-14T16:49:09.317 に答える