for l=1:10
[xdata,ydata]=tichu(datath ,lab3,l);
[xdata,ydata]=elim(xdata,ydata);
xxdata=snv(xdata);
[Xtrain,Xtest,Ytrain,Ytest]=divided(xxdata,ydata);
model=pls(Xtrain,Ytrain,5);
y1=Xtest*model(5,:)';
r(l)=corrcoef(y1,Ytest);
R(l)=r(01,2)^2;
RMSE(l)=sqrt(mse(y1 - Ytest));
MRE(l)=mean(abs(Ytest-y1)./Ytest);
RSD(l)=std(y1);
end