Analyze the results from computations on the uk2006 graph

The uk2006 graph is a large ... This code will be documented and explained soon...

Contents

Setup the experiment

This experiment should be run from the rapr/experiments/comparison directory

cwd = pwd;
dirtail = 'experiments/comparison';
if strcmp(cwd(end-length(dirtail)+1:end),dirtail) == 0
    warning('rapr:dir','%s should be executed from rapr/%s\n',mfilename,dirtail);
end
addpath('../tools');
addpath('../../matlab');

Load the results

results = [];
results(end+1).x = binload('../webspam/uk-2006-05.50.pr');
results(end).name = 'pagerank 0.5';
results(end+1).x = binload('../webspam/uk-2006-05.85.pr');
results(end).name = 'pagerank 0.85';
results(end+1).x = binload('../webspam/uk-2006-05.95.pr');
results(end).name = 'pagerank 0.95';
results(end+1).x = binload('../uk-spam-2006/uk-2006-05-beta-2-16-21.mean');
results(end).name = 'ex beta(2,16)';
results(end+1).x = binload('../uk-spam-2006/uk-2006-05-beta-1-1-10.mean');
results(end).name = 'ex beta(1,1)';
results(end+1).x = binload('../uk-spam-2006/uk-2006-05-beta-2-16-21.std');
results(end).name = 'stdx beta(2,16)';
results(end+1).x = binload('../uk-spam-2006/uk-2006-05-beta-1-1-10.std');
results(end).name = 'stdx beta(1,1)';
nx = length(results);
ndiff = zeros(nx,nx);
ntau = zeros(nx,nx);
nttau = zeros(nx,nx);
ttol = 1e-10;

Compute all the pairwise comparisons between these vectors

for i=1:nx
    for j=i+1:nx
        x = results(i).x; y = results(j).x;
        tx = round(x/ttol)*ttol; ty = round(u/ttol)*ttol;
        ndiff(i,j) = norm(x-y,1);
        ntau(i,j) = ktau(x,y);
        nttau(i,j) = ktau(tx,ty);
        save 'uk-2006-results.mat' ndiff ntau nttau;
    end
end
Error using ==> evalin
Undefined function or variable 'u'.

Load and plot the output. For the saved scripts, this code wasn't complete.