calculatey.asv

一个简单的pls程序吧

function yy=calculatey(x,y,xt)
[inds, dists] = findneib(x,xt,6);
% [inds, dists] = findneibkpls(x,xt);
x1=x(inds,:);y1=y(inds,:);
[nx,mx] =size(x1);
avgx=mean(x1);
varx=std(x1);
centerx = x1 - avgx(ones(nx,1),:);
centerx=(x1- avgx(ones(nx,1),:))./ varx(ones(nx,1),:);
avgy=mean(y1);
vary=std(y1);
centery = (y1 - avgy(ones(nx,1),:))./vary(ones(nx,1),:);
% [Y_hat,Yt_hat,B]=PCR_cent(x1,xt,y1,[],4);
% yy=Yt_hat;

[B,T,U,R]= PLS_SIMPLS(centerx,centery,5);
yy=(((xt-avgx)./varx)*B).*vary+avgy;