lm_robust_fit function

Internal method that creates linear fits