
以下使用scikit-learn中數據集進行分享。
如果選用隨機森林作為最終的模型,那么找出它的最佳參數可能有1000多種組合的可能,你可以使用使用窮盡的網格搜索(Exhaustive Grid Seaarch)方法,但時間成本將會很高(運行很久...),或者使用隨機搜索(Randomized Search)方法,僅分析超參數集合中的子集合。
該例子以手寫數據集為例,使用支持向量機的方法對數據進行建模,然后調用scikit-learn中validation_surve方法將模型交叉驗證的結果進行可視化。需要注意的是,在使用validation_curve方法時,只能驗證一個超參數與模型訓練集和驗證集得分的關系(即二維的可視化),而不能實現多參數與得分間關系的可視化。以下搜索的參數是gamma,需要給定參數范圍,用param_range進行傳遞,評分策略用scoring參數進行傳遞。其代碼示例如下所示:
print(__doc__) import matplotlib.pyplot as plt import numpy as np from sklearn.datasets import load_digits from sklearn.svm import SVC from sklearn.model_selection import validation_curve X, y = load_digits(return_X_y=True) param_range = np.logspace(-6, -1, 5) train_scores, test_scores = validation_curve( SVC(), X, y, param_name="gamma", param_range=param_range, scoring="accuracy", n_jobs=1) train_scores_mean = np.mean(train_scores, axis=1) train_scores_std = np.std(train_scores, axis=1) test_scores_mean = np.mean(test_scores, axis=1) test_scores_std = np.std(test_scores, axis=1) plt.title("Validation Curve with SVM") plt.xlabel(r"$\gamma$") plt.ylabel("Score") plt.ylim(0.0, 1.1) lw = 2 plt.semilogx(param_range, train_scores_mean, label="Training score", color="darkorange", lw=lw) plt.fill_between(param_range, train_scores_mean - train_scores_std, train_scores_mean + train_scores_std, alpha=0.2, color="darkorange", lw=lw) plt.semilogx(param_range, test_scores_mean, label="Cross-validation score", color="navy", lw=lw) plt.fill_between(param_range, test_scores_mean - test_scores_std, test_scores_mean + test_scores_std, alpha=0.2, color="navy", lw=lw) plt.legend(loc="best") plt.show();
代碼中:
X, y = load_digits(return_X_y=True) # 等價于 digits = load_digits() X_digits = digits.data y_digits = digits.target
以下是支持向量機的驗證曲線,調節的超參數gamma共有5個值,每一個點的分數是五折交叉驗證(cv=5)的均值。
當想看模型多個超參數與模型評分之間的關系時,使用scikit-learn中validation curve就難以實現,因此可以考慮繪制三維坐標圖。
主要用plotly的庫繪制3D Scatter(3d散點圖)。以下的例子使用scikit-learn中的鶯尾花的數據集(iris)。以下例子選用隨機森林模型(RandomForestRegressor),利用scikit-learn中的GridSearchCV方法調試最佳超參(tuning hyper-parameters),分別設置超參數"n_estimators","max_features","min_samples_split"的參數范圍,詳見代碼如下:
import numpy as np from sklearn.model_selection import validation_curve from sklearn.datasets import load_iris from sklearn.ensemble import RandomForestRegressor from plotly.offline import iplot from plotly.graph_objs as go model = RandomForestRegressor(n_jobs=-1, random_state=2, verbose=2) grid = {'n_estimators': [10,110,200], 'max_features': [0.05, 0.07, 0.09, 0.11, 0.13], 'min_samples_split': [2, 3, 5, 8]} rf_gridsearch = GridSearchCV(estimator=model, param_grid=grid, n_jobs=4, cv=5, verbose=2, return_train_score=True) rf_gridsearch.fit(X, y) # and after some hours... df_gridsearch = pd.DataFrame(rf_gridsearch.cv_results_) trace = go.Scatter3d( x=df_gridsearch['param_max_features'], y=df_gridsearch['param_n_estimators'], z=df_gridsearch['param_min_samples_split'], mode='markers', marker=dict( # size=df_gridsearch.mean_fit_time ** (1 / 3), size = 10, color=df_gridsearch.mean_test_score, opacity=0.99, colorscale='Viridis', colorbar=dict(title = 'Test score'), line=dict(color='rgb(140, 140, 170)'), ), text=df_gridsearch.Text, hoverinfo='text' ) data = [trace] layout = go.Layout( title='3D visualization of the grid search results', margin=dict( l=30, r=30, b=30, t=30 ), scene = dict( xaxis = dict( title='max_features', nticks=10 ), yaxis = dict( title='n_estimators', ), zaxis = dict( title='min_samples_split', ), ), ) fig = go.Figure(data=data, layout=layout) iplot(fig)
其運行結果如果,是一個三維散點圖(3D Scatter)。
可以看到顏色越淺,分數越高。n_estimators(子估計器)越多,分數越高,max_features的變化對模型分數的影響較小,在圖中看不到變化,min_samples_split的個數并不是越高越好,但與模型分數并不呈單調關系,在min_samples_split取2時(此時,其它條件不變),模型分數最高。
除了使用scikit-learn中validation curve繪制超參數與得分的可視化,還可以利用seaborn庫中heatmap方法來實現兩個超參數之間的關系圖,如下代碼示例:
import seaborn as sns title = '''Maximum R2 score on test set VS max_features, min_samples_split''' sns.heatmap(max_scores.mean_test_score, annot=True, fmt='.4g'); plt.title(title); plt.savefig("heatmap_test.png", dpi = 300);
import seaborn as sns title = '''Maximum R2 score on train set VS max_features, min_samples_split''' sns.heatmap(max_scores.mean_train_score, annot=True, fmt='.4g'); plt.title(title); plt.savefig("heatmap_train.png", dpi = 300);
max_features和min_samples與模型得分關系的可視化如下圖所示(分別為網格搜索中測試集和訓練集的得分):
由于一般人很難迅速的在大量數據中找到隱藏的關系,因此,可以考慮繪圖,將數據關系以圖表的形式,清晰的顯現出來。
綜上,當關注單個超參數的學習曲線時,可以使用scikit-learn中validation curve,找到拐點,作為模型的最佳參數。
當關注兩個超參數的共同變化對模型分數的影響時,可以使用seaborn庫中的heatmap方法,制作“熱圖”,以找到超參數協同變化對分數影響的趨勢。
當關注三個參數的協同變化與模型得分的關系時,可以使用poltly庫中的iplot和go方法,繪制3d散點圖(3D Scatter),將其協同變化對模型分數的影響展現在高維圖中。
數據分析咨詢請掃描二維碼
若不方便掃碼,搜微信號:CDAshujufenxi
在本文中,我們將探討 AI 為何能夠加速數據分析、如何在每個步驟中實現數據分析自動化以及使用哪些工具。 數據分析中的AI是什么 ...
2025-05-20當數據遇見人生:我的第一個分析項目 記得三年前接手第一個數據分析項目時,我面對Excel里密密麻麻的銷售數據手足無措。那些跳動 ...
2025-05-20在數字化運營的時代,企業每天都在產生海量數據:用戶點擊行為、商品銷售記錄、廣告投放反饋…… 這些數據就像散落的拼圖,而相 ...
2025-05-19在當今數字化營銷時代,小紅書作為國內領先的社交電商平臺,其銷售數據蘊含著巨大的商業價值。通過對小紅書銷售數據的深入分析, ...
2025-05-16Excel作為最常用的數據分析工具,有沒有什么工具可以幫助我們快速地使用excel表格,只要輕松幾步甚至輸入幾項指令就能搞定呢? ...
2025-05-15數據,如同無形的燃料,驅動著現代社會的運轉。從全球互聯網用戶每天產生的2.5億TB數據,到制造業的傳感器、金融交易 ...
2025-05-15大數據是什么_數據分析師培訓 其實,現在的大數據指的并不僅僅是海量數據,更準確而言是對大數據分析的方法。傳統的數 ...
2025-05-14CDA持證人簡介: 萬木,CDA L1持證人,某電商中廠BI工程師 ,5年數據經驗1年BI內訓師,高級數據分析師,擁有豐富的行業經驗。 ...
2025-05-13CDA持證人簡介: 王明月 ,CDA 數據分析師二級持證人,2年數據產品工作經驗,管理學博士在讀。 學習入口:https://edu.cda.cn/g ...
2025-05-12CDA持證人簡介: 楊貞璽 ,CDA一級持證人,鄭州大學情報學碩士研究生,某上市公司數據分析師。 學習入口:https://edu.cda.cn/g ...
2025-05-09CDA持證人簡介 程靖 CDA會員大咖,暢銷書《小白學產品》作者,13年頂級互聯網公司產品經理相關經驗,曾在百度、美團、阿里等 ...
2025-05-07相信很多做數據分析的小伙伴,都接到過一些高階的數據分析需求,實現的過程需要用到一些數據獲取,數據清洗轉換,建模方法等,這 ...
2025-05-06以下的文章內容來源于劉靜老師的專欄,如果您想閱讀專欄《10大業務分析模型突破業務瓶頸》,點擊下方鏈接 https://edu.cda.cn/g ...
2025-04-30CDA持證人簡介: 邱立峰 CDA 數據分析師二級持證人,數字化轉型專家,數據治理專家,高級數據分析師,擁有豐富的行業經驗。 ...
2025-04-29CDA持證人簡介: 程靖 CDA會員大咖,暢銷書《小白學產品》作者,13年頂級互聯網公司產品經理相關經驗,曾在百度,美團,阿里等 ...
2025-04-28CDA持證人簡介: 居瑜 ,CDA一級持證人國企財務經理,13年財務管理運營經驗,在數據分析就業和實踐經驗方面有著豐富的積累和經 ...
2025-04-27數據分析在當今信息時代發揮著重要作用。單因素方差分析(One-Way ANOVA)是一種關鍵的統計方法,用于比較三個或更多獨立樣本組 ...
2025-04-25CDA持證人簡介: 居瑜 ,CDA一級持證人國企財務經理,13年財務管理運營經驗,在數據分析就業和實踐經驗方面有著豐富的積累和經 ...
2025-04-25在當今數字化時代,數據分析師的重要性與日俱增。但許多人在踏上這條職業道路時,往往充滿疑惑: 如何成為一名數據分析師?成為 ...
2025-04-24以下的文章內容來源于劉靜老師的專欄,如果您想閱讀專欄《劉靜:10大業務分析模型突破業務瓶頸》,點擊下方鏈接 https://edu.cda ...
2025-04-23