如何从R中的mapply或Map函数中的数据列表中选择2个变量？

Question

我有一个数据框列表，其中有两个变量：utmX和utmY将用于此分析。

在下面的代码中，我试图：1）写L.fun（）2）将此L.fun（）应用于数据帧列表3）分别绘制每个数据帧的结果。

L.fun <- function(x, y){
  window <- ripras(x, y)
  p.patt <- ppp(x, y, window=window)
  L <- Lest(p.patt, correction="Ripley")
  return(L)
}

 data.list <- list(data1,data2,data3...data30)

 L.res <- mapply(L.fun, data.list$utmX, data.list$utmY)  # empty
 L.res <- mapply(L.fun, x$utmX, x$utmY) # x object not found
 L.res <- Map(L.fun, data.list)         # arg y is missing 
 plot(L.res)

此外，我通过在函数内部包含x和y来尝试这个，仍然无法正常工作。

 L.res <- lapply(data.list, function(x) {
    for(i in 1:nrow(x)){
       window <- ripras(x$utmX, x$utmY)
       p.patt <- ppp(x$utmX, x$utmY, window=window)
       x$L <- Lest(p.patt, correction="Ripley")
   }
   x
 })

所以我的问题是如何从数据列表中选择2个用于mappy和Map函数的变量？

非常感谢你！

Answer 1

稍微改动你的第二次尝试应该有效。

library(spatstat)

首先生成2个假数据集的列表：

data.list <- replicate(2, data.frame(utmX=runif(20), utmY=runif(20)), simplify = FALSE)

运行列表并将L函数应用于列表中的每个数据集：

L.res <- lapply(data.list, function(x) {
  window <- ripras(x$utmX, x$utmY)
  p.patt <- ppp(x$utmX, x$utmY, window=window)
  return(Lest(p.patt, correction="Ripley"))
})

结果是两个L函数的列表，如图所示：

plot.anylist(L.res)

由reprex package创建于2019-03-14（v0.2.1）

Answer 2

另一种选择：

coor <- lapply(data.list, function(x) cbind(x$utmX, x$utmY))
pats <- lapply(coor, as.ppp, W=ripras)
Lfuns <- anylapply(pats, Lest, correction="Ripley")
plot(Lfuns)