レベル3のランダム切片をR

Question

に指定する

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のランダム切片を追加する方法を誰もが知っていますか？

Answer 1

(1|SubjectID)という用語を使用すると、異なったベースラインがSubjectIDの異なるインスタンスに対してのみ期待されるようにモデルに指示しています。モデルに固定効果のX1に対する異なる応答を期待するには、(1+X1|SubjectID)を使用します。したがって、ご使用のモデルには

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

が必要です。

ところで、これに関する多くの良い情報がCross Validatedにあります。