如何在循环中为R中的autoKrige {automap}编写公式

Question

我有一个很大的矩阵，列名如下：

colid=vector(length = 60)
for(i in 1"60) {
  colid[i]=paste0("V",i)
}

当我在autoKrige中使用automap函数时，必须用这样的格式编写公式: V1~1。当我试图使用一个循环来完成这个操作时，会出现一个错误：

library(automap)
value=list()
for(i in 1:60) {
  value[[i]]=autoKrige(colid[i]~1,Mydata,new_data = newgrid)
}

Error: too many spatial dimensions: 3068 In addition: Warning message:
NAs introduced by coercion 

然后我试着解决这个问题，我测试公式：

> colid[10]~1
colid[10] ~ 1

因此，问题在于公式格式不正确。我想知道我怎么能解决这个问题？非常感谢。

Answer 1

使用函数as.formula将字符转换为公式：

colid=vector(length = 60)
value=list()
for(i in 1:60) {
  colid[i] <- paste0("V",i, " ~ 1")
  value[[i]]=autoKrige(as.formula(colid[i]),Mydata,new_data = newgrid)
}