SS_IdentifyBestOps.m

function SS_OutputBestOpsTxtFile

load('run_options.mat');
load('HCTSA_N.mat');
load('linkage_clusters_with_member_corrs.mat');

fprintf('Saving clustering information to %s \n',outTxtFileName);

fID = fopen(outTxtFileName,'w');
fprintf(fID,['Output from kmedoids clustering (k = %i) followed by ',...
    'linkage clustering (corr_dist_threshold = %f)\n\n'],kmedoidsClusters.k,corr_dist_threshold);

autoChosenOps = [];
autoChosenIdxs = [];

for i = 1:length(linkageClusters)
    dists = linkageClusters(i).memDists;
    mems = linkageClusters(i).memIdxs;
    kmedCentres = linkageClusters(i).kmedCentres;
    ops = Operations(mems);
    
    % Sort by increasing distance
    [sortedDists sortIdx] = sort(dists);
    sortedMems = mems(sortIdx);
    sortedOps = ops(sortIdx);
    
    kmCentreOps = Operations(kmedoidsClusters.CCi(kmedCentres));
    kmCentreNames = {kmCentreOps.Name};
    kmCentreKeys = {kmCentreOps.Keywords};
    
    % Print all centres associated with cluster
    fprintf(fID,'CLUSTER %i: ',i);
    for j = 1:length(kmCentreNames)
        fprintf(fID,'%s (%s), ',kmCentreNames{j},kmCentreKeys{j}); 
    end
    fprintf(fID,'\n');
    
    % Print top 10 operations highly correlated with cluster centre
    for j = 1:min(length(sortedOps),10)
        fprintf(fID,'(%f) %s - %s\n',1-sortedDists(j),sortedOps(j).Name,...
            sortedOps(j).Keywords);
    end
    
    autoChosenOps = [autoChosenOps sortedOps(1)];
    autoChosenIdxs = [autoChosenIdxs sortedMems(1)];

    fprintf(fID,'\n');
end

fclose(fID);

save('auto_chosen_ops.mat','autoChosenOps','autoChosenIdxs');

end