Add data_analysis.R - benchamark regression models