R中的接近分析

Question

我在R中有一组GPS坐标，我想将其视为对应于感兴趣的患者的另一组GPS坐标的“曝光”。如何在R中进行接近度分析，将患者分为两组：在暴露坐标的x米范围内和远处的患者？

Answer 1

我认为这样的事情应该有效：

library(sp)
exp <- data.frame(lat= 40.741895,long = -73.989308)
patients <- data.frame(lat = rnorm(10,exp$lat,0.1),long = rnorm(10,exp$long,0.1))
coordinates(patients) <- ~ long + lat
coordinates(exp) <- ~ long + lat
d <- spDistsN1(coordinates(patients),coordinates(exp),longlat = TRUE)

强烈启发：http://handlebarsjs.com/block_helpers.html