从coeftest对象列表中提取列
|
是否存在可以从“ 0”对象中提取两列或更多列的函数？一次很容易实现一个“ 0”对象，但是我可以对一个列表做同样的事情（一个“ 2”循环除外）吗？
> # 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
[[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
没有找到相关结果
已邀请:
2 个回复
唤副埂侧壬
是错误的子集功能。请注意，当您在列表上使用时，要使用的内容是列表的组成部分，使用将获得的位，其中是第i个组成部分。 这样，您只需在调用中将的位。由于的参数，它有点棘手，但with8ѭ可以轻松实现：> 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荆怖赡
> 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