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