topic Distance Calculation in SAS Programming

Distance Calculation

mauri0623 — Wed, 14 Mar 2018 13:09:51 GMT

SAS_M SAS_F
myzip facility_zip
22180 33180
22181 33181
22182 33182
22183 33183
22184 .....
22185 numobs=270 zip codes
22186
22187
22188
22189
22190
.....
.....
.....
numobs=3000 zip codes

SAS_M and SAS_F are two SAS datasets each containing only a column of zip. SAS_M has 3000 rows and SAS_F has 270 rows.
ZIPCITYDISTANCE function works fine for a pair of zips. I need some sort of looping that can do the followings,

ZIPCITYDISTANCE(22180,33180) then
ZIPCITYDISTANCE(22180,33181) then
ZIPCITYDISTANCE(22180,33182) then
ZIPCITYDISTANCE(22180,33183) then
ZIPCITYDISTANCE(22181,33180) then
ZIPCITYDISTANCE(22181,33181) and soforth.
............................

As said the ZIPCITYDISTANCE function works fine when you only have a pair of zips. In my case I need the distance between
first zip from SAS_M to every row of zip in the SAS_F dataset and then the second zip from the SAS_M to every zip rows of
SAS_F dataset and soforth.

Thank you for your help.
* ZIPCITYDISTANCE function;
d1 = zipcitydistance(12203,27513);

Re: Distance Calculation

Kurt_Bremser — Wed, 14 Mar 2018 13:26:05 GMT

SQL should do it:

proc sql;
create table want as
select
  a.myzip,
  b.facility_zip,
  zipcitydistance(a.myzip,b.facility_zip) as distance
from sas_m a, sas_f b;
quit;

as it automatically creates a cartesian join.

Re: Distance Calculation

RW9 — Wed, 14 Mar 2018 13:28:53 GMT

If you need all rows of the second calculated against the rows from the first, them merge 1 to many:

data want;
  merge myzip (in=a) facility_zip (in=b;
  if a;
  diff=zipdistance(myzip,facility_zip);
run;

Note, not tested, post test data in the form of a datastep, and show required output. We are not here to type this in for you!

Re: Distance Calculation

mauri0623 — Wed, 14 Mar 2018 13:49:42 GMT

Thank you. I will try.