我通过使用相同的自变量集来运行近40个逻辑回归,但是有40个不同的因变量。我想从40条回归线中提取系数,并创建一个数据框来绘制点图。我认为最好将其可视化,而不是将所有回归线放在同一个表中。
我无法弄清楚如何提取自变量的效果(每个方程的相同类型的系数)
如果你可以帮助我,我会感激不尽!
这是一个可重复的例子。
set.seed(10)
y <- matrix(rnorm(10000 * 14), ncol = 14)
x <- matrix(rnorm(10000 * 2), ncol = 2)
res <- lapply(1:14, function(i) lm(y[, i] ~ x))
答案 0 :(得分:1)
尝试使用列表do.call
rbind
和res
> do.call(rbind, lapply(res, function(x) coef(x)))
(Intercept) x1 x2
[1,] 1.440115e-03 -0.0198232209 -0.0005720764
[2,] -2.227644e-02 -0.0134155339 -0.0092420757
[3,] 3.535811e-03 -0.0284229117 0.0140198529
[4,] -2.031279e-02 0.0032004789 -0.0036719760
[5,] -1.127532e-02 -0.0004463859 -0.0116754425
[6,] -1.369851e-02 0.0174797415 0.0112791379
这是一个矩阵,使用data.frame
获取data.frame