我想通过for循环做批量线性回归，然后提取线性回归中关于系数部分的结果，老师给的是一个因变量和不同的自变量分别做线性回归，代码如下：
data <- read.csv("data.csv")
result <- data.frame(estimate=0,se=0,p=0)
O3_lag0 <- data[,1]
for (i in 1:5) {
  linear <- lm(O3_lag0~data[,i+4])
  result[i,1]<-summary(linear)$coefficients[2,1]
  result[i,2]<-summary(linear)$coefficients[2,3]
  result[i,3]<-summary(linear)$coefficients[2,4]
  }
result


我的问题，我想根据这个代码作进一步的改动，就是多个因变量和多个自变量分别做线性回归，目前这个我通过2个for循环实现了，可是问题是关于提取系数部分的结果不会写w(ﾟДﾟ)w，不知道怎么才能依次把得到的结果提取出来，下面是我写的渣渣代码，就是标红部分有问题，不会写o(╥﹏╥)o：

data <- read.csv("data.csv")
result <- data.frame(estimate=0,se=0,p=0)
for (j in 1:4) {
  for (i in 1:5) {
  linear <- lm(data[,j]~data[,i+4])
  result[i,1]<-summary(linear)$coefficients[2,1]
  result[i,2]<-summary(linear)$coefficients[2,3]
  result[i,3]<-summary(linear)$coefficients[2,4]
}
}
result


附上我的数据~