2017-06-20 64 views
0

我在R中有一组GPS坐标,我想将其视为对应于感兴趣的患者的另一组GPS坐标的“曝光”。我如何在R中进行接近性分析以将患者分成两组:在接触坐标x米内的患者和距离较远的患者?R中的接近度分析

回答

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) 

强烈灵感来自:Calculate distance from GPS data

+0

这完美的作品。谢谢!唯一值得注意的是,sDistsN1计算矩阵和一个点之间的距离,而spDists计算城镇矩阵中每个点之间的距离。 – novice