DaveM@9: function linkList = treeLinkFeatures(data)
DaveM@9: 
DaveM@9: linkList = aglomCluster(data);
DaveM@9: 
DaveM@9: linkList(:,4) = 0;
DaveM@9: listSize = size(data,1);
DaveM@9: currentRow = [size(linkList,1)];
DaveM@9: 
DaveM@9: 
DaveM@9: while (linkList(currentRow(1),1) > listSize && linkList(currentRow(1),2) > listSize)
DaveM@9:     classList = traceLinkageToBinary(linkList, currentRow);
DaveM@9:     linkList(:,4) = rfFeatureSelection(data(classList>0,:), classList(classList>0));
DaveM@9:     currentRow = [currentRow(2:end); currentRow(1),1; currentRow(1),2];
DaveM@9: end
DaveM@9: 
DaveM@9: end