lnoxvalsompls.m

有监督自组织映射-偏最小二乘算法(A supervised self-organising map–partial least squares algorithm),可以用语多变量数据的回归分析

function [MeanRmseVal, StdRmseVal, RmseVal, Ytest, Ypred] = LNOXValSOMPLS(BestMLK,MLKP,OPT,Ntest);

Tbegin=cputime;

XSav=[BestMLK.Xtr; BestMLK.Xtest];
YSav=[BestMLK.Ytr; BestMLK.Ytest];
[Nobj,NvarX]=size(XSav);
[Nobj,NvarY]=size(YSav);
if (upper(MLKP.ProblemType) == 'REG')
    RmseVal=zeros(1,NvarY);
else
    RmseVal=0;
end
Nval=floor(Nobj/Ntest);
Nrest=Nobj-Nval*Ntest;
MLKP.SilentMode='Y';
Ytest=[];
Ypred=[];

for ival=1:Nval
    X=XSav;
    BestMLK.Xtest=X(1+(ival-1)*Ntest:ival*Ntest,:);
    X(1+(ival-1)*Ntest:ival*Ntest,:)=[];
    BestMLK.Xtr=X;
    Y=YSav;
    BestMLK.Ytest=Y(1+(ival-1)*Ntest:ival*Ntest,:);
    Y(1+(ival-1)*Ntest:ival*Ntest,:)=[];
    BestMLK.Ytr=Y;
    [PCTrain,PCTest,NBestMLK,NMLKP] = MLZValidate(BestMLK,MLKP,OPT);
    Ytest=[Ytest; NBestMLK.Ytest];
    Ypred=[Ypred; NBestMLK.PredTest];
    if (upper(MLKP.ProblemType) == 'REG')
        RmseVal(ival,:)=NBestMLK.RmseTest;
    else
        RmseVal(ival)=NBestMLK.PercCorrTest;
    end
end
if (Nrest > 0)
    ival=Nval+1;
    X=XSav;
    BestMLK.Xtest=X(1+(ival-1)*Ntest:end,:);
    X(1+(ival-1)*Ntest:end,:)=[];
    BestMLK.Xtr=X;
    Y=YSav;
    BestMLK.Ytest=Y(1+(ival-1)*Ntest:end,:);
    Y(1+(ival-1)*Ntest:end,:)=[];
    BestMLK.Ytr=Y;
    [PCTrain,PCTest,NBestMLK,NMLKP] = MLZValidate(BestMLK,MLKP,OPT);
    Ytest=[Ytest; NBestMLK.Ytest];
    Ypred=[Ypred; NBestMLK.PredTest];
    if (upper(MLKP.ProblemType) == 'REG')
        RmseVal(ival,:)=NBestMLK.RmseTest;
    else
        RmseVal(ival)=NBestMLK.PercCorrTest;
    end
    Nval=Nval+1;
end
MeanRmseVal=mean(RmseVal);
StdRmseVal=std(RmseVal);

Tend=cputime;
Message=sprintf('Total time %g seconds',Tend-Tbegin);
disp(Message);