Choisir le meilleur modèle

Modélisation avec tidymodels en R

David Svancer

Data Scientist

Résultats détaillés de l'ajustement

La fonction collect_metrics() donne un résumé par défaut

  • Avec summarize = FALSE, vous obtenez tous les résultats d'ajustement des hyperparamètres
dt_tuning %>% 
  collect_metrics(summarize = FALSE)
# A tibble: 150 x 8
 id     cost_complexity tree_depth min_n .metric  ...  .estimate  .config
<chr>        <dbl>         <int>   <int>  <chr>   ...    <dbl>      <chr>  
Fold01    0.0000000758     14       39    sens    ...     0.75     Model1 
Fold01    0.0000000758     14       39    spec    ...     0.906    Model1 
Fold01    0.0000000758     14       39    roc_auc ...     0.888    Model1 
.....     ............     ..       ..    ......  ...     .....    ......
Fold10    0.00380          5        36    roc_auc ...     0.789    Model5
Modélisation avec tidymodels en R

Explorer les résultats d'ajustement

Avec summarise = FALSE, collect_metrics() renvoie un tibble

  • Exploration facile avec dplyr
  • Explorer l'AUC ROC
    • Sélectionner la métrique roc_auc
    • Regrouper par la colonne id
    • Calculer les stats de synthèse de .estimate
dt_tuning %>% 
  collect_metrics(summarize = FALSE) %>% 

filter(.metric == 'roc_auc') %>%
group_by(id) %>%
summarize(min_roc_auc = min(.estimate), median_roc_auc = median(.estimate), max_roc_auc = max(.estimate))
# A tibble: 10 x 4
 id     min_roc_auc  median_roc_auc  max_roc_auc
<chr>      <dbl>          <dbl>       <dbl>
Fold01     0.830          0.885       0.888
Fold02     0.857          0.882       0.885
Fold03     0.818          0.836       0.836
......     ....           ....        ....
Fold10     0.762          0.790       0.813
Modélisation avec tidymodels en R

Voir les meilleurs modèles

La fonction show_best()

  • Affiche les n meilleurs modèles selon la moyenne de metric
  • Model1 est le gagnant
dt_tuning %>% 
  show_best(metric = 'roc_auc', n = 5)
# A tibble: 5 x 9
cost_complexity  tree_depth  min_n  .metric .estimator   mean    n    std_err  .config
    <dbl>           <int>    <int>    <chr>   <chr>      <dbl>  <int>  <dbl>    <chr>
0.0000000758         14       39     roc_auc  binary     0.827   10   0.0147   Model1 
0.00380               5       36     roc_auc  binary     0.825   10   0.0146   Model5 
0.0243                5       34     roc_auc  binary     0.823   10   0.0147   Model2 
0.00000443           11       8      roc_auc  binary     0.816   10   0.00786  Model3 
0.000000600           3       5      roc_auc  binary     0.814   10   0.0131   Model4
Modélisation avec tidymodels en R

Sélection d'un modèle

La fonction select_best()

  • Transmettre les résultats dt_tuning à select_best()
  • Choisir la metric pour évaluer la performance

 

Retourne un tibble avec le meilleur modèle et ses hyperparamètres

best_dt_model <- dt_tuning %>% 
  select_best(metric = 'roc_auc')

best_dt_model

 

# A tibble: 1 x 4
cost_complexity tree_depth  min_n  .config
     <dbl>         <int>    <int>   <chr>  
0.0000000758        14       39     Model1
Modélisation avec tidymodels en R

Finaliser le flux de travail

La fonction finalize_workflow() finalise un workflow contenant un modèle avec hyperparamètres à ajuster

  • Transmettre l'objet workflow
  • Un tibble à une ligne avec les valeurs finales des hyperparamètres
    • Les noms de colonnes doivent correspondre aux hyperparamètres du modèle

 

Retourne un objet workflow avec hyperparamètres fixés

final_leads_wkfl <- leads_tune_wkfl %>% 
  finalize_workflow(best_dt_model)

final_leads_wkfl
== Workflow ========================================
Preprocessor: Recipe
Model: decision_tree()
-- Preprocessor ------------------------------------
3 Recipe Steps
* step_corr()
* step_normalize()
* step_dummy()
-- Model --------------------------------------------
Decision Tree Model Specification (classification)
Main Arguments:
  cost_complexity = 0.0000000758
  tree_depth = 14
  min_n = 39
Computational engine: rpart
Modélisation avec tidymodels en R

Ajustement du modèle

L'objet workflow finalisé peut être entraîné avec last_fit() et l'objet de découpage des données d'origine, leads_split

 

En arrière-plan

  • Création des ensembles d'entraînement et de test
  • recipe entraînée et appliquée
  • Arbre de décision ajusté entraîné sur tout l'ensemble d'entraînement
  • Prédictions et métriques sur les données de test
leads_final_fit <- final_leads_wkfl %>% 
  last_fit(split = leads_split)

leads_final_fit %>% collect_metrics()

 

# A tibble: 2 x 3
  .metric  .estimator .estimate
  <chr>    <chr>          <dbl>
1 accuracy binary         0.771
2 roc_auc  binary         0.793
Modélisation avec tidymodels en R

Passons à la pratique !

Modélisation avec tidymodels en R

Preparing Video For Download...