Many thanks to Jon Peck. This is now working.
The extendedTransforms.zip was difficult to find but the resources on
devcentral are worth browsing.
*The purpose of this is crosscheck 3 sources for the location of hospitals.
data list list/ ID(f3) acc1 (f1) acc2(f1) latitude1 longitude1 latitude2
longitude2 latitude3 longitude3 (6f11.7).
begin data.
2 6 8 41.6674894 -70.5158194 41.667489 -70.515819
41.6697540 -70.5155370
20 8 6 42.7115620 -71.0061448 42.711562 -71.006145
42.7110741 -71.0034720
41 9 6 42.3902597 -72.5216735 42.390260 -72.521674
42.3899722 -72.5212817
70 8 9 42.2562632 -71.0792062 42.256263 -71.079206
42.2563830 -71.0792450
104 9 8 42.3936720 -71.0333220 42.393672 -71.033322
42.3936720 -71.0333220
166 8 9 42.2754785 -71.0936511 42.275478 -71.093651
42.2753140 -71.0937040
233 6 8 42.3695410 -71.0366920 42.369541 -71.036692
42.3695327 -71.0366994
308 9 8 42.3085910 -71.1149198 42.308591 -71.114920
42.3085910 -71.1149198
326 5 6 42.3398361 -71.0933014 42.339836 -71.093301
42.3398106 -71.0930371
379 9 8 41.7280625 -70.5854771 41.728063 -70.585477
41.7280625 -70.5854771
421 9 8 41.9373140 -70.0207360 41.937314 -70.020736
41.9373140 -70.0207360
end data.
var labels
latitude1 'latitude from source 1'/
longitude1 'longitude from source 1'/
latitude2 'latitude from source 2'/
longitude2 'longitude from source 2'/
latitude3 'latitude from source 3'/
longitude3 'longitude from source 3'.
numeric miles1_2 miles1_3 miles2_3 (f10.5).
spssinc trans result=miles1_2
/formula "extendedTransforms.ellipseDist(latitude1, longitude1,
latitude2, longitude2, inradians=False)".
spssinc trans result=miles1_3
/formula "extendedTransforms.ellipseDist(latitude1, longitude1,
latitude3, longitude3, inradians=False)".
spssinc trans result=miles2_3
/formula "extendedTransforms.ellipseDist(latitude2, longitude2,
latitude3, longitude3, inradians=False)".
numeric feet1_2 feet1_3 feet2_3 (f6).
compute feet1_2 = miles1_2 *5280.
compute feet1_3 = miles1_3 *5280.
compute feet2_3 = miles2_3 *5280.
list /variables = ID miles1_2 miles1_3 miles2_3 feet1_2 feet1_3 feet2_3 .
--
Art Kendall
Social Research Consultants
=====================
To manage your subscription to SPSSX-L, send a message to
[hidden email] (not to SPSSX-L), with no body text except the
command. To leave the list, send the command
SIGNOFF SPSSX-L
For a list of commands to manage subscriptions, send the command
INFO REFCARD
Art Kendall
Social Research Consultants
Locked
