plspred.m

来自「wang xiao ping 版遗传算法」· M 代码 · 共 20 行

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