Rick,

proc iml;
use maps.uscity;
read all var {City StateCode Pop} where(POP=0);
close;
print city;
call tabulate(levels,freq,StateCode);
idx=loc(freq=freq[<>]);
max=freq[idx];
l=levels[idx];
print max[r=l] ;
print freq[c=levels];
quit;