实战项目

第29课:房价预测

数据探索、特征工程、模型选择、调参、集成、提交结果

1. 项目概述

基于房屋特征预测售价。评估:RMSLE(对数误差,房价跨度大)。

2. 数据探索

GrLivArea(地上面积)与售价强相关。Neighborhood(区位)影响巨大。OverallQual是最强特征。面积类特征有异常值。

3. 特征工程

4. 模型训练

基线(Ridge)→强力(XGBoost/LightGBM)→调参(Optuna/GridSearch)→融合

线性+树模型互补,权重由验证集确定。

5. 完整流程

数据加载→EDA→缺失值→特征构造→编码→变换→训练→CV→调参→融合→提交

Python代码实现与验证

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

完整代码


import numpy as np
import pandas as pd
from sklearn.datasets import fetch_california_housing
from sklearn.model_selection import cross_val_score, train_test_split
from sklearn.preprocessing import StandardScaler
from sklearn.linear_model import Ridge, Lasso
from sklearn.ensemble import RandomForestRegressor, GradientBoostingRegressor
from sklearn.pipeline import Pipeline
from sklearn.metrics import mean_squared_error, r2_score
np.random.seed(42)

housing = fetch_california_housing()
X = pd.DataFrame(housing.data, columns=housing.feature_names)
y = np.log1p(housing.target)

print("数据: {} 目标: mean={:.4f} std={:.4f}".format(X.shape, y.mean(), y.std()))

Xtr,Xte,ytr,yte = train_test_split(X,y,test_size=0.2,random_state=42)
print("\n=== 模型对比 ===")
for name,mdl in [('Ridge',Pipeline([('s',StandardScaler()),('m',Ridge(1))])),
                 ('Lasso',Pipeline([('s',StandardScaler()),('m',Lasso(0.001))])),
                 ('RF',RandomForestRegressor(n_estimators=100,random_state=42,n_jobs=-1)),
                 ('GBDT',GradientBoostingRegressor(n_estimators=200,learning_rate=0.1,max_depth=4,random_state=42))]:
    mdl.fit(Xtr,ytr)
    r2 = r2_score(yte,mdl.predict(Xte))
    rmse = np.sqrt(mean_squared_error(yte,mdl.predict(Xte)))
    print("{}: R2={:.4f} RMSE={:.4f}".format(name,r2,rmse))

print("\n=== 模型融合 ===")
rf = RandomForestRegressor(n_estimators=100,random_state=42,n_jobs=-1).fit(Xtr,ytr)
gb = GradientBoostingRegressor(n_estimators=200,learning_rate=0.1,max_depth=4,random_state=42).fit(Xtr,ytr)
for w in [0.3,0.5,0.7]:
    yb = w*rf.predict(Xte)+(1-w)*gb.predict(Xte)
    print("RF权重={}: R2={:.4f}".format(w, r2_score(yte,yb)))

运行结果

数据: (20640, 8) 目标: mean=1.0561 std=0.3562 === 模型对比 === Ridge: R2=0.6006 RMSE=0.2244 Lasso: R2=0.6031 RMSE=0.2237 RF: R2=0.8249 RMSE=0.1486 GBDT: R2=0.8302 RMSE=0.1463 === 模型融合 === RF权重=0.3: R2=0.8363 RF权重=0.5: R2=0.8367 RF权重=0.7: R2=0.8342
✅ 验证通过

深入理解:从项目到产品的关键步骤

ML项目成功的关键因素

  1. 定义清晰的目标:业务指标 > ML指标(提升AUC不是目的,提升收入才是)
  2. 数据质量:Garbage in, garbage out。花80%时间在数据上是正常的
  3. 简单优先:简单基线 > 复杂模型(如果简单就够了)
  4. 持续监控:模型上线不是终点,漂移检测和重训练是必须的
  5. 沟通能力:能向非技术人解释模型为什么做某个预测

项目时间分配建议

阶段时间占比关键产出
问题定义10%明确的评估指标和成功标准
数据获取/清洗30%干净可用的数据集
特征工程25%强特征集
模型训练/调优20%最优模型
评估/部署15%上线模型+监控方案

从竞赛到生产的差距

Kaggle竞赛和实际生产的差异:

ML系统技术债

D. Sculley等人在2015年的经典论文中指出,ML系统中只有很小一部分是实际的ML代码,其余都是胶水代码、配置、数据管道、监控等。隐藏的技术债包括:

项目复盘与进阶路线

5个项目的核心收获

  1. 房价预测:回归全流程——EDA→特征工程(对数变换+交叉特征)→多模型比较→融合。关键教训:特征工程比模型选择更重要,对数变换处理右偏目标是标准操作。
  2. 欺诈检测:极度不平衡——SMOTE/类别权重/异常检测框架。关键教训:AUC-PR比AUC-ROC更有意义,监督方法远优于无监督方法。
  3. 客户分群:无监督业务应用——RFM+聚类+画像解读。关键教训:业务可解释性比聚类指标重要,K值选择要考虑运营成本。
  4. 推荐基础:推荐系统经典算法——协同过滤+矩阵分解。关键教训:冷启动是最难的问题,简单方法(SVD)常常足够好。
  5. 文本分类:NLP传统方法——TF-IDF+NB/SVM。关键教训:预处理影响巨大,线性SVM在文本上通常最好。

从竞赛到生产的差距

维度竞赛生产
数据已清洗CSV多源异构、实时流
目标单一指标多目标(精度+延迟+成本)
模型复杂集成简单高效、可维护
评估离线CV在线A/B测试
部署不涉及Docker/K8s/服务化
监控不涉及漂移检测+重训练

继续学习路线

推荐学习资源

类别资源特点
教材ISLR (Python版)统计学习入门经典
教材ESL统计学习圣经(进阶)
课程Andrew Ng ML最佳入门课
课程Stanford CS229更数学化
实践Kaggle Learn动手学ML
文档Scikit-learn最优秀的ML库文档

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

核心直觉

常见面试问题

  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. 完成房价预测完整EDA
  2. 实现特征工程(对数变换、交叉特征)
  3. 比较多个回归模型并调参
  4. 实现模型融合
  5. 分析预测误差分布
🏠
房价预言家
完成房价预测端到端项目