7

是否有可以从对象中提取两列或多列的coeftest函数？这很容易一次一个coeftest对象，但是我可以对列表（for()循环除外）做同样的事情吗？

> # meaningless data
> temp <- data.frame(a = rnorm(100, mean = 5), b = rnorm(100, mean = 1),
+                    c = 1:100)
> formulas <- list(a ~ b, a ~ c)
> models <- lapply(formulas, lm, data = temp)
> library(lmtest)
> cts <- lapply(models, coeftest)

> # easy to extract columns one object at a time
> cts[[1]][, 1:2]
              Estimate Std. Error
(Intercept)  5.0314196  0.1333705
b           -0.1039264  0.0987044

> # but more difficult algorithmically
> # either one column
> lapply(cts, "[[", 1)
[[1]]
[1] 5.03142

[[2]]
[1] 5.312007

> # or two
> lapply(cts, "[[", 1:2)
Error in FUN(X[[1L]], ...) : attempt to select more than one element

也许更基本的问题是是否有办法将coeftest对象的肉变成数据框，这将允许我单独提取列，然后使用mapply(). 谢谢！

编辑：我想以第一列和第二列的矩阵（或数据框）结束。

    [[1]]
              Estimate Std. Error
(Intercept)  5.0314196  0.1333705
b           -0.1039264  0.0987044

[[2]]
                Estimate  Std. Error
(Intercept)  5.312007153 0.199485363
c           -0.007378529 0.003429477
4

2 回答 2

12

[[在这种情况下是错误的子集函数。请注意，当您lapply()遍历列表时，您正在操作的是列表的组件，您将获得的位list[[i]]是i第 i 个组件。

因此，您只需要调用中的[, 1:2]一点。由于 的论点，这有点棘手，但很容易通过以下方式实现：cts[[1]][, 1:2]lapply()[lapply()

> lapply(cts, `[`, , 1:2)
[[1]]
                Estimate Std. Error
(Intercept)  4.926679544  0.1549482
b           -0.001967657  0.1062437

[[2]]
               Estimate  Std. Error
(Intercept) 4.849041327 0.204342067
c           0.001494454 0.003512972

注意<space>,前面1:2; 这相当于[ , 1:2].

于 2011-06-15T15:43:48.897 回答
5

我不确定这是否是您想要的，但是如何：

> do.call("rbind", cts)[, 1:2]
                 Estimate  Std. Error
(Intercept)  4.8200993881 0.142381642
b           -0.0421189130 0.092620363
(Intercept)  4.7459340076 0.206372906
c            0.0005770324 0.003547885
于 2011-06-15T15:13:32.207 回答