У меня есть список фреймов данных, из которых 2 переменные: 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?
Большое спасибо!