plspred.m

来自「偏最小二乘算法在MATLAB中的实现」· M 代码 · 共 32 行

function ypred = plspred(x,b,p,q,w,lv)
%PLSPRED Predictions based on existing PLS model.
%  The inputs are the matrix of predictor variables (x),
%  the PLS model inner-relation coefficients (b), the
%  x-block loadings (p), the y-block loadings (q),
%  the x-block weights (w), and the number of  latent 
%  variables to use in prediction (lv). The output is
%  the vector or matrix of the predicted values (ypred).
%  I/O format is: ypred = plspred(x,b,p,q,w,lv);

%  Copyright
%  Barry M. Wise
%  1991
%  Modified by B.M. Wise, November 1993

[mx,nx] = size(x);
[mq,nq] = size(q);
[mw,nw] = size(w);
that = zeros(mx,lv);
ypred = zeros(mx,mq);
if lv > nw
  s = sprintf('Maximum number of latent variables exceeded (Max = %g)',nw);
  error(s)
end
for i = 1:lv
  that(:,i) = x*w(:,i);
  x = x - that(:,i)*p(:,i)';
end
for i = 1:lv
  ypred = ypred + b(1,i)*that(:,i)*q(:,i)';
end