返回中的R的观察给定半径内的属性

Question

新手ř程序员奋力空间组件添加到大的数据集（> 1500个OBS）

鉴于列表：观察，对应的x/y点（未LAT /长）和属性;我怎样才能写出一个算法来返回给定半径（10）内所有观测值的变量“Z”的总和。我需要为每个观察做到这一点。

简化的数据报头的内容如下：

OBS  X  Y   Z** 
A  56.55  -289.65 
B  52.59  -287.82 
C  58.34  -284.58  

任何援助将不胜感激。谢谢。

Answer 1

的关键是在同一时间使用mapply()到申请一个函数米 ultiple载体 - 在这里：X和您的data.frame的Y分量。哪个功能？为data.frame中的所有观测结果添加Z变量，该变量最多不超过“定位点”（我们运行所有条目）的radius。

> set.seed(1) 
> foo <- data.frame(X=runif(100),Y=runif(100),Z=runif(100)) 
> radius <- 0.5 
> result <- mapply(function(X,Y,Z) 
    sum(foo$Z[(foo$X-X)^2+(foo$Y-Y)^2<=radius^2]),foo$X,foo$Y) 
> head(cbind(foo,result)) 
      X   Y   Z result 
1 0.2655087 0.6547239 0.2675082 24.99153 
2 0.3721239 0.3531973 0.2186453 30.51512 
3 0.5728534 0.2702601 0.5167968 28.15519 
4 0.9082078 0.9926841 0.2689506 12.10058 
5 0.2016819 0.6334933 0.1811683 22.50695 
6 0.8983897 0.2132081 0.5185761 19.05273 

请注意，当然每个点的Z本身都计入其半径。如果你不想这样做，你需要减去foo$Z。

看看?mapply。