最近我们被客户要求撰写关于主成分回归的研究报告。
此示例显示如何在matlab中应用偏最小二乘回归（PLSR）和主成分回归（PCR），并讨论这两种方法的有效性。当存在大量预测变量时，PLSR和PCR都是对因变量建模的方法，并且这些预测变量高度相关或甚至共线性。
在MATLAB中，可以使用Statistics and Machine Learning Toolbox中的函数来执行偏最小二乘回归（PLSR）和主成分回归（PCR）。下面，我将通过一个简单的例子来说明如何在MATLAB中应用这两种方法，并讨论它们的有效性。
首先，确保你已经安装了MATLAB的Statistics and Machine Learning Toolbox。
偏最小二乘回归（PLSR）
偏最小二乘回归是一种用于解决预测变量之间共线性和多重共线性问题的回归方法。它结合了多元线性回归、主成分分析和典型相关分析的思想。
可下载资源
作者
两种方法都将新的预测变量（称为成分）构建为原始预测变量的线性组合，但它们以不同的方式构造这些成分。PCR创建成分来解释预测变量中观察到的变异性，而根本不考虑因变量。另一方面，PLSR确实将因变量考虑在内，因此通常会导致模型能够使用更少的成分来适应因变量。
偏最小二乘法与主成分回归，都是对主成分进行回归，解决多元线性回归中的多重共线性的问题，但是它们之间还是有很多不同的点。
主成分回归，在这里也不细说了，学统计的同学肯定都知道主成分分析，提取其中K个主成分对y进行回归，就是主成分回归了。主成分个数的确定，可以通过信息量占比来确定，也可以通过交叉验证来确定最优的主成分个数。
偏最小二乘法则是：多元线性回归分析+典型相关分析+主成分分析
偏最小二乘法不仅能够解决变量多重共线性的问题，还能解决多因变量与自变量的回归问题。现只举单因变量回归的步骤：
Step1：找出与y最相关的主成分t1
Step2:对求得y对t1回归后的残差向量，再重复Step1，得到主成分t2
Step3：依此类推，根据你要的主成分数进行相应的回归
到底选多少个主成分，可以根据交叉检验选择最优的主成分数
对于多因变量的回归问题，偏最小二乘法利用了典型相关分析，在这里也就不细说了。想要了解偏最小二乘法更多的相关理论知识，大家可以看看《偏最小二乘回归方法及应用》，非常详细。
到这里，应该就可以看出两者的区别了：
1.主成分回归只利用了x的信息，根据协方差矩阵求得主成分，而偏最小二乘法利用了y和x的信息
2.偏最小二乘法要比主成分回归效果要好
3.主成分回归是对数据做了一个正交变换，因此主成分之间都是正交的，而偏最小二乘法不是
4.在确定主成分个数的时候，两者都可以通过交叉检验确定，但是主成分回归一般使用信息占比值来确定，一般来说大于0.8就行了
加载数据
加载包括401个波长的60个汽油样品的光谱强度及其辛烷值的数据集。
set(gcf,'DefaultAxesColorOrder',jet(60));
plot3(repmat(1:401,60,1)',repmat(octane(h),1,401)',NIR(h,:)');
set(gcf,'DefaultAxesColorOrder',oldorder);
xlabel('Wavelength Index'); ylabel('Octane'); axis('tight');
grid on
使用两个拟合数据
使PLSR模型拟合10个PLS成分和一个因变量。
为了充分拟合数据，可能需要十个成分，但可以使用此拟合的诊断来选择具有更少成分的更简单模型。例如，选择成分数量的一种快速方法是将因变量中解释的方差百分比绘制为成分数量的函数。
在实践中，在选择成分数量时可能需要更加谨慎。例如，交叉验证是一种广泛使用的方法，稍后将在本示例中进行说明。目前，上图显示具有两个成分的PLSR解释了观察到的大部分方差y。计算双组分模型的拟合因变量。
plot(1:10,cumsum(100*PLSPctVar(2,:)),'-bo');
xlabel('Number of PLS components');
ylabel('Percent Variance Explained in Y');接下来，拟合具有两个主要成分的PCR模型。第一步是X使用该pca函数执行主成分分析，并保留两个主成分。
然后，PCR只是这两个成分的因变量的线性回归。当变量具有非常不同的可变性时，通常首先通过其标准偏差来规范每个变量。
plot(y,yfitPLS,'bo',y,yfitPCR,'r^');xlabel('Observed Response');ylabel('Fitted Response');legend({'PLSR with 2 Components' 'PCR with 2 Components'}, ... 'location','NW');从某种意义上说，上图中的比较并不合理 – 通过观察双组分PLSR模型预测因变量的程度来选择成分数（两个），并且没有说明为什么PCR模型应该限制相同数量的成分。然而，使用相同数量的成分，PLSR做得更好。实际上，观察上图中拟合值的水平分布，使用两个分量的PCR几乎不比使用常数模型好。回归的r方值证实了这一点。
比较两种模型的预测能力的另一种方法是在两种情况下将因变量绘制成两个预测变量。
plot3(Xscores(:,1),Xscores(:,2),y-mean(y),'bo');
legend('PLSR');
grid on; view(-30,30);
如果不能以交互方式旋转图形，有点难以看到，但上面的PLSR图显示了紧密分散在平面上的点。另一方面，下面的PCR图显示点几乎没有线性关系。
plot3(PCAScores(:,1),PCAScores(:,2),y-mean(y),'r^');
legend('PCR');
grid on; view(-30,30);
请注意，尽管两个PLS成分是观察到的更好的预测因子，但下图显示它们解释的方差比例比PCR中使用的前两个主成分少。
随时关注您喜欢的主题
PCR曲线一致性较高的事实表明，为什么使用两种成分的PCR相对于PLSR在拟合时表现很差。PCR构建成分以便最好地解释X，因此，前两个成分忽略了数据拟合中观察到的重要信息y。
拟合更多成分
随着在PCR中添加更多成分，它必然会更好地拟合原始数据y，这仅仅是因为在某些时候，大多数重要的预测信息X将存在于主要成分中。例如，使用10个成分时，两种方法的残差远小于两个成分的残差。
交叉验证
在预测未来变量的观察结果时，选择成分数量以减少预期误差通常很有用。简单地使用大量成分将很好地拟合当前观察到的数据，但这是一种导致过度拟合的策略。过于拟合当前数据会导致模型不能很好地推广到其他数据，并对预期误差给出过度乐观的估计。
交叉验证是一种统计上更加合理的方法，用于选择PLSR或PCR中的成分数量。它通过不重复使用相同的数据来拟合模型和估计预测误差来避免过度拟合数据。因此，预测误差的估计不会乐观地向下偏差。
pls可以选择通过交叉验证来估计均方预测误差（MSEP），在这种情况下使用10倍CV。
[X1，Y1，Xs，Ys，β，pctVar，PLSmsep] = plsregress（X，y，10，'CV'，10）;对于PCR，crossval结合用于计算PCR的平方误差之和，可以再次使用10倍交叉验证来估计MSEP。
PCRmsep = sum（crossval（@ pcrsse，X，y，'KFold'，10），1）/ n;PLSR的MSEP曲线表明两个或三个成分好。另一方面，PCR需要四个成分才能获得相同的预测精度。
事实上，PCR中的第二个成分会增加模型的预测误差，这表明该成分中包含的预测变量的组合与其没有很强的相关性y。再次，这是因为PCR构建成分来解释X，而不是y。
模型简约
因此，如果PCR需要四个成分来获得与具有三个成分的PLSR相同的预测精度，那么PLSR模型是否更加简约？这取决于您考虑的模型的哪个方面。
PLS权重是定义PLS分量的原始变量的线性组合，即，它们描述了PLSR中的每个分量依赖于原始变量的权重。
类似地，PCA载荷描述了PCR中每个成分依赖于原始变量的强度。
对于PLSR或PCR，可以通过检查每个成分最重要的变量来为每个成分提供有意义的解释。例如，利用这些光谱数据，可以根据汽油中存在的化合物解释强度峰值，然后观察特定成分的权重挑选出少量这些化合物。从这个角度来看，更少的成分更易于解释，并且由于PLSR通常需要更少的成分来充分预测因变量，因此会导致更简约的模型。
另一方面，PLSR和PCR都导致每个原始预测变量的一个回归系数加上截距。从这个意义上讲，两者都不是更简约，因为无论使用多少成分，两种模型都依赖于所有预测变量。更具体地，对于这些数据，两个模型都需要401个光谱强度值以进行预测。
然而，最终目标可能是将原始变量集减少到仍然能够准确预测因变量的较小子集。例如，可以使用PLS权重或PCA载荷来仅选择对每个成分贡献最大的那些变量。如前所示，来自PCR模型拟合的一些成分可主要用于描述预测变量的变化，并且可包括与因变量不强相关的变量的权重。因此，PCR会导致保留预测不必要的变量。
对于本例中使用的数据，PLSR和PCR所需的成分数量之间的差异不是很大，PLS权重和PCA载荷选择了相同的变量。其他数据可能并非如此。
可下载资源
关于作者
Kaizong Ye是拓端研究室（TRL）的研究员。Kaizong Ye是拓端研究室（TRL）的研究员。在此对他对本文所作的贡献表示诚挚感谢，他在上海财经大学完成了统计学专业的硕士学位，专注人工智能领域。擅长Python.Matlab仿真、视觉处理、神经网络、数据分析。
本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。
非常感谢您阅读本文，如需帮助请联系我们！
R语言混合模型回归GBTM群组轨迹模型绘图可视化研究
MATLAB用CNN-LSTM神经网络的语音情感分类深度学习研究
R语言非参数回归预测摩托车事故、收入数据：局部回归、核回归、LOESS可视化
【视频讲解】非参数重采样bootstrap的逻辑回归Logistic应用及模型差异Python实现
1 comment on “matlab中的偏最小二乘回归（PLSR）和主成分回归（PCR）”