ベイズ最適化法を使用してxgboostのパラメーターを最適化します

通常行うグリッド検索とランダム検索に加えて、ベイズ最適化手法が非常に優れていることがわかりました。それから試してみたところ、効果が非常に優れていることがわかりました。ここでコードを共有します。

ベイズ最適化は、目的関数の過去の評価結果に基づいて置換関数(確率モデル)を確立することにより、目的関数を最小化する値を見つけます。ベイジアン法とランダム検索またはグリッド検索の違いは、次のハイパーパラメータのセットを試行するときに前の評価結果を参照するため、無駄な作業を大幅に節約できることです。

ハイパーパラメータの評価は、モデルをトレーニングするために評価するハイパーパラメータを使用する必要があるため、非常にコストがかかります。また、多くの深層学習モデルは、トレーニングを完了してモデルを評価するのに数時間と数日かかるため、コストがかかります。ベイジアンチューニングは、常に更新される確率モデルを使用して、過去の結果を推測することにより、有望なハイパーパラメータに「焦点を合わせ」ます。

1ライブラリパッケージをインポートします

from skopt import BayesSearchCV
import xgboost as xgb
from sklearn.model_selection import train_test_split
import pandas as pd
from sklearn.model_selection import StratifiedKFold
import numpy as np
from sklearn.utils import shuffle

2データをロードします


train_path='ads_train.csv'
train_data=pd.read_csv(train_path)

3データセット機能の処理

train_data = shuffle(train_data)
X=train_data[['isbuyer', 'buy_freq', 'visit_freq', 'buy_interval',
       'sv_interval', 'expected_time_buy', 'expected_time_visit', 'last_buy', 'multiple_buy', 'multiple_visit', 'uniq_urls',
       'num_checkins']]
Y=train_data[['y_buy']]
X_train,X_test,y_train,y_test=train_test_split(X,Y,test_size=0.2)

最適化されたコード

ITERATIONS=100
# Classifier
bayes_cv_tuner = BayesSearchCV(
    estimator = xgb.XGBClassifier(
        n_jobs = 1,
        objective = 'binary:logistic',
        eval_metric = 'auc',
        silent=1,
        tree_method='approx'
    ),
    search_spaces = {
        'learning_rate': (0.01, 1.0, 'log-uniform'),
        'min_child_weight': (0, 10),
        'max_depth': (0, 50),
        'max_delta_step': (0, 20),
        'subsample': (0.01, 1.0, 'uniform'),
        'colsample_bytree': (0.01, 1.0, 'uniform'),
        'colsample_bylevel': (0.01, 1.0, 'uniform'),
        'reg_lambda': (1e-9, 1000, 'log-uniform'),
        'reg_alpha': (1e-9, 1.0, 'log-uniform'),
        'gamma': (1e-9, 0.5, 'log-uniform'),
        'min_child_weight': (0, 5),
        'n_estimators': (50, 100),
        'scale_pos_weight': (1e-6, 500, 'log-uniform')
    },    
    scoring = 'roc_auc',
    cv = StratifiedKFold(
        n_splits=5,
        shuffle=True,
        random_state=42
    ),
    n_jobs = 6,
    n_iter = ITERATIONS,   
    verbose = 0,
    refit = True,
    random_state = 42
)

def status_print(optim_result):
    """Status callback durring bayesian hyperparameter search"""
    
    # Get all the models tested so far in DataFrame format
    all_models = pd.DataFrame(bayes_cv_tuner.cv_results_)    
    
    # Get current parameters and the best parameters    
    best_params = pd.Series(bayes_cv_tuner.best_params_)
    print('Model #{}\nBest ROC-AUC: {}\nBest params: {}\n'.format(
        len(all_models),
        np.round(bayes_cv_tuner.best_score_, 4),
        bayes_cv_tuner.best_params_
    ))
    print(dict(bayes_cv_tuner.best_params_))
    
    
    # Save all model results
    clf_name = bayes_cv_tuner.estimator.__class__.__name__
    all_models.to_csv(clf_name+"_cv_results.csv")

result = bayes_cv_tuner.fit(X.values, Y.values, callback=status_print)

参照

xgBoost自動機械学習ハイパーパラメーター調整のベイズハイパーパラメーター調整
(ベイズ最適化)

おすすめ

転載: blog.csdn.net/w5688414/article/details/113484064
おすすめ