在现有数据帧中存储有效的回归系数

Question

我正在尝试提高脚本的效率，基本上，我运行了一些线性回归，对于每个拟合模型，我将估计系数和标准误差结果存储在先前创建的数据框中，比如说{{1} }。

因此，在存储任何回归系数之前，数据框results已经具有所需的尺寸。

此外，对于每个results - 回归我做：

i

然后我跑：

mod.fit <- plm(y ~ x1 + x2, index="group", sample)

有没有办法让上述存储步骤更快？

Answer 1

您可以使用矩阵索引：

results[i,1:4] <- summary(m.fit)$coefficients[matrix(c(1,2,1,2,1,1,2,2),ncol=2)]

如果results只有4列宽，则可以消除左侧的1:4。

另外

results[i,] <- summary(m.fit)$coefficients[1:2,1:2]

应该有效，因为R以列优先的顺序存储矩阵。

我建议您使用coef()访问者而非$coefficients，如果它是为summary.plm类定义的...