cross_fit_robust function

Cross fit robust linear models