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.