title 'Modeclus Analysis of 10 American Cities';
title2 'Based on Flying Mileages';

data mileages(type=distance);
   input (Atlanta Chicago Denver Houston LosAngeles
   Miami NewYork SanFrancisco Seattle DC) (5.)
   @53 City $15.;
   datalines;
   0                                                Atlanta
 587    0                                           Chicago
1212  920    0                                      Denver
 701  940  879    0                                 Houston
1936 1745  831 1374    0                            Los Angeles
 604 1188 1726  968 2339    0                       Miami
 748  713 1631 1420 2451 1092    0                  New York
2139 1858  949 1645  347 2594 2571    0             San Francisco
2182 1737 1021 1891  959 2734 2408  678    0        Seattle
 543  597 1494 1220 2300  923  205 2442 2329    0   Washington D.C.
;

*-----Fill in Upper Triangle of Distance Matrix---------------;
proc transpose out=tran;
   copy city;
run;

data mileages(type=distance);
   merge mileages tran;
   array var[*] atlanta--dc;
   array col[*] col1-col10;
   do i = 1 to 10;
      var[i] = sum(var[i], col[i]);
   end;
   drop col1-col10 _name_ i;
run;

*-----Clustering with K-Nearest-Neighbor Density Estimates-----;
proc modeclus data=mileages all m=1 k=3;
   id CITY;
run;

*------Clustering with Uniform-Kernel Density Estimates--------;
proc modeclus data=mileages all m=1 r=600 800;
   id CITY;
run;

*------Clustering Neighborhoods Extended to Nearest Neighbor--------;
proc modeclus data=mileages list m=1 ck=2 r=600 800;
   id CITY;
run;