reg.fit function

Regression modelling