DefinedCentroids = importdata('centroids_out.txt');
SampleCentroids = importdata('out.txt');

Threshold=0.15;

clear Corr;
clear ST;

UncenteredData=SampleCentroids;

for i=1:numel(SampleCentroids(:,1))
    SampleCentroids(i,:)=SampleCentroids(i,:)-mean(SampleCentroids(i,:));
end

%for MAD centering, use:
%SampleCentroids=manorm(SampleCentroids, 'Method', 'MAD');

for i=1:numel(SampleCentroids(1,:))
    MaxScore=-1;
    ChosenSubtype(i)=99;
    
    for j=1:5
        Corr(i,j)=corr(SampleCentroids(:,i), DefinedCentroids(:,j), 'rows', 'pairwise');
        if(Corr(i,j)>Threshold && Corr(i,j)>MaxScore)
            ChosenSubtype(i)=j;
            MaxScore=Corr(i,j);
        end
    end
end

fid=fopen('ChosenSubtype', 'w');
fprintf(fid, '%i\t', ChosenSubtype(:));
fclose(fid);