#!/usr/bin/env Rscript

get.estimate <- function(data, MB.value){
  my.lm = lm(minutes ~ MB, data=data)
  coeffs = coefficients(my.lm); 
  est.time = coeffs[["(Intercept)"]] + coeffs[["MB"]]*MB.value
  return(data.frame(est.time=est.time, intercept.coeff=coeffs[["(Intercept)"]],
                    MB.coeff=coeffs[["MB"]]))
}

MB.value = 34000
results <- list()
for(file in list.files(pattern=".*\\.csv")){
  d <- read.csv(file)
  est <- get.estimate(d, MB.value)
  est.ext <- cbind(data.frame(file=file), est)
  results <- append(results, list(est.ext))
}
results.d <- do.call(rbind, results)
results.d