实战项目

第35课:毕业项目:端到端ML Pipeline

完整ML项目流程、数据版本控制、模型部署、监控、MLOps实践

1. 端到端ML Pipeline

数据获取→数据验证→特征工程→模型训练→模型验证→模型部署→性能监控→反馈循环

2. 数据管理

DVC(数据版本)、Great Expectations(数据质量)、Feature Store(共享特征)

3. 模型训练

实验跟踪(MLflow/W&B)、超参数优化(Optuna/Bayesian)、模型选择策略(快基线→增复杂度→集成→微调)

4. 模型部署

模式延迟适用
批处理分钟-小时离线推荐
实时API毫秒在线预测
流式实时监控

序列化:Pickle/Joblib/ONNX。容器化:Docker+K8s。

5. 模型监控

性能漂移(准确率下降)、数据漂移(PSI>0.25=显著)、概念漂移(P(Y|X)变化)

重训练触发:定期/性能低/漂移显著

6. MLOps成熟度

L0(手动)→L1(Pipeline自动化)→L2(CI/CD)→L3(持续训练全自动)

7. 课程总结

阶段课次核心
统计基础1-7概率→分布→检验→信息论
监督学习8-14线性→树→集成→SVM→KNN
无监督15-21聚类→降维→关联→异常
特征工程22-28选择→构造→编码→评估
实战项目29-355个完整项目

继续学习: 深度学习(CNN/RNN/Transformer)、强化学习、MLOps生产化

Python代码实现与验证

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

完整代码


import numpy as np
import pandas as pd
import time
from sklearn.datasets import fetch_california_housing
from sklearn.model_selection import train_test_split, GridSearchCV
from sklearn.preprocessing import StandardScaler
from sklearn.compose import ColumnTransformer
from sklearn.pipeline import Pipeline
from sklearn.impute import SimpleImputer
from sklearn.ensemble import GradientBoostingRegressor, RandomForestRegressor
from sklearn.linear_model import Ridge
from sklearn.metrics import mean_squared_error, r2_score, mean_absolute_error
import joblib
from scipy import stats
np.random.seed(42)

housing = fetch_california_housing()
X = pd.DataFrame(housing.data, columns=housing.feature_names)
y = housing.target
Xtr, Xte, ytr, yte = train_test_split(X, y, test_size=0.2, random_state=42)

preproc = ColumnTransformer([('num', Pipeline([('imp',SimpleImputer(strategy='median')),('s',StandardScaler())]), X.columns)])

print("=== 端到端ML Pipeline ===")
print("1.数据获取: {} {}".format(X.shape, '无缺失' if X.isnull().sum().sum()==0 else '有缺失'))

results = {}
for name,mdl in [('Ridge',Ridge(1)),('RF',RandomForestRegressor(n_estimators=100,random_state=42,n_jobs=-1)),
                 ('GBDT',GradientBoostingRegressor(n_estimators=200,max_depth=4,learning_rate=0.1,random_state=42))]:
    pipe = Pipeline([('p',preproc),('m',mdl)])
    t0=time.time(); pipe.fit(Xtr,ytr); t1=time.time()
    yp=pipe.predict(Xte)
    results[name]={'R2':r2_score(yte,yp),'RMSE':np.sqrt(mean_squared_error(yte,yp)),'MAE':mean_absolute_error(yte,yp),'time':t1-t0}
    print("{}: R2={:.4f} RMSE={:.4f} MAE={:.4f} ({:.1f}s)".format(name,results[name]['R2'],results[name]['RMSE'],results[name]['MAE'],results[name]['time']))

print("\n5.超参数优化")
pipe_gb = Pipeline([('p',preproc),('m',GradientBoostingRegressor(random_state=42))])
grid = GridSearchCV(pipe_gb,{'m__n_estimators':[100,200],'m__max_depth':[3,5],'m__learning_rate':[0.05,0.1]},cv=3,scoring='r2',n_jobs=-1)
grid.fit(Xtr,ytr)
print("最佳: {} CV={:.4f} Test={:.4f}".format(grid.best_params_, grid.best_score_, grid.score(Xte,yte)))

joblib.dump(grid.best_estimator_, '/tmp/best_ml_model.pkl')
loaded = joblib.load('/tmp/best_ml_model.pkl')
print("\n6.模型保存/加载: 预测一致={}".format(np.allclose(loaded.predict(Xte),grid.predict(Xte))))

print("\n7.数据漂移检测(KS检验)")
for col_idx, col in enumerate(X.columns):
    ks_p = stats.ks_2samp(Xtr[col].values if hasattr(Xtr,'columns') else Xtr[:,col_idx], Xte[col].values if hasattr(Xte,'columns') else Xte[:,col_idx]).pvalue
    if ks_p < 0.05:
        print("  {} 可能漂移 (p={:.4f})".format(col, ks_p))

print("\n=== Pipeline完成 ===")
print("全流程: 数据获取->验证->划分->预处理->训练->优化->序列化->监控")

运行结果

⚠️ Command '['python3', '-c', '\nimport numpy as np\nimport pandas as pd\nimport time\nfrom sklearn.datasets import fetch_california_housing\nfrom sklearn.model_selection import train_test_split, GridSearchCV\nfrom sklearn.preprocessing import StandardScaler\nfrom sklearn.compose import ColumnTransformer\nfrom sklearn.pipeline import Pipeline\nfrom sklearn.impute import SimpleImputer\nfrom sklearn.ensemble import GradientBoostingRegressor, RandomForestRegressor\nfrom sklearn.linear_model import Ridge\nfrom sklearn.metrics import mean_squared_error, r2_score, mean_absolute_error\nimport joblib\nfrom scipy import stats\nnp.random.seed(42)\n\nhousing = fetch_california_housing()\nX = pd.DataFrame(housing.data, columns=housing.feature_names)\ny = housing.target\nXtr, Xte, ytr, yte = train_test_split(X, y, test_size=0.2, random_state=42)\n\npreproc = ColumnTransformer([(\'num\', Pipeline([(\'imp\',SimpleImputer(strategy=\'median\')),(\'s\',StandardScaler())]), X.columns)])\n\nprint("=== 端到端ML Pipeline ===")\nprint("1.数据获取: {} {}".format(X.shape, \'无缺失\' if X.isnull().sum().sum()==0 else \'有缺失\'))\n\nresults = {}\nfor name,mdl in [(\'Ridge\',Ridge(1)),(\'RF\',RandomForestRegressor(n_estimators=100,random_state=42,n_jobs=-1)),\n (\'GBDT\',GradientBoostingRegressor(n_estimators=200,max_depth=4,learning_rate=0.1,random_state=42))]:\n pipe = Pipeline([(\'p\',preproc),(\'m\',mdl)])\n t0=time.time(); pipe.fit(Xtr,ytr); t1=time.time()\n yp=pipe.predict(Xte)\n results[name]={\'R2\':r2_score(yte,yp),\'RMSE\':np.sqrt(mean_squared_error(yte,yp)),\'MAE\':mean_absolute_error(yte,yp),\'time\':t1-t0}\n print("{}: R2={:.4f} RMSE={:.4f} MAE={:.4f} ({:.1f}s)".format(name,results[name][\'R2\'],results[name][\'RMSE\'],results[name][\'MAE\'],results[name][\'time\']))\n\nprint("\\n5.超参数优化")\npipe_gb = Pipeline([(\'p\',preproc),(\'m\',GradientBoostingRegressor(random_state=42))])\ngrid = GridSearchCV(pipe_gb,{\'m__n_estimators\':[100,200],\'m__max_depth\':[3,5],\'m__learning_rate\':[0.05,0.1]},cv=3,scoring=\'r2\',n_jobs=-1)\ngrid.fit(Xtr,ytr)\nprint("最佳: {} CV={:.4f} Test={:.4f}".format(grid.best_params_, grid.best_score_, grid.score(Xte,yte)))\n\njoblib.dump(grid.best_estimator_, \'/tmp/best_ml_model.pkl\')\nloaded = joblib.load(\'/tmp/best_ml_model.pkl\')\nprint("\\n6.模型保存/加载: 预测一致={}".format(np.allclose(loaded.predict(Xte),grid.predict(Xte))))\n\nprint("\\n7.数据漂移检测(KS检验)")\nfor col_idx, col in enumerate(X.columns):\n ks_p = stats.ks_2samp(Xtr[col].values if hasattr(Xtr,\'columns\') else Xtr[:,col_idx], Xte[col].values if hasattr(Xte,\'columns\') else Xte[:,col_idx]).pvalue\n if ks_p < 0.05:\n print(" {} 可能漂移 (p={:.4f})".format(col, ks_p))\n\nprint("\\n=== Pipeline完成 ===")\nprint("全流程: 数据获取->验证->划分->预处理->训练->优化->序列化->监控")\n']' timed out after 120 seconds
⚠️ 验证部分失败

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

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. 实现完整数据Pipeline
  2. 用MLflow跟踪实验
  3. 实现超参数优化
  4. 模型序列化并部署为API
  5. 实现数据漂移检测
🎓
ML工程师
完成端到端ML Pipeline毕业项目!