## Programming the statistical procedures from SAS

Solved
Occasional Contributor
Posts: 16

Accepted Solutions
Solution
‎07-17-2014 10:22 AM
Super User
Posts: 10,207

SQL + geodis()  -- sorry I have forgetton. You need check documentation. Here is dummy code.

proc sql;

create table distance as

select a.* ,b.id as _id, geodis(.......)  as distance

from have a, have b;

quit;

after that you need transpose your table distance to get final output.

Xia Keshan

All Replies
Solution
‎07-17-2014 10:22 AM
Super User
Posts: 10,207

SQL + geodis()  -- sorry I have forgetton. You need check documentation. Here is dummy code.

proc sql;

create table distance as

select a.* ,b.id as _id, geodis(.......)  as distance

from have a, have b;

quit;

after that you need transpose your table distance to get final output.

Xia Keshan

Super User
Posts: 11,810

Syntaxt for

GEODIST(latitude-1, longitude-1, latitude-2, longitude-2 <,options>)

so use Ksharps code and bring in the b version of lat and long :

b.longitude as blong, b.latitude as blat, geodist (latitude, longitude, blat, blong) as distance

you may likely want the store from the b data set as well so you know which distance which would bring up b.id as bstore in the select as well.

Occasional Contributor
Posts: 16