2

lme 出力を正しく解釈していることを確認したかっただけです。この方法で個人内SDを見つけることができると言われたため、ヌルモデルを設定しました(Nelzek、2012 MLMの章)。しかし、この章は R コードではなく、数式で書かれています。

使用したモデルはこちら。SubjectID は、人物のグループ化変数です。

nullmodel <- lme(dv1 ~ 1, data=dataset, random = ~1|subjectID,  na.action=na.exclude, method = "ML")
summary(nullmodel)

そして、これが出力の要約です。従属変数 (dv1) の個人内 SD は切片 StdDev であると想定していますが、確認したかっただけです。または、これを間違って解釈している場合はお知らせください:)

Linear mixed-effects model fit by maximum likelihood
Data: dataset 
AIC      BIC    logLik
34404.33 34423.05 -17199.16

Random effects:
Formula: ~1 | subjectID
          (Intercept) Residual
StdDev:    13.39356   21.844

Fixed effects: dv1 ~ 1 
               Value Std.Error   DF  t-value p-value
(Intercept) 37.54255  1.336242 3677 28.09562       0

Standardized Within-Group Residuals:
       Min         Q1        Med         Q3        Max 
-2.7019919 -0.6909552 -0.1045376  0.6288470  3.4520256
4

0 に答える 0