compare_models function

Compare different models fitted by est_multi_poly