网格搜索

对给定参数进行组合,用某标准进行评价,只适合小数据集

class sklearn.model_selection.GridSearchCV(estimator, param_grid, scoring=None, fit_params=None, n_jobs=None, iid=’warn’, refit=True, cv=’warn’, verbose=0, 
pre_dispatch=‘2*n_jobs’, error_score=’raise-deprecating’, return_train_score=’warn’)

参数

estimator:一个学习器对象,它必须有.fit方法用于学习,.predict方法用于预测,.score方法用于评分

param_grid:字典或者字典的列表,每个字典是学习器的一个参数,key是参数名,value是备选参数序列

scoring:评分函数,可以是一个字符串,或一个可调用对象,或者None

    如果是None,使用学习器的.score方法,

    如果是字符串,可以为如下值:

      accuracy:用 metrics.accuracy_score 评分函数

      average_precision:用 metrics.average_precision_score 评分函数

      f1系列:用 metrics.f1_score 评分函数

        f1:用于二分类

        f1_micro:micro_averaged

        f1_macro:macro_averaged

        f1_weighted:weighted average

        f1_samples:by multilabel sample

      log_loss:用 metrics.log_loss 评分函数

      precision:用 metrics.precision_score 评分函数,具体形式类似f1系列

      recall:用 metrics.recall_score 评分函数,具体形式类似f1系列

      roc_auc:用的是 metrics.roc_auc_score 评分函数

      adjusted_rand_score:

      mean_absolute_error:

      mean_squared_error:

      median_absolute_error:

      r2:

fit_params:字典形式,用于给学习器的fit传递参数

n_jobs:并行性,默认为-1或None,表示派发任务到所有计算机的cpu上

iid:样本独立同分布,bool值

refit:bool值,如果为True,表示在参数优化之后,使用整个数据集来重新训练该最优的estimator

cv:一个整数,k折交叉生成器,一个迭代器,或者None

  如果为None,使用3折交叉验证

  如果为整数k,使用k这交叉验证

  如果为k折交叉验证生成器

  迭代器

verbose:整数,控制输出日志的内容,该值越大,输出越多的内容

pre_dispatch:整数或者字符串,控制并行执行时分发的总的任务量

error_score:

return_train_score:

 

属性

grid_scores_:命名元组组成的列表,列表中每个元素代表参数组合的得分

best_estimator_:一个学习器对象,代表最后的最优学习器

best_score_:最佳学习器的得分

best_params_:最佳参数组合

 

方法

fit:执行参数优化  

predict:使用学习到的最佳学习器预测数据

predict_log_proba(x):使用最佳学习器来预测属于每个类别的概率的对数值

predict_proba(x):使用最佳学习器来预测属于每个类别的概率

score:最佳学习器的评分

 

示例代码

from sklearn.datasets import load_digits
from sklearn.linear_model import LogisticRegression
from sklearn.metrics import classification_report
from sklearn.model_selection import train_test_split, GridSearchCV

###
digits = load_digits()
X_train, X_test, y_train, y_test = train_test_split(digits.data, digits.target, test_size=0.25, random_state=0, stratify=digits.target)

###
tuned_parameters = [{'penalty':['l1', 'l2'], 'C':[0.01, 0.05, 0.1, 0.5, 1, 5, 10, 50, 100], 'solver':['liblinear'], 'multi_class':['ovr']},
                    {'penalty':['l2'], 'C':[0.01, 0.05, 0.1, 0.5, 1, 5, 10, 50, 100], 'solver':['lbfgs'], 'multi_class':['ovr', 'multinomial']}]

clf = GridSearchCV(LogisticRegression(tol=1e-6), tuned_parameters, cv=10)
clf.fit(X_train, y_train)
print(clf.best_params_)                     # {'penalty': 'l2', 'multi_class': 'multinomial', 'C': 0.01, 'solver': 'lbfgs'}

for i in clf.grid_scores_:
    print(i)

# mean: 0.93022, std: 0.01571, params: {'penalty': 'l1', 'multi_class': 'ovr', 'C': 0.01, 'solver': 'liblinear'}
# mean: 0.96511, std: 0.01095, params: {'penalty': 'l2', 'multi_class': 'ovr', 'C': 0.01, 'solver': 'liblinear'}
# mean: 0.96437, std: 0.01082, params: {'penalty': 'l1', 'multi_class': 'ovr', 'C': 0.05, 'solver': 'liblinear'}
# mean: 0.96511, std: 0.01164, params: {'penalty': 'l2', 'multi_class': 'ovr', 'C': 0.05, 'solver': 'liblinear'}
# ...
# mean: 0.95100, std: 0.01262, params: {'penalty': 'l1', 'multi_class': 'ovr', 'C': 1, 'solver': 'liblinear'}
# mean: 0.95323, std: 0.01390, params: {'penalty': 'l2', 'multi_class': 'ovr', 'C': 1, 'solver': 'liblinear'}
# mean: 0.95249, std: 0.01096, params: {'penalty': 'l1', 'multi_class': 'ovr', 'C': 5, 'solver': 'liblinear'}
# mean: 0.95100, std: 0.01237, params: {'penalty': 'l2', 'multi_class': 'ovr', 'C': 5, 'solver': 'liblinear'}
# ...
# mean: 0.94358, std: 0.01260, params: {'penalty': 'l2', 'multi_class': 'ovr', 'C': 50, 'solver': 'lbfgs'}
# mean: 0.95768, std: 0.01487, params: {'penalty': 'l2', 'multi_class': 'multinomial', 'C': 50, 'solver': 'lbfgs'}
# mean: 0.93690, std: 0.01520, params: {'penalty': 'l2', 'multi_class': 'ovr', 'C': 100, 'solver': 'lbfgs'}
# mean: 0.95620, std: 0.01314, params: {'penalty': 'l2', 'multi_class': 'multinomial', 'C': 100, 'solver': 'lbfgs'}

print(clf.score(X_test, y_test))            # 0.973333333333

pred = clf.predict(X_test)
print(classification_report(y_test, pred))
#               precision    recall  f1-score   support
#
#           0       1.00      1.00      1.00        45
#           1       0.87      0.98      0.92        46
#           2       1.00      0.98      0.99        44
#           3       1.00      1.00      1.00        46
#           4       1.00      0.96      0.98        45
#           5       0.96      0.98      0.97        46
#           6       1.00      0.98      0.99        45
#           7       1.00      1.00      1.00        45
#           8       0.93      0.88      0.90        43
#           9       1.00      0.98      0.99        45
#
# avg / total       0.97      0.97      0.97       450

 

随机搜索

随机搜索的原理和网格搜索是一致的,只是随机搜索在取参数值时是随机取值,其主要解决以下问题

1. 参数很多

2. 参数取值是连续的,或者参数取值范围很大,或者参数可取值非常多

随机取值其实不是完全随机,如根据参数取值的分布进行采样,其实也可以自己定义

class sklearn.model_selection.RandomizedSearchCV(estimator, param_distributions, n_iter=10, scoring=None, fit_params=None, n_jobs=None, iid=’warn’, refit=True, 
cv=’warn’, verbose=0, pre_dispatch=‘2*n_jobs’, random_state=None, error_score=’raise-deprecating’, return_train_score=’warn’)[source]

参数

参数与gridsearchCV大致相同,简单解释下

param_distributions:相当于上面的param_grid,只是字典的值是个分布,通常使用scipy.stats模块提供的分布类,如 scipy.expon 指数分布,scipy.gamma gamma分布,scipy.uniform 均匀分布等

n_iter:整数,指定每个参数采样的数量,通常该值越大,参数优化效果越好,但计算量也越大

{'C': scipy.stats.expon(scale=100), 'gamma': scipy.stats.expon(scale=.1),'kernel': ['rbf'], 'class_weight':['balanced', None]}

 

属性

cv_results_:不同参数组合的得分

best_params_,best_score_ 等

 

示例代码

from sklearn.datasets import load_iris
from sklearn.neighbors import KNeighborsClassifier  # 要估计的是knn里面的参数,包括k的取值和样本权重分布方式 
from sklearn.model_selection import GridSearchCV,RandomizedSearchCV 

iris = load_iris()
X = iris.data  
y = iris.target 

k_range = range(1, 31)  # 优化参数k的取值范围
weight_options = ['uniform', 'distance']  # 代估参数权重的取值范围。uniform为统一取权值,distance表示距离倒数取权值
param_grid = {'n_neighbors':k_range,'weights':weight_options}  # 字典中的key值必须是分类算法的函数的参数名
print(param_grid)

knn = KNeighborsClassifier(n_neighbors=5)  

# ================================网格搜索=======================================
# 这里GridSearchCV的参数形式和cross_val_score的形式差不多,其中param_grid是parameter grid所对应的参数
# GridSearchCV中的n_jobs设置为-1时,可以实现并行计算(如果你的电脑支持的情况下)
grid = GridSearchCV(estimator = knn, param_grid = param_grid, cv=10, scoring='accuracy') #针对每个参数对进行了10次交叉验证。scoring='accuracy'使用准确率为结果的度量指标。可以添加多个度量指标
grid.fit(X, y)

print('网格搜索-度量记录:',grid.cv_results_)  # 包含每次训练的相关信息
print('网格搜索-最佳度量值:',grid.best_score_)  # 获取最佳度量值
print('网格搜索-最佳参数:',grid.best_params_)  # 获取最佳度量值时的代定参数的值。是一个字典
print('网格搜索-最佳模型:',grid.best_estimator_)  # 获取最佳度量时的分类器模型

# 使用获取的最佳参数生成模型,预测数据
knn = KNeighborsClassifier(n_neighbors=grid.best_params_['n_neighbors'], weights=grid.best_params_['weights'])  
knn.fit(X, y)  
print(knn.predict([[3, 5, 4, 2]]))  

# =====================================随机搜索===========================================
rand = RandomizedSearchCV(knn, param_grid, cv=10, scoring='accuracy', n_iter=10, random_state=5)  #
rand.fit(X, y)

print('随机搜索-度量记录:',grid.cv_results_)  # 包含每次训练的相关信息
print('随机搜索-最佳度量值:',grid.best_score_)  # 获取最佳度量值
print('随机搜索-最佳参数:',grid.best_params_)  # 获取最佳度量值时的代定参数的值。是一个字典
print('随机搜索-最佳模型:',grid.best_estimator_)  # 获取最佳度量时的分类器模型

# 使用获取的最佳参数生成模型,预测数据
knn = KNeighborsClassifier(n_neighbors=grid.best_params_['n_neighbors'], weights=grid.best_params_['weights'])  
knn.fit(X, y)  
print(knn.predict([[3, 5, 4, 2]]))  

# =====================================自定义度量===========================================
from sklearn import metrics
# 自定义度量函数
def scorerfun(estimator, X, y):
    y_pred = estimator.predict(X)
    return metrics.accuracy_score(y, y_pred)

rand = RandomizedSearchCV(knn, param_grid, cv=10, scoring=scorerfun, n_iter=10, random_state=5)  #
rand.fit(X, y)

print('随机搜索-最佳度量值:',grid.best_score_)  # 获取最佳度量值

 

 

参考资料:

https://scikit-learn.org/stable/modules/generated/sklearn.model_selection.GridSearchCV.html#sklearn.model_selection.GridSearchCV

https://scikit-learn.org/stable/modules/generated/sklearn.model_selection.RandomizedSearchCV.html

https://blog.csdn.net/luanpeng825485697/article/details/79831703