需要构建基于距离的空间权重矩阵,下面图一为我本人按照下面的stata命令做的各省的地理坐标结果,但发现我自己做的结果(如图1)跟网上做的结果x_stub 和 y_stub所对应的数值刚好反了。不知道x_stub 和 y_stub是不是各自对应的纬度和经度?如果本来应该是x_stub对应纬度,y_stub对应经度的话,那我的结果是不是做错了?应该怎么修改呢?求助。。
cd c:\w
shp2dta using w, database(data_db) coordinates(data_xy) genid(weightid) gencentroids(stub) replace
use data_db, clear
rename 省 province
list province x_stub y_stub in 1/30
图一(一部分) :
--------------------------------------+
| province x_stub y_stub |
|--------------------------------------|
1. | 北京市 116.4123 40.18554 |
2. | 天津市 117.3341 39.29343 |
3. | 河北省 116.1374 39.54777 |
4. | 山西省 112.2895 37.57175 |
5. | 内蒙古自治区 113.9221 44.08926 |
|--------------------------------------|
6. | ?? 125.5647 40.68582 |
7. | 辽宁省 122.6026 41.28811 |
8. | 吉林省 126.1884 43.66769 |
9. | 黑龙江省 127.7662 47.86168 |
10. | 上海市 121.4692 31.24011 |
|--------------------------------------|
11. | 江苏省 119.4755 32.96331 |
12. | 浙江省 120.1108 29.19609 |
13. | 安徽省 117.2265 31.82603 |
14. | 福建省 117.99 26.07406 |
15. | 江西省 115.7222 27.61407 |
|--------------------------------------|
16. | 山东省 118.1532 36.34917 |
17. | 河南省 113.6146 33.88233 |
18. | 湖北省 112.2715 30.97546 |
19. | 湖南省 111.7087 27.60982 |
20. | 广东省 113.4216 23.33963 |
|--------------------------------------|