Ottimizzazione degli iperparametri in caret

Ottimizzazione degli iperparametri in R

Dr. Shirin Elsinghorst

Senior Data Scientist

Dataset elettori USA elezioni 2016

Dividi in training e test set

library(tidyverse)
glimpse(voters_train_data)

Observations: 6,692
Variables: 42
$ turnout16_2016       <chr> "Did not vote", "Did not vote", "Did not vote", "Did not vote", ...
$ RIGGED_SYSTEM_1_2016 <int> 2, 2, 3, 2, 2, 3, 3, 1, 2, 3, 4, 4, 4, 3, 1, 2, 2, 2, 3, 2, 1, 2, 3, 2, 1, ...
$ RIGGED_SYSTEM_2_2016 <int> 3, 3, 2, 2, 3, 3, 2, 2, 1, 2, 4, 2, 3, 2, 3, 4, 3, 2, 2, 2, 4, 1, 2, 2, 3, ...
$ RIGGED_SYSTEM_3_2016 <int> 1, 1, 3, 1, 1, 1, 2, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 2, 3, 1, 3, 2, 1, 1, ...
$ RIGGED_SYSTEM_4_2016 <int> 2, 1, 2, 2, 2, 2, 2, 2, 1, 3, 3, 1, 3, 3, 1, 3, 3, 2, 1, 1, 1, 2, 1, 2, 2, ...
$ RIGGED_SYSTEM_5_2016 <int> 1, 2, 2, 2, 2, 3, 1, 1, 2, 3, 2, 2, 1, 3, 1, 1, 2, 2, 1, 2, 1, 2, 2, 2, 1, ...
$ RIGGED_SYSTEM_6_2016 <int> 1, 1, 2, 1, 2, 2, 2, 1, 2, 2, 1, 3, 1, 3, 1, 1, 1, 2, 1, 1, 1, 2, 2, 2, 1, ...
$ track_2016           <int> 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 2, 1, 2, 1, 1, 2, 2, 3, 2, 2, 2, 2, 3, 2, 2, ...
...

Alleniamo un altro modello con caret

Stochastic Gradient Boosting

library(caret)
library(tictoc)
fitControl <- trainControl(method = "repeatedcv", number = 3, repeats = 5)
tic()
set.seed(42)
gbm_model_voters <- train(turnout16_2016 ~ ., 
                   data = voters_train_data, 
                   method = "gbm", 
                   trControl = fitControl,
                   verbose = FALSE)
toc()

32.934 sec elapsed

Alleniamo un altro modello con caret

gbm_model_voters

Stochastic Gradient Boosting 
...
Risultati del resampling sui parametri di tuning:
  interaction.depth  n.trees  Accuracy   Kappa        
  1                   50      0.9604603  -0.0001774346
 ...
Il parametro di tuning 'shrinkage' è stato mantenuto costante a 0.1
Il parametro di tuning 'n.minobsinnode' è stato mantenuto costante a 10

L'Accuracy è stata usata per selezionare il modello ottimale con il valore più alto.
I valori finali usati per il modello sono stati n.trees = 50,
interaction.depth = 1, shrinkage = 0.1 e n.minobsinnode = 10.

Ricerca a griglia cartesiana con caret

Definisci una griglia cartesiana di iperparametri:

man_grid <-  expand.grid(n.trees = c(100, 200, 250), interaction.depth = c(1, 4, 6), 
                         shrinkage = 0.1, n.minobsinnode = 10)

fitControl <- trainControl(method = "repeatedcv", number = 3, repeats = 5)
tic()
set.seed(42)
gbm_model_voters_grid <- train(turnout16_2016 ~ ., 
                   data = voters_train_data, 
                   method = "gbm", 
                   trControl = fitControl,
                   verbose = FALSE,
                   tuneGrid = man_grid)
toc()

85.745 sec elapsed

Ricerca a griglia cartesiana con caret

gbm_model_voters_grid

Stochastic Gradient Boosting 
...
Risultati del resampling sui parametri di tuning:
  interaction.depth  n.trees  Accuracy   Kappa      
  1                  100      0.9603108  0.000912769
 ...
Il parametro di tuning 'shrinkage' è stato mantenuto costante a 0.1
Il parametro di tuning 'n.minobsinnode' è stato mantenuto costante a 10

L'Accuracy è stata usata per selezionare il modello ottimale con il valore più alto.
I valori finali usati per il modello sono stati n.trees = 100, 
interaction.depth = 1, shrinkage = 0.1 e n.minobsinnode = 10.

Grafici dei modelli per iperparametri

plot(gbm_model_voters_grid)

plot(gbm_model_voters_grid, 
     metric = "Kappa", 
     plotType = "level")

Provalo tu!

Ottimizzazione degli iperparametri in R