从 R 中的 lme4 摘要中提取一列

Question

我想知道从下面的vc对象中提取(而不是像这里那样打印(Std.Dev.列作为向量的最有效方法是什么？

library(lme4)
library(nlme)
data(Orthodont, package = "nlme")
fm1 <- lmer(distance ~ age + (age|Subject), data = Orthodont)
vc <- VarCorr(fm1) ## extract only the `Std.Dev.` column as a vector

Answer 1

"vc"的str表示它是一个具有单元素"Subject"的list，而"stddev"是一个属性

str(vc)
#List of 1
# $ Subject: num [1:2, 1:2] 6.3334 -0.3929 -0.3929 0.0569
#  ..- attr(*, "dimnames")=List of 2
#  .. ..$ : chr [1:2] "(Intercept)" "age"
#  .. ..$ : chr [1:2] "(Intercept)" "age"
#  ..- attr(*, "stddev")= Named num [1:2] 2.517 0.239  ####

因此，直接提取属性

attr(vc$Subject, "stddev")

残差标准差是外部属性

attr(vc, "sc")
#[1] 1.297364

---

如果我们将它们与c结合起来，我们会得到一个vector

c(attr(vc$Subject, "stddev"), attr(vc, "sc"))
#   (Intercept)         age             
#  2.5166317   0.2385853   1.2973640 

用as.numeric/as.vector换行以删除名称，因为它是named向量

或使用attributes

c(attributes(vc)$sc, attributes(vc$Subject)$stddev)

Answer 2

如果需要列中的三个元素，可以使用：

as.numeric(c(attr(vc[[1]], "stddev"), attr(vc, "sc")))