fit_linear_step function

Function to fit linear model for variable of interest