load risk_rate.dat

prox = risk_rate;

n=18;

best_vaf = 0.0;

store_vaf = zeros(100,1);

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

store_vaf'

sorted_vafs = sort(store_vaf')

figure(1)


    z = seqneighjoin(squareform(find),'equivar')
    
    plot(z,'type','angular','Orientation','top')



