特征工程

第27课:交叉验证

K折交叉验证、分层K折、留一法、时间序列CV、嵌套CV、Group CV

1. K折交叉验证

K=5或10是折中。K小→高偏差;K大(n=LOO)→高方差。

2. 分层K折

每个折中类别比例与整体一致→不平衡数据必须用。

3. 留一交叉验证(LOO)

K=n→几乎无偏但方差高、计算量大。适用:样本极小+训练快。

4. 时间序列CV

前向链式:Train[1]→Test[2], Train[1,2]→Test[3], ... 不能随机划分!

5. 嵌套CV

外层估计泛化性能,内层选择超参数→无偏估计。非嵌套→过于乐观。

6. GroupKFold

同一组(同用户/同患者)不同时出现在训练和测试→防止泄露

Python代码实现与验证

💡 代码说明:以下代码使用 Python / Scikit-learn 实现,已实机运行验证。

完整代码


import numpy as np
from sklearn.datasets import load_iris, make_classification
from sklearn.model_selection import (KFold, StratifiedKFold, LeaveOneOut,
    cross_val_score, GridSearchCV, TimeSeriesSplit, GroupKFold)
from sklearn.linear_model import LogisticRegression
from sklearn.preprocessing import StandardScaler
from sklearn.pipeline import Pipeline
np.random.seed(42)

iris = load_iris()
pipe = Pipeline([('s',StandardScaler()),('lr',LogisticRegression(max_iter=200))])

print("=== CV对比 ===")
for name,cv in [('5-Fold',KFold(5)),('10-Fold',KFold(10)),('Stratified5',StratifiedKFold(5))]:
    s = cross_val_score(pipe, iris.data, iris.target, cv=cv)
    print("{}: mean={:.4f} std={:.4f}".format(name, s.mean(), s.std()))

s_loo = cross_val_score(pipe, iris.data, iris.target, cv=LeaveOneOut())
print("LOO: mean={:.4f}".format(s_loo.mean()))

print("\n=== 时间序列CV ===")
Xts = np.random.randn(200,3)
yts = (Xts[:,0]+Xts[:,1]>0).astype(int)
s = cross_val_score(LogisticRegression(), Xts, yts, cv=TimeSeriesSplit(5))
print("TimeSeriesSplit: mean={:.4f}".format(s.mean()))

print("\n=== 嵌套CV ===")
inner = StratifiedKFold(3); outer = StratifiedKFold(5)
grid = GridSearchCV(pipe, {'lr__C':[0.01,0.1,1,10]}, cv=inner)
ns = cross_val_score(grid, iris.data, iris.target, cv=outer)
print("嵌套CV: mean={:.4f} std={:.4f}".format(ns.mean(),ns.std()))

non_nest = GridSearchCV(pipe, {'lr__C':[0.01,0.1,1,10]}, cv=5).fit(iris.data,iris.target)
print("非嵌套(可能过估): {:.4f}".format(non_nest.best_score_))

groups = np.repeat(np.arange(50),4)
Xg,yg = make_classification(n_samples=200,n_features=5,random_state=42)
s = cross_val_score(LogisticRegression(max_iter=200),Xg,yg,cv=GroupKFold(5),groups=groups)
print("\nGroupKFold: mean={:.4f}".format(s.mean()))

运行结果

=== CV对比 === 5-Fold: mean=0.9067 std=0.0929 10-Fold: mean=0.9467 std=0.0653 Stratified5: mean=0.9600 std=0.0389 LOO: mean=0.9533 === 时间序列CV === TimeSeriesSplit: mean=0.9636 === 嵌套CV === 嵌套CV: mean=0.9733 std=0.0327 非嵌套(可能过估): 0.9733 GroupKFold: mean=0.8450
✅ 验证通过

深入理解:特征工程的精髓

"特征工程是决定模型上限的关键"

业界共识:好的特征工程比好的模型更重要。一个线性模型配上优秀特征可以超越复杂模型配普通特征。

特征工程流程

  1. 理解业务:领域知识是最好的特征来源
  2. EDA发现:从数据中发现模式→构造特征
  3. 变换:对数/Box-Cox处理偏态,分箱处理非线性
  4. 编码:类别→数值,选择合适编码方式
  5. 缩放:数值特征标准化/归一化
  6. 选择:删除无关/冗余特征
  7. 验证:交叉验证评估特征贡献

Kaggle竞赛中的特征工程技巧

常见特征工程错误

  1. 数据泄露:在特征构造中使用了未来信息
  2. 目标泄露:特征中包含了目标变量的信息
  3. 过拟合特征:在训练集上有效但测试集失效
  4. 忽略缺失模式:缺失本身可能是有价值的特征
  5. 不验证:没有用交叉验证确认特征的实际贡献

进阶特征工程技巧

Kaggle竞赛中的特征工程

在Kaggle竞赛中,特征工程是区分Top 1%和Top 10%的关键。以下是从顶级选手那里学到的技巧:

聚合特征:对时序数据,按实体ID聚合历史行为的统计量——均值、标准差、最小最大、趋势斜率、最近N次的均值变化。这些"在用户层面的历史摘要"往往是最强特征。
目标编码:高基数类别特征(如城市名2000+类)的杀手锏。关键:必须在5折CV内计算,使用平滑参数防止小样本类别不稳定。公式:enc = (count*mean + m*global_mean)/(count+m),m通常取10-100。
比率特征:两个相关特征的比值常常比原始特征更有信息量。贷款金额/年收入、月供/月收入、信用额度使用率等。这些比率捕捉了"相对"而非"绝对"信息。

数据泄露防范完整清单

  1. 目标泄露:特征包含目标信息——检查每个特征与目标的互信息
  2. 时间泄露:用了未来数据——特征构造时只用截止时间前的数据
  3. 组泄露:同组数据分布在训练和测试集——使用GroupKFold
  4. 预处理泄露:全数据fit scaler——用Pipeline自动处理
  5. 特征选择泄露:在CV外选特征——特征选择必须在CV内

缺失值处理策略决策树

交叉验证的常见错误

错误的交叉验证比没有交叉验证更危险——它给你虚假的信心。

补充:关键概念与面试要点

核心直觉

常见面试问题

  1. L1正则化为什么产生稀疏解?(几何:菱形与等高线在角上相交)
  2. 逻辑回归为什么用交叉熵而不用MSE?(MSE导致梯度消失)
  3. 随机森林为什么不容易过拟合?(Bagging降方差,特征随机降相关性)
  4. K-Means为什么收敛到局部最优?(目标函数非凸)
  5. PCA为什么需要中心化?(否则第一主成分指向均值方向)

补充:关键概念与面试要点

核心直觉

常见面试问题

  1. L1正则化为什么产生稀疏解?(几何:菱形与等高线在角上相交)
  2. 逻辑回归为什么用交叉熵而不用MSE?(MSE导致梯度消失)
  3. 随机森林为什么不容易过拟合?(Bagging降方差,特征随机降相关性)
  4. K-Means为什么收敛到局部最优?(目标函数非凸)
  5. PCA为什么需要中心化?(否则第一主成分指向均值方向)

补充:关键概念与面试要点

核心直觉

常见面试问题

  1. L1正则化为什么产生稀疏解?(几何:菱形与等高线在角上相交)
  2. 逻辑回归为什么用交叉熵而不用MSE?(MSE导致梯度消失)
  3. 随机森林为什么不容易过拟合?(Bagging降方差,特征随机降相关性)
  4. K-Means为什么收敛到局部最优?(目标函数非凸)
  5. PCA为什么需要中心化?(否则第一主成分指向均值方向)

补充:关键概念与面试要点

核心直觉

常见面试问题

  1. L1正则化为什么产生稀疏解?(几何:菱形与等高线在角上相交)
  2. 逻辑回归为什么用交叉熵而不用MSE?(MSE导致梯度消失)
  3. 随机森林为什么不容易过拟合?(Bagging降方差,特征随机降相关性)
  4. K-Means为什么收敛到局部最优?(目标函数非凸)
  5. PCA为什么需要中心化?(否则第一主成分指向均值方向)

补充:关键概念与面试要点

核心直觉

常见面试问题

  1. L1正则化为什么产生稀疏解?(几何:菱形与等高线在角上相交)
  2. 逻辑回归为什么用交叉熵而不用MSE?(MSE导致梯度消失)
  3. 随机森林为什么不容易过拟合?(Bagging降方差,特征随机降相关性)
  4. K-Means为什么收敛到局部最优?(目标函数非凸)
  5. PCA为什么需要中心化?(否则第一主成分指向均值方向)

补充:关键概念与面试要点

核心直觉

常见面试问题

  1. L1正则化为什么产生稀疏解?(几何:菱形与等高线在角上相交)
  2. 逻辑回归为什么用交叉熵而不用MSE?(MSE导致梯度消失)
  3. 随机森林为什么不容易过拟合?(Bagging降方差,特征随机降相关性)
  4. K-Means为什么收敛到局部最优?(目标函数非凸)
  5. PCA为什么需要中心化?(否则第一主成分指向均值方向)

补充:关键概念与面试要点

核心直觉

常见面试问题

  1. L1正则化为什么产生稀疏解?(几何:菱形与等高线在角上相交)
  2. 逻辑回归为什么用交叉熵而不用MSE?(MSE导致梯度消失)
  3. 随机森林为什么不容易过拟合?(Bagging降方差,特征随机降相关性)
  4. K-Means为什么收敛到局部最优?(目标函数非凸)
  5. PCA为什么需要中心化?(否则第一主成分指向均值方向)

补充:关键概念与面试要点

核心直觉

常见面试问题

  1. L1正则化为什么产生稀疏解?(几何:菱形与等高线在角上相交)
  2. 逻辑回归为什么用交叉熵而不用MSE?(MSE导致梯度消失)
  3. 随机森林为什么不容易过拟合?(Bagging降方差,特征随机降相关性)
  4. K-Means为什么收敛到局部最优?(目标函数非凸)
  5. PCA为什么需要中心化?(否则第一主成分指向均值方向)

📝 课后练习

  1. 实现本课核心算法并用Scikit-learn验证
  2. 调优关键超参数并记录性能变化
  3. 用交叉验证评估模型稳定性
  4. 与之前学过的方法对比分析
  5. 分析本课方法的失效条件
🔄
验证专家
掌握多种交叉验证策略