ols function

Fit linear regression model