stata geodist 命令算两两之间距离

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

查了一下stata的命令是geodist。我一共400多个地区，有400多个经纬度
我编的是geodist lati loni latj lonj, gen (v)但出来的是
geodist lati loni latj lonj, gen (v)
Was expecting a variable, scalar, or number instead of lati
r(198);
请问一下这是什么意思？我刚刚接触软件这一方面还不是太熟悉，希望能给我些指点，谢谢。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏10 回帖

关键词：Stata ODIS tata IST GEO instead number 经纬度 软件

Stata的自带例子是可以

1. clear
   
2. set obs 100
   
3. set seed 1234
   
4. gen double lat1 = 37 + (41 - 37) * uniform()
   
5. gen double lon1 = -109 + (109 - 102) * uniform()
   
6. gen double lat2 = 37 + (41 - 37) * uniform()
   
7. gen double lon2 =-109 + (109 - 102) * uniform()
   
8. sum lat1 lon1 lat2 lon2
   
9. 
   
10. geodist lat1 lon1 lat2 lon2, gen(v)

复制代码

你看下help geodist文件  里边有个实例就是关于你这个问题的介绍，就是那个cross using 什么的

我仔细读了stata关于geodist的help文件，楼上坛友回答的代码就是里面示例的截取。但是不必限于那个例子，楼主的问题核心和我一样，是求“l两两”的距离，核心其实在cross命令的应用，利用cross可以把楼主400条数据变成400*400=160000条数据，然后分别用geodist求出距离。虽然是很多年前的问题，但好像论坛里相关问题说得不多，既然我也遇到了就和大家分享一下。

其实，这句话的意思就是，你的经纬度必须要interger或double或其他数值型变量才行，不能是string型

Andrewkings 发表于 2018-3-19 01:07
我仔细读了stata关于geodist的help文件，楼上坛友回答的代码就是里面示例的截取。但是不必限于那个例子，楼 ...

麻烦问下cross命令怎么使用啊  怎么把数据拆分   跪谢！！！

cheng123cc 发表于 2018-7-11 17:40
麻烦问下cross命令怎么使用啊  怎么把数据拆分   跪谢！！！

cross命令你可以查一下stata 的help，原理就是两两引用，我觉得不存在拆分的问题，因为cross语句就可以实现了。我现在也说不太清楚，不好意思哈

继续求问如何实现一组几百个点之间距离的测算

求问，求一个点到几百个点的距离