Allowing for higher-order fits
Make it work so that color and column can both represent the same factor
adding func_name to regplot for more flexibility
let corr_func return a single value
adding linear models tutorial