偏最小二乘回归分析实例plsregress命令

偏最小二乘回归分析实例

原文：http://blog.csdn.net/qq_34861102/article/details/77102683

---

plsregress命令

• 步骤
  
  • 数据标准化
  • 求相关系数矩阵
  • 分别提出自变量组合因变量组的成分
  • 求两个成分对标准化指标变量和成分变量之间的回归方程
  • 求因变量组与自变量组之间的回归方程
  • 模型的解释与检验

• 实例代码

  • ab0 = load('data.txt');%均值和方差mu = mean(ab0);sig = std(ab0);%标准化以后的自变量和因变量数据a = ab(:,[1:end-1]);b = ab(:,end);[XL YL XS YS BETA PCTVAR MSE status] = plsregress(a,b);%求自变量的主成分系数xw = a\XS;%求因变量的主成分系数yw = b\YS;ncomp = input('请根据PCTVAR的值确定提出成分对的个数');[XL2 YL2 XS2 YS2 BETA2 PCTVAR2 MSE2 status2] = plsregress(a,b,ncomp);%自变量的个数n = size(a,2);%原始数据回归方程的常数项beta3(1) = mu(end) - mu(1:n)./sig(1:n)*BETA2([2:end]).*sig(end);beta3([2:n+1]) = (1./sig(1:n))'*sig(n+1:end).*BETA2([2:end]);%直方图bar(BETA2','k');%y的预测值yhat = repmat(beta3(1),[size(a,1),1]) + ab0(:,[1:n])*beta3([2:end])';%预测值的最大值ymax = max([yhat;ab0(:,end)]);