model_selection function

Compare fitted models