load lincoln_mapquest_time.dat

objectlabels{1} = 'Urbana';
objectlabels{2} = 'Danville';
objectlabels{3} = 'Paris';
objectlabels{4} = 'Shelbyville';
objectlabels{5} = 'Sullivan';
objectlabels{6} = 'Decatur';
objectlabels{7} = 'Taylorville';
objectlabels{8} = 'Springfield';
objectlabels{9} = 'Pekin';
objectlabels{10} = 'Metamora';
objectlabels{11} = 'Bloomington';
objectlabels{12} = 'Mt. Pulaski';
objectlabels{13} = 'Clinton';
objectlabels{14} = 'Monticello';




PhyloTree = seqlinkage(squareform(lincoln_mapquest_time),'complete',...
    objectlabels);

plot(PhyloTree)

PhyloAdditiveTree = seqneighjoin(squareform(lincoln_mapquest_time), ...
   'equivar', objectlabels);

plot(PhyloAdditiveTree)

