Is there a faster lm function

Question

I would like to get the slope of a linear regression fit for 1M separate data sets (1M * 50 rows for data.frame, or 1M * 50 for array). Now I am using the lm()

Answer 1

speedlm from speedglm should do it as it works on large data sets.