function [linkList, featureList]= treeLinkFeatures(data, depthThresh)
%% linkList = treeLinkFeatures(data)
% given a dataset, a hierarchical cluster of the data is produced, and then
% the data is traversed, such that, for each split in the data, a set of
% features are produced, which are the ranked features that can be used to
% separate the given dataset at that point.



if(nargin < 2)
    depthThresh = 9999999;
end
linkList = aglomCluster(data);
linkList = depthCheck(linkList);
listSize = size(data,1);

% linkList(:,4) = 0;
featureList = cell(listSize-1,1);
currentRow = [2*listSize-1];

%%
while (~isempty(currentRow))
    if(currentRow(1) > listSize)
        row = currentRow(1) - listSize
        if(linkList(row,4) < depthThresh)
            classList = traceLinkageToBinary(linkList, row);
            featureList{row} = rfFeatureSelection(data(classList>0,:), classList(classList>0));
            currentRow = [currentRow(2:end); linkList(row,1); linkList(row,2)];
        end
    else
        currentRow = currentRow(2:end);
    end
end

end