r - 如何在 Rstudio 的 georob 包中插入更密集采样的辅助变量

Question

我需要根据 georob 包中的辅助变量来预测土壤变量。

我的个人数据集有 200 个观察值，我的辅助变量集有 19940 个数据，但是在代码中，我无法输入辅助变量的坐标作为预测点。

dat= read.csv("malhas amostrais/solo_200.csv", sep = ",")

covar = read.csv("../dados/csv/variaveis_auxiliares.csv", sep = ";")

ku_georob_cpeso <- georob(argila ~ CV + CH + dist_bebedouros + Eca_0.5m + Eca_1m + elevacao + IH_0.5m + sd_ndvi_01 + sd_ndvi_02 + twi + S_P_T + sd_b4 +sd_b5 + sd_b6+ sd_b7,
   data= dat,
   locations= ~ x + y,
   variogram.model="RMexp",
   param=c(variance=200, nugget=600, scale=150),
   verbose = 3,
   psi.func = "huber") 

ku_georob_cpeso <- georob(argila ~ CV + CH + dist_bebedouros + Eca_0.5m + Eca_1m + elevacao + IH_0.5m + sd_ndvi_01 + sd_ndvi_02 + twi + S_P_T + sd_b4 +sd_b5 + sd_b6+ sd_b7,
   data= dat1,
   subset = cova,
   locations= ~ x + y,
   variogram.model="RMexp",
   param=c(variance=200, nugget=600, scale=150),+ verbose = 3,
   psi.func = "huber") 

我收到错误：

xj[i] 中的错误：无效的下标类型“列表”

标签： rrstudio