R - 子集列表data.frames由矢量值

Posted 2021-03-27

我有一个SpatialLinesDataFrames列表，并希望基于将值与数值向量中的值进行比较来对列表进行子集化。

具体来说，我想删除具有data.frame的特定列（'lineID'）中的向量中包含的值之一的列表元素。可重复的例子：

#create list of single-feature SpatialLineDataFrame
library(raster)
l1 <- cbind(c(0,3), c(0,3))
l2 <- cbind(c(0, 13), c(0, 1))
l3 <- cbind(c(0, 24), c(0,22.5))
l4 <- cbind(c(0, 1), c(0,13))
l5 <- cbind(c(0, 6), c(0,6))
Sldf <- spLines(l1, l2, l3, l4, l5, attr=data.frame(lineID=1:5))

sldfl <- list()
sldfl[[1]] <- Sldf[1,]
sldfl[[2]] <- Sldf[2,]
sldfl[[3]] <- Sldf[3,]
sldfl[[4]] <- Sldf[4,]
sldfl[[5]] <- Sldf[5,]

#create numeric vector
m <- c(1,3,5,7,10)

#attempt to keep only list elements that are not equal to any 
#of the values contained in vector
final <- list()
for (i in 1:length(sldfl)) {
  for (j in 1:length(m)) {
    if (factor(sldfl[[i]]@data$lineID) != m[j]) {
      final[[i]] <- sldfl[[i]]
    }}}

循环的结果返回整个原始列表。我的循环出了什么问题？

答案

你基本上有两个向量，ids和m：

> ids
[1] 1 2 3 4 5
> m
[1]  1  3  5  7 10

并且基本上运行：

for(i in 1:length(ids)){
 for(j in 1:length(m)){
  if(i != m[j]){
    message("add ",i,j)
  }else{
    message("Not adding ",i,j)
  }
 }
}

运行它，你会看到它添加了许多元素，因为你正在使用m中的每个元素测试每个ID，并且至少有一个m元素不在ID中，因此添加了一个（和更多）。

你真正想要的是：

for(i in 1:length(ids)){

  if(!(i %in% m)){
    message("add ",i,j)
  }else{
    message("Not adding ",i,j)
  }
 }

打印：

Not adding 15
add 25
Not adding 35
add 45
Not adding 55

这增加了ID为2和4的元素，这些元素不在m中。

或者，使用基本R Filter函数，该函数通过列表元素上的函数减少列表：

> Filter(function(L){!(L@data$lineID  %in% m)}, sldfl)
[[1]]
class       : SpatialLinesDataFrame 
features    : 1 
extent      : 0, 13, 0, 1  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
variables   : 1
names       : lineID 
value       :      2 

[[2]]
class       : SpatialLinesDataFrame 
features    : 1 
extent      : 0, 1, 0, 13  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
variables   : 1
names       : lineID 
value       :      4