QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2778|回复: 0
打印 上一主题 下一主题

偏最小二乘回归分析

[复制链接]
字体大小: 正常 放大

1189

主题

4

听众

2934

积分

该用户从未签到

跳转到指定楼层
1#
发表于 2023-11-30 16:59 |只看该作者 |正序浏览
|招呼Ta 关注Ta
偏小二乘回归提供一种多对多线性回归建模的方法,特别当两组变量的个数很多,且都存在多重相关性,而观测数据的数量(样本量)又较少时,用偏小二乘回归建立的模型具有传统的经典回归分析等方法所没有的优点。
5 N3 K9 {# T( T+ D0 f9 x& M偏小二乘回归分析在建模过程中集中了主成分分析,典型相关分析和线性回归分析方法的特点,因此在分析结果中,除了可以提供一个更为合理的回归模型外,还可以同时完成一些类似于主成分分析和典型相关分析的研究内容,提供更丰富、深入的一些信息。
  1. import numpy as np) ~, M* {7 f/ ^
  2. from sklearn.cross_decomposition import PLSRegression
  3. 5 R) M) I, o* E* _) T1 G
  4. from sklearn.preprocessing import StandardScaler
  5. 1 U( b, \: ~) c5 C6 q\\" F  t) I
  6. import matplotlib.pyplot as plt# b; v\\" @8 ], J+ T
  7. - |/ X2 K- \- D
  8. # 原始数据,自变量X和因变量y3 I& o  `$ x0 E7 u
  9. X = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9], [10, 11, 12]])  v\\" c7 l3 ?+ S$ o, E1 }+ T
  10. y = np.array([1, 2, 3, 4])
  11. & E0 J\\" G7 F& ~: F& ?; x

  12. ( q, ^\\" O3 r8 G6 i* P& ?
  13. # 对自变量进行标准化处理% t; {+ h( l- `' \5 m  l8 F( `$ g
  14. scaler = StandardScaler()1 g$ Z& i$ k( F+ K, y
  15. X_scaled = scaler.fit_transform(X); p6 ^( M8 D- z0 i1 ]' A
  16. - Y4 j9 n\\" |5 d) Z
  17. # 创建PLSRegression对象,并指定主成分个数为2
  18. 7 M3 a% A; g2 i$ B3 w
  19. pls = PLSRegression(n_components=2)3 P. w) e6 ]- Z/ z- c7 u
  20. 8 h0 F+ m) t+ T$ x/ b# u; L' R\\" C! h8 }
  21. # 进行偏最小二乘回归分析
  22. , k% ^# X, T\\" I( l( a1 Y\\" U& t3 D
  23. pls.fit(X_scaled, y)
  24. $ _2 e0 p' \5 k& b! m+ ^# M
  25. : \; [* w+ W! O* p9 ~! v6 z  a
  26. # 获得预测值. U& K( S- w1 D\\" h
  27. y_pred = pls.predict(X_scaled)
  28. 2 U' e) @# P1 Z$ x# D\\" ~. V

  29. 0 c! h: F  t& O\\" [) K/ x\\" C
  30. # 绘制原始数据和预测值
  31. / y( \$ ]: j) i# d8 ^( Z, U\\" p- k0 A
  32. plt.figure(figsize=(8, 6))) _/ m( {- S9 J' P4 x$ [\\" [
  33. plt.scatter(y, y_pred, c='b', label='Predicted vs. Observed')
  34. . L) T\\" T. L6 m! [$ U  s0 S8 F
  35. plt.plot([min(y), max(y)], [min(y), max(y)], 'k--', lw=2, label='Perfect Fit')
  36. ( B1 x' x2 Y4 S! d. Y+ ]$ k
  37. plt.xlabel('Observed')
  38. , r8 t/ W! l5 ?) o
  39. plt.ylabel('Predicted')
  40. $ p3 z5 C- Z\\" t1 J0 }
  41. plt.title('PLS Regression'): v\\" [7 ?+ D2 e. g' j
  42. plt.legend()+ F. Z3 W3 k6 B. K) c7 J  m
  43. plt.grid(True)8 j! \9 w$ v1 @' u
  44. plt.show()
  45. , E* y* P9 H/ }; j: H/ E

  46. 4 q& q8 ^5 S6 Z! s8 U5 R' m* D
在上述代码中,我们首先定义了原始数据矩阵X和因变量向量y,然后创建了一个PLSRegression对象并指定保留的主成分个数为2。接下来,使用fit()方法进行偏最小二乘回归分析,并使用coef_属性获取回归系数。最后,使用predict()方法对原始数据进行预测,并获得预测值y_pred。请根据实际情况调整原始数据和保留的主成分个数,并根据需要对结果进行解释和分析。2 o- r$ n: l- v! Q8 p5 A
运行上述代码,将绘制散点图来展示观察值(原始数据)与预测值之间的关系。如果模型拟合良好,散点图中的点应该基本落在对角线上。如果点主要集中在对角线上方,则预测值高估了观察值;如果点主要集中在对角线下方,则预测值低估了观察值。在理想情况下,所有点都应该在对角线上。
- B# K6 y/ J1 B VeryCapture_20231130164832.jpg
1 P1 ?& [2 l$ q; J/ l0 p, {( X8 @. d+ h+ I/ E! z, f4 {# Z
6 z% d) E, D% V  F6 a
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-7-10 09:50 , Processed in 0.453631 second(s), 54 queries .

回顶部