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