Hyperparameters van een RF afstellen

Machine Learning met boomgebaseerde modellen in Python

Elie Kawerk

Data Scientist

Hyperparameters van Random Forests

CART-hyperparameters
aantal estimators
bootstrap
....

Tuning is duur

Hyperparameter-tuning:

rekentechnisch duur,
soms maar kleine verbetering,

Weeg de impact van tuning op het hele project.

RF-hyperparameters in sklearn bekijken

# Importeer RandomForestRegressor 
from sklearn.ensemble import RandomForestRegressor

# Zet seed voor reproduceerbaarheid
SEED = 1

# Instantieer een Random Forests-regressor 'rf' 
rf = RandomForestRegressor(random_state= SEED)

# Bekijk de hyperparameters van rf
rf.get_params()

{'bootstrap': True,
 'criterion': 'mse',
 'max_depth': None,
 'max_features': 'auto',
 'max_leaf_nodes': None,
 'min_impurity_decrease': 0.0,
 'min_impurity_split': None,
 'min_samples_leaf': 1,
 'min_samples_split': 2,
 'min_weight_fraction_leaf': 0.0,
 'n_estimators': 10,
 'n_jobs': -1,
 'oob_score': False,
 'random_state': 1,
 'verbose': 0,
 'warm_start': False}

# Basisimports
from sklearn.metrics import mean_squared_error as MSE
from sklearn.model_selection import GridSearchCV

# Definieer een hyperparameterrooster 'params_rf'
params_rf = {
              'n_estimators': [300, 400, 500],
              'max_depth': [4, 6, 8],
              'min_samples_leaf': [0.1, 0.2],
              'max_features': ['log2', 'sqrt']
             }

# Instantieer 'grid_rf'
grid_rf = GridSearchCV(estimator=rf,
                       param_grid=params_rf, 
                       cv=3,
                       scoring='neg_mean_squared_error',
                       verbose=1,
                       n_jobs=-1)

Zoeken naar de beste hyperparameters

# Fit 'grid_rf' op de trainingset
grid_rf.fit(X_train, y_train)

Fitting 3 folds for each of 36 candidates, totalling 108 fits
[Parallel(n_jobs=-1)]: Done  42 tasks      | elapsed:   10.0s
[Parallel(n_jobs=-1)]: Done 108 out of 108 | elapsed:   24.3s finished
RandomForestRegressor(bootstrap=True, criterion='mse', max_depth=4,
           max_features='log2', max_leaf_nodes=None,
           min_impurity_decrease=0.0, min_impurity_split=None,
           min_samples_leaf=0.1, min_samples_split=2,
           min_weight_fraction_leaf=0.0, n_estimators=400, n_jobs=1,
           oob_score=False, random_state=1, verbose=0, warm_start=False)

Beste hyperparameters ophalen

# Haal de beste hyperparameters uit 'grid_rf'
best_hyperparams = grid_rf.best_params_

print('Beste hyperparameters:\n', best_hyperparams)

Beste hyperparameters:
        {'max_depth': 4,
         'max_features': 'log2', 
         'min_samples_leaf': 0.1,
         'n_estimators': 400}

Prestatie van het beste model evalueren

# Haal het beste model uit 'grid_rf'
best_model = grid_rf.best_estimator_
# Voorspel de testlabels
y_pred = best_model.predict(X_test)
# Evalueer de RMSE op de testset
rmse_test = MSE(y_test, y_pred)**(1/2)
# Print de RMSE van de testset
print('Test set RMSE of rf: {:.2f}'.format(rmse_test))

Test set RMSE of rf: 3.89

Laten we oefenen!

Machine Learning met boomgebaseerde modellen in Python