final_model function

The function calculates all models and find final model according Gilmour