已知每个农户的经纬度，求距离每个农户最近的五个农户，hhcode代表农户

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

hhdata-311-new-2.dta (362.66 KB, 需要: 8 个论坛币)

2015-6-4 15:13:29 上传

已知每个农户的经纬度，求距离每个农户最近的五个农户。hhcode代表农户
需要: 8 个论坛币  [购买]

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：code 经纬度 ODE COD 经纬度

哪个变量是经纬度？

andruw 发表于 2015-6-4 15:40
哪个变量是经纬度？

latitude
longitude

andruw 发表于 2015-6-4 15:40
哪个变量是经纬度？

longitude
latitude

1. gen distance1 = .
   
2. gen distance2 = .
   
3. gen distance3 = .
   
4. gen distance4 = .
   
5. gen distance5 = .
   
6. 
   
7. gen close1 = .
   
8. gen close2 = .
   
9. gen close3 = .
   
10. gen close4 = .
    
11. gen close5 = .
    
12. 
    
13. 
    
14. 
    
15. local n = _N
    
16. quietly{
    
17. forval i = 1/`n' {
    
18.         forval j = 1/`n' {
    
19.                 if hhcode[`i'] & hhcode[`j'] & (`i' != `j'){
    
20.                         scalar d = (longitude[`i'] - longitude[`j'])^2 + (latitude[`i'] - latitude[`j'])^2
    
21.                          if d < distance1[`i'] {
    
22.                           replace distance5 = distance4 in `i'
    
23.                           replace close5 = close4 in `i'
    
24.                           replace distance4 = distance3 in `i'
    
25.                           replace close4 = close3 in `i'
    
26.                           replace distance3 = distance2 in `i'
    
27.                           replace close3 = close2 in `i'
    
28.                           replace distance2 = distance1 in `i'
    
29.                           replace close2 = close1 in `i'
    
30.                           replace distance1 = d in `i'
    
31.                           replace close1 = hhcode[`j'] in `i'
    
32.                           }
    
33.                          if d < distance2[`i'] & d > distance1[`i'] {
    
34.                           replace distance5 = distance4 in `i'
    
35.                           replace close5 = close4 in `i'
    
36.                           replace distance4 = distance3 in `i'
    
37.                           replace close4 = close3 in `i'
    
38.                           replace distance3 = distance2 in `i'
    
39.                           replace close3 = close2 in `i'
    
40.                           replace distance2 = d in `i'
    
41.                           replace close2 = hhcode[`j'] in `i'
    
42.                           }
    
43.                          if d < distance3[`i'] & d > distance1[`i']  & d > distance2[`i'] {
    
44.                           replace distance5 = distance4 in `i'
    
45.                           replace close5 = close4 in `i'
    
46.                           replace distance4 = distance3 in `i'
    
47.                           replace close4 = close3 in `i'
    
48.                           replace distance3 = d in `i'
    
49.                           replace close3 = hhcode[`j'] in `i'
    
50.                           }
    
51.                          if d < distance4[`i'] & d > distance3[`i'] & d > distance1[`i']  & d > distance2[`i'] {
    
52.                          replace distance5 = distance4 in `i'
    
53.                          replace close5 = close4 in `i'
    
54.                          replace distance4 = d in `i'
    
55.                          replace close4 = hhcode[`j'] in `i'
    
56.                          }
    
57.                          if d < distance5[`i'] & d > distance4[`i'] & d > distance3[`i'] & d > distance1[`i']  & d > distance2[`i'] {
    
58.                          replace distance5 = distance4 in `i'
    
59.                          replace close5 = hhcode[`j'] in `i'
    
60.                          }
    
61.                 }
    
62.         }
    
63. }
    
64. }
    
65. 
    
66. 
    
67. drop distance1 distance2 distance3 distance4 distance5

复制代码

andruw 发表于 2015-6-4 16:57

抛砖引玉，ssc应该有简单的方法。