test <- function() {
  vec1 <- c(1,2,3)
  vec2 <- c(4,5,6)
  
  line <- list()
  for (i in 1:2) {
    line$i <- lm(vec2 ~ vec1)
    names(line)[i] <- paste0("reg", i)
  }
  # return(line)
  
  coef <- list()
  for (i in 1:2) {
    
    coef$i <- line[i]$coefficient[[1]]
    return(coef$i)
    names(coef)[i] <- paste0("coe", i)
  }
  return(coef)
}

test()