r - 使用 glmer 提取随机效应分量的标准偏差

Question

我正在使用 glmer 并且我希望提取随机效应(截距和斜率)的方差分量的标准偏差。

我试过使用:

VarCorr(model)

它返回两个标准偏差值(加上相关性)，但我只想提取截距和斜率 SD 值。

我尝试使用:

VarrCorr(model)[1]

提取随机截距 SD，这让我知道:

attr(,"stddev")
(Intercept)        year 
      0.075       0.011 

但我不知道如何将这些提取为单独的元素。

Answer 1

有两种方法可以做到这一点。

## make up a model
library(lme4)
(gm <- glmer(incidence ~ period + (size | herd),
              family = poisson, data = cbpp))

方法一

当前版本 lme4允许你强制一个 VarCorr对象到数据框:

as.data.frame(VarCorr(gm))

然后您可以选择第 1:2 行和第 5 列来提取随机截距和斜率的标准偏差。

方法二

如果你想以老式的方式提取值，你可以使用 attributes :

attributes(VarCorr(gm)$herd)$stddev
(Intercept)        size 
 1.18970662  0.08826278 

如果你想去掉名字(即 (intercept) 、 size )，那么你可以使用 as.numeric或 unname .