Есть ли способ перебрать переменные внутри объекта SpatialPointsDataFrame в R, используя интерполяцию IDW пакета gstat?
Я имею в виду что-то вроде этого
for (day in list_days) {
P.idw <- gstat::idw(day~1, P, newdata=grd, idp=2.0)
r <- raster(P.idw)
plot(r)
}
Изменить: потому что эта структура дает мне ошибку:
Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim, :
too many spatial dimensions: 58
In addition: Warning message:
In predict.gstat(g, newdata = newdata, block = block, nsim = nsim, :
NAs introduced by coercion
Но если я пишу напрямую имя переменной, все работает нормально:
P.idw <- gstat::idw(X2018.01.14~1, P, newdata=grd, idp=2.0)
P.idw <- gstat::idw(get(day)~1, P, newdata=grd, idp=2.0)
- person Marco Sandri   schedule 22.01.2018