plspred.m

王小平《遗传算法——理论、应用与软件实现》随书光盘

function ypred = plspred(x,b,p,q,w,lv)


[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