calccorrvalues.m

来自「有监督自组织映射-偏最小二乘算法(A supervised self-organ」· M 代码 · 共 22 行

function [CorrValues] = CalcCorrValues(X,Map);

[Nvar,Nunit]=size(Map);

X=X-mean(X);
XX=X*X';
if (XX == 0)
    CorrValues=zeros(1,Nunit);
else
    MapMean=mean(Map);
    Map=Map-MapMean(ones(Nvar,1),:);
    YY=sum(Map.*Map);
    XY=X*Map;
    for iunit=1:Nunit
        if (YY(iunit) > 0)
            CorrValues(iunit)=XY(iunit)/sqrt(XX*YY(iunit));
        else
            CorrValues(iunit)=0;
        end
    end
end