COMPUTE GEODIST FOR EACH ROW MANUAL

Now that you have a way to replicate the distances calculated by Stata's sp* commands, you can repeat what you were trying to do in #1 using great-circle distances. * redo with the radius used by spdistance * spdistance uses a different earth radius Great-circle distance (haversine formula, radius of 6371km) = 24.300728 km * redo using geodist (from SSC) using haversine formula (data currently use latitude and longitude) Linked shapefile: NUTS_RG_01M_2013_shp.dta spset, modify coordsys(latlong, kilometers)Д¬oordinates: _CY, _CX (latitude-and-longitude, kilometers) If you may want to have a look at the files, this is the link to the dataset I am using (already spset): I have been struggling all the day, but in fact I obtain different results. Now if everything works, I should observe the same entries in each ij position. Insheet using 'dW2.txt', delim(' ') clear Spmatrix export dW2 using dW2.txt, replace

COMPUTE GEODIST FOR EACH ROW MANUAL

I am using an example that I found into the sp manual (p.

Nonetheless, to be sure that I proceed correctly, I firstly try to replicate the inverse distance I have done previously. I generate an inverse distance function using spmatrix idistance.ДЇor some specific reasons, I need to use mata to create my own Matrix. I am using Stata 15 to carry out my spatial analysis.

Compute geodist for each row