特征工程

第25课:缺失值处理

缺失机制、删除法、单值填充、模型填充、多重插补、Scikit-learn实现

1. 缺失机制

类型定义处理
MCAR完全随机删除/简单填充
MAR与已观测有关模型填充
MNAR与未观测有关领域知识

2. 方法

删除法: MCAR+低缺失率可用,MAR/MNAR下有偏差

单值填充: 均值/中位数/众数(低估方差)、常数+缺失指示列

KNN填充: 用K个近邻均值。迭代填充(MICE): 用其他特征训练模型预测缺失值

多重插补: M个完整数据集→分别分析→Rubin合并。M=5-20。

3. 流程

分析缺失机制→>50%缺失考虑删除→MCAR低缺失率→删除→MAR→MICE/KNN→MNAR→指示列+敏感性分析

Python代码实现与验证

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

完整代码


import numpy as np
from sklearn.impute import SimpleImputer, KNNImputer
from sklearn.experimental import enable_iterative_imputer
from sklearn.impute import IterativeImputer
from sklearn.datasets import load_iris
from sklearn.model_selection import cross_val_score
from sklearn.linear_model import LogisticRegression
from sklearn.preprocessing import StandardScaler
from sklearn.pipeline import Pipeline
np.random.seed(42)

iris = load_iris()
X = iris.data.copy()
mask = np.random.random(X.shape) < 0.1
Xm = X.copy(); Xm[mask] = np.nan
print("缺失率: {:.2%}".format(np.mean(mask)))

print("\n=== 填充方法对比 ===")
for name,imp in [('均值',SimpleImputer(strategy='mean')),('中位数',SimpleImputer(strategy='median')),
                 ('KNN',KNNImputer(n_neighbors=5)),('迭代',IterativeImputer(max_iter=10,random_state=42))]:
    Xf = imp.fit_transform(Xm)
    mse = np.mean((Xf[mask]-X[mask])**2)
    print("{}: MSE={:.4f}".format(name, mse))

print("\n=== 填充后模型性能 ===")
for name,imp in [('均值',SimpleImputer(strategy='mean')),('KNN',KNNImputer(n_neighbors=5)),
                 ('迭代',IterativeImputer(max_iter=10,random_state=42))]:
    pipe = Pipeline([('imp',imp),('s',StandardScaler()),('lr',LogisticRegression(max_iter=200))])
    s = cross_val_score(pipe, Xm, iris.target, cv=5)
    print("{}: CV准确率={:.4f}".format(name, s.mean()))

运行结果

缺失率: 12.00% === 填充方法对比 === 均值: MSE=1.2993 中位数: MSE=1.3093 KNN: MSE=0.2245 迭代: MSE=0.1393 === 填充后模型性能 === 均值: CV准确率=0.9200 KNN: CV准确率=0.9533 迭代: CV准确率=0.9667
✅ 验证通过

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

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

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

特征工程流程

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

📝 课后练习

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