load risk_rate.dat

tic;

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

for i = 1:100
    
    [find,vaf] = ...
        atreefnd(risk_rate,randperm(18));
    
    store_vaf(i) = vaf;
    
    if (vaf > best_vaf)
        
        best_vaf = vaf;
        best_find = find;
        
       
    end
    
end

sorted_vafs = sort(store_vaf');

sorted_vafs;

best_vaf;
best_find;

[ulmetric,ctmetric] = atreedec(best_find,0.0);

[orderprox,orderperm] = ultraorder(ulmetric);

best_risk_additivetree = best_find;




toc;
