先交代一句《RF-RFE-BP基于隨機(jī)森林遞歸特征消除與BP神經(jīng)網(wǎng)絡(luò)回歸預(yù)測》這事我已經(jīng)在回歸預(yù)測的坑里反復(fù)折騰過很多輪了。早些年拿原始維度直接懟進(jìn)BP網(wǎng)絡(luò)結(jié)果不是過擬合就是訓(xùn)練時間感人后來學(xué)了聰明辦法先用隨機(jī)森林做個遞歸特征消除把冗余變量砍掉再送給BP做回歸。今天把這套流程里我實測過的東西全部寫透從特征選擇原理到MATLAB代碼逐步照抄級別的細(xì)節(jié)以及我踩過的坑和調(diào)參心法都放在這里。1. 為什么偏偏是RF-RFE加BP這套組合解決的問題與適用人群1.1 先聊痛點維度一高BP網(wǎng)絡(luò)就翻車做回歸預(yù)測的人尤其是拿MATLAB吭哧吭哧寫代碼的那種大概率遇到過這樣的局面手里的數(shù)據(jù)表格一打開列數(shù)比行數(shù)還多或者幾十上百個特征里真正有用的就那幾個剩下的全是噪聲。此時如果直接敲fitnet建BP神經(jīng)網(wǎng)絡(luò)會立刻撞上兩個老熟人——過擬合和維數(shù)災(zāi)難。BP網(wǎng)絡(luò)的本質(zhì)是一個萬能函數(shù)逼近器輸入維度越高需要擬合的參數(shù)空間就越大而你的樣本量往往撐不起這個復(fù)雜度。我拿一組80維特征、300個樣本的工業(yè)數(shù)據(jù)試過直接BP訓(xùn)練訓(xùn)練集R2能到0.98驗證集直接掉到0.6以下典型的泛化崩盤。另一個麻煩是神經(jīng)網(wǎng)絡(luò)對無關(guān)特征極其敏感一堆純噪聲輸入會把隱藏層的權(quán)重攪得亂七八糟你根本分不清模型到底學(xué)到了規(guī)律還是記住了噪聲。所以問題的核心不是“把BP調(diào)多好”而是“送什么進(jìn)BP”。特征選擇這一步價值比調(diào)參大多了。1.2 RF-RFE的原理優(yōu)勢為什么不是單純的隨機(jī)森林重要性排序隨機(jī)森林Random Forest本身就能輸出特征重要性分?jǐn)?shù)很多人拿到重要性排名后直接取Top K說實話這種做法能用但不嚴(yán)謹(jǐn)。因為隨機(jī)森林的特征重要性是有偏的——它偏向于那些取值水平多、數(shù)值范圍大的連續(xù)特征并且在特征之間存在強相關(guān)時重要性會被分散到相關(guān)特征組里導(dǎo)致排名失真。遞歸特征消除Recursive Feature Elimination, RFE解決的就是這個問題。它在每輪迭代里訓(xùn)練一次模型根據(jù)特征重要性或模型系數(shù)消除最不重要或貢獻(xiàn)最小的一個或一批特征然后在剩余特征上繼續(xù)訓(xùn)練、繼續(xù)消除直到達(dá)到目標(biāo)特征數(shù)量或模型性能開始下降。這個“遞歸”過程的關(guān)鍵在于它考慮的是每個特征在當(dāng)前特征子集下的邊際貢獻(xiàn)而不是一次性用全局重要性來拍板。把隨機(jī)森林和RFE結(jié)合就是讓隨機(jī)森林這個對非線性關(guān)系捕捉能力強、自帶OOB評估的模型來扮演“重要性裁判”每輪用它的特征重要性結(jié)果做排序再逐步刪減。這種組合在學(xué)術(shù)圈和工業(yè)軟測量里都被大量驗證效果穩(wěn)定的原因在于隨機(jī)森林對異常值和噪聲相對魯棒而且每輪迭代用OOB誤差或交叉驗證誤差來監(jiān)測性能變化可以科學(xué)地找到“特征數(shù)量-預(yù)測精度”的最優(yōu)點。1.3 這套方案適合誰和什么場景如果手上是這幾類項目我覺得這套代碼你基本可以無腦嘗試高光譜遙感數(shù)據(jù)建模波段數(shù)量幾百個樣本量又稀有必須做波段選擇后再進(jìn)模型工業(yè)過程軟測量傳感器變量幾十上百個目標(biāo)變量如產(chǎn)品質(zhì)量指標(biāo)難以在線測量生物醫(yī)學(xué)、組學(xué)數(shù)據(jù)基因/蛋白/代謝物特征動輒上千有效特征極少氣象、電力、負(fù)荷預(yù)測候選特征包括不同時段、不同站點的觀測值相關(guān)性和冗余度很高任何想用BP網(wǎng)絡(luò)但擔(dān)心輸入太雜導(dǎo)致過擬合的場景注意這套方案解決的是“特征多、樣本中等”的問題不是“樣本太少”的銀彈。如果樣本量只有三五十個再怎么特征選擇BP也救不回來這種情況先去考慮小樣本方法或數(shù)據(jù)增廣。2. RF-RFE特征選擇原理解密與MATLAB逐步實操2.1 核心流程到底長什么樣理清概念后我們要把RF-RFE落在代碼上。在MATLAB里沒有現(xiàn)成的rf_rfe函數(shù)你需要組合幾個官方工具箱函數(shù)來實現(xiàn)核心包括TreeBagger或fitensemble裝袋回歸樹、predictorImportance特征重要性提取以及自己寫的外層循環(huán)或直接用sequentialfs。我按照原理把標(biāo)準(zhǔn)流程拆成六個步驟輸入標(biāo)準(zhǔn)化后的訓(xùn)練數(shù)據(jù)X和標(biāo)簽y定義特征全集和候選集用當(dāng)前特征子集訓(xùn)練隨機(jī)森林回歸模型建議用袋裝回歸樹200~500棵樹提取特征重要性向量OOBPermutedPredictorDeltaError對應(yīng)的就是predictorImportance的輸出按重要性從小到大排序消除最不重要當(dāng)前貢獻(xiàn)最低的一個特征用交叉驗證或OOB誤差評估當(dāng)前特征子集對應(yīng)的模型性能記錄下特征數(shù)-誤差軌跡循環(huán)執(zhí)行直到特征全集為空對比各輪誤差曲線選取誤差最低或兼顧簡潔度的特征子集這里有非常關(guān)鍵的一點解釋為什么要記錄每一輪的性能而不是直接干到只剩一個特征因為特征數(shù)量和模型性能不是簡單的單調(diào)關(guān)系——刪除噪聲特征會提升性能但刪到某個臨界點后繼續(xù)刪就會丟掉有效信息性能會掉頭向下。所以整個RFE過程的“誤差-特征數(shù)量”曲線一般是一個U型曲線最低點對應(yīng)的特征數(shù)量就是我們要找的最佳子集大小。2.2 MATLAB代碼實現(xiàn)手寫RFE循環(huán)版先給大家看我最常用的一套手寫版代碼邏輯這對理解RFE的本質(zhì)更有幫助。function [selectedIdx, featureHistory] rf_rfe_hand(X, y, minFeatures, nTrees, kFold) % X: nSample x nFeature 的訓(xùn)練輸入 % y: nSample x 1 的訓(xùn)練目標(biāo) % minFeatures: 最少保留的特征數(shù) % nTrees: 隨機(jī)森林中樹的棵數(shù) % kFold: 交叉驗證折數(shù) allIdx 1:size(X, 2); currentIdx allIdx; featureHistory []; rng(42); % 固定隨機(jī)種子保證可復(fù)現(xiàn) while length(currentIdx) minFeatures % 1. 在當(dāng)前特征子集上訓(xùn)練隨機(jī)森林 rf TreeBagger(nTrees, X(:, currentIdx), y, ... Method, regression, ... OOBPrediction, on, ... OOBPredictorImportance, on, ... MinLeafSize, 5); % 2. 提取特征重要性 imp rf.OOBPermutedPredictorDeltaError; % 每個特征的重要性分?jǐn)?shù) % 3. 找到重要性最小的特征從當(dāng)前索引集中剔除 [~, minPos] min(imp); removedIdx currentIdx(minPos); currentIdx(minPos) []; % 4. 用交叉驗證評估剔除后的特征子集 cvMse crossValMse(X(:, currentIdx), y, kFold, nTrees); featureHistory(end1, :) [length(currentIdx), cvMse]; %#okAGROW fprintf(剩余特征數(shù)%dCV-MSE%.4f剔除特征%d\n, ... length(currentIdx), cvMse, removedIdx); end % 5. 根據(jù)歷史記錄選擇CV誤差最小的特征子集 [~, bestPos] min(featureHistory(:, 2)); bestNum featureHistory(bestPos, 1); % 從歷史記錄反推最佳特征集合 % 注意這里需要在循環(huán)中記錄每個階段的具體特征索引才能準(zhǔn)確回退 selectedIdx recoverIdxAtFeatureCount(featureHistory, bestNum); end特別注意最后一步“反推最佳特征集合”——因為循環(huán)過程中特征是一個一個被剔除的你不僅要記錄誤差歷史還要記錄每次剔除的是哪個特征。我實際的代碼通常會把currentIdx的整個快照都存進(jìn)一個cell數(shù)組上例中我用recoverIdxAtFeatureCount這個函數(shù)示意每輪結(jié)束后保存{length(currentIdx), currentIdx, cvMse}三元組最后從歷史里挑 CV誤差最小時的索引快照。crossValMse是我封裝的一個小函數(shù)里面用cvpartition做K折劃分每折用TreeBagger訓(xùn)練計算驗證集MSEfunction mse crossValMse(Xsub, y, kFold, nTrees) rng(1); cvp cvpartition(length(y), KFold, kFold); foldMse zeros(cvp.NumTestSets, 1); for i 1:cvp.NumTestSets trIdx cvp.training(i); teIdx cvp.test(i); mdl TreeBagger(nTrees, Xsub(trIdx, :), y(trIdx), ... Method, regression, MinLeafSize, 5); yhat predict(mdl, Xsub(teIdx, :)); foldMse(i) mean((yhat - y(teIdx)).^2); end mse mean(foldMse); end這套手寫方案的好處是每一步都透明可控你能看到每輪剔除哪個特征、誤差怎么變化。適合學(xué)習(xí)和調(diào)試。2.3 快速版實現(xiàn)用 sequentialfs 優(yōu)雅搞定RFE如果想少寫點代碼MATLAB的sequentialfs函數(shù)可以承擔(dān)外層RFE循環(huán)配合函數(shù)句柄來實現(xiàn)隨機(jī)森林的評估。整體代碼短很多適合已經(jīng)理解原理、想快速出結(jié)果的場景。rng(42); % 定義RFE評估函數(shù)訓(xùn)練隨機(jī)森林并返回MSE fun (XTrain, yTrain, XTest, yTest) ... mse(predict(TreeBagger(200, XTrain, yTrain, Method, regression), XTest) - yTest); opts statset(Display, iter, MaxIter, 200); % 使用Holdout交叉驗證按比例留出驗證集 cv cvpartition(length(y), Holdout, 0.3); % backward方式從全特征開始逐個消除 [fs, history] sequentialfs(fun, X, y, ... cv, cv, direction, backward, options, opts); % fs是邏輯向量true表示保留的特征 selectedIdx find(fs); fprintf(RFE最終選擇特征%d 個\n, length(selectedIdx));sequentialfs內(nèi)部做的事和我上面手寫循環(huán)基本一致它默認(rèn)用1個特征淘汰粒度也可以設(shè)置nfeatures控制最小保留數(shù)量。注意direction, backward指從全集開始刪除對應(yīng)RFE語義如果設(shè)成forward則是逐個添加特征這種是SFS速度雖然更快但效果通常不如RFE。我的個人建議是正式跑項目用sequentialfs快速確認(rèn)最優(yōu)特征數(shù)量范圍然后再用手寫版跑一輪詳細(xì)的迭代軌跡方便畫誤差曲線圖發(fā)論文或做匯報。2.4 特征數(shù)-誤差曲線怎么看才科學(xué)跑完RFE之后你會得到一條“特征數(shù)-CV誤差”的曲線。有兩點判斷經(jīng)驗供參考第一U型曲線的谷底并不總是唯一的最優(yōu)解。如果谷底特征數(shù)量是30但在特征數(shù)20~35這個區(qū)間內(nèi)的誤差差異很小比如CV-MSE變化小于3%~5%我通常會選更小的特征數(shù)那一端。道理很簡單特征數(shù)量減半模型復(fù)雜度降低泛化能力更強而精度損失可以忽略。這不叫妥協(xié)這叫工程理性。第二RFE的結(jié)果受隨機(jī)種子影響。隨機(jī)森林本身有抽樣隨機(jī)性每輪排序可能波動導(dǎo)致最終選出的特征集合不完全一致。多跑幾次RFE比如跑5次固定不同種子看哪些特征被反復(fù)選中。那些在多次運行中穩(wěn)定登場的特征才是真正信號偶爾閃現(xiàn)的多半是運氣好。這個叫“穩(wěn)定性分析”在提供報告時可以加分不少。3. BP神經(jīng)網(wǎng)絡(luò)回歸模型結(jié)構(gòu)選擇、參數(shù)設(shè)定與MATLAB實現(xiàn)3.1 特征選好了BP的結(jié)構(gòu)怎么定RFE選出的特征子集接下來就要喂給我們親愛的BP神經(jīng)網(wǎng)絡(luò)。在MATLAB里構(gòu)建BP網(wǎng)的主要方式有feedforwardnet前饋網(wǎng)絡(luò)、fitnet函數(shù)擬合網(wǎng)絡(luò)本質(zhì)也是前饋BP、newff老版本接口。強烈建議用fitnet它內(nèi)部實現(xiàn)更規(guī)范自帶數(shù)據(jù)劃分和訓(xùn)練配置。網(wǎng)絡(luò)結(jié)構(gòu)上需要確定的參數(shù)有三個輸入層節(jié)點數(shù)由RFE選出的特征數(shù)量決定不需要你糾結(jié)隱藏層節(jié)點數(shù)最常被問的參數(shù)。經(jīng)驗公式有sqrt(nIn nOut) 1~10、2*nIn1、(nInnOut)/2等但我的實測經(jīng)驗是在樣本量幾百到幾千的回歸任務(wù)里隱藏層節(jié)點數(shù)取輸入特征數(shù)的0.5~1.5倍通常夠用。比如篩選后特征數(shù)是15隱藏層取10~20都可以具體用交叉驗證微調(diào)輸出層節(jié)點數(shù)單輸出回歸就是1個節(jié)點無需激活函數(shù)偏置多輸出就是對應(yīng)維度還有隱藏層層數(shù)。絕大多數(shù)回歸任務(wù)單隱藏層就夠因為單隱層BP已經(jīng)能逼近任意連續(xù)函數(shù)這是BP的萬能逼近定理。只有數(shù)據(jù)極其非線性、單隱層怎么調(diào)都不收斂時才考慮雙隱層。我的建議是一次建模先從單隱層開始不要一上來就堆深度深度越大越容易過擬合調(diào)試難度也指數(shù)上升。3.2 訓(xùn)練參數(shù)設(shè)置的邏輯和實戰(zhàn)取值fitnet創(chuàng)建網(wǎng)絡(luò)后需要設(shè)置的幾個關(guān)鍵參數(shù)hiddenLayerSize 12; % 隱藏層節(jié)點數(shù)結(jié)合RFE后的特征數(shù)調(diào)整 net fitnet(hiddenLayerSize, trainlm); % trainlm是Levenberg-Marquardt優(yōu)化器 % 關(guān)鍵訓(xùn)練參數(shù) net.trainParam.epochs 1000; % 最大迭代次數(shù) net.trainParam.goal 1e-6; % 均方誤差目標(biāo) net.trainParam.min_grad 1e-7; % 最小梯度閾值 net.trainParam.max_fail 20; % 驗證集連續(xù)不改善的最大次數(shù)早停用 net.trainParam.lr 0.01; % 學(xué)習(xí)率(針對梯度下降類; trainlm下影響較小) net.trainParam.mc 0.9; % 動量因子這里重點解釋為什么推薦trainlm。Levenberg-Marquardt算法是高斯牛頓法和梯度下降法的結(jié)合在小規(guī)模網(wǎng)絡(luò)中收斂速度極快精度也高。注意它需要計算雅可比矩陣內(nèi)存開銷與參數(shù)數(shù)量平方相關(guān)因此只適合幾百個權(quán)重的輕量網(wǎng)絡(luò)。如果隱藏層節(jié)點數(shù)超過50且特征很多導(dǎo)致參數(shù)爆炸就要考慮換成trainscg縮放共軛梯度或trainbr貝葉斯正則化后兩者內(nèi)存壓力小其中trainbr自帶正則化對過擬合有天生抑制能力。數(shù)據(jù)劃分也是個常被忽略但影響巨大的設(shè)置net.divideFcn divideblock; % 按順序分塊劃分適合時序數(shù)據(jù) net.divideParam.trainRatio 0.7; net.divideParam.valRatio 0.15; net.divideParam.testRatio 0.15;如果是普通截面數(shù)據(jù)用默認(rèn)的dividerand隨機(jī)劃分沒問題但涉及時間序列比如負(fù)荷預(yù)測、氣象預(yù)測絕對不要隨機(jī)劃分必須用divideblock否則未來信息會泄漏進(jìn)訓(xùn)練集測試效果虛高上了實際應(yīng)用瞬間翻車。這個問題我見過太多人栽過。3.3 輸入輸出歸一化99%的人都會踩的坑BP神經(jīng)網(wǎng)絡(luò)的激活函數(shù)tansig、logsig對輸入范圍敏感輸入數(shù)據(jù)范圍差距過大的時候比如一個特征范圍是0.01~0.99另一個是500~5000訓(xùn)練會非常不穩(wěn)定梯度要么爆炸要么消失。所以訓(xùn)練前必須歸一化。MATLAB里最方便的方式可能是用mapminmax但我更推薦自己手工做Z-score標(biāo)準(zhǔn)化% 用訓(xùn)練集的均值和標(biāo)準(zhǔn)差做標(biāo)準(zhǔn)化注意測試集也要用訓(xùn)練集的統(tǒng)計量 mu_X mean(X_train); sigma_X std(X_train); X_train_std (X_train - mu_X) ./ sigma_X; X_test_std (X_test - mu_X) ./ sigma_X;這句“測試集也要用訓(xùn)練集的統(tǒng)計量”是無數(shù)教程里沒寫清楚的關(guān)鍵。如果你拿測試集自己的均值標(biāo)準(zhǔn)差去做相當(dāng)于讓模型偷偷看到了測試集的分布信息這叫數(shù)據(jù)泄漏會導(dǎo)致驗證指標(biāo)虛高。正確的做法就是上面代碼那樣訓(xùn)練集算出來的mu_X, sigma_X保存變量留在工作區(qū)就行測試時直接用。標(biāo)簽y同樣需要標(biāo)準(zhǔn)化mu_y mean(y_train); sigma_y std(y_train); y_train_std (y_train - mu_y) ./ sigma_y; % 預(yù)測得到y(tǒng)_pred_std后反標(biāo)準(zhǔn)化 y_pred y_pred_std * sigma_y mu_y;更省事的是用mapminmax但那個函數(shù)有個坑它默認(rèn)針對行向量處理需要轉(zhuǎn)置來轉(zhuǎn)回去。我前幾年在這上面浪費過不少時間后來干脆全部手寫標(biāo)準(zhǔn)化代碼清晰也不容易出錯。3.4 完整BP訓(xùn)練代碼可直接套用function [net, info, y_pred_denorm, y_test_denorm] train_bp_regression(X_train, y_train, X_test, y_test, hiddenSize) % 1. 標(biāo)準(zhǔn)化 [Xtr, mu_x, sigma_x] zscore(X_train); ytr_mean mean(y_train); ytr_std std(y_train); ytr (y_train - ytr_mean) / ytr_std; Xte (X_test - mu_x) ./ sigma_x; yte_mean mean(y_test); yte_std std(y_test); yte (y_test - yte_mean) / yte_std; % 2. 構(gòu)建網(wǎng)絡(luò) net fitnet(hiddenSize, trainlm); net.divideFcn divideblock; net.divideParam.trainRatio 0.7; net.divideParam.valRatio 0.15; net.divideParam.testRatio 0.15; net.trainParam.epochs 1000; net.trainParam.goal 1e-6; net.trainParam.min_grad 1e-7; net.trainParam.max_fail 20; % 3. 訓(xùn)練 rng(42); % 固定種子 [net, ~] train(net, Xtr, ytr); % 4. 測試集預(yù)測與反標(biāo)準(zhǔn)化 y_pred_std net(Xte); y_pred_denorm y_pred_std * ytr_std ytr_mean; y_test_denorm y_test; % 原始測試標(biāo)簽 end注意MATLAB神經(jīng)網(wǎng)絡(luò)工具箱默認(rèn)要求輸入數(shù)據(jù)是樣本數(shù)×特征數(shù)的矩陣但train函數(shù)內(nèi)部經(jīng)常要用轉(zhuǎn)置形式Xtr特征數(shù)×樣本數(shù)這個轉(zhuǎn)置操作容易搞暈人。我的經(jīng)驗是統(tǒng)一約定在進(jìn)入train之前就把數(shù)據(jù)轉(zhuǎn)置好出來后預(yù)測也注意轉(zhuǎn)置多用size()檢查維度別偷懶。4. RF-RFE-BP全流程串聯(lián)數(shù)據(jù)準(zhǔn)備、聯(lián)調(diào)與效果評測4.1 完整的項目級流程編排從原始數(shù)據(jù)到預(yù)測結(jié)果整個流程串聯(lián)起來其實是固定的幾個環(huán)節(jié)我這里直接給出一套我自己項目上反復(fù)使用的編排邏輯數(shù)據(jù)清洗處理缺失值均值/中位數(shù)填充或刪除樣本、異常值3σ原則或分位數(shù)截斷、重復(fù)樣本合并樣本劃分按時間順序時序數(shù)據(jù)或隨機(jī)分層截面數(shù)據(jù)劃分訓(xùn)練集和測試集。測試集原則上只在這個流程的最后出現(xiàn)一次中間的特征選擇、模型調(diào)參都不能碰它RFE特征選擇只用訓(xùn)練集做RF-RFE確定最優(yōu)特征數(shù)量輸出特征索引子集提取把訓(xùn)練集和測試集都按選出的特征索引取出對應(yīng)列BP模型訓(xùn)練在篩選后的訓(xùn)練集上訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)內(nèi)部再分訓(xùn)練/驗證/測試子集模型評價在保留的測試集上做最終預(yù)測計算R2、RMSE、MAE、MAPE等指標(biāo)畫回歸散點圖和誤差圖這里要特別強調(diào)特征選擇只能利用訓(xùn)練集很多新手把整個數(shù)據(jù)集拿去跑RFE選好特征之后才劃分訓(xùn)練測試這樣做其實是在特征選擇環(huán)節(jié)就偷看了測試集的信息結(jié)果就是測試指標(biāo)好看得不得了一到新數(shù)據(jù)就崩。科學(xué)做法是劃分訓(xùn)練/測試集第一步就完成之后所有特征選擇操作都限定在訓(xùn)練集內(nèi)。4.2 核心評測指標(biāo)的計算代碼評價回歸模型只用MSE或R2其實不夠全面。我一般至少算4個指標(biāo)function [R2, RMSE, MAE, MAPE] regression_metrics(y_true, y_pred) % R2決定系數(shù) SS_res sum((y_true - y_pred).^2); SS_tot sum((y_true - mean(y_true)).^2); R2 1 - SS_res / SS_tot; % RMSE均方根誤差 RMSE sqrt(mean((y_true - y_pred).^2)); % MAE平均絕對誤差 MAE mean(abs(y_true - y_pred)); % MAPE平均絕對百分比誤差注意y中有0時要處理 nonzeroIdx y_true ~ 0; MAPE mean(abs((y_true(nonzeroIdx) - y_pred(nonzeroIdx)) ./ y_true(nonzeroIdx))) * 100; endMAPE 在標(biāo)簽存在0值或接近0的值時要小心分母趨近于0會導(dǎo)致指標(biāo)爆炸這種情況建議改用 sMAPE 或直接省略。我遇到過負(fù)荷預(yù)測的測試集里剛好有幾個0點結(jié)果MAPE算出個幾百的怪值差點以為是模型崩了其實是分母問題。4.3 一個完整示例從80維特征到精準(zhǔn)預(yù)測為了讓這個流程更具體我用之前做過的風(fēng)速預(yù)測項目來演示效果。原始特征80個歷史風(fēng)速、溫度、氣壓、濕度、風(fēng)向正弦余弦、時間特征等一堆候選變量樣本量500其中前350個做訓(xùn)練集后150個做測試集時序數(shù)據(jù)按時間切分。RFE跑完后得到的特征數(shù)量-誤差曲線顯示特征數(shù)量從80削減到20左右時CV-MSE開始顯著下降降到10~15時處于谷底平臺低于6后誤差急劇上升。最終我選了9個特征的穩(wěn)定子集占了原始維度的約11%但保留了解釋風(fēng)速變化的全部核心變量。然后把這些特征送進(jìn)BP隱藏層節(jié)點數(shù)設(shè)為10trainlm訓(xùn)練迭代大概150步就收斂了。測試集結(jié)果R20.913RMSE0.62m/sMAE0.47m/s。而對比直接拿80維特征跑BP測試集R2只有0.74RMSE高達(dá)1.05m/s對比只用RF重要性Top10特征跑BPR20.86RMSE0.81m/s。RF-RFE比純RF篩選多出的這0.05個R2提升基本就是“遞歸考慮特征組合效應(yīng)”換來的優(yōu)勢。上面這張表是我個人經(jīng)驗里統(tǒng)計出來的典型對比不同數(shù)據(jù)表現(xiàn)會有差異但規(guī)律是一致的直接BP最差純重要性排序BP中等RF-RFE-BP最好。方案測試集R2RMSE特征數(shù)原始80維直接BP0.741.0580RF重要性Top10 BP0.860.8110RF-RFE(9特征) BP0.910.6294.4 歸一化和標(biāo)準(zhǔn)化如何與特征選擇聯(lián)動在RFE階段隨機(jī)森林不需要歸一化樹模型對特征尺度不敏感所以在做RFE時直接拿原始X就行。但進(jìn)入BP階段前必須歸一化。這就帶來一個銜接細(xì)節(jié)你需要在RFE確定特征索引后用選出的特征列重新做標(biāo)準(zhǔn)化。順序不能反——如果先標(biāo)準(zhǔn)化再做RFE其實也不會有問題因為特征選擇只依賴排序不依賴scale但最規(guī)范的做法還是全流程在特征選擇完成后再統(tǒng)一做標(biāo)準(zhǔn)化避免無謂的計算誤差。我個人流程會保存一個selected_idx向量和一個mu_x/sigma_x結(jié)構(gòu)體這樣將來新數(shù)據(jù)進(jìn)入模型時只需要取對應(yīng)特征列、套同樣的標(biāo)準(zhǔn)化參數(shù)再進(jìn)BP預(yù)測即可整體封裝得像一個完整的預(yù)測管線。5. 高頻踩坑記錄與調(diào)參經(jīng)驗速查表價值超高的部分5.1 特征選擇結(jié)果不穩(wěn)定的排查癥狀同一份數(shù)據(jù)不加rng固定種子跑兩次RF-RFE選出的特征集合差好多。原因很直接隨機(jī)森林本身基于bootstrap抽樣和隨機(jī)特征子集構(gòu)建重要性分?jǐn)?shù)帶有隨機(jī)波動當(dāng)多個特征相關(guān)性高時重要性會在它們之間隨機(jī)挪移導(dǎo)致排序不穩(wěn)定。解決方案按優(yōu)先級排列增加樹的數(shù)量200顆不行就500甚至1000樹多則重要性估計方差降低固定隨機(jī)種子rng(42)讓實驗可復(fù)現(xiàn)這是學(xué)術(shù)報告的基本要求多次運行取交集跑5次RFE保留至少出現(xiàn)3次的特征這個“穩(wěn)定特征集合”更可靠檢查多重共線性如果業(yè)務(wù)上允許先做相關(guān)矩陣剔除高度相關(guān)的特征對|r|0.8可以顯著提升RFE穩(wěn)定性5.2 BP網(wǎng)絡(luò)訓(xùn)練不收斂或陷入局部極小癥狀1訓(xùn)練MSE一直在高位怎么都壓不下來。癥狀2loss曲線劇烈震蕩不往下降癥狀3訓(xùn)練集收斂但驗證集糟糕過擬合。排查順序如下檢查數(shù)據(jù)標(biāo)準(zhǔn)化是否所有輸入和輸出都在合理范圍均值為0方差為1少了這步神仙難救減少隱藏層節(jié)點數(shù)節(jié)點太多會讓網(wǎng)絡(luò)過于自由容易繞進(jìn)過擬合更換訓(xùn)練算法trainlm不收斂就換trainscg或者從較小的學(xué)習(xí)率開始逐步加大降低學(xué)習(xí)率如果震蕩把lr從0.01降到0.001甚至更低配合動量因子0.8~0.9檢查標(biāo)簽是否異常y里有沒有極端離群點一個大離群點就能把誤差函數(shù)拉出天際建議先做Winsorize分位數(shù)截斷或剔除我這里還要單獨提醒BP訓(xùn)練前必須處理數(shù)據(jù)的排列順序。如果訓(xùn)練集的前70%全是低值、后30%全是高值尤其在截面數(shù)據(jù)隨機(jī)劃分不當(dāng)時可能出現(xiàn)trainlm會在局部數(shù)據(jù)范圍內(nèi)過度擬合泛化能力極差。用divideblock時也要確保訓(xùn)練塊內(nèi)部樣本足夠多樣。最穩(wěn)妥的做法是在劃分前對樣本做隨機(jī)打亂時序數(shù)據(jù)除外然后用dividerand劃分。5.3 驗證集誤差一直高過訓(xùn)練集一倍以上這個現(xiàn)象十有八九是數(shù)據(jù)泄漏的反面——特征選擇選得太激進(jìn)選出了一組在訓(xùn)練集上表現(xiàn)好但結(jié)構(gòu)脆弱的特征組合。還有一種常見原因測試集和訓(xùn)練集分布不一致比如時序數(shù)據(jù)前半年和后半年環(huán)境完全變了。解決辦法退回RFE曲線選擇略大特征數(shù)量那端的方案犧牲一點點訓(xùn)練精度換泛化穩(wěn)健用分層抽樣保證訓(xùn)練測試分布一致OLS那種寫法cvpartition指定Stratify不過這要求y是分類變量回歸中常用的是對特征做分箱后分層增加訓(xùn)練集比例到80%讓模型見到更多模式5.4 超參數(shù)選擇的經(jīng)驗表給出我個人常用的網(wǎng)格搜索范圍直接套用可以省不少事參數(shù)推薦范圍我的首選RF樹數(shù)量100~500200RF最小葉節(jié)點數(shù)3~105RF最大特征數(shù)自動(默認(rèn)1/3回歸)默認(rèn)BP隱藏層節(jié)點數(shù)4~20特征數(shù)的0.8~1.2倍BP訓(xùn)練算法trainlm/trainscgtrainlm(小數(shù)據(jù))學(xué)習(xí)率0.001~0.10.01動量因子0.8~0.950.9訓(xùn)練/驗證/測試比例7:1.5:1.57:1.5:1.5這些參數(shù)不是拍腦袋定的背后原理我在前面各節(jié)都解釋過了——樹數(shù)量影響穩(wěn)定性、葉節(jié)點影響平滑度、隱藏層節(jié)點影響擬合容量、學(xué)習(xí)率和動量影響優(yōu)化軌跡。調(diào)參的本質(zhì)是理解每個旋鈕控制的偏差-方差權(quán)衡位置。5.5 畫圖輸出的細(xì)節(jié)拉滿論文質(zhì)感的技巧MATLAB里出圖時有幾個讓人眼前一亮的細(xì)節(jié)% 訓(xùn)練集預(yù)測對比 figure(Color, w, Position, [100 100 900 400]); plot(y_train_denorm, b-, LineWidth, 1.2); hold on; plot(y_train_pred, r--, LineWidth, 1.2); xlabel(樣本序號); ylabel(目標(biāo)值); legend({真實值, 訓(xùn)練集預(yù)測}, Location, best); grid on; set(gca, FontSize, 12); title(訓(xùn)練集預(yù)測效果對比);測試集預(yù)測對比圖、散點Q-Q圖預(yù)測值vs真實值散點加對角線、誤差分布直方圖這三張圖一張都別省。散點圖上那條45度對角線是判斷模型是否系統(tǒng)性偏差的利器如果散點都落在線下方說明預(yù)測偏低都在線上方偏高越貼合對角線模型越好。6. 進(jìn)階擴(kuò)展這套組合的更多玩法與體會分享先說一個我最近在嘗試的方向把RF-RFE選出的特征不只是喂給BP還同時喂給多個模型支持向量回歸、極端梯度提升、高斯過程回歸做對比。為什么值得做因為RF-RFE是一個與最終回歸器無關(guān)的包裹式特征選擇器它的特征集對任何非線性模型都適用。你可以用同一組特征訓(xùn)練多個回歸器再做加權(quán)集成。這比用不同特征集分別訓(xùn)練再融合要簡潔得多也比單模型更穩(wěn)健。當(dāng)然這屬于進(jìn)階玩法先把RF-RFE-BP基礎(chǔ)鏈路跑通再考慮。另一個很有價值的擴(kuò)展是找出特征選擇的物理解釋。RFE給出了特征重要性的遞進(jìn)排名這些信息在撰寫實驗報告、論文或向業(yè)務(wù)方匯報時非常有用。舉個例子在軟測量項目中我發(fā)現(xiàn)某個與目標(biāo)變量間接相關(guān)的滯后變量被RFE保留了而不是那個直覺上應(yīng)該更重要的直接變量——這往往意味著數(shù)據(jù)中存在動態(tài)特性或隱藏的時滯效應(yīng)。向工藝工程師解釋時這張?zhí)卣鞅Ａ襞c否的表比任何R2指標(biāo)都能打動他們。還有一個我經(jīng)常用的變體把RFE的評價指標(biāo)從MSE換成MAE或者在RF-RFE選特征時把樣本加權(quán)策略加入進(jìn)去對異常值降權(quán)這樣得到的特征子集對極值預(yù)測更友好。如果你的目標(biāo)變量存在長尾分布這個變體值得一試。關(guān)于運行效率RF-RFE在特征維度較高時比較耗時。一次迭代訓(xùn)練200棵樹的隨機(jī)森林80維特征跑80輪單核大約需要幾分鐘到十幾分鐘。我的經(jīng)驗是先粗跑一輪樹100每輪刪2個特征快速確定特征數(shù)量的大致范圍再在候選區(qū)間內(nèi)精跑樹300每輪刪1個。這套“粗篩精選”的策略能把總體耗時壓縮到三分之一結(jié)果基本一致。最后想提醒一個我在項目交付時反復(fù)強調(diào)的原則所謂“最優(yōu)特征數(shù)量”并不是煉金術(shù)般的定數(shù)每個數(shù)據(jù)集都有它自己的最優(yōu)區(qū)間。RFE曲線只是一個工具、一個參考真的要做上線決策時別只盯著曲線谷底還要綜合考慮采集成本某個傳感器特征需要額外設(shè)備成本的話只要誤差增加在5%以內(nèi)我會果斷刪掉它、計算開銷和可解釋性。把工程約束放進(jìn)模型設(shè)計里你做的就不再是一個漂亮的實驗室玩具而是一個真正可落地的預(yù)測系統(tǒng)。這套RF-RFE-BP流程我前前后后在風(fēng)電功率預(yù)測、水質(zhì)軟測量、設(shè)備健康管理等項目上反復(fù)用過每一次測試集誤差都明顯優(yōu)于直接建模。它的優(yōu)雅之處在于兩個算法都很成熟不需要黑科技關(guān)鍵價值全在流程編排和細(xì)節(jié)控制上——早一點劃分?jǐn)?shù)據(jù)、認(rèn)真做穩(wěn)定性分析、堅持歸一化放最后、測試集絕不參與選特征。把這些細(xì)節(jié)做到位這套方法就是你回歸預(yù)測工具箱里最穩(wěn)的武器。