我在R中有一组GPS坐标,我想将其视为对应于感兴趣的患者的另一组GPS坐标的“曝光”。如何在R中进行接近度分析,将患者分为两组:在暴露坐标的x米范围内和远处的患者?
答案 0 :(得分:0)
我认为这样的事情应该有效:
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)