Announcement

Collapse
No announcement yet.
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • working with clusters and creating buffers to match 2 geocoded locations

    Dear Statalist Users,

    I have a panel data (let's call it dataset1) on refugee camps from 2000 to 2016 for 15 countries, where I have longitude and latitude of the camps. In this dataset, I also have total number of refugees prsent in the given longitudes and latitudes.

    Then, I have another dataset (let's call it dataset2) which is at individual level and an appended crossectional form. The time interval is still between 2000 and 2016 and the observations come from the same 15 countries. I have each individuals' longitude and latitude information as well.
    By using "dataex", I provided a small example from both datasets. My question here is that, I would like to merge my individual dataset with refugee information (i.e., grandtotal) by creating buffers at the cluster level.

    I first created cluster_id for dataset1 and then dataset2 using the following command: egen cluster_id=group(longitude latitude). For example, to proxy refugee presence, I would like to create a buffer in my individual dataset (lets say it is within 20 km) and I would like to know whether there is a refugee presence within 20km of an individual. I have used the following command but it did not worked at all. The code I was using after creating buffers: forvalues i = 1(100)`max_cluster_id' {
    use "C:\Users\cansu\Dropbox\PC\Desktop\AB_All_Roun ds.d ta", clear
    keep if cluster_id >= `i' & cluster_id < `i' + `chunk_size'
    gen clat = latitude
    gen clong = longitude
    cross using `camps'
    geodist clat clong camplat camplong, gen(distance)
    gen within_buffer = distance <= `buffer_radius'
    keep if within_buffer
    tempfile temp_result
    save `temp_result'
    append using final_results, force
    save final_results, replace
    } )

    Hıge thanks in advance.

    Example from datase1:


    Code:
    * Example generated by -dataex-. To install: ssc install dataex
    clear
    input float year str32 ctr_n double(Long Lat grandtotal) float cluster_id
    2002 "Benin"   2.021   6.461    2 61
    2002 "Benin"   2.021   6.461    2 61
    2002 "Benin"   2.021   6.461   38 61
    2002 "Benin"   2.021   6.461    1 61
    2002 "Benin"   2.021   6.461   86 61
    2002 "Benin"   2.021   6.461  121 61
    2002 "Benin"   2.021   6.461    4 61
    2002 "Benin"   2.021   6.461    3 61
    2002 "Benin"   2.021   6.461  229 61
    2002 "Benin"   2.021   6.461   33 61
    2002 "Benin"   2.021   6.461    1 61
    2002 "Benin"   2.021   6.461    2 61
    2002 "Benin"   2.021   6.461    2 61
    2002 "Benin"   2.021   6.461  645 61
    2003 "Benin"   2.021   6.461    2 61
    2003 "Benin"   2.021   6.461    2 61
    2003 "Benin"   2.021   6.461   38 61
    2003 "Benin"   2.021   6.461    1 61
    2003 "Benin"   2.021   6.461   85 61
    2003 "Benin"   2.021   6.461  106 61
    2003 "Benin"   2.021   6.461    4 61
    2003 "Benin"   2.021   6.461    3 61
    2003 "Benin"   2.021   6.461  250 61
    2003 "Benin"   2.021   6.461   30 61
    2003 "Benin"   2.021   6.461    1 61
    2003 "Benin"   2.021   6.461    2 61
    2003 "Benin"   2.021   6.461    2 61
    2003 "Benin"   2.021   6.461  633 61
    2004 "Benin"   2.021   6.461    2 61
    2004 "Benin"   2.021   6.461    2 61
    2004 "Benin"   2.021   6.461   38 61
    2004 "Benin"   2.021   6.461   81 61
    2004 "Benin"   2.021   6.461   99 61
    2004 "Benin"   2.021   6.461    4 61
    2004 "Benin"   2.021   6.461    3 61
    2004 "Benin"   2.021   6.461  242 61
    2004 "Benin"   2.021   6.461   26 61
    2004 "Benin"   2.021   6.461    1 61
    2004 "Benin"   2.021   6.461    1 61
    2004 "Benin"   2.021   6.461    2 61
    2004 "Benin"   2.021   6.461  634 61
    2005 "Benin"   2.021   6.461    2 61
    2005 "Benin"   2.021   6.461    2 61
    2005 "Benin"   2.021   6.461   21 61
    2005 "Benin"   2.021   6.461   84 61
    2005 "Benin"   2.021   6.461   97 61
    2005 "Benin"   2.021   6.461    4 61
    2005 "Benin"   2.021   6.461    3 61
    2005 "Benin"   2.021   6.461  238 61
    2005 "Benin"   2.021   6.461   26 61
    2005 "Benin"   2.021   6.461    1 61
    2005 "Benin"   2.021   6.461    1 61
    2005 "Benin"   2.021   6.461    2 61
    2005 "Benin"   2.021   6.461  514 61
    2006 "Benin" 1.76048 6.72748 9248 58
    2006 "Benin" 1.88198 6.40764 1141 59
    2006 "Benin"   2.021   6.461    2 61
    2006 "Benin"   2.021   6.461    1 61
    2006 "Benin"   2.021   6.461   22 61
    2006 "Benin"   2.021   6.461   84 61
    2006 "Benin"   2.021   6.461   97 61
    2006 "Benin"   2.021   6.461    4 61
    2006 "Benin"   2.021   6.461    3 61
    2006 "Benin"   2.021   6.461  244 61
    2006 "Benin"   2.021   6.461   27 61
    2006 "Benin"   2.021   6.461    1 61
    2006 "Benin"   2.021   6.461    1 61
    2006 "Benin"   2.021   6.461    2 61
    2006 "Benin"   2.021   6.461  516 61
    2007 "Benin" 1.76048 6.72748 5303 58
    2007 "Benin"   2.021   6.461    2 61
    2007 "Benin"   2.021   6.461   19 61
    2007 "Benin"   2.021   6.461   31 61
    2007 "Benin"   2.021   6.461   35 61
    2007 "Benin"   2.021   6.461    2 61
    2007 "Benin"   2.021   6.461  180 61
    2007 "Benin"   2.021   6.461    9 61
    2007 "Benin"   2.021   6.461    1 61
    2007 "Benin"   2.021   6.461  437 61
    2008 "Benin" 1.76048 6.72748 3241 58
    2008 "Benin"   2.021   6.461   13 61
    2008 "Benin"   2.021   6.461   29 61
    2008 "Benin"   2.021   6.461   15 61
    2008 "Benin"   2.021   6.461    2 61
    2008 "Benin"   2.021   6.461  146 61
    2008 "Benin"   2.021   6.461    4 61
    2008 "Benin"   2.021   6.461    1 61
    2008 "Benin"   2.021   6.461  431 61
    2009 "Benin" 1.76048 6.72748 2976 58
    2009 "Benin"   2.021   6.461   12 61
    2009 "Benin"   2.021   6.461   15 61
    2009 "Benin"   2.021   6.461   14 61
    2009 "Benin"   2.021   6.461   94 61
    2009 "Benin"   2.021   6.461    6 61
    2009 "Benin"   2.021   6.461    1 61
    2009 "Benin"   2.021   6.461  234 61
    2010 "Benin" 1.76048 6.72748 2783 58
    2010 "Benin"   2.021   6.461   12 61
    2010 "Benin"   2.021   6.461   13 61
    2010 "Benin"   2.021   6.461   13 61
    end

    Example from dataset2:

    Code:
    * Example generated by -dataex-. To install: ssc install dataex
    clear
    input str8 respno str32 ctr_n float yr double(longitude latitude) float cluster_id
    "BEN0995" "Benin" 2005 3.1050999  12.10646 2236
    "BEN1078" "Benin" 2005   3.38327  11.86819 2257
    "BEN1001" "Benin" 2005   2.72104  10.81281 2199
    "BEN0985" "Benin" 2005 2.2929499  11.18131 2063
    "BEN1189" "Benin" 2005 3.1050999  12.10646 2236
    "BEN0984" "Benin" 2005 2.2929499  11.18131 2063
    "BEN1186" "Benin" 2005   2.72104  10.81281 2199
    "BEN1140" "Benin" 2005   3.22174  11.72937 2246
    "BEN1130" "Benin" 2005   2.77878  10.72015 2207
    "BEN1143" "Benin" 2005   2.72104  10.81281 2199
    "BEN1184" "Benin" 2005   3.38327  11.86819 2257
    "BEN1071" "Benin" 2005    2.9061  11.08049 2217
    "BEN1083" "Benin" 2005   2.77878  10.72015 2207
    "BEN0998" "Benin" 2005   3.22174  11.72937 2246
    "BEN1131" "Benin" 2005 2.9386101  11.13417 2220
    "BEN1069" "Benin" 2005 2.9386101  11.13417 2220
    "BEN0987" "Benin" 2005   2.43856  11.29845 2110
    "BEN1000" "Benin" 2005   3.38327  11.86819 2257
    "BEN1190" "Benin" 2005   3.58459  10.78888 2272
    "BEN1072" "Benin" 2005   3.58459  10.78888 2272
    "BEN1110" "Benin" 2005   2.41646  11.23617 2097
    "BEN1146" "Benin" 2005 2.2929499  11.18131 2063
    "BEN1180" "Benin" 2005   3.22174  11.72937 2246
    "BEN1133" "Benin" 2005    2.9061  11.08049 2217
    "BEN1193" "Benin" 2005 2.9386101  11.13417 2220
    "BEN1129" "Benin" 2005   2.77878  10.72015 2207
    "BEN0997" "Benin" 2005   3.22174  11.72937 2246
    "BEN1138" "Benin" 2005 3.1050999  12.10646 2236
    "BEN1074" "Benin" 2005 3.1050999  12.10646 2236
    "BEN1141" "Benin" 2005   3.38327  11.86819 2257
    "BEN0991" "Benin" 2005    2.9061  11.08049 2217
    "BEN0994" "Benin" 2005 3.1050999  12.10646 2236
    "BEN1112" "Benin" 2005 2.2929499  11.18131 2063
    "BEN1075" "Benin" 2005 3.1050999  12.10646 2236
    "BEN1077" "Benin" 2005   3.22174  11.72937 2246
    "BEN1194" "Benin" 2005    2.9061  11.08049 2217
    "BEN0989" "Benin" 2005 2.9386101  11.13417 2220
    "BEN1062" "Benin" 2005 2.2929499  11.18131 2063
    "BEN0986" "Benin" 2005   2.43856  11.29845 2110
    "BEN1142" "Benin" 2005   3.38327  11.86819 2257
    "BEN1073" "Benin" 2005   3.58459  10.78888 2272
    "BEN1113" "Benin" 2005   2.43856  11.29845 2110
    "BEN1068" "Benin" 2005 2.9386101  11.13417 2220
    "BEN1144" "Benin" 2005   2.41646  11.23617 2097
    "BEN1065" "Benin" 2005   2.41646  11.23617 2097
    "BEN1134" "Benin" 2005    2.9061  11.08049 2217
    "BEN0999" "Benin" 2005   3.38327  11.86819 2257
    "BEN1183" "Benin" 2005   2.77878  10.72015 2207
    "BEN1076" "Benin" 2005   3.22174  11.72937 2246
    "BEN1067" "Benin" 2005   2.43856  11.29845 2110
    "BEN1066" "Benin" 2005   2.43856  11.29845 2110
    "BEN1188" "Benin" 2005 3.1050999  12.10646 2236
    "BEN1191" "Benin" 2005   3.58459  10.78888 2272
    "BEN1132" "Benin" 2005 2.9386101  11.13417 2220
    "BEN1145" "Benin" 2005 2.2929499  11.18131 2063
    "BEN1111" "Benin" 2005 2.2929499  11.18131 2063
    "BEN1147" "Benin" 2005   2.43856  11.29845 2110
    "BEN1195" "Benin" 2005 2.9386101  11.13417 2220
    "BEN0996" "Benin" 2005   2.41646  11.23617 2097
    "BEN1136" "Benin" 2005   3.58459  10.78888 2272
    "BEN1002" "Benin" 2005   2.72104  10.81281 2199
    "BEN0988" "Benin" 2005 2.9386101  11.13417 2220
    "BEN1192" "Benin" 2005    2.9061  11.08049 2217
    "BEN1185" "Benin" 2005   2.72104  10.81281 2199
    "BEN0982" "Benin" 2005   2.41646  11.23617 2097
    "BEN1182" "Benin" 2005   2.77878  10.72015 2207
    "BEN1003" "Benin" 2005   2.77878  10.72015 2207
    "BEN1064" "Benin" 2005   2.41646  11.23617 2097
    "BEN1181" "Benin" 2005   3.38327  11.86819 2257
    "BEN0990" "Benin" 2005    2.9061  11.08049 2217
    "BEN0983" "Benin" 2005   2.41646  11.23617 2097
    "BEN1079" "Benin" 2005   3.38327  11.86819 2257
    "BEN1081" "Benin" 2005   2.72104  10.81281 2199
    "BEN1179" "Benin" 2005   2.43856  11.29845 2110
    "BEN1139" "Benin" 2005   3.22174  11.72937 2246
    "BEN1109" "Benin" 2005   2.41646  11.23617 2097
    "BEN1004" "Benin" 2005   2.77878  10.72015 2207
    "BEN1080" "Benin" 2005   2.72104  10.81281 2199
    "BEN1070" "Benin" 2005    2.9061  11.08049 2217
    "BEN1135" "Benin" 2005   3.58459  10.78888 2272
    "BEN1128" "Benin" 2005   2.72104  10.81281 2199
    "BEN1063" "Benin" 2005 2.2929499  11.18131 2063
    "BEN0993" "Benin" 2005   3.58459  10.78888 2272
    "BEN1082" "Benin" 2005   2.77878  10.72015 2207
    "BEN1187" "Benin" 2005   3.22174  11.72937 2246
    "BEN1114" "Benin" 2005   2.43856  11.29845 2110
    "BEN0992" "Benin" 2005   3.58459  10.78888 2272
    "BEN1137" "Benin" 2005 3.1050999  12.10646 2236
    "BEN1121" "Benin" 2005   1.37962  10.30416 1874
    "BEN1156" "Benin" 2005 1.0630701  10.69623 1856
    "BEN1122" "Benin" 2005    1.4368  10.12554 1879
    "BEN1023" "Benin" 2005   1.10621  10.17711 1857
    "BEN1021" "Benin" 2005   1.37962  10.30416 1874
    "BEN1059" "Benin" 2005   1.73333      10.3 1923
    "BEN1100" "Benin" 2005    .94198  10.40492 1844
    "BEN1160" "Benin" 2005   1.26651  10.62118 1867
    "BEN1028" "Benin" 2005   1.52147 9.9956503 1886
    "BEN1047" "Benin" 2005   1.10782  10.35972 1858
    "BEN0973" "Benin" 2005   1.03711  10.84793 1853
    "BEN1060" "Benin" 2005   2.10479  10.82578 2015
    end

  • #2
    One thing that Forum FAQ makes clear is that it is useless to say that your code "did not work," because there are so many ways that code can fail. So with "did not work," somebody who wants to help you won't even know where to start troubleshooting. You need to say what went wrong. Did Stata crash? Or hang? Did the code halt with error messages? (If so, show the output you got including those error messages?) Did the code run without error messages but produce incorrect results? (If so, in what way were the results different from what was required? Show example results and, unless it is blatantly obvious even to people unfamiliar with the general topic, explain why and how they are wrong.)

    Comment


    • #3
      I named the first dataset as locations.

      Code:
      geonear respno latitude longitude using locations, neighbors(cluster_id Lat Long)
      It gives km distance to nearest camp.

      Comment

      Working...
      X