load risk_rate.dat

tic;

best_vaf = 0.0;
store_vaf = zeros(1000,1);

for i = 1:1000
    
    [find,vaf,targone,targtwo] = ...
        ms_biatreefnd(risk_rate,randperm(18));
    
    store_vaf(i) = vaf;
    
    if (vaf > best_vaf)
        
        best_vaf = vaf;
        best_find = find;
        best_targone = targone;
        best_targtwo = targtwo;
       
    end
    
end

sorted_vafs = sort(store_vaf');

sorted_vafs;

best_vaf
best_find
best_targone
best_targtwo

[ulmetricone,ctmetricone] = atreedec(best_targone,0.0)
[ulmetrictwo,ctmetrictwo] = atreedec(best_targtwo,0.0)

[orderproxone,orderpermone] = ultraorder(ulmetricone)
[orderproxtwo,orderpermtwo] = ultraorder(ulmetrictwo)


toc;
