This function will take a data frame and the dependent variable and fit all possible combinations of models.
The result will be a data frame of models and test statistics for all the models possible. The test statistics are current set and contain all
the following: R-squared, Adjusted R-squared, Degree of freedom, Residual standard error, AIC, BIC, Durbin-Watson statistic.