0

我正在处理sf对象，我有一个关于将用户定义的函数应用于数据框的问题。这是一个愚蠢的例子，但它类似于我试图解决一个更复杂的问题。我有一个名为的数据框names，它有两列，每列都有一个县的名称。sf我还阅读了包裹中包含的北卡罗来纳州数据。

require(sf)
name_1 <- c('Ashe','Alleghany','Surry')
name_2 <- c('Currituck','Northampton','Hertford')
names <- data.frame(cbind(name_1,name_2))
nc <- st_read(system.file("shape/nc.shp", package="sf"))

我要做的是创建一个函数，该函数沿着我的数据框的每一行，names从 name_1 列中取出名称，从 name_2 列中取出名称，并查看它们是否使用中的几何数据相交nc。我有：

check_intersection <- function(x){
  return(st_intersects(nc[nc$NAME== x$name_1,],nc[nc$NAME==x$name_2,],sparse = FALSE))
}
apply(names,1,check_intersection)

但这会产生一个错误Error in x$name_1 : $ operator is invalid for atomic vectors。对于 中的每一行，我如何告诉函数从适当的列中获取字符值names？

4

1 回答 1

1

尝试

apply(names, 1, 
  function(x) 
    st_intersects(nc[nc$NAME == x[1],], nc[nc$NAME == x[2],], 
       sparse = FALSE)
)

在匿名函数中，x是一个向量，而不是一个data.frame，因此$不起作用。

于 2017-09-26T20:00:26.627 回答